Showing 1–50 of 75 results for author: Rothkopf, A

Exact symmetry conservation and automatic mesh refinement in discrete initial boundary value problems

Authors: Alexander Rothkopf, W. A. Horowitz, Jan Nordström

Abstract: We present a novel solution procedure for initial boundary value problems. The procedure is based on an action principle, in which coordinate maps are included as dynamical degrees of freedom. This reparametrization invariant action is formulated in an abstract parameter space and an energy density scale associated with the space-time coordinates separates the dynamics of the coordinate maps and o… ▽ More We present a novel solution procedure for initial boundary value problems. The procedure is based on an action principle, in which coordinate maps are included as dynamical degrees of freedom. This reparametrization invariant action is formulated in an abstract parameter space and an energy density scale associated with the space-time coordinates separates the dynamics of the coordinate maps and of the propagating fields. Treating coordinates as dependent, i.e. dynamical quantities, offers the opportunity to discretize the action while retaining all space-time symmetries and also provides the basis for automatic adaptive mesh refinement (AMR). The presence of unbroken space-time symmetries after discretization also ensures that the associated continuum Noether charges remain exactly conserved. The presence of coordinate maps in addition provides new freedom in the choice of boundary conditions. An explicit numerical example for wave propagation in $1+1$ dimensions is provided, using recently developed regularized summation-by-parts finite difference operators. △ Less

Submitted 29 April, 2024; originally announced April 2024.

Comments: 50 pages, 12 figures

In-medium static inter-quark potential on high resolution quenched lattices

Authors: Rasmus N. Larsen, Gaurang Parkar, Alexander Rothkopf, Johannes Heinrich Weber

Abstract: We re-investigate the interactions between static color sources in a finite temperature gluonic medium using both high resolution isotropic and anisotropic quenched lattice QCD ensembles. The underlying ill-posed inverse problem, related to the extraction of spectral functions, is attacked with a range of different methods, including Bayesian inference, Padé interpolation and model fits. Among the… ▽ More We re-investigate the interactions between static color sources in a finite temperature gluonic medium using both high resolution isotropic and anisotropic quenched lattice QCD ensembles. The underlying ill-posed inverse problem, related to the extraction of spectral functions, is attacked with a range of different methods, including Bayesian inference, Padé interpolation and model fits. Among the latter we include a tail amended Gaussian ansatz and a HTL-inspired fit ansatz. We reconfirm the presence of a dominant low-lying spectral feature that supports the existence of a potential picture for the in-medium evolution of the static charges at late real times. Using the raw unmodified lattice data, all applicable methods show clear signs of screening of the real-part of the potential. After applying a subtraction procedure featured in a previous study we find however that screening disappears from the extracted potential. Paths towards the resolution of this puzzle are discussed. △ Less

Submitted 16 February, 2024; originally announced February 2024.

Comments: 54 pages, 40 figures

Report number: HU-EP-24/06-RTG

The crucial role of Lagrange multipliers in a space-time symmetry preserving discretization scheme for IVPs

Authors: Alexander Rothkopf, Jan Nordström

Abstract: In a recently developed variational discretization scheme for second order initial value problems ( J. Comput. Phys. 498, 112652 (2024) ), it was shown that the Noether charge associated with time translation symmetry is exactly preserved in the interior of the simulated domain. The obtained solution also fulfils the naively discretized equations of motions inside the domain, except for the last t… ▽ More In a recently developed variational discretization scheme for second order initial value problems ( J. Comput. Phys. 498, 112652 (2024) ), it was shown that the Noether charge associated with time translation symmetry is exactly preserved in the interior of the simulated domain. The obtained solution also fulfils the naively discretized equations of motions inside the domain, except for the last two grid points. Here we provide an explanation for the deviations at the boundary as stemming from the Lagrange multipliers used to implement initial and connection conditions. We show explicitly that the Noether charge including the boundary corrections is exactly preserved at its continuum value over the whole simulation domain, including the boundary points. △ Less

Submitted 15 December, 2023; originally announced December 2023.

Comments: 9 pages, 1 figure. arXiv admin note: text overlap with arXiv:2307.04490

Lattice real-time simulations with learned optimal kernels

Authors: Daniel Alvestad, Alexander Rothkopf, Dénes Sexty

Abstract: We present a simulation strategy for the real-time dynamics of quantum fields, inspired by reinforcement learning. It builds on the complex Langevin approach, which it amends with system specific prior information, a necessary prerequisite to overcome this exceptionally severe sign problem. The optimization process underlying our machine learning approach is made possible by deploying inherently s… ▽ More We present a simulation strategy for the real-time dynamics of quantum fields, inspired by reinforcement learning. It builds on the complex Langevin approach, which it amends with system specific prior information, a necessary prerequisite to overcome this exceptionally severe sign problem. The optimization process underlying our machine learning approach is made possible by deploying inherently stable solvers of the complex Langevin stochastic process and a novel optimality criterion derived from insight into so-called boundary terms. This conceptual and technical progress allows us to both significantly extend the range of real-time simulations in 1+1d scalar field theory beyond the state-of-the-art and to avoid discretization artifacts that plagued previous real-time field theory simulations. Limitations of and promising future directions are discussed. △ Less

Submitted 12 October, 2023; originally announced October 2023.

Comments: 5 pages, 5 figures

Un-screened forces in Quark-Gluon Plasma?

Authors: Alexei Bazavov, Daniel Hoying, Rasmus N. Larsen, Swagato Mukherjee, Peter Petreczky, Alexander Rothkopf, Johannes Heinrich Weber

Abstract: We study the correlator of temporal Wilson lines at non-zero temperature in 2+1 flavor lattice QCD with the aim to define the heavy quark-antiquark potential at non-zero temperature. For temperatures $153~{\rm MeV} \leq T \leq 352~{\rm MeV}$ the spectral representation of this correlator is consistent with a broadened peak in the spectral function, position or width of which then defines the real… ▽ More We study the correlator of temporal Wilson lines at non-zero temperature in 2+1 flavor lattice QCD with the aim to define the heavy quark-antiquark potential at non-zero temperature. For temperatures $153~{\rm MeV} \leq T \leq 352~{\rm MeV}$ the spectral representation of this correlator is consistent with a broadened peak in the spectral function, position or width of which then defines the real or imaginary parts of the heavy quark-antiquark potential at non-zero temperature, respectively. We find that the potential's real part is not screened contrary to the widely-held expectations. We comment on how this fact may modify the picture of quarkonium melting in the quark-gluon plasma. △ Less

Submitted 9 April, 2024; v1 submitted 31 August, 2023; originally announced August 2023.

Journal ref: Physical Review D 109, 074504 (2024)

A symmetry and Noether charge preserving discretization of initial value problems

Authors: Alexander Rothkopf, Jan Nordström

Abstract: Taking insight from the theory of general relativity, where space and time are treated on the same footing, we develop a novel geometric variational discretization for second order initial value problems (IVPs). By discretizing the dynamics along a world-line parameter, instead of physical time directly, we retain manifest translation symmetry and conservation of the associated continuum Noether c… ▽ More Taking insight from the theory of general relativity, where space and time are treated on the same footing, we develop a novel geometric variational discretization for second order initial value problems (IVPs). By discretizing the dynamics along a world-line parameter, instead of physical time directly, we retain manifest translation symmetry and conservation of the associated continuum Noether charge. A non-equidistant time discretization emerges dynamically, realizing a form of automatic adaptive mesh refinement (AMR), guided by the system symmetries. Using appropriately regularized summation by parts finite difference operators, the continuum Noether charge, defined via the Killing vector associated with translation symmetry, is shown to be exactly preserved in the interior of the simulated time interval. The convergence properties of the approach are demonstrated with two explicit examples. △ Less

Submitted 10 July, 2023; originally announced July 2023.

Comments: 35 pages, 12 figures

Kernel controlled real-time Complex Langevin simulation

Authors: Daniel Alvestad, Rasmus Larsen, Alexander Rothkopf

Abstract: This study explores the utility of a kernel in complex Langevin simulations of quantum real-time dynamics on the Schwinger-Keldysh contour. We give several examples where we use a systematic scheme to find kernels that restore correct convergence of complex Langevin. The schemes combine prior information we know about the system and the correctness of convergence of complex Langevin to construct a… ▽ More This study explores the utility of a kernel in complex Langevin simulations of quantum real-time dynamics on the Schwinger-Keldysh contour. We give several examples where we use a systematic scheme to find kernels that restore correct convergence of complex Langevin. The schemes combine prior information we know about the system and the correctness of convergence of complex Langevin to construct a kernel. This allows us to simulate up to $1.5β$ on the real-time Schwinger-Keldysh contour with the 0+1 dimensional anharmonic oscillator using $m=1$, $λ=24$, which was previously unattainable using the complex Langevin equation. △ Less

Submitted 15 December, 2022; originally announced December 2022.

