BZPEERPapers, people, evidence and research paths
Menu

Discover

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

Research tools

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

Network

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

Opportunities

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

Account

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

Source author record

Alexander Rothkopf

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

ResearcherUnclaimed source record
hep-lathep-phnucl-thphysics.comp-phnucl-exphysics.data-anquant-phastro-phhep-exMachine Learning

Catalog footprint

What is connected

44works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

44 published item(s)

preprint2022arXiv

Inverse problems, real-time dynamics and lattice simulations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The effective complex heavy-quark potential in an anisotropic quark-gluon plasma

We introduce a method for reducing anisotropic heavy-quark potentials to isotropic potentials by using an effective screening mass that depends on the quantum numbers $l$ and $m$ of a given state. We demonstrate that, using the resulting 1D effective potential model, one can solve a 1D Schrödinger equation and reproduce the full 3D results for the energies and binding energies of low-lying heavy-quarkonium bound states to relatively high accuracy. This includes the splitting of different p-wave polarizations. The resulting 1D effective model provides a way to include momentum anisotropy effects in open quantum system simulations of heavy-quarkonium dynamics in the quark-gluon plasma.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Conserving Lattice Gauge Theory for Finite Systems

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Spectral Reconstruction with Deep Neural Networks

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Quarkonium production and suppression: Theory

Heavy quarkonium theory has seen significant progress over the past years. Not only has the community steadily improved the understanding of fully equilibrated quarkonium using first principles lattice QCD studies, but it has taken major steps towards developing a genuine dynamical understanding of quark-antiquark pairs immersed in a hot medium by embracing the paradigm of open quantum systems. We review recent developments, placing an emphasis on quarkonium real-time dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Real-Time-Evolution of Heavy Quarks in the Glasma

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Thermalization and dynamical spectral properties in the quark-meson model

We investigate the nonequilibrium evolution of the quark-meson model using two-particle irreducible effective action techniques. Our numerical simulations, which include the full dynamics of the order parameter of chiral symmetry, show how the model thermalizes into different regions of its phase diagram. In particular, by studying quark and meson spectral functions, we shed light on the real-time dynamics approaching the crossover transition, revealing e.g. the emergence of light effective fermionic degrees of freedom in the infrared. At late times in the evolution, the fluctuation-dissipation relation emerges naturally among both meson and quark degrees of freedom, confirming that the simulation successfully reaches thermal equilibrium.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Trace preserving quantum dynamics using a novel reparametrization-neutral summation-by-parts difference operator

We develop a novel numerical scheme for the simulation of dissipative quantum dynamics following from two-body Lindblad master equations. All defining continuum properties of the Lindblad dynamics, hermiticity, positivity and in particular trace conservation of the evolved density matrix are preserved. The central ingredient is a new spatial difference operator, which not only fulfils the summation by parts (SBP) property but also implements a continuum reparametrization property. Using the time evolution of a heavy-quark anti-quark bound state in a hot thermal medium as an explicit example, we show how the reparametrization neutral summation-by-parts (RN-SBP) operator preserves the continuum properties of the theory.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Heavy Quarkonium in Extreme Conditions

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Quantum Brownian motion of a heavy quark pair in the quark-gluon plasma

In this paper we study the real-time evolution of heavy quarkonium in the quark-gluon plasma (QGP) on the basis of the open quantum systems approach. In particular, we shed light on how quantum dissipation affects the dynamics of the relative motion of the quarkonium state over time. To this end we present a novel non-equilibrium master equation for the relative motion of quarkonium in a medium, starting from Lindblad operators derived systematically from quantum field theory. In order to implement the corresponding dynamics, we deploy the well established quantum state diffusion method. In turn we reveal how the full quantum evolution can be cast in the form of a stochastic non-linear Schrödinger equation. This for the first time provides a direct link from quantum chromodynamics (QCD) to phenomenological models based on non-linear Schrödinger equations. Proof of principle simulations in one-dimension show that dissipative effects indeed allow the relative motion of the constituent quarks in a quarkonium at rest to thermalize. Dissipation turns out to be relevant already at early times well within the QGP lifetime in relativistic heavy ion collisions.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Towards the nucleon hadronic tensor from lattice QCD

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A first look at Bottomonium melting via a stochastic potential