Comments: 8 pages, 10 figures, talk given at the 39th International Symposium on Lattice Field Theory, 8th-13th August, 2022, Bonn, Germany

Towards learning optimized kernels for complex Langevin

Authors: Daniel Alvestad, Rasmus Larsen, Alexander Rothkopf

Abstract: We present a novel strategy aimed at restoring correct convergence in complex Langevin simulations. The central idea is to incorporate system-specific prior knowledge into the simulations, in order to circumvent the NP-hard sign problem. In order to do so, we modify complex Langevin using kernels and propose the use of modern auto-differentiation methods to learn optimal kernel values. The optimiz… ▽ More We present a novel strategy aimed at restoring correct convergence in complex Langevin simulations. The central idea is to incorporate system-specific prior knowledge into the simulations, in order to circumvent the NP-hard sign problem. In order to do so, we modify complex Langevin using kernels and propose the use of modern auto-differentiation methods to learn optimal kernel values. The optimization process is guided by functionals encoding relevant prior information, such as symmetries or Euclidean correlator data. Our approach recovers correct convergence in the non-interacting theory on the Schwinger-Keldysh contour for any real-time extent. For the strongly coupled quantum anharmonic oscillator we achieve correct convergence up to three-times the real-time extent of the previous benchmark study. An appendix sheds light on the fact that for correct convergence not only the absence of boundary terms, but in addition the correct Fokker-Plank spectrum is crucial. △ Less

Submitted 18 April, 2023; v1 submitted 28 November, 2022; originally announced November 2022.

Journal ref: J. High Energ. Phys. 2023, 57 (2023)

Static quark anti-quark interactions at non-zero temperature from lattice QCD

Authors: Gaurang Parkar, Dibyendu Bala, Olaf Kaczmarek, Rasmus Larsen, Swagato Mukherjee, Peter Petreczky, Alexander Rothkopf, Johannes Heinrich Weber

Abstract: We present results on the in-medium interactions of static quark anti-quark pairs using realistic 2+1 HISQ flavor lattice QCD. Focus is put on the extraction of spectral information from Wilson line correlators in Coulomb gauge using four complementary methods. Our results indicate that on HISQ lattices, the position of the dominant spectral peak associated with the real-part of the interquark pot… ▽ More We present results on the in-medium interactions of static quark anti-quark pairs using realistic 2+1 HISQ flavor lattice QCD. Focus is put on the extraction of spectral information from Wilson line correlators in Coulomb gauge using four complementary methods. Our results indicate that on HISQ lattices, the position of the dominant spectral peak associated with the real-part of the interquark potential remains unaffected by temperature. This is in contrast to prior work in quenched QCD and we present follow up comparisons to newly generated quenched ensembles. △ Less

Submitted 30 November, 2022; v1 submitted 23 November, 2022; originally announced November 2022.

Comments: 7 pages, 4 figures, talk given at the XVth Quark confinement and the Hadron spectrum conference, Aug. 1st - 6th, 2022, Stavanger, Norway. Date of presentation Aug. 5th 2022

Kernel controlled real-time Complex Langevin simulation

Authors: Daniel Alvestad, Rasmus Larsen, Alexander Rothkopf

Abstract: This study explores the utility of a kernel in complex Langevin simulations of quantum real-time dynamics on the Schwinger-Keldysh contour. We give several examples where we use a systematic scheme to find kernels that restore correct convergence of complex Langevin. The schemes combine prior information we know about the system and the correctness of convergence of complex Langevin to construct a… ▽ More This study explores the utility of a kernel in complex Langevin simulations of quantum real-time dynamics on the Schwinger-Keldysh contour. We give several examples where we use a systematic scheme to find kernels that restore correct convergence of complex Langevin. The schemes combine prior information we know about the system and the correctness of convergence of complex Langevin to construct a kernel. This allows us to simulate up to $1.5β$ on the real-time Schwinger-Keldysh contour with the $0+1$ dimensional anharmonic oscillator using $m=1,λ=24$, which was previously unattainable using the complex Langevin equation. △ Less

Submitted 19 November, 2022; originally announced November 2022.

Comments: 6 pages, 3 figures, talk given at the XVth Quark confinement and the Hadron spectrum conference, Aug. 1st - 6th, 2022, Stavanger, Norway. Date of presentation Aug. 1st 2022

Inverse problems, real-time dynamics and lattice simulations

Authors: Alexander Rothkopf

Abstract: The determination of real-time dynamics of strongly coupled quantum fields is a central goal of modern nuclear and particle physics, which requires insight into quantum field theory beyond the weak-coupling approximation. While lattice QCD has provided vital insights into the non-perturbative static properties of quarks and gluons it hides their real-time dynamics behind an ill-posed inverse probl… ▽ More The determination of real-time dynamics of strongly coupled quantum fields is a central goal of modern nuclear and particle physics, which requires insight into quantum field theory beyond the weak-coupling approximation. While lattice QCD has provided vital insights into the non-perturbative static properties of quarks and gluons it hides their real-time dynamics behind an ill-posed inverse problem. In this proceeding I will discuss developments in tackling the inverse problem on the lattice and touch upon progress in the direct simualtion of real-time dynamics. △ Less

Submitted 19 November, 2022; originally announced November 2022.

Comments: 11 pages, 4 figures, plenary talk given at the XVth Quark confinement and the Hadron spectrum conference, August 1st - 6th, 2022, Stavanger, Norway

Towards symmetric discretization schemes via weak boundary conditions

Authors: Alexander Rothkopf

Abstract: The Szymanzik improvement program for gauge theories is most commonly implemented using forward finite difference corrections to the Wilson action. Central symmetric schemes naively applied, suffer from a doubling of degrees of freedom, identical to the well known fermion doubling phenomenon. And while adding a complex Wilson term remedies the problem for fermions, it does not easily transfer to r… ▽ More The Szymanzik improvement program for gauge theories is most commonly implemented using forward finite difference corrections to the Wilson action. Central symmetric schemes naively applied, suffer from a doubling of degrees of freedom, identical to the well known fermion doubling phenomenon. And while adding a complex Wilson term remedies the problem for fermions, it does not easily transfer to real-valued gauge fields. In this talk I report on recent progress in formulating symmetric discretization schemes for classical actions of simple one-dimensional problems. They avoid doubling by exploiting the weak imposition of initial/boundary conditions. Inspired by recent work in the field of numerical analysis of partial differential equations, I construct a regularized summation-by-parts finite difference operator using boundary data based on affine coordinates. Application to a classical initial value problems with second order derivatives are presented. △ Less

Submitted 19 November, 2022; originally announced November 2022.

Comments: 9 pages, 5 figures, talk given at the 39th International Symposium on Lattice Field Theory, 8th-13th August, 2022, Bonn, Germany

Bayesian inference of real-time dynamics from lattice QCD

Authors: Alexander Rothkopf

Abstract: The computation of dynamical properties of nuclear matter, ranging from parton distribution functions of nucleons and nuclei to transport properties in the quark-gluon plasma, constitutes a central goal of modern theoretical physics. This real-time physics often defies a perturbative treatment and the most successful strategy so far is to deploy lattice QCD simulations. These numerical computation… ▽ More The computation of dynamical properties of nuclear matter, ranging from parton distribution functions of nucleons and nuclei to transport properties in the quark-gluon plasma, constitutes a central goal of modern theoretical physics. This real-time physics often defies a perturbative treatment and the most successful strategy so far is to deploy lattice QCD simulations. These numerical computations are based on Monte-Carlo sampling and formulated in an artificial Euclidean time. Real-time physics is most conveniently formulated in terms of spectral functions, which are hidden in lattice QCD behind an ill-posed inverse problem. I will discuss the current methods state-of-the art in the extraction of spectral functions from lattice QCD simulations, based on Bayesian inference and emphasize the importance of prior domain knowledge, vital to regularizing the otherwise ill-posed extraction task. With Bayesian inference allowing us to make explicit the uncertainty in both observations and in our prior knowledge, a systematic estimation of the total uncertainties in the extracted spectral functions is nowadays possible. Two implementations of the Bayesian Reconstruction (BR) method for spectral function extraction, one for MAP point estimates and one based on an open access Monte-Carlo sampler are provided.I will briefly touch on the use of machine learning for spectral function reconstruction and discuss some new insight it has brought to the Bayesian community. △ Less

Submitted 29 August, 2022; originally announced August 2022.

Comments: 33 pages, 4 figures, invited review for the Frontiers in Physics collection "Uncertainty Quantification in Nuclear Physics"

Journal ref: Front. Phys. 10:1028995 (2022)

A new variational discretization technique for initial value problems bypassing governing equations

Authors: Alexander Rothkopf, Jan Nordström