We investigate the phenomenon of Bottomonium melting in a thermal quark-gluon plasma using three-dimensional stochastic simulations based on the concept of open-quantum systems. In this non-relativistic framework, introduced in [Phys.Rev. D85 (2012) 105011], which makes close contact to the potentials derived in effective field theory, the $b\bar{b}$ system evolves unitarily under the incessant kicks by the constituents of the surrounding heat bath. In particular thermal fluctuations and the presence of a complex potential in the EFT are naturally related. An intricate interplay between state mixing and thermal excitations emerges as we show how non-thermal initial conditions of Bottomonium states evolve over time. We emphasize that the dynamics of these states gives us access to information beyond what is encoded in the thermal Bottomonium spectral functions. Assumptions underlying our approach and their limitations, as well as the refinements necessary to connect to experimental measurements under more realistic conditions are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

In-medium heavy-quarkonium from lattice QCD spectral functions

We discuss recent progress in the study of in-medium heavy quarkonium using first-principles lattice QCD calculations and effective field theory. In particular our focus lies on real-time information carried by QCD spectral functions and we report on a novel strategy for their extraction from Euclidean-time simulations. Combining an effective description for heavy quarks with finite temperature lattice QCD we compute in-medium spectra providing insight in heavy quarkonium melting in a static thermal medium. On the other hand we use spectral information to extract the complex in-medium potential which provides the basis for a dynamical description of quarkonium real-time evolution via the Schrödinger equation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

In-medium quarkonium properties from a lattice QCD based effective field theory

In order to understand the experimental data on heavy quarkonium production in heavy ion collisions at RHIC and LHC it is necessary (though not sufficient) to pinpoint the properties of heavy $Q\bar{Q}$ bound states in the deconfined quark-gluon plasma, including their dissolution. Here we present recent results on the temperature dependence of bottomonium and charmonium correlators, as well as their spectral functions in a lattice QCD based effective field theory called NRQCD, surveying temperatures close to the crossover transition $140 {\rm MeV} < T< 249 {\rm MeV}$. The spectra are reconstructed based on a novel Bayesian prescription, whose systematic uncertainties are assessed. We present indications for sequential melting of different quarkonium species with respect to their vacuum binding energies and give estimates on the survival of S-wave and P-wave ground states.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Lattice NRQCD study of quarkonium at non-zero temperature