Abstract: Motivated by the fact that both the classical and quantum description of nature rest on causality and a variational principle, we develop a novel and highly versatile discretization prescription for classical initial value problems (IVPs). It is based on an optimization (action) functional with doubled degrees of freedom, which is discretized using a single regularized summation-by-parts (SBP) ope… ▽ More Motivated by the fact that both the classical and quantum description of nature rest on causality and a variational principle, we develop a novel and highly versatile discretization prescription for classical initial value problems (IVPs). It is based on an optimization (action) functional with doubled degrees of freedom, which is discretized using a single regularized summation-by-parts (SBP) operator. Formulated as optimization task it allows us to obtain classical trajectories without the need to derive an equation of motion. The novel regularization we develop in this context is inspired by the weak imposition of initial data, often deployed in the modern treatment of IVPs and is implemented using affine coordinates. We demonstrate numerically the stability, accuracy and convergence properties of our approach in systems with classical equations of motion featuring both first and second order derivatives in time. onvergence properties of our approach in systems with classical equations of motion featuring both first and second order derivatives in time. △ Less

Submitted 24 January, 2023; v1 submitted 27 May, 2022; originally announced May 2022.

Comments: 28 pages, 15 figures, accepted in Journal of Computational Physics

The complex potential from 2+1 flavor QCD using HTL inspired approach

Authors: Dibyendu Bala, Olaf Kaczmarek, Rasmus Larsen, Swagato Mukherjee, Gaurang Parkar, Peter Petreczky, Alexander Rothkopf, Johannes Heinrich Weber

Abstract: We have studied finite temperature complex static quark-antiquark potentials for 2+1 flavor QCD using highly improved staggered action with physical strange quark masses and light quark masses corresponding to a pion mass of 161 MeV. We calculated the potential using Wilson line correlators fixed in Coulomb gauge. For the extraction, we have used HTL motivated parametrization of the correlators. W… ▽ More We have studied finite temperature complex static quark-antiquark potentials for 2+1 flavor QCD using highly improved staggered action with physical strange quark masses and light quark masses corresponding to a pion mass of 161 MeV. We calculated the potential using Wilson line correlators fixed in Coulomb gauge. For the extraction, we have used HTL motivated parametrization of the correlators. We found that the real part of the potential is screened above the crossover temperature and it's close to singlet free energies, whereas the imaginary part is increasing with both distance and temperature. △ Less

Submitted 1 December, 2021; originally announced December 2021.

Comments: 9 pages, 4 figures, Contribution to the 38th International Symposium on Lattice Field Theory, 26th-30th July 2021, Zoom/Gather@Massachusetts Institute of Technology

In-medium static quark potential from spectral functions on realistic HISQ ensembles

Authors: Gaurang Parkar, Dibyendu Bala, Olaf Kaczmarek, Rasmus Larsen, Swagato Mukherjee, Peter Petreczky, Alexander Rothkopf, Johannes Heinrich Weber

Abstract: We explore the interactions between a quark anti-quark pair in a thermal medium based on lattice QCD ensembles with $N_f = 2+1$ dynamical HISQ flavors. Our dataset spans the phenomenologically relevant temperature range between T=140MeV-2GeV based on lattice sizes $N_τ=10,12$ and $16$, with an aspect ratio of $N_σ/N_τ=4$. The peak position $Ω$ and the width $Γ$ of the spectral function of Wilson-l… ▽ More We explore the interactions between a quark anti-quark pair in a thermal medium based on lattice QCD ensembles with $N_f = 2+1$ dynamical HISQ flavors. Our dataset spans the phenomenologically relevant temperature range between T=140MeV-2GeV based on lattice sizes $N_τ=10,12$ and $16$, with an aspect ratio of $N_σ/N_τ=4$. The peak position $Ω$ and the width $Γ$ of the spectral function of Wilson-line correlators in Coulomb gauge is computed. We assess the information content in the correlation functions and deploy three complementary strategies to reconstruct spectral information: model fits, Padé approximation and the Bayesian BR method. Limitations of each approach are carefully assessed. △ Less

Submitted 30 November, 2021; originally announced November 2021.

Comments: Contribution to The 38th International Symposium on Lattice Field Theory, LATTICE2021 26th-30th July, 2021 Zoom/Gather@Massachusetts Institute of Technology, 9 pages, 7 figure, LaTeX, uses PoS class

Journal ref: PoS LATTICE2021 (2022) 239

Static quark anti-quark interactions at non-zero temperature from lattice QCD

Authors: Dibyendu Bala, Olaf Kaczmarek, Rasmus Larsen, Swagato Mukherjee, Gaurang Parkar, Peter Petreczky, Alexander Rothkopf, Johannes Heinrich Weber

Abstract: We study the interactions of a static quark antiquark pair at non-zero temperature using realistic 2+1 flavor lattice QCD calculations. The study consists of two parts: the first investigates the properties of Wilson line correlators in Coulomb gauge and compares to predictions of hard-thermal loop perturbation theory. As a second step we extract the spectral functions underlying the correlators u… ▽ More We study the interactions of a static quark antiquark pair at non-zero temperature using realistic 2+1 flavor lattice QCD calculations. The study consists of two parts: the first investigates the properties of Wilson line correlators in Coulomb gauge and compares to predictions of hard-thermal loop perturbation theory. As a second step we extract the spectral functions underlying the correlators using four conceptually different methods: spectral function fits, a HTL inspired fit for the correlation function, Padé rational approximation and the Bayesian BR spectral reconstruction. We find that our high statistics Euclidean lattice data are amenable to different hypotheses for the shapes of the spectral function and we compare the implications of each analysis method for the existence and properties of a well defined ground state spectral peak. △ Less

Submitted 25 March, 2022; v1 submitted 22 October, 2021; originally announced October 2021.

Journal ref: Phys. Rev. D 105, 054513 , Published 24 March 2022

Static Potential At Non-zero Temperatures From Fine Lattices

Authors: D. Hoying, A. Bazavov, D. Bala, G. Parkar, O. Kaczmarek, R. Larsen, Swagato Mukherjee, P. Petreczky, A. Rothkopf, J. H. Weber

Abstract: We report on a preliminary study of static quark anti-quark potential at non-zero temperature in $2+1$ flavor QCD using $96^3\times N_τ$ lattices with lattice spacing $a=0.028$fm, physical strange quark mass and light quark masses corresponding to pion mass of about $300$ MeV. We use $N_τ=32,~24,~20$ and $16$ that correspond to temperature range $T=220-441$ MeV. The in order to obtain the potentia… ▽ More We report on a preliminary study of static quark anti-quark potential at non-zero temperature in $2+1$ flavor QCD using $96^3\times N_τ$ lattices with lattice spacing $a=0.028$fm, physical strange quark mass and light quark masses corresponding to pion mass of about $300$ MeV. We use $N_τ=32,~24,~20$ and $16$ that correspond to temperature range $T=220-441$ MeV. The in order to obtain the potential we calculate the Wilson line correlator in Coulomb gauge with additional HYP smearing to reduce the noise at large quark anti-quark separations. We apply $0$, $5$ and $10$ steps of HYP smearing to ensure that there is no physical effect from over-smearing. At the two highest temperatures we also consider a noise reduction technique that is based on an interpolation in the spatial separation between the static quark and anti-quark. △ Less

Submitted 5 October, 2021; v1 submitted 1 October, 2021; originally announced October 2021.

Comments: Contribution to The 38th International Symposium on Lattice Field Theory, LATTICE2021 26th-30th July, 2021 Zoom/Gather@Massachusetts Institute of Technology, 7 pages, 4 figure, LaTeX, uses PoS class, typos corrected

Journal ref: PoS LATTICE2021 (2022) 178

The QCD Equation of State in Small Systems

Authors: W. A. Horowitz, Alexander Rothkopf

Abstract: We present first results on finite system size corrections to the equation of state, trace anomaly, and speed of sound for two model systems: 1) free, massless scalar theory and 2) quenched QCD with periodic boundary conditions (PBC). We further present work-in-progress results for quenched QCD with Dirichlet boundary conditions. We present first results on finite system size corrections to the equation of state, trace anomaly, and speed of sound for two model systems: 1) free, massless scalar theory and 2) quenched QCD with periodic boundary conditions (PBC). We further present work-in-progress results for quenched QCD with Dirichlet boundary conditions. △ Less

Submitted 1 March, 2022; v1 submitted 3 September, 2021; originally announced September 2021.

Comments: Proceedings for the 50th International Symposium on Multiparticle Dynamics (ISMD2021). 4 pages, 2 figures. Minor revisions from first version

Conserving Lattice Gauge Theory for Finite Systems

Authors: Alexander Rothkopf

Abstract: In this contribution I discuss a recent proposal of a novel action for lattice gauge theory for finite systems, which accommodates non-periodic spatial boundary conditions. Drawing on the summation-by-parts formulation of finite differences and finite volume strategies of computational electrodynamics, an action is constructed that implements the proper integral form of Gauss' law and exhibits an… ▽ More In this contribution I discuss a recent proposal of a novel action for lattice gauge theory for finite systems, which accommodates non-periodic spatial boundary conditions. Drawing on the summation-by-parts formulation of finite differences and finite volume strategies of computational electrodynamics, an action is constructed that implements the proper integral form of Gauss' law and exhibits an inherently symmetric energy momentum tensor, all while realizing automatic ${\cal O}(a)$ improvement. Its central ingredients are illustrated using Abelian gauge theory as example. △ Less

Submitted 31 August, 2021; originally announced August 2021.

Comments: 9 pages, 4 figures, Talk given at the 38th International Symposium on Lattice Field Theory, 26-30 July, 2021, Zoom/Gather@MIT, virtual

Stable solvers for real-time Complex Langevin

Authors: Daniel Alvestad, Rasmus Larsen, Alexander Rothkopf

Abstract: This study explores the potential of modern implicit solvers for stochastic partial differential equations in the simulation of real-time complex Langevin dynamics. Not only do these methods offer asymptotic stability, rendering the issue of runaway solution moot, but they also allow us to simulate at comparatively largeLangevin time steps, leading to lower computational cost. We compare different… ▽ More This study explores the potential of modern implicit solvers for stochastic partial differential equations in the simulation of real-time complex Langevin dynamics. Not only do these methods offer asymptotic stability, rendering the issue of runaway solution moot, but they also allow us to simulate at comparatively largeLangevin time steps, leading to lower computational cost. We compare different ways of regularizing the underlying path integral and estimate the errors introduced due to the finite Langevin time. Based on that insight, we implement benchmark (non-)thermal simulations of the quantum anharmonic oscillator on the canonical Schwinger-Keldysh contour of short real-time extent. △ Less

Submitted 6 May, 2021; originally announced May 2021.

Conserving Lattice Gauge Theory for Finite Systems

Authors: Alexander Rothkopf

Abstract: In this study I develop a novel action for lattice gauge theory for finite systems, which accommodates non-periodic boundary conditions, implements the proper integral form of Gauss' law and exhibits an inherently symmetric energy momentum tensor, all while realizing automatic ${\cal O}(a)$ improvement. Taking the modern summation-by-parts formulation for finite differences as starting point and c… ▽ More In this study I develop a novel action for lattice gauge theory for finite systems, which accommodates non-periodic boundary conditions, implements the proper integral form of Gauss' law and exhibits an inherently symmetric energy momentum tensor, all while realizing automatic ${\cal O}(a)$ improvement. Taking the modern summation-by-parts formulation for finite differences as starting point and combining it with insight from the finite volume strategies of computational electrodynamics I show how the concept of a conserving discretization can be realized for non-Abelian lattice gauge theory. Major steps in the derivation are illustrated using Abelian gauge theory as example. △ Less

Submitted 17 February, 2021; originally announced February 2021.

Comments: 9 pages, 3 figures; computations in the Abelian setting provided as ancillary Mathematica notebooks

Proper static potential in classical lattice gauge theory at finite T

Authors: Alexander Lehmann, Alexander Rothkopf

Abstract: We compute the proper real-time interaction potential between a static quark and antiquark in classical lattice gauge theory at finite temperature. Our central result is the determination of the screened real-part of this potential, and we reconfirm the presence of an imaginary part. The real part is intimately related to the back-reaction of the static sources onto the gauge fields, incorporated… ▽ More We compute the proper real-time interaction potential between a static quark and antiquark in classical lattice gauge theory at finite temperature. Our central result is the determination of the screened real-part of this potential, and we reconfirm the presence of an imaginary part. The real part is intimately related to the back-reaction of the static sources onto the gauge fields, incorporated via Gauss's law. Differences in the treatment of static sources in quantum and classical lattice gauge theory are discussed. △ Less

Submitted 18 December, 2020; originally announced December 2020.

Comments: 24 pages, 10 figures, supplementary material: 2 Mathematica scripts, source code available on Zenodo

Real-Time-Evolution of Heavy Quarks in the Glasma

Authors: Alexander Lehmann, Alexander Rothkopf

Abstract: We introduce a novel real-time formulation of lattice NRQCD designed for simulations in the background of an highly occupied gluon field. By evolving quarks in the background of a dynamically evolving gluon field we computed the time-evolution of heavy-quarkonium spectral functions as well as the static and for finitely heavy quarks generalised potential. We conclude that the back reaction of the… ▽ More We introduce a novel real-time formulation of lattice NRQCD designed for simulations in the background of an highly occupied gluon field. By evolving quarks in the background of a dynamically evolving gluon field we computed the time-evolution of heavy-quarkonium spectral functions as well as the static and for finitely heavy quarks generalised potential. We conclude that the back reaction of the quarks is necessary for any binding process. Here we discuss the methodology, our results and the origin of the absence of a binding process. △ Less

Submitted 5 March, 2020; originally announced March 2020.

Comments: 37th International Symposium on Lattice Field Theory - Lattice2019. Pre-published: 2019 December 05. To be published in: POS, Vol. 363 (2020)

Bryan's Maximum Entropy Method -- diagnosis of a flawed argument and its remedy

Authors: Alexander Rothkopf

Abstract: The Maximum Entropy Method (MEM) is a popular data analysis technique based on Bayesian inference, which has found various applications in the research literature. While the MEM itself is well-grounded in statistics, I argue that its state-of-the-art implementation, suggested originally by Bryan, artificially restricts its solution space. This restriction leads to a systematic error often unaccoun… ▽ More The Maximum Entropy Method (MEM) is a popular data analysis technique based on Bayesian inference, which has found various applications in the research literature. While the MEM itself is well-grounded in statistics, I argue that its state-of-the-art implementation, suggested originally by Bryan, artificially restricts its solution space. This restriction leads to a systematic error often unaccounted for in contemporary MEM studies. The goal of this paper is to carefully revisit Bryan's train of thought, point out its flaw in applying linear algebra arguments to an inherently nonlinear problem, and suggest possible ways to overcome it. △ Less

Submitted 2 November, 2020; v1 submitted 23 February, 2020; originally announced February 2020.

Comments: 16 pages, 3 figures, peer reviewed and published version

Journal ref: Data 2020, 5(3), 85

Heavy Quarkonium in Extreme Conditions

Authors: Alexander Rothkopf

Abstract: In this report we review recent progress achieved in the understanding of heavy quarkonium under extreme conditions from a theory perspective. Its focus lies both on quarkonium properties in thermal equilibrium, as well as recent developments towards a genuine real-time description, valid also out-of-equilibrium. We will give an overview of the theory tools developed and deployed over the last dec… ▽ More In this report we review recent progress achieved in the understanding of heavy quarkonium under extreme conditions from a theory perspective. Its focus lies both on quarkonium properties in thermal equilibrium, as well as recent developments towards a genuine real-time description, valid also out-of-equilibrium. We will give an overview of the theory tools developed and deployed over the last decade, including effective field theories, lattice field theory simulations, modern methods for spectral reconstructions and the the open-quantum systems paradigm. The report will discuss in detail the concept of quarkonium melting, providing the reader with a contemporary perspective. In order to judge where future progress is needed we will also discuss recent results from experiments and phenomenological modeling of quarkonium in relativistic heavy-ion collisions. △ Less

Submitted 4 December, 2019; originally announced December 2019.

Comments: 129 pages, 55 figures, submitted to Physics Reports on July 23rd 2019

Towards the nucleon hadronic tensor from lattice QCD

Authors: Jian Liang, Terrence Draper, Keh-Fei Liu, Alexander Rothkopf, Yi-Bo Yang

Abstract: We present the first calculation of the hadronic tensor on the lattice for the nucleon. The hadronic tensor can be used to extract the structure functions in deep inelastic scatterings and also provide information for the neutrino-nucleon scattering which is crucial to the neutrino-nucleus scattering experiments at low energies. The most challenging part in the calculation is to solve an inverse p… ▽ More We present the first calculation of the hadronic tensor on the lattice for the nucleon. The hadronic tensor can be used to extract the structure functions in deep inelastic scatterings and also provide information for the neutrino-nucleon scattering which is crucial to the neutrino-nucleus scattering experiments at low energies. The most challenging part in the calculation is to solve an inverse problem. We have implemented and tested three algorithms using mock data, showing that the Bayesian Reconstruction method has the best resolution in extracting peak structures while the Backus-Gilbert and Maximum Entropy methods are somewhat more stable for the flat spectral function. Numerical results are presented for both the elastic case (clover fermions on domain wall configuration with $m_π\sim$ 370 MeV and $a\sim$ 0.06 fm) and a case (anisotropic clover lattice with $m_π\sim$ 380 MeV and $a_t\sim$ 0.035 fm) with large momentum transfer. For the former case, the reconstructed Minkowski hadronic tensor gives precisely the vector charge which proves the feasibility of the approach. While for the latter case, the nucleon resonances and possibly shallow inelastic scattering contributions around $ν=1$ GeV are clearly observed but no information is obtained for higher excited states with $ν>2$ GeV. A check of the effective masses of $ρ$ meson with different lattice setups indicates that, in order to reach higher energy transfers, using lattices with smaller lattice spacings is essential. △ Less