To study the in-medium modification of quarkonium properties, charmonium correlators at $140.4 (β=6.664) \le T \le 221 (β= 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 ($β=6.664, 6.740, 6.800, 6.880, 6.950, 7.030, 7.150$ and $7.280$). 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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Lattice NRQCD study of S- and P-wave bottomonium states in a thermal medium with $N_f=2+1$ light flavors

We investigate the properties of S-and P-wave bottomonium states in the vicinity of the deconfinement transition temperature. The light degrees of freedom are represented by dynamical lattice QCD configurations of the HotQCD collaboration with $N_f = 2 + 1$ flavors. Bottomonium correlators are obtained from bottom quark propagators, computed in nonrelativistic quantum chromodynamics (NRQCD) under the background of these gauge field configurations. The spectral functions for the $^3S_1$ ($Υ$) and $^3P_1$ ($χ_{b1}$) channel are extracted from the Euclidean time correlators using a novel Bayesian approach in the temperature region $140 {\rm MeV} \le T \le 249 {\rm MeV}$ and the results are contrasted to those from the standard maximum entropy method. We find that the new Bayesian approach is far superior to the maximum entropy method. It enables us to study reliably the presence or absence of the lowest state signal in the spectral function of a certain channel, even under the limitations present in the finite temperature setup. We find that $χ_{b1}$ survives up to $T=249 {\rm MeV}$, the highest temperature considered in our study and put stringent constraints on the size of the medium modification of $Υ$ and $χ_{b1}$ states.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Non-Abelian chiral instabilities at high temperature on the lattice

We report on an exploratory lattice study on the phenomenon of chiral instabilities in non-Abelian gauge theories at high temperature. It is based on a recently constructed anomalous Langevin-type effective theory of classical soft gauge fields in the presence of a chiral number density $n_5=n_{\rm R}-n_{\rm L}$. Evaluated in thermal equilibrium using classical lattice techniques it reveals that the fluctuating soft fields indeed exhibit a rapid energy increase at early times and we observe a clear dependence of the diffusion rate of topological charge (sphaleron rate) on the the initial $n_5$, relevant in both early universe baryogenesis and relativistic heavy-ion collisions. The topological charge furthermore shows a drift among distinct vacuum sectors, roughly proportional to the initial $n_5$ and in turn the chiral imbalance is monotonously reduced as required by helicity conservation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Quarkonium at finite temperature: Towards realistic phenomenology from first principles

We present the finite temperature spectra of both bottomonium and charmonium, obtained from a consistent lattice QCD based potential picture. Starting point is the complex in-medium potential extracted on full QCD lattices with dynamical u,d and s quarks, generated by the HotQCD collaboration. Using the generalized Gauss law approach, vetted in a previous study on quenched QCD, we fit ${\rm Re}[V]$ with a single temperature dependent parameter $m_D$, the Debye screening mass, and confirm the up to now tentative values of ${\rm Im}[V]$. The obtained analytic expression for the complex potential allows us to compute quarkonium spectral functions by solving an appropriate Schrödinger equation. These spectra exhibit thermal widths, which are free from the resolution artifacts that plague direct reconstructions from Euclidean correlators using Bayesian methods. In the present adiabatic setting, we find clear evidence for sequential melting and derive melting temperatures for the different bound states. Quarkonium is gradually weakened by both screening (${\rm Re}[V]$) and scattering (${\rm Im}[V]$) effects that in combination lead to a shift of their in-medium spectral features to smaller frequencies, contrary to the mass gain of elementary particles at finite temperature.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Lattice NRQCD study on in-medium bottomonium spectra using a novel Bayesian reconstruction approach

We present recent results on the in-medium modification of S- and P-wave bottomonium states around the deconfinement transition. Our study uses lattice QCD with $N_f=2+1$ light quark flavors to describe the non-perturbative thermal QCD medium between $140$MeV$<T<249$MeV and deploys lattice regularized non-relativistic QCD (NRQCD) effective field theory to capture the physics of heavy quark bound states immersed therein. The spectral functions of the $^3S_1$ $(Υ)$ and $^3P_1$ $(χ_{b1})$ bottomonium states are extracted from Euclidean time Monte Carlo simulations using a novel Bayesian prescription, which provides higher accuracy than the Maximum Entropy Method. Based on a systematic comparison of interacting and free spectral functions we conclude that the ground states of both the S-wave $(Υ)$ and P-wave $(χ_{b1})$ channel survive up to $T=249$MeV. Stringent upper limits on the size of the in-medium modification of bottomonium masses and widths are provided.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

NRQCD based S- and P-wave bottomonium spectra at finite temperature from $48^3 \times 12$ lattices with $N_f=2+1$ light HISQ flavors

We study S-wave and P-wave bottomonium spectral functions at non-zero temperature in 2+1 flavor QCD using the NRQCD formulation for bottom quarks. We use a novel Bayesian approach to reconstruct the spectral functions and find that $χ_{b1}$ survives up to $T=249 {\rm MeV}$. We also study the effect of different temporal discretizations of the NRQCD formalism on the bottomonium correlation functions

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Static quark-antiquark potential in the quark-gluon plasma from lattice QCD

We present a state-of-the-art determination of the complex valued static quark-antiquark potential at phenomenologically relevant temperatures around the deconfinement phase transition. Its values are obtained from non-perturbative lattice QCD simulations using spectral functions extracted via a novel Bayesian inference prescription. We find that the real part, both in a gluonic medium as well as in realistic QCD with light $u$, $d$ and $s$ quarks, lies close to the color singlet free energies in Coulomb gauge and shows Debye screening above the (pseudo) critical temperature $T_c$. The imaginary part is estimated in the gluonic medium, where we find that it is of the same order of magnitude as in hard-thermal loop resummed perturbation theory in the deconfined phase.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Bayesian reconstruction of the in-medium heavy quark potential from lattice QCD and its stability

We report recent results of a non-perturbative determination of the static heavy-quark potential in quenched and dynamical lattice QCD at finite temperature. The real and imaginary part of this complex quantity are extracted from the spectral function of Wilson line correlators in Coulomb gauge. To obtain spectral information from Euclidean time numerical data, our study relies on a novel Bayesian prescription that differs from the Maximum Entropy Method. We perform simulations on quenched $32^3\times N_τ$ $(β=7.0,ξ=3.5)$ lattices with $N_τ=24,...,96$, which cover $839{\rm MeV} \geq T\geq 210 {\rm MeV}$. To investigate the potential in a quark-gluon plasma with light u,d and s quarks we utilize $N_f=2+1$ ASQTAD lattices with $m_l=m_s/20$ by the HotQCD collaboration, giving access to temperatures between $286 {\rm MeV} \geq T\geq 148{\rm MeV}$. The real part of the potential exhibits a clean transition from a linear, confining behavior in the hadronic phase to a Debye screened form above deconfinement. Interestingly its values lie close to the color singlet free energies in Coulomb gauge at all temperatures. We estimate the imaginary part on quenched lattices and find that it is of the same order of magnitude as in hard-thermal loop perturbation theory. From among all the systematic checks carried out in our study, we discuss explicitly the dependence of the result on the default model and the number of datapoints.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The in-medium heavy quark potential from quenched and dynamical lattice QCD

We present our latest results for the the complex valued static heavy-quark potential at finite temperature from lattice QCD. The real and imaginary part of the potential are obtained from the position and width of the lowest lying peak in the spectral function of the Wilson line correlator in Coulomb gauge. Spectral information is extracted from Euclidean time data using a novel Bayesian approach different from the Maximum Entropy Method. In order to extract both the real and imaginary part, we generated anisotropic quenched lattices $32^3\times N_τ$ $(β=7.0,ξ=3.5)$ with $N_τ=24,\ldots,96$, corresponding to $839{\rm MeV} \geq T\geq 210 {\rm MeV}$. For the case of a realistic QCD medium with light u, d and s quarks we use isotropic $48^3\times12$ ASQTAD lattices with $m_l=m_s/20$ provided by the HotQCD collaboration, which span $286 {\rm MeV} \geq T\geq 148{\rm MeV}$. We find a clean transition from a confining to a Debye screened real part and observe that its values lie close to the color singlet free energies in Coulomb gauge. The imaginary part, estimated on quenched lattices, is found to be of the same order of magnitude as in hard-thermal loop (HTL) perturbation theory.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A hard thermal loop benchmark for the extraction of the nonperturbative $Q\bar{Q}$ potential

The extraction of the finite temperature heavy quark potential from lattice QCD relies on a spectral analysis of the Wilson loop. General arguments tell us that the lowest lying spectral peak encodes, through its position and shape, the real and imaginary part of this complex potential. Here we benchmark this extraction strategy using leading order hard-thermal loop (HTL) calculations. I.e. we analytically calculate the Wilson loop and determine the corresponding spectrum. By fitting its lowest lying peak we obtain the real- and imaginary part and confirm that the knowledge of the lowest peak alone is sufficient for obtaining the potential. Access to the full spectrum allows an investigation of spectral features that do not contribute to the potential but can pose a challenge to numerical attempts of an analytic continuation from imaginary time data. Differences in these contributions between the Wilson loop and gauge fixed Wilson line correlators are discussed. To better understand the difficulties in a numerical extraction we deploy the Maximum Entropy method with extended search space to HTL correlators in Euclidean time and observe how well the known spectral function and values for the real and imaginary part are reproduced. Possible venues for improvement of the extraction strategy are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A new Bayesian approach to the reconstruction of spectral functions

We present a novel approach for the reconstruction of spectra from Euclidean correlator data that makes close contact to modern Bayesian concepts. It is based upon an axiomatically justified dimensionless prior distribution, which in the case of constant prior function $m(ω)$ only imprints smoothness on the reconstructed spectrum. In addition we are able to analytically integrate out the only relevant overall hyper-parameter $α$ in the prior, removing the necessity for Gaussian approximations found e.g. in the Maximum Entropy Method. Using a quasi-Newton minimizer and high-precision arithmetic, we are then able to find the unique global extremum of $P[ρ|D]$ in the full $N_ω\gg N_τ$ dimensional search space. The method actually yields gradually improving reconstruction results if the quality of the supplied input data increases, without introducing artificial peak structures, often encountered in the MEM. To support these statements we present mock data analyses for the case of zero width delta peaks and more realistic scenarios, based on the perturbative Euclidean Wilson Loop as well as the Wilson Line correlator in Coulomb gauge.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A novel Bayesian approach to spectral function reconstruction

We present a novel approach to the inference of spectral functions from Euclidean time correlator data that makes close contact with modern Bayesian concepts. Our method differs significantly from the maximum entropy method (MEM). A new set of axioms is postulated for the prior probability, leading to an improved expression, which is devoid of the asymptotically flat directions present in the Shanon-Jaynes entropy. Hyperparameters are integrated out explicitly, liberating us from the Gaussian approximations underlying the evidence approach of the MEM. We present a realistic test of our method in the context of the non-perturbative extraction of the heavy quark potential. Based on hard-thermal-loop correlator mock data, we establish firm requirements in the number of data points and their accuracy for a successful extraction of the potential from lattice QCD. An improved potential estimation from previously investigated quenched lattice QCD correlators is provided.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Benchmarking the Bayesian reconstruction of the non-perturbative heavy $Q\bar{Q}$ potential

The extraction of the finite temperature heavy quark potential from lattice QCD relies on a spectral analysis of the real-time Wilson loop. Through its position and shape, the lowest lying spectral peak encodes the real and imaginary part of this complex potential. We benchmark this extraction strategy using leading order hard-thermal loop (HTL) calculations. I.e. we analytically calculate the Wilson loop and determine the corresponding spectrum. By fitting its lowest lying peak we obtain the real- and imaginary part and confirm that the knowledge of the lowest peak alone is sufficient for obtaining the potential. We deploy a novel Bayesian approach to the reconstruction of spectral functions to HTL correlators in Euclidean time and observe how well the known spectral function and values for the real and imaginary part are reproduced. Finally we apply the method to quenched lattice QCD data and perform an improved estimate of both real and imaginary part of the non-perturbative heavy $Q\bar{Q}$ potential.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

From Complex to Stochastic Potential: Heavy Quarkonia in the Quark-Gluon Plasma

The in-medium physics of heavy quarkonium is an ideal proving ground for our ability to connect knowledge about the fundamental laws of physics to phenomenological predictions. One possible route to take is to attempt a description of heavy quark bound states at finite temperature through a Schroedinger equation with an instantaneous potential. Here we review recent progress in devising a comprehensive approach to define such a potential from first principles QCD and extract its, in general complex, values from non-perturbative lattice QCD simulations. Based on the theory of open quantum systems we will show how to interpret the role of the imaginary part in terms of spatial decoherence by introducing the concept of a stochastic potential. Shortcomings as well as possible paths for improvement are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Improved Maximum Entropy Analysis with an Extended Search Space

The standard implementation of the Maximum Entropy Method (MEM) follows Bryan and deploys a Singular Value Decomposition (SVD) to limit the dimensionality of the underlying solution space apriori. Here we present arguments based on the shape of the SVD basis functions and numerical evidence from a mock data analysis, which show that the correct Bayesian solution is not in general recovered with this approach. As a remedy we propose to extend the search basis systematically, which will eventually recover the full solution space and the correct solution. In order to adequately approach problems where an exponentially damped kernel is used, we provide an open-source implementation, using the C/C++ language that utilizes high precision arithmetic adjustable at run-time. The LBFGS algorithm is included in the code in order to attack problems without the need to resort to a particular search space restriction.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Lattice NRQCD study of in-medium bottomonium states using $N_f = 2+1, 48^3 \times 12$ HotQCD configurations

The behavior of bottomonium state correlators at non-zero temperature, $140.4 (β= 6.664) \le T \le 221 (β= 7.280)$ (MeV), where the transition temperature is $154(9)$ (MeV), is studied, using lattice NRQCD on $48^3 \times 12$ HotQCD HiSQ action configurations with light dynamical $N_f = 2 + 1$ ($m_{u,s}/m_s = 0.05$) staggered quarks. In order to understand finite temperature effects on quarkonium states, zero temperature behavior of bottomonium correlators is compared based on $32^4$ ($β= 6.664, 6.800$ and $6.950$) and $48^3 \times 64$ ($β= 7.280)$ lattices. We find that temperature effects on S-wave bottomonium states are small but P-wave bottomonium states show a noticeable temperature dependence above the transition temperature.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Complex Heavy-Quark Potential at Finite Temperature from Lattice QCD

We calculate for the first time the complex potential between a heavy quark and antiquark at finite temperature across the deconfinement transition in lattice QCD. The real and imaginary part of the potential at each separation distance $r$ is obtained from the spectral function of the thermal Wilson loop. We confirm the existence of an imaginary part above the critical temperature $T_C$, which grows as a function of $r$ and underscores the importance of collisions with the gluonic environment for the melting of heavy quarkonia in the quark-gluon-plasma.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Disentangling the timescales behind the non-perturbative heavy quark potential

The static part of the heavy quark potential has been shown to be closely related to the spectrum of the rectangular Wilson loop. In particular the lowest lying positive frequency peak encodes the late time evolution of the two-body system, characterized by a complex potential. While initial studies assumed a perfect separation of early and late time physics, where a simple Lorentian (Breit-Wigner) shape suffices to describe the spectral peak, we argue that scale decoupling in general is not complete. Thus early time, i.e. non-potential effects, significantly modify the shape of the lowest peak. We derive on general grounds an improved peak distribution that reflects this fact. Application of the improved fit to non-perturbative lattice QCD spectra now yields a potential that is compatible with a transition to a deconfined screening plasma.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Heavy Quarkonium in the Quark Gluon Plasma from Effective Field Theories and Potentials

The measurements of heavy quarkonium suppression at RHIC and LHC urge theory to develop intuitive as well as quantitative methods for the description of $Q\bar{Q}$ melting in the quark-gluon plasma. Here I will present a brief sketch on the effective field theory strategies underlying the definition of the heavy quark static potential and report on two recent advances in the extraction and interpretation of such a potential. On the one side, progress has been made in obtaining its values from lattice QCD, which promises to make possible investigating its real and imaginary part non-perturbatively. On the other side, the existence of an imaginary part emphasizes the dynamical nature of the melting process and invites us to make a direct connection to the framework of open quantum systems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Improved Maximum Entropy Method with an Extended Search Space

We report on an improvement to the implementation of the Maximum Entropy Method (MEM). It amounts to departing from the search space obtained through a singular value decomposition (SVD) of the Kernel. Based on the shape of the SVD basis functions we argue that the MEM spectrum for given $N_τ$ data-points $D(τ)$ and prior information $m(ω)$ does not in general lie in this $N_τ$ dimensional singular subspace. Systematically extending the search basis will eventually recover the full search space and the correct extremum. We illustrate this idea through a mock data analysis inspired by actual lattice spectra, to show where our improvement becomes essential for the success of the MEM. To remedy the shortcomings of Bryan's SVD prescription we propose to use the real Fourier basis, which consists of trigonometric functions. Not only does our approach lead to more stable numerical behavior, as the SVD is not required for the determination of the basis functions, but also the resolution of the MEM becomes independent from the position of the reconstructed peaks.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Stochastic potential and quantum decoherence of heavy quarkonium in the quark-gluon plasma

We propose an open quantum systems approach to the physics of heavy quarkonia in a thermal medium, based on stochastic quantum evolution. This description emphasizes the importance of collisions with the environment and focuses on the concept of spatial decoherence of the heavy quarkonium wave function. It is shown how to determine the parameters of the dynamical evolution, i.e. the real potential and the noise strength, from a comparison with quantities to be obtained from lattice QCD. Furthermore the imaginary part of the lattice QCD heavy quark potential is found to be naturally related to the strength of the noise correlations. We discuss the time evolution of $Q\bar{Q}$ analytically in a limiting scenario for the spatial decoherence and provide a qualitative 1-dimensional numerical simulation of the real-time dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

A close look into the carbon disk at the core of the planetary nebula CPD-568032

We present high spatial resolution observations of the dusty core of the Planetary Nebula with Wolf-Rayet central star CPD-568032. These observations were taken with the mid-infrared interferometer VLTI/MIDI in imaging mode providing a typical 300 mas resolution and in interferometric mode using UT2-UT3 47m baseline providing a typical spatial resolution of 20 mas. The visible HST images exhibit a complex multilobal geometry dominated by faint lobes. The farthest structures are located at 7" from the star. The mid-IR environment of CPD-568032 is dominated by a compact source, barely resolved by a single UT telescope in a 8.7 micron filter. The infrared core is almost fully resolved with the three 40-45m projected baselines ranging from -5 to 51 degree but smooth oscillating fringes at low level have been detected in spectrally dispersed visibilities. This clear signal is interpreted in terms of a ring structure which would define the bright inner rim of the equatorial disk. Geometric models allowed us to derive the main geometrical parameters of the disk. For instance, a reasonably good fit is reached with an achromatic and elliptical truncated Gaussian with a radius of 97+/-11 AU, an inclination of 28+/-7 degree and a PA for the major axis at 345+/-7 degree. Furthermore, we performed some radiative transfer modeling aimed at further constraining the geometry and mass content of the disk, by taking into account the MIDI dispersed visibilities, spectra, and the large aperture SED of the source. These models show that the disk is mostly optically thin in the N band and highly flared.

0 citations0 reviewsNo code linkedOpen access