Submitted 12 June, 2019; originally announced June 2019.

Comments: 22 pages, 8 figures

Journal ref: Phys. Rev. D 101, 114503 (2020)

Quarkonium Phenomenology from a Generalised Gauss Law

Authors: David Lafferty, Alexander Rothkopf

Abstract: We present an improved analytic parametrisation of the complex in-medium heavy quark potential derived rigorously from the generalised Gauss law. To this end we combine in a self-consistent manner a non-perturbative vacuum potential with a weak-coupling description of the QCD medium. The resulting Gauss-law parametrisation is able to reproduce full lattice QCD data by using only a single temperatu… ▽ More We present an improved analytic parametrisation of the complex in-medium heavy quark potential derived rigorously from the generalised Gauss law. To this end we combine in a self-consistent manner a non-perturbative vacuum potential with a weak-coupling description of the QCD medium. The resulting Gauss-law parametrisation is able to reproduce full lattice QCD data by using only a single temperature dependent parameter, the Debye mass $m_D$. Using this parametrisation we model the in-medium potential at finite baryo-chemical potential, which allows us to estimate the $Ψ^\prime$/$J/Ψ$ ratio in heavy-ion collisions at different beam energies. △ Less

Submitted 6 June, 2019; originally announced June 2019.

Comments: 10 pages, 4 figures. arXiv admin note: text overlap with arXiv:1906.00035

Journal ref: Universe 2019, 5 (2019) 119

Improved Gauss law model and in-medium heavy quarkonium at finite density and velocity

Authors: David Lafferty, Alexander Rothkopf

Abstract: We explore the in-medium properties of heavy-quarkonium states at finite baryo-chemical potential and finite transverse momentum based on a modern complex-valued potential model. Our starting point is a novel, rigorous derivation of the generalized Gauss law for in-medium quarkonium, combining the non-perturbative physics of the vacuum bound state with a weak coupling description of the medium deg… ▽ More We explore the in-medium properties of heavy-quarkonium states at finite baryo-chemical potential and finite transverse momentum based on a modern complex-valued potential model. Our starting point is a novel, rigorous derivation of the generalized Gauss law for in-medium quarkonium, combining the non-perturbative physics of the vacuum bound state with a weak coupling description of the medium degrees of freedom. Its relation to previous models in the literature is discussed. We show that our approach is able to reproduce the complex lattice QCD heavy quark potential even in the non-perturbative regime, using a single temperature dependent parameter, the Debye mass $m_D$. After vetting the Gauss-law potential with state-of-the-art lattice QCD data, we extend it to the regime of finite baryon density and finite velocity, currently inaccessible to first principles simulations. In-medium spectral functions computed from the Gauss-law potential are subsequently used to estimate the $ψ'/J/ψ$ ratio in heavy-ion collisions at different beam energies and transverse momenta. We find qualitative agreement with the predictions from the statistical model of hadronization for the $\sqrt{s_{NN}}$ dependence and a mild dependence on the transverse momentum. △ Less

Submitted 11 March, 2020; v1 submitted 31 May, 2019; originally announced June 2019.

Comments: 26 pages, 11 figures

Journal ref: Phys. Rev. D 101, 056010 (2020)

Spectral Reconstruction with Deep Neural Networks

Authors: Lukas Kades, Jan M. Pawlowski, Alexander Rothkopf, Manuel Scherzer, Julian M. Urban, Sebastian J. Wetzel, Nicolas Wink, Felix P. G. Ziegler

Abstract: We explore artificial neural networks as a tool for the reconstruction of spectral functions from imaginary time Green's functions, a classic ill-conditioned inverse problem. Our ansatz is based on a supervised learning framework in which prior knowledge is encoded in the training data and the inverse transformation manifold is explicitly parametrised through a neural network. We systematically in… ▽ More We explore artificial neural networks as a tool for the reconstruction of spectral functions from imaginary time Green's functions, a classic ill-conditioned inverse problem. Our ansatz is based on a supervised learning framework in which prior knowledge is encoded in the training data and the inverse transformation manifold is explicitly parametrised through a neural network. We systematically investigate this novel reconstruction approach, providing a detailed analysis of its performance on physically motivated mock data, and compare it to established methods of Bayesian inference. The reconstruction accuracy is found to be at least comparable, and potentially superior in particular at larger noise levels. We argue that the use of labelled training data in a supervised setting and the freedom in defining an optimisation objective are inherent advantages of the present approach and may lead to significant improvements over state-of-the-art methods in the future. Potential directions for further research are discussed in detail. △ Less

Submitted 1 February, 2021; v1 submitted 10 May, 2019; originally announced May 2019.

Comments: 20 pages, 16 figures

Journal ref: Phys. Rev. D 102, 096001 (2020)

Spectral properties of light and charm mesons from $N_f=2+1$ anisotropic lattice QCD

Authors: Ryan Quinn, Jonas Glesaaen, Alexander Rothkopf, Jon-Ivar Skullerud

Abstract: We compute temporal correlators and spectral functions for light, open charm and charmonium mesons in the pseudoscalar and vector channel for a range of temperatures below and above the deconfinement transition. The study is carried out using anisotropic lattice QCD with 2+1 dynamical flavours, $a_s=0.123\,$fm and $a_s/a_τ=3.5$. The high-temperature results are benchmarked by comparing them to rec… ▽ More We compute temporal correlators and spectral functions for light, open charm and charmonium mesons in the pseudoscalar and vector channel for a range of temperatures below and above the deconfinement transition. The study is carried out using anisotropic lattice QCD with 2+1 dynamical flavours, $a_s=0.123\,$fm and $a_s/a_τ=3.5$. The high-temperature results are benchmarked by comparing them to reconstructed correlators obtained by direct summation of the zero temperature correlator. We use two Bayesian methods to reconstruct the spectral functions: the maximum entropy method and the more recent BR method. △ Less

Submitted 26 March, 2019; originally announced March 2019.

Comments: 6 pages, 7 figures, based on poster presented at the XIIIth Quark Confinement and the Hadron Spectrum Conference 2018, Maynooth, Ireland

In-medium heavy quarkonium from lattice NRQCD

Authors: Seyong Kim, Peter Petreczky, Alexander Rothkopf

Abstract: We present the final results from a multi-year study of the in-medium spectral properties of heavy quarkonium bound states on the lattice. In this work we combine high statistics $N_f=2+1$ ensembles from the HotQCD collaboration with the effective theory NRQCD and improved Bayesian spectral reconstruction methods. We corroborate with high precision the hierarchical in-medium modification of quarko… ▽ More We present the final results from a multi-year study of the in-medium spectral properties of heavy quarkonium bound states on the lattice. In this work we combine high statistics $N_f=2+1$ ensembles from the HotQCD collaboration with the effective theory NRQCD and improved Bayesian spectral reconstruction methods. We corroborate with high precision the hierarchical in-medium modification of quarkonium states with respect to their vacuum binding energy and provide updated values on melting temperatures. In particular we are able to understand previous disagreements between different Bayesian methods as resulting from underestimated systematic uncertainties. The main quantitative result is a robust determination of the in-medium mass shifts of quarkonium ground states, which we find are negative, consistent with the behavior observed in strongly coupled pNRQCD potential based computations. △ Less

Submitted 6 March, 2019; originally announced March 2019.

Comments: 7 pages, 9 figures, talk given at the XIIIth Quark Confinement and the Hadron Spectrum Conference 2018, Maynooth, Ireland

Reconstructing parton distribution functions from Ioffe time data: from Bayesian methods to Neural Networks

Authors: Joseph Karpie, Kostas Orginos, Alexander Rothkopf, Savvas Zafeiropoulos

Abstract: The computation of the parton distribution functions (PDF) or distribution amplitudes (DA) of hadrons from first principles lattice QCD constitutes a central open problem. In this study, we present and evaluate the efficiency of a selection of methods for inverse problems to reconstruct the full $x$-dependence of PDFs. Our starting point are the so called Ioffe time PDFs, which are accessible from… ▽ More The computation of the parton distribution functions (PDF) or distribution amplitudes (DA) of hadrons from first principles lattice QCD constitutes a central open problem. In this study, we present and evaluate the efficiency of a selection of methods for inverse problems to reconstruct the full $x$-dependence of PDFs. Our starting point are the so called Ioffe time PDFs, which are accessible from Euclidean time calculations in conjunction with a matching procedure. Using realistic mock data tests, we find that the ill-posed incomplete Fourier transform underlying the reconstruction requires careful regularization, for which both the Bayesian approach as well as neural networks are efficient and flexible choices. △ Less

Submitted 16 January, 2019; originally announced January 2019.

Comments: 1+41 pages, 20 figures

Realistic in-medium heavy-quark potential from high statistics lattice QCD simulations

Authors: Peter Petreczky, Alexander Rothkopf, Johannes Weber

Abstract: We present our first results on a direct computation of the complex in-medium heavy quark potential from realistic lattice QCD simulations. Ensembles with $N_τ=12$ from the HotQCD and TUMQCD collaboration offer unprecedented high statistics, those with $N_τ=16$ unprecedented time resolution, making possible a robust extraction of the real part from the spectral functions of Wilson line correlators… ▽ More We present our first results on a direct computation of the complex in-medium heavy quark potential from realistic lattice QCD simulations. Ensembles with $N_τ=12$ from the HotQCD and TUMQCD collaboration offer unprecedented high statistics, those with $N_τ=16$ unprecedented time resolution, making possible a robust extraction of the real part from the spectral functions of Wilson line correlators. To this end we deploy a combination of a Bayesian reconstruction (BR method), as well as a Padé-like approximation. We corroborate findings made on less realistic lattices that ${\rm Re}[V]$ smoothly transitions from a confining to a screened behavior at high temperatures and its values lie close to the color singlet free energies. A finite value of the ${\rm Im}[V]$ is observed in the quark-gluon-plasma phase. △ Less

Submitted 4 October, 2018; originally announced October 2018.

Comments: 4 pages, 4 figures, talk given at the XXVIIth International Conference on Ultra-relativistic Nucleus-Nucleus Collisions (Quark Matter 2018), Venice, Italy

Quarkonium in-medium properties from realistic lattice NRQCD

Authors: Seyong Kim, Peter Petreczky, Alexander Rothkopf

Abstract: We present the final results of our high statistics study on the properties of bottomonium and charmonium at finite temperature. We focus on the temperature range around the crossover transition $150\leq T\leq 410$MeV, relevant for current heavy ion collision experiments. The QCD medium degrees of freedom which consist of dynamical u,d, and s quarks and gluons are captured by realistic state-of-th… ▽ More We present the final results of our high statistics study on the properties of bottomonium and charmonium at finite temperature. We focus on the temperature range around the crossover transition $150\leq T\leq 410$MeV, relevant for current heavy ion collision experiments. The QCD medium degrees of freedom which consist of dynamical u,d, and s quarks and gluons are captured by realistic state-of-the art ($m_π\approx 161$MeV) lattice QCD simulations of the HotQCD collaboration. For the heavy quarks we deploy the non-relativistic effective field theory of QCD, NRQCD. The in-medium properties of quarkonium are deduced from their spectral functions, which are reconstructed using improved and novel Bayesian approaches. Through a systematic analysis we shed light on the origin of the discrepancies in melting temperatures previously reported in the literature, showing that they are owed to underestimated methods uncertainties of the deployed spectral reconstructions. Our simulations corroborate a picture of sequential in-medium modification, ordered according to the vacuum binding energy of the states. As a central quantitative result, our study reveals how the mass of the heavy quarkonium ground state reduces as temperature increases. The observed spectral modifications are interpreted in the light of, and compared to previous studies based on the complex lattice potential for heavy quarkonium. Thus for the first time we provide a robust picture of in-medium heavy quarkonium modification in the quark-gluon plasma consistent among different non-relativistic methods. We also critically discuss the perspectives for improving on these results. △ Less

Submitted 6 November, 2018; v1 submitted 27 August, 2018; originally announced August 2018.

Comments: 84 pages, 45 figures (t.b.p. in JHEP, submission includes ancillary datafiles for in-medium mass shifts and tentative widths)

Lattice study of static quark-antiquark interactions in dense quark matter

Authors: N. Yu. Astrakhantsev, V. G. Bornyakov, V. V. Braguta, E. -M. Ilgenfritz, A. Yu. Kotov, A. A. Nikolaev, A. Rothkopf

Abstract: In this paper we study the interactions among a static quark-antiquark pair in the presence of dense two-color quark matter with lattice simulation. To this end we compute Polyakov line correlation functions and determine the renormalized color averaged, color singlet and color triplet grand potentials. The color singlet grand potential allows us to elucidate the number of quarks induced by a stat… ▽ More In this paper we study the interactions among a static quark-antiquark pair in the presence of dense two-color quark matter with lattice simulation. To this end we compute Polyakov line correlation functions and determine the renormalized color averaged, color singlet and color triplet grand potentials. The color singlet grand potential allows us to elucidate the number of quarks induced by a static quark antiquark source, as well as the internal energy of such a pair in dense quark matter. We furthermore determine the screening length, which in the confinement phase is synonymous with the string breaking distance. The screening length is a decreasing function of baryon density, due to the possibility to break the interquark string via a scalar diquark condensate at high density. We also study the large distance properties of the color singlet grand potential in a dense medium and find that it is well described by a simple Debye screening formula, parameterized by a Debye mass and an effective coupling constant. The latter is of order of unity, i.e. even at large density two-color quark matter is a strongly correlated system. △ Less

Submitted 9 February, 2019; v1 submitted 20 August, 2018; originally announced August 2018.

Comments: 14 figures, 1 table

Journal ref: JHEP05(2019)171

Reconstructing the gluon

Authors: Anton K. Cyrol, Jan M. Pawlowski, Alexander Rothkopf, Nicolas Wink

Abstract: We reconstruct the gluon spectral function in Landau gauge QCD from numerical data for the gluon propagator. The reconstruction relies on two novel ingredients: Firstly we derive analytically the low frequency asymptotics of the spectral function. Secondly we construct a functional basis from a careful consideration of the analytic properties of the gluon propagator in Landau gauge. This allows us… ▽ More We reconstruct the gluon spectral function in Landau gauge QCD from numerical data for the gluon propagator. The reconstruction relies on two novel ingredients: Firstly we derive analytically the low frequency asymptotics of the spectral function. Secondly we construct a functional basis from a careful consideration of the analytic properties of the gluon propagator in Landau gauge. This allows us to reliably capture the non-perturbative regime of the gluon spectrum. We also compare different reconstruction methods and discuss the respective systematic errors. △ Less

Submitted 7 November, 2018; v1 submitted 3 April, 2018; originally announced April 2018.

Comments: 12 pages, 10 figures

Journal ref: SciPost Phys. 5, 065 (2018)

Bayesian study of relativistic open and hidden charm in anisotropic lattice QCD

Authors: Aoife Kelly, Alexander Rothkopf, Jon-Ivar Skullerud

Abstract: We present the first combined study of correlators and spectral properties of charmonium and open-charm-mesons at finite temperature using a fully relativistic lattice QCD approach. The QCD medium is captured by second generation anisotropic 24^3xN_t lattices from the FASTSUM collaboration, including 2+1 flavors of clover discretized quarks with m_pi~380 MeV. Two Bayesian methods are deployed to r… ▽ More We present the first combined study of correlators and spectral properties of charmonium and open-charm-mesons at finite temperature using a fully relativistic lattice QCD approach. The QCD medium is captured by second generation anisotropic 24^3xN_t lattices from the FASTSUM collaboration, including 2+1 flavors of clover discretized quarks with m_pi~380 MeV. Two Bayesian methods are deployed to reconstruct the spectral functions, the recent BR method as well as the Maximum Entropy Method with Fourier basis. We take particular care to disentangle genuine in-medium effects from method artifacts with the help of the reconstructed correlator. Consistent with the direct inspection of correlators, we observe no significant in-medium modification for J/Psi and eta_c around the crossover, while the chi_c states on the other hand show clear changes around the transition. At the highest temperature, T=352 MeV, J/Psi and eta_c also exhibit discernible changes compared to the vacuum. For D mesons around $T_c$, no significant modifications are observed, but we find clear indications that no bound state survives at the highest temperature of T=352 MeV. Above T_c we discover a significant difference between D and D* mesons, the latter being much more strongly affected by the medium. △ Less

Submitted 14 September, 2018; v1 submitted 2 February, 2018; originally announced February 2018.

Comments: 17 pages, many pdf figures. Matches version published in PRD

Journal ref: Phys. Rev. D 97, 114509 (2018)

Transversal and longitudinal gluon spectral functions from twisted mass lattice QCD with $N_f=2+1+1$ flavors

Authors: E. -M. Ilgenfritz, J. M. Pawlowski, A. Rothkopf, A. M. Trunin

Abstract: I report on the first application of a novel, generalized Bayesian reconstruction (BR) method for spectral functions to the characterization of QCD constituents. These spectral functions find applications in off-shell kinetics of the quark-gluon plasma and in calculations of transport coefficients. The new BR method is applied to Euclidean propagator data, obtained in Landau gauge on lattices with… ▽ More I report on the first application of a novel, generalized Bayesian reconstruction (BR) method for spectral functions to the characterization of QCD constituents. These spectral functions find applications in off-shell kinetics of the quark-gluon plasma and in calculations of transport coefficients. The new BR method is applied to Euclidean propagator data, obtained in Landau gauge on lattices with $N_f=2+1+1$ dynamical flavors by the "twisted mass at finite temperature" (tmfT) collaboration. The deployed reconstruction method is designed for spectral functions that can exhibit positivity violation (opposed to that of hadronic bound states). The transversal and longitudinal gluon spectral functions show a robust structure composed of quasiparticle peak and a negative trough. Characteristic differences between the hadronic and the plasma phase and between the two channels become visible. We obtain the temperature dependence of the transversal and longitudinal gluon masses. △ Less

Submitted 12 October, 2017; originally announced October 2017.

Comments: talk given by E.-M. Ilgenfritz at the 11th APCTP-BLTP JINR-PNPI-SPbSU Joint Workshop , Peterhof, 24 -- 28 July 2017, appears in the Proceedings (to be published in PEPAN Letters)

Journal ref: Physics of Particles and Nuclei Letters (2018) Vol. 15, No. 4, pp 362-366

Improved real-time dynamics from imaginary frequency lattice simulations

Authors: Jan M. Pawlowski, Alexander Rothkopf

Abstract: The computation of real-time properties, such as transport coefficients or bound state spectra of strongly interacting quantum fields in thermal equilibrium is a pressing matter. Since the sign problem prevents a direct evaluation of these quantities, lattice data needs to be analytically continued from the Euclidean domain of the simulation to Minkowski time, in general an ill-posed inverse probl… ▽ More The computation of real-time properties, such as transport coefficients or bound state spectra of strongly interacting quantum fields in thermal equilibrium is a pressing matter. Since the sign problem prevents a direct evaluation of these quantities, lattice data needs to be analytically continued from the Euclidean domain of the simulation to Minkowski time, in general an ill-posed inverse problem. Here we report on a novel approach to improve the determination of real-time information in the form of spectral functions by setting up a simulation prescription in imaginary frequencies. By carefully distinguishing between initial conditions and quantum dynamics one obtains access to correlation functions also outside the conventional Matsubara frequencies. In particular the range between $ω_0$ and $ω_1=2πT$, which is most relevant for the inverse problem may be more highly resolved. In combination with the fact that in imaginary frequencies the kernel of the inverse problem is not an exponential but only a rational function we observe significant improvements in the reconstruction of spectral functions, demonstrated in a simple 0+1 dimensional scalar field theory toy model. △ Less

Submitted 7 October, 2017; originally announced October 2017.

Comments: 8 pages, 5 figures, Talk given at the XXXVth International Symposium on Lattice Field Theory, June 18-24, 2017, Granada, Spain

High statistics study of in-medium S- and P-wave quarkonium states in lattice Non-relativistic QCD

Authors: Seyong Kim, Peter Petreczky, Alexander Rothkopf

Abstract: Many measurements of quarkonium suppression at the LHC, e.g. the nuclear modification factor $R_{AA}$ of $J/Ψ$, are well described by a multitude of different models. Thus pinpointing the underlying physics aspects is difficult and guidance based on first principles is needed. Here we present the current status of our ongoing high precision study of in-medium spectral properties of both bottomoniu… ▽ More Many measurements of quarkonium suppression at the LHC, e.g. the nuclear modification factor $R_{AA}$ of $J/Ψ$, are well described by a multitude of different models. Thus pinpointing the underlying physics aspects is difficult and guidance based on first principles is needed. Here we present the current status of our ongoing high precision study of in-medium spectral properties of both bottomonium and charmonium based on NRQCD on the lattice. This effective field theory allows us to capture the physics of quarkonium without modeling assumptions in a thermal QCD medium. In our study a first principles and realistic description of the QCD medium is provided by state-of-the-art lattices of the HotQCD collaboration at almost physical pion mass. Our updated results corroborate a picture of sequential modification of states with respect to their vacuum binding energy. Using a novel low-gain variant of the Bayesian BR method for reconstructing spectral functions we find that remnant features of the Upsilon may survive up to $T\sim400$MeV, while the $χ_b$ signal disappears around $T\sim270$MeV. The $c\bar{c}$ analysis hints at melting of $χ_c$ below $T\sim190$MeV while some $J/Ψ$ remnant feature might survive up to $T\sim245$MeV. An improved understanding of the numerical artifacts in the Bayesian approach and the availability of increased statistics have made possible a first quantitative study of the in-medium ground state masses, which tend to lower values as $T$ increases, consistent with lattice potential based studies. △ Less

Submitted 18 April, 2017; originally announced April 2017.

Comments: 4 pages, 3 figures, oral contribution to the XXVI International Conference on Ultrarelativistic Nucleus-Nucleus Collisions (Quark Matter 2017), Chicago, USA, February 6-11, 2017

Finite temperature gluon spectral functions from $N_f=2+1+1$ lattice QCD

Authors: Ernst-Michael Ilgenfritz, Jan M. Pawlowski, Alexander Rothkopf, Anton Trunin

Abstract: We investigate gluon spectral functions at finite temperature in Landau gauge, based on a subset of lattice QCD ensembles with $N_f=2+1+1$ dynamical twisted mass quarks flavors, generated by the tmfT collaboration. Our study uses a novel Bayesian approach for the extraction of non-positive definite spectral functions, which for each binned spatial momentum takes into account the gluon correlation… ▽ More We investigate gluon spectral functions at finite temperature in Landau gauge, based on a subset of lattice QCD ensembles with $N_f=2+1+1$ dynamical twisted mass quarks flavors, generated by the tmfT collaboration. Our study uses a novel Bayesian approach for the extraction of non-positive definite spectral functions, which for each binned spatial momentum takes into account the gluon correlation functions at all available discrete imaginary frequencies. The spectral functions are extracted at three different lattice spacing, where for each of them, a scan of temperatures around the crossover transition is carried out at fixed scale. We find indications for the existence of a well defined quasi-particle peak. Due to a relatively small number of imaginary frequencies available, we focus on the momentum and temperature dependence of the position of this spectral feature. This dispersion relation reveals different in-medium masses for longitudinal and transversal gluons at high temperatures, qualitatively consistent with weak coupling expectations. △ Less

Submitted 30 January, 2017; originally announced January 2017.

Comments: 23 pages, 30 figures

Heavy-flavor production and medium properties in high-energy nuclear collisions - What next?

Authors: G. Aarts, J. Aichelin, C. Allton, R. Arnaldi, S. A. Bass, C. Bedda, N. Brambilla, E. Bratkovskaya, P. Braun-Munzinger, G. E. Bruno, T. Dahms, S. K. Das, H. Dembinski, M. Djordjevic, E. G. Ferreiro, A. Frawley, P. -B. Gossiaux, R. Granier de Cassagnac, A. Grelli, M. He, W. Horowitz, G. M. Innocenti, M. Jo, O. Kaczmarek, P. G. Kuijer , et al. (24 additional authors not shown)

Abstract: Open and hidden heavy-flavor physics in high-energy nuclear collisions are entering a new and exciting stage towards reaching a clearer understanding of the new experimental results with the possibility to link them directly to the advancement in lattice Quantum Chromo-dynamics (QCD). Recent results from experiments and theoretical developments regarding open and hidden heavy-flavor dynamics have… ▽ More Open and hidden heavy-flavor physics in high-energy nuclear collisions are entering a new and exciting stage towards reaching a clearer understanding of the new experimental results with the possibility to link them directly to the advancement in lattice Quantum Chromo-dynamics (QCD). Recent results from experiments and theoretical developments regarding open and hidden heavy-flavor dynamics have been debated at the Lorentz Workshop "Tomography of the quark-gluon plasma with heavy quarks}, which was held in October 2016 in Leiden, the Netherlands. In this contribution, we summarize identified common understandings and developed strategies for the upcoming five years, which aim at achieving a profound knowledge of the dynamical properties of the quark-gluon plasma. △ Less

Submitted 3 April, 2017; v1 submitted 23 December, 2016; originally announced December 2016.

Comments: 18 pages, 4 figures, Lorentz workshop document, submitted to EPJ A

Quarkonia at $T>0$ and lattice QCD

Authors: Alexander Rothkopf

Abstract: We report here on recent progress in the determination of S-wave and P-wave heavy-quarkonium states at finite temperature. Our results are based on the combination of effective field theories with numerical lattice QCD simulations. These non-perturbative tools allow us to compute the heavy-quarkonium in-medium spectral functions, from which we in turn determine the melting temperatures of individu… ▽ More We report here on recent progress in the determination of S-wave and P-wave heavy-quarkonium states at finite temperature. Our results are based on the combination of effective field theories with numerical lattice QCD simulations. These non-perturbative tools allow us to compute the heavy-quarkonium in-medium spectral functions, from which we in turn determine the melting temperatures of individual states and estimate phenomenologically relevant observables, such as the $ψ^\prime$ to J/$ψ$ ratio in heavy-ion collisions. △ Less

Submitted 20 November, 2016; originally announced November 2016.

Comments: 11 pages, 9 figures, talk given at the XIIth Quark Confinement and the Hadron Spectrum conference, Thessaloniki, Greece

What lattice QCD spectral functions can tell us about heavy quarkonium in the QGP

Authors: Alexander Rothkopf

Abstract: The bound states of a heavy quark and antiquark ($c\bar{c}, b\bar{b}$) are ideal probes to explore the quark-gluon plasma created in relativistic heavy-ion collisions at the RHIC and LHC. Not only have they become experimentally accessible with high precision but also efficient tools, so called effective field theories (EFT) have been developed to treat them theoretically. Here we present recent p… ▽ More The bound states of a heavy quark and antiquark ($c\bar{c}, b\bar{b}$) are ideal probes to explore the quark-gluon plasma created in relativistic heavy-ion collisions at the RHIC and LHC. Not only have they become experimentally accessible with high precision but also efficient tools, so called effective field theories (EFT) have been developed to treat them theoretically. Here we present recent progress in understanding the in-medium behavior of heavy-quarkonium with the help of EFT's combined with non-perturbative and first principles simulations in lattice QCD. In particular we discuss computations of heavy quarkonium spectral functions with the help of Bayesian unfolding methods and the physics we can extract from them. Limitations and the underlying assumptions of the used approaches are pointed out. △ Less

Submitted 20 November, 2016; originally announced November 2016.

Comments: 6 pages, 5 figures, talk given at the XXXVIIIth International Conference on High Energy Physics, Chicago, USA

A gauge invariant Debye mass for the complex heavy-quark potential

Authors: Yannis Burnier, Alexander Rothkopf

Abstract: The concept of a screening mass is a powerful tool to simplify the intricate physics of in-medium test charges surrounded by light charge carriers. While it has been successfully used to describe electromagnetic properties, its definition and computation in QCD is plagued by questions of gauge invariance and the presence of non-perturbative contributions from the magnetic sector. Here we present a… ▽ More The concept of a screening mass is a powerful tool to simplify the intricate physics of in-medium test charges surrounded by light charge carriers. While it has been successfully used to describe electromagnetic properties, its definition and computation in QCD is plagued by questions of gauge invariance and the presence of non-perturbative contributions from the magnetic sector. Here we present a recent alternative definition of a gauge invariant Debye mass parameter following closely the original idea of Debye and Hueckel. Our test charges are a static heavy quark-antiquark pair whose complex potential and its in-medium modification can be extracted using lattice QCD. By combining in a generalized Gauss-Law the non-perturbative aspects of quark binding with a perturbative ansatz for the medium effects, we succeed to describe the lattice values of the potential with a single temperature dependent parameter, in turn identified with a Debye mass. We find that its behavior, as evaluated in a recent quenched lattice QCD study, deviates from that in other approaches, such as hard-thermal-loop perturbation theory or from electric field correlators on the lattice. In particular around the phase transition its values tend to zero significantly faster than at weak-coupling. △ Less

Submitted 20 November, 2016; originally announced November 2016.

Comments: 7 pages, 2 figures, talk given at the XXXIVth International Symposium on Lattice Field Theory, 24-30 July 2016, Southampton, UK

Bayesian inference of non-positive spectral functions in quantum field theory

Authors: Alexander Rothkopf

Abstract: We present the generalization to non positive definite spectral functions of a recently proposed Bayesian deconvolution approach (BR method). The novel prior used here retains many of the beneficial analytic properties of the original method, in particular it allows us to integrate out the hyperparameter $α$ directly. To preserve the underlying axiom of scale invariance, we introduce a second defa… ▽ More We present the generalization to non positive definite spectral functions of a recently proposed Bayesian deconvolution approach (BR method). The novel prior used here retains many of the beneficial analytic properties of the original method, in particular it allows us to integrate out the hyperparameter $α$ directly. To preserve the underlying axiom of scale invariance, we introduce a second default-model related function, whose role is discussed. Our reconstruction prescription is contrasted with existing direct methods, as well as with an approach where shift functions are introduced to compensate for negative spectral features. A mock spectrum analysis inspired by the study of gluon spectral functions in QCD illustrates the capabilities of this new approach. △ Less

Submitted 29 March, 2017; v1 submitted 2 November, 2016; originally announced November 2016.

Comments: 14 pages, 15 figures

Journal ref: Phys. Rev. D 95, 056016 (2017)

Thermal dynamics on the lattice with exponentially improved accuracy

Authors: Jan Pawlowski, Alexander Rothkopf

Abstract: We present a novel simulation prescription for thermal quantum fields on a lattice that operates directly in imaginary frequency space. By distinguishing initial conditions from quantum dynamics it provides access to correlation functions also outside of the conventional Matsubara frequencies $ω_n=2πn T$. In particular it resolves their frequency dependence between $ω=0$ and $ω_1=2πT$, where the t… ▽ More We present a novel simulation prescription for thermal quantum fields on a lattice that operates directly in imaginary frequency space. By distinguishing initial conditions from quantum dynamics it provides access to correlation functions also outside of the conventional Matsubara frequencies $ω_n=2πn T$. In particular it resolves their frequency dependence between $ω=0$ and $ω_1=2πT$, where the thermal physics $ω\sim T$ of e.g.~transport phenomena is dominantly encoded. Real-time spectral functions are related to these correlators via an integral transform with rational kernel, so their unfolding is exponentially improved compared to Euclidean simulations. We demonstrate this improvement within a $0+1$-dimensional scalar field theory and show that spectral features inaccessible in standard Euclidean simulations are quantitatively captured. △ Less

Submitted 28 March, 2017; v1 submitted 29 October, 2016; originally announced October 2016.

Comments: 5 pages, 4 figures

Complex heavy-quark potential and Debye mass in a gluonic medium from lattice QCD

Authors: Yannis Burnier, Alexander Rothkopf

Abstract: We improve and extend our study of the complex in-medium heavy quark potential and its Debye mass $m_D$ in a gluonic medium with a finer scan around the deconfinement transition and newly generated ensembles closer to the thermodynamic limit. On the lattices with larger physical volume, Re[V] shows signs of screening, i.e. a finite $m_D$, only in the deconfined phase, reminiscent of a genuine phas… ▽ More We improve and extend our study of the complex in-medium heavy quark potential and its Debye mass $m_D$ in a gluonic medium with a finer scan around the deconfinement transition and newly generated ensembles closer to the thermodynamic limit. On the lattices with larger physical volume, Re[V] shows signs of screening, i.e. a finite $m_D$, only in the deconfined phase, reminiscent of a genuine phase transition. Consistently Im[V] exhibits nonzero values also only above $T_C$. We compare the behavior of Re[V] with the color singlet free-energies that have been used historically to extract the Debye mass. An effective coupling constant is computed to assess the residual influence of the confining part of the potential at $T>0$. Our previous finding of a gradual screening of Re[V] around $T_C$ on finer lattices is critically reassessed and interpreted to originate from finite volume artifacts that affect the deployed $β=7$, $ξ_b=3.5$ parameter set at $N_s=32$. △ Less

Submitted 18 April, 2017; v1 submitted 14 July, 2016; originally announced July 2016.

Comments: 13 pages, 10 figures

Journal ref: Phys. Rev. D 95, 054511 (2017)

In-medium P-wave quarkonium from the complex lattice QCD potential

Authors: Yannis Burnier, Olaf Kaczmarek, Alexander Rothkopf

Abstract: We extend our lattice QCD potential based study arXiv:1509.07366 (JHEP 1512 (2015) 101) of the in-medium properties of heavy quark bound states to P-wave bottomonium and charmonium. Similar to the behavior found in the S-wave channel their spectra show a characteristic broadening, as well as mass shifts to lower energy with increasing temperature. In contrast to the S-wave states, finite angular m… ▽ More We extend our lattice QCD potential based study arXiv:1509.07366 (JHEP 1512 (2015) 101) of the in-medium properties of heavy quark bound states to P-wave bottomonium and charmonium. Similar to the behavior found in the S-wave channel their spectra show a characteristic broadening, as well as mass shifts to lower energy with increasing temperature. In contrast to the S-wave states, finite angular momentum leads to the survival of spectral peaks even at temperatures, where the continuum threshold reaches below the bound state remnant mass. We elaborate on the ensuing challenges in defining quarkonium dissolution, present estimates of melting temperatures for the spin averaged $χ_b$ and $χ_c$ states and contrast the findings to recent direct lattice NRQCD studies of P-wave quarkonium. As an application to heavy-ion collisions we estimate the contribution of feed down to S-wave quarkonium through the P-wave states after freezeout. △ Less

Submitted 20 June, 2016; originally announced June 2016.

Comments: 26 pages, 9 figures

Report number: BI-TP 2016/06