Source author record

Christian Schmidt

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

50works
22topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

50 published item(s)

preprint2024arXiv

Overcoming the compression limit of the individualsequence (zero order empirical entropy) using the Set Shaping Theory

Given the importance of the claim, we want to start by exposing the following consideration: this claim comes out more than a year after the article "Practical applications of Set Shaping Theory in Huffman coding" which reports the program that carried out an experiment of data compression in which the coding limit NH0(S) of a single sequence was questioned. We waited so long because, before making a claim of this type, we wanted to be sure of the consistency of the result. All this time the program has always been public; anyone could download it, modify it and independently obtain the reported results. In this period there have been many information theory experts who have tested the program and agreed to help us, we thank these people for the time dedicated to us and their precious advice. Given a sequence S of random variables i.i.d. with symbols belonging to an alphabet A; the parameter NH0(S) (the zero-order empirical entropy multiplied by the length of the sequence) is considered the average coding limit of the symbols of the sequence S through a uniquely decipherable and instantaneous code. Our experiment that calls into question this limit is the following: a sequence S is generated in a random and uniform way, the value NH0(S) is calculated, the sequence S is transformed into a new sequence f(S), longer but with the symbols belonging to the same alphabet, finally we code f(S) using Huffman coding. By generating a statistically significant number of sequences we obtain that the average value of the length of the encoded sequence f(S) is less than the average value of NH0(S). In this way, a result is obtained which is incompatible with the meaning given to NH0(S).

preprint2023arXiv

Determination of Lee-Yang edge singularities in QCD by rational approximations

We report updated results on the determination of Lee-Yang edge (LYE) singularities in $N_f = 2+1$ QCD using highly improved staggered quarks (HISQ) with physical masses on $N_τ= 4, 6, 8$ lattices. The singularity structure of QCD in the complex $μ_B$ plane is probed using conserved charges calculated at imaginary $μ_B$. The location of the singularities is determined by studying the (uncancelled) poles of multi-point Padé approximants. We show that close to the Roberge-Weiss (RW) transition, the location of the LYE singularities scales according to the $3$-$d$ $Z(2)$ universality class. By combining the new $N_τ= 6$ data with the $N_τ= 4$ data from our previous analysis we extract a rough estimate for the RW temperature in the continuum limit. We also discuss some preliminary results for the singularities close to the chiral phase transition obtained from simulations on $N_τ= 6, 8$ lattices.

preprint2023arXiv

Fourier coefficients of the net-baryon number density

We calculate Fourier coefficients of the net-baryon number as a function of a purely imaginary chemical potential. The asymptotic behavior of these coefficients is governed by the singularity structure of the QCD partition function and thus encodes information on phase transitions. For the calculation of the Fourier coefficients from lattice data of the Bielefeld-Parma collaboration we use a novel Filon-type quadrature, designed for highly oscillatory integrals. We find sensitivity to chiral scaling in a narrow temperature interval below the Roberge-Weiss transition temperature. Scaling fits yield reasonable values for the position of the Lee-Yang edge singularity in the complex chemical potential plane. Our lattice data has been obtained from simulations with (2+1)-flavors of highly improved staggered quarks (HISQ) at imaginary chemical potential on $N_τ=4, 6$ and $8$ lattices at physical quark masses.

preprint2022arXiv

Detecting critical points from Lee-Yang edge singularities in lattice QCD

A new approach is presented to explore the singularity structure of lattice QCD in the complex chemical potential plane. Our method can be seen as a combination of the Taylor expansion and analytic continuation approaches. Its novelty lies in using rational (Padé) approximants for studying Lee-Yang edge singularities. We present a calculation of the cumulants of the net-baryon number as a function of a purely imaginary baryon number chemical potential, obtained with highly improved staggered quarks at temporal lattice extent of $N_τ=4,6$. We construct various rational function approximations of the lattice data and determine their poles (and roots) in the complex plane. We compare the position of the closest pole to the theoretically expected position of the Lee-Yang edge singularity. At high temperature, we find scaling that is in accordance with the expected power law behavior of the Roberge-Weiss transition while a different behavior is found for $T\lesssim 170$ MeV.

preprint2022arXiv

Model predictive control for retinal laser treatment at 1 kHz

Laser photocoagulation is a technique applied in the treatment of retinal diseases. While this is often done manually or using simple control schemes, we pursue an optimization-based approach, namely Model Predictive Control (MPC), to enforce bounds on the peak temperature and, thus, safety during the medical treatment procedure - despite the spot-dependent absorption of the tissue. To this end, a repetition rate of 1 kHz is desirable rendering the real-time requirements a major challenge. We present a tailored MPC scheme using parametric model reduction, an extended Kalman filter for the parameter and state estimation, and suitably constructed stage costs and verify its applicability both in simulation and experiments with porcine eyes. Moreover, we give some insight on the implementation specifically tailored for fast numerical computations.

preprint2022arXiv

New Way to Resum the Lattice QCD Taylor Series Equation of State at Finite Chemical Potential

Taylor expansion of the thermodynamic potential in powers of the (baryo)chemical potential $μ_B$ is a well-known method to bypass the Sign Problem of Lattice QCD. Due to the difficulty in calculating the higher order Taylor coefficients, various alternative expansion schemes as well as resummation techniques have been suggested to extend the Taylor series to larger values of $μ_B$. Recently, a way to resum the contribution of the first $N$ charge density correlation functions $D_1,\dots,D_N$ to the Taylor series to all orders in $μ_B$ was proposed in Phys. Rev. Lett. 128, 2, 022001 (2022). The resummation takes the form of an exponential factor. Since the correlation functions are calculated stochastically, the exponential factor contains a bias which can be significant for large $N$ and $μ_B$. In this paper, we present a new method to calculate the QCD equation of state based on the well-known cumulant expansion from statistics. By truncating the expansion at a maximum order $M$, we end up with only finite products of the correlation functions which can be evaluated in an unbiased manner. Although our formalism is also applicable for $μ_B\ne0$, here we present it for the simpler case of a finite isospin chemical potential $μ_I$ for which there is no Sign Problem. We present and compare results for the pressure and the isospin density obtained using Taylor expansion, exponential resummation and cumulant expansion, and provide evidence that the absence of bias in the latter actually improves the convergence.

preprint2022arXiv

Practical applications of Set Shaping Theory in Huffman coding

One of the biggest criticisms of the Set Shaping Theory is the lack of a practical application. This is due to the difficulty of its application. In fact, to apply this technique from an experimental point of view we must use a table that defines the correspondences between two sets. However, this approach is not usable in practice, because the table has A^N elements, with A number of symbols and N length of the message to be encoded. Consequently, these tables can be implemented in a program only when A and N have a low value. Unfortunately, in these cases, there are no compression algorithms with such efficiency as to detect the improvement introduced by this method. In this article, we use a function capable of performing the transform without using the correspondence table; this allows us to apply this theory to a wide range of values of A and N. The results obtained confirm the theoretical predictions.

preprint2022arXiv

RF regulation with superconducting cavities and beam operation using a frequency shifted cavity

The free-electron laser FLASH at DESY and the European XFEL are operated with superconducting radio frequency cavities and supply beam to several user experiments. The switching time between experiments is limited to dozens of microseconds. This contribution will show a regulation with a frequency shifted superconducting cavity to manipulate and change the accelerating properties of electron bunches with 250 kHz.The main challenge of the concept presented in this contribution can be summarized in this statement: finding a way to modulate the energy of individual bunches in a single-source multiple-cavity scheme for a potential CW upgrade of the EuXFEL.

preprint2022arXiv

Searching for the QCD critical point along the pseudo-critical/freeze-out line using Padé-resummed Taylor expansions of cumulants of conserved charge fluctuations

Using high-statistics datasets generated in (2+1)-flavor QCD calculations at finite temperature we construct estimators for the radius of convergence from an eighth order series expansion of the pressure as well as the number density. We show that the estimator for pressure and number density will be identical in the asymptotic limit. In the vicinity of the pseudo-critical temperature, $T_{pc}\simeq 156.5$~MeV, we find the estimator of the radius of convergence to be $μ_B/T \gtrsim\ 3$ for strangeness-neutral matter. We also present results for the pole structure of the Padé approximants for the pressure at non-zero values of the baryon chemical potential and show that the pole structure of the [4,4] Padé is consistent with not having a critical point at temperatures larger than $135~$MeV and a baryon chemical potential smaller than $μ_B/T \sim \ 2.5$.

preprint2022arXiv

The Chiral Phase Transition in three-flavor QCD from Lattice QCD

We analyze the pseudo-critical behavior of three-flavor QCD using highly improved staggered quarks (HISQ) on lattices with temporal extent $N_τ=8$ and for quark masses corresponding to a pseudoscalar Goldstone mass in the range $80 ~ {\rm MeV} ~ \lesssim ~ m_π~ \lesssim ~ 140 ~ {\rm MeV}$. Our findings are consistent with the occurrence of a second order chiral phase transition at vanishing values of the quark masses. The chiral phase transition temperature at this finite value of the lattice spacing is determined to be $T_c = 98_{-6}^{+3}~{\rm MeV}$. A comparison with a corresponding analysis performed in (2+1)-flavor QCD suggests that the continuum limit extrapolated chiral phase transition temperature in 3-flavor QCD will turn out to be below $90 ~ {\rm MeV}$.

preprint2022arXiv

Towards the chiral phase transition in the Roberge-Weiss plane

We discuss the interplay between chiral and center sector phase transitions that occur in QCD with an imaginary quark chemical potential $μ=i(2n+1) πT/3$. Based on a finite size scaling analysis in (2+1)-flavor QCD using HISQ fermions with a physical strange quark mass and a range of light quark masses, we show that the endpoint of the line of first-order Roberge-Weiss (RW) transitions between center sectors is second order for light quark masses $m_l\ge m_s/320$, and that it belongs to the $3$-d, $Z(2)$ universality class. The operator for the chiral condensate behaves like an energy-like operator in an effective spin model for the RW phase transition. As a consequence, for any non-zero value of the quark mass, the chiral condensate will have an infinite slope at the RW phase transition temperature, $T_{RW}$. Its fluctuation, the disconnected chiral susceptibility, behaves like the specific heat in $Z(2)$ symmetric models and diverges in the infinite volume limit at the RW phase transition temperature $T_{RW}$ for any non-zero value of the light quark masses. Our analysis suggests the critical temperatures for the RW phase transition and the chiral phase transition coincide in the RW plane. On lattices with temporal extent $N_τ=4$, we find in the chiral limit $T_χ=T_{RW}=195(1)~$MeV.

preprint2021arXiv

Conserved charge fluctuations at vanishing net-baryon density from Lattice QCD

We present here continuum extrapolated results for all 2nd order cumulants using the most recent results obtained by the HotQCD collaboration in (2+1)-flavor QCD. We constrain the applicability of various HRG models by presenting a detailed comparison of our results based on different sets of hadron spectra as well as with virial expansion based model calculations. A comparison with our lattice QCD results for 2nd order cumulants with models that parametrize repulsive interactions among baryons and anti-baryons in a hadron resonance gas through a single excluded volume parameter (EVHRG) is also shown.

preprint2021arXiv

Critical behavior towards the chiral limit at vanishing and non-vanishing chemical potentials

We study the scaling behavior of the (2+1)-flavor QCD crossover region towards the chiral limit with smaller-than-physical light quark mass gauge ensembles, generated using the HISQ fermion discretization. At zero chemical potential, we study the fluctuations of conserved charges and their correlations with the chiral condensate, towards the chiral limit. We analyse the role of universal and regular contributions to the above quantities. We find a preliminary estimate of the leading curvature coefficient of the chiral phase transition line using scaling arguments.

preprint2021arXiv

HotQCD on Multi-GPU Systems

We present $\texttt{SIMULATeQCD}$, HotQCD's software for performing lattice QCD calculations on GPUs. Started in late 2017 and intended as a full replacement of the previous single GPU lattice QCD code used by the HotQCD collaboration, our software has been developed into an extensive framework for lattice QCD calculations distributed on multiple GPUs over many compute nodes. The code is built on C++, CUDA, and MPI and leverages modern C++ language features to provide high-level data structures, objects, and algorithms that allow users to express lattice QCD calculations in an intuitive way without sacrificing performance. Implemented algorithms range from gradient flow, correlator measurements, and mixed precision conjugate gradient solvers all the way to full HISQ gauge field configuration generation using RHMC. After successful deployment in large-scale computing projects, we want to share the result of our efforts with the lattice QCD community by making it publicly available. In these proceedings, we will present some of the key features of our code, demonstrate its ease of use, and show benchmarks of performance critical kernels on state-of-the-art supercomputers.

preprint2020arXiv

Chiral phase structure of three flavor QCD in a background magnetic field

We investigate the chiral phase structure of three flavor QCD in a background $U(1)$ magnetic field using the standard staggered action and the Wilson plaquette gauge action. We perform simulations on lattices with a temporal extent of $N_τ=4$ and four spatial extents of $N_σ= 8,16, 20$ and 24. We choose a smaller-than-physical quark mass in lattice spacing as $am = 0.030$ such that there exists a crossover transition at vanishing magnetic fields, and adopt two values of magnetic field strength in lattice spacing $a \sqrt{ e{B}}\simeq 1.5$ and 2. We find that the transition becomes stronger in the presence of a background magnetic field, and turns into a first order as seen from the volume scaling of the order parameter susceptibility as well as the metastable states in the time history of the chiral condensate. On the other hand, the chiral condensate and transition temperature always increase with $B$ even within the regime of a first order phase transition. This suggests that the discrepancy in the behavior of chiral condensates and transition temperature as a function of $B$ between earlier lattice studies using larger-than-physical pion masses with standard staggered fermions and those using physical pions with improved staggered fermions is mainly due to lattice cutoff effects.

preprint2020arXiv

Deep Deterministic Portfolio Optimization

Can deep reinforcement learning algorithms be exploited as solvers for optimal trading strategies? The aim of this work is to test reinforcement learning algorithms on conceptually simple, but mathematically non-trivial, trading environments. The environments are chosen such that an optimal or close-to-optimal trading strategy is known. We study the deep deterministic policy gradient algorithm and show that such a reinforcement learning agent can successfully recover the essential features of the optimal trading strategies and achieve close-to-optimal rewards.

preprint2020arXiv

Higher order cumulants of net baryon-number distributions at non-zero $μ_B$

Using recent results on higher order cumulants of conserved charge fluctuations from lattice QCD, we construct mean, variance, skewness, kurtosis, hyper-skewness and hyper-kurtosis of net-baryon number distributions for small baryon chemical potentials $μ_B$. For the strangeness neutral case ($μ_S=0$) at fixed ratio of electric charge to baryon number density ($\frac{n_Q}{n_B}=0.4$), which is appropriate for a comparison with heavy ion collisions, we present results for $κ_B σ_B^2$, $S_B σ_B^3/M_B$, $κ^{H}_{B}σ_{B}^4$ and $S^{H}_{B}σ^5_{B}/M_{B}$ on the crossover line for the chiral transition, $T_{pc}(μ_B)$. Continuum extrapolations for this pseudo-critical transition line have recently been reported by HotQCD up to baryon chemical potentials $μ_B\simeq 300$ MeV [arXiv:1812.08235]. These cumulant ratios are of direct relevance for comparisons with corresponding ratios measured by STAR in the BES-I and II runs at beam energies $\sqrt{s_{NN}}\ge 20$ GeV. In particular, we point out that recent high statistics results on skewness and kurtosis of net-baryon number distributions obtained by STAR at $\sqrt{s_{NN}} = 54.4$ GeV put strong constraints on freeze-out parameters and are consistent with predictions from thermal QCD.

preprint2020arXiv

QCD phase transition in the chiral limit

We present a lattice QCD based determination of the chiral phase transition temperature in QCD with two massless (up and down) and one strange quark having its physical mass. We propose and calculate two novel estimators for the chiral transition temperature for several values of the light quark masses, corresponding to Goldstone pion masses in the range of $58~{\rm MeV}\lesssim m_π\lesssim 163~{\rm MeV}$. The chiral phase transition temperature is determined by extrapolating to vanishing pion mass using universal scaling analysis. After thermodynamic, continuum and chiral extrapolations we find the chiral phase transition temperature $T_c^0=132^{+3}_{-6}$ MeV. We also present some preliminary calculations on interplay of effective $U_A(1)$ restoration and chiral phase transition towards chiral limit.

preprint2020arXiv

Simulating gauge theories on Lefschetz thimbles

Lefschetz thimbles have been proposed recently as a possible solution to the complex action problem (sign problem) in Monte Carlo simulations. Here we discuss pure abelian gauge theory with a complex coupling $β$ and apply the concept of Generalized Lefschetz thimbles. We propose to simulate the theory on the union of the tangential manifolds to the thimbles. We construct a local Metropolis-type algorithm, that is constrained to a specific tangential manifold. We also discuss how, starting from this result, successive subleading tangential manifolds can be taken into account via a reweighting approach. We demonstrate the algorithm on $U(1)$ gauge theory in 1+1 dimensions and investigate the residual sign problem.

preprint2019arXiv

Critical end points in (2+1)-flavor QCD with imaginary chemical potential

We present here the results from an ongoing determination of the critical quark mass in simulations of (2+1)-flavor QCD with an imaginary chemical potential. Studies with unimproved actions found the existence of a critical quark mass value at which the crossover transition ends on a second order phase transition and becomes first order for smaller values of the quark mass for the case of both vanishing and imaginary chemical potential. We use the Highly Improved Staggered Quark (HISQ) action and perform calculations in the Roberge-Weiss (RW) plane, where the value of the critical mass is expected to be largest. The lowest quark mass value used in our simulation corresponds to the pion mass $m_π$, down to $40$ MeV. Contrary to calculations performed with unimproved actions we find no evidence for the occurrence of first order transitions at the smallest quark mass values explored so far. Moreover we also show that the chiral observables are sensitive to the RW transition. Our results also indicate that the RW transition and chiral transition could coincide in the chiral limit.

preprint2019arXiv

Meson Screening Masses in (2+1)-Flavor QCD

We present lattice QCD results for mesonic screening masses in the temperature range 140 MeV $\lesssim T \lesssim$ 2500 MeV. Our calculations were carried out using (2+1)-flavors of the Highly Improved Staggered Quark (HISQ) action, with a physical value for the strange quark mass and two values of the light quark mass corresponding to pion masses of 160 MeV and 140 MeV. Continuum-extrapolated results were obtained using calculations with a variety of lattice spacings corresponding to temporal lattice extents $N_τ= 6 - 16$. We discuss the implications of these results for the effective restoration of various symmetries in the high temperature phase of QCD, as well as the approach toward the perturbative limit.

preprint2019arXiv

Uncertainty Modeling and Analysis of the European X-ray Free Electron Laser Cavities Manufacturing Process

This paper reports on comprehensive efforts on uncertainty quantification and global sensitivity analysis for accelerator cavity design. As a case study object the TESLA shaped superconducting cavities, as produced for the European X-ray Free Electron Laser (EXFEL), are selected. The choice for these cavities is explained by the available measurement data that can be leveraged to substantiate the simulation model. Each step of the manufacturing chain is documented together with the involved uncertainties. Several of these steps are mimicked on the simulation side, e.g. by introducing a random eigenvalue problem. The uncertainties are then quantified numerically and in particular the sensitivities give valuable insight into the systems behavior. We also compare these findings to purely statistical studies carried out for the manufactured cavities. More advanced, adaptive, surrogate modeling techniques are adopted, which are crucial to incorporate a large number of uncertain parameters. The main contribution is the detailed comparison and fusion of measurement results for the EXFEL cavities on the one hand and simulation based uncertainty studies on the other hand. After introducing the quantities of physical interest for accelerator cavities and the Maxwell eigenvalue problem, the details on the manufacturing of the EXFEL cavities and measurements are reported. This is followed by uncertainty modeling with quantification studies.

preprint2018arXiv

Dense Limit of the Dawid-Skene Model for Crowdsourcing and Regions of Sub-optimality of Message Passing Algorithms

Crowdsourcing is a strategy to categorize data through the contribution of many individuals. A wide range of theoretical and algorithmic contributions are based on the model of Dawid and Skene [1]. Recently it was shown in [2,3] that, in certain regimes, belief propagation is asymptotically optimal for data generated from the Dawid-Skene model. This paper is motivated by this recent progress. We analyze the dense limit of the Dawid-Skene model. It is shown that it belongs to a larger class of low-rank matrix estimation problems for which it is possible to express the asymptotic, Bayes-optimal, performance in a simple closed form. In the dense limit the mapping to a low-rank matrix estimation problem provides an approximate message passing algorithm that solves the problem algorithmically. We identify the regions where the algorithm efficiently computes the Bayes-optimal estimates. Our analysis refines the results of [2,3] about optimality of message passing algorithms by characterizing regions of parameters where these algorithms do not match the Bayes-optimal performance. We further study numerically the performance of approximate message passing, derived in the dense limit, on sparse instances and carry out experiments on a real world dataset.

preprint2016arXiv

Circular Coloring of Random Graphs: Statistical Physics Investigation

Circular coloring is a constraints satisfaction problem where colors are assigned to nodes in a graph in such a way that every pair of connected nodes has two consecutive colors (the first color being consecutive to the last). We study circular coloring of random graphs using the cavity method. We identify two very interesting properties of this problem. For sufficiently many color and sufficiently low temperature there is a spontaneous breaking of the circular symmetry between colors and a phase transition forwards a ferromagnet-like phase. Our second main result concerns 5-circular coloring of random 3-regular graphs. While this case is found colorable, we conclude that the description via one-step replica symmetry breaking is not sufficient. We observe that simulated annealing is very efficient to find proper colorings for this case. The 5-circular coloring of 3-regular random graphs thus provides a first known example of a problem where the ground state energy is known to be exactly zero yet the space of solutions probably requires a full-step replica symmetry breaking treatment.

preprint2016arXiv

On unitary reconstruction of linear optical networks

Linear optical elements are pivotal instruments in the manipulation of classical and quantum states of light. The vast progress in integrated quantum photonic technology enables the implementation of large numbers of such elements on chip while providing interferometric stability. As a trade-off these structures face the intrinsic challenge of characterizing their optical transformation as individual optical elements are not directly accessible. Thus the unitary transformation needs to be reconstructed from a dataset generated with having access to the input and output ports of the device only. Here we present a novel approach to unitary reconstruction that significantly improves upon existing approaches. We compare its performance to several approaches via numerical simulations for networks up to 14 modes. We show that an adapted version of our approach allows to recover all mode-dependent losses and to obtain highest reconstruction fidelities under such conditions.

preprint2016arXiv

Optical Activation of Germanium Plasmonic Antennas in the Mid Infrared

Impulsive interband excitation with femtosecond near-infrared pulses establishes a plasma response in intrinsic germanium structures fabricated on a silicon substrate. This direct approach activates the plasmonic resonance of the Ge structures and enables their use as optical antennas up to the mid-infrared spectral range. The optical switching lasts for hundreds of picoseconds until charge recombination red-shifts the plasma frequency. The full behavior of the structures is modeled by the electrodynamic response established by an electron-hole plasma in a regular array of antennas.

preprint2015arXiv

Controlled Polar Asymmetry of Few-Cycle and Intense Mid-Infrared Pulses

We demonstrate synthesis of super-octave-spanning and phase-locked transients in the multi-terahertz frequency range with amplitudes exceeding 13 MV/cm. Sub-cycle polar asymmetry of the electric field is adjusted by changing the relative phase between superposed fundamental and second harmonic components. The resultant broken symmetry of the field profile is directly resolved via ultrabroadband electro-optic sampling. Access to such waveforms provides a direct route for control of low-energy degrees of freedom in condensed matter as well as non-perturbative light-matter interactions.

preprint2015arXiv

On U- and V-statistics for discontinuous Ito semimartingales

In this paper we examine the asymptotic theory for U-statistics and V-statistics of discontinuous Ito semimartingales that are observed at high frequency. For different types of kernel functions we show laws of large numbers and associated stable central limit theorems. In most of the cases the limiting process will be conditionally centered Gaussian. The structure of the kernel function determines whether the jump and/or the continuous part of the semimartingale contribute to the limit.

preprint2014arXiv

Equation of state at finite net-baryon density using Taylor coefficients up to sixth order

We employ the lattice QCD data on Taylor expansion coefficients up to sixth order to construct an equation of state at finite net-baryon density. When we take into account how hadron masses depend on lattice spacing and quark mass, the coefficients evaluated using the p4 action are equal to those of hadron resonance gas at low temperature. Thus the parametrised equation of state can be smoothly connected to the hadron resonance gas equation of state. We see that the equation of state using Taylor coefficients up to second order is realistic only at low densities, and that at densities corresponding to s/n_B > 40, the expansion converges by the sixth order term.

preprint2014arXiv

Limit theorems for nondegenerate U-statistics of continuous semimartingales

This paper presents the asymptotic theory for nondegenerate $U$-statistics of high frequency observations of continuous Itô semimartingales. We prove uniform convergence in probability and show a functional stable central limit theorem for the standardized version of the $U$-statistic. The limiting process in the central limit theorem turns out to be conditionally Gaussian with mean zero. Finally, we indicate potential statistical applications of our probabilistic results.

preprint2013arXiv

Chiral restoration and deconfinement in two-color QCD with two flavors of staggered quarks

In preparation of lattice studies of the two-color QCD phase diagram we study chiral restoration and deconfinement at finite temperature with two flavors of staggered quarks using an RHMC algorithm on GPUs. We first study unquenching effects in local Polyakov loop distributions, and the Polyakov loop potential obtained via Legendre transformation, in a fixed-scale approach for heavier quarks. We also present the chiral condensate and the corresponding susceptibility over the lattice coupling across the chiral transition for lighter quarks. Using Ferrenberg-Swendsen reweighting we extract the maxima of the chiral susceptibility in order to determine pseudo-critical couplings on various lattices suitable for chiral extrapolations. These are then used to fix the relation between coupling and temperature in the chiral limit.

preprint2013arXiv

Conserved charge fluctuations with HISQ fermions

We calculate cumulants of fluctuations of net-baryon number, net-electric charge and net-strangeness, in the framework of lattice regularized QCD. We use a highly improved staggered quark (HISQ) action on lattices with temporal extent of N_tau=6,8 and 12 and almost physical quark masses. By means of a Taylor expansion in various chemical potentials and under demanding both strangeness neutrality as well as the correct isospin asymmetry, we evaluate these fluctuations at conditions met in heavy ion collisions. Cumulants of net-electric charge fluctuations can, in principle, also be measured in heavy ion experiments. We therefore propose a method to extract freeze-out parameters, such as the freeze-out temperature and baryon chemical potential, based on a comparison of lattice results and experimental measurements of two different ratios of net-electric charge cumulants. As this method involves only (lattice) QCD results and experimental measurements it is model independent and, for the purpose of extracting freeze-out parameters, does not require any input from hadron resonance gas model calculations.

preprint2013arXiv

Strangeness at high temperatures

We use up to fourth order cumulants of net strangeness fluctuations and their correlations with net baryon number fluctuations to extract information on the strange meson and baryon contribution to the low temperature hadron resonance gas, the dissolution of strange hadronic states in the crossover region of the QCD transition and the quasi-particle nature of strange quark contributions to the high temperature quark-gluon plasma phase.

preprint2013arXiv

The strange degrees of freedom in QCD at high temperature

We discuss recent results on fluctuations of conserved charges, and their approach to the hadron resonance gas at temperatures below the chiral crossover (Tc) as well as to a gas of free quarks at very high temperatures. We will focus on the strange degrees of freedom and verify that they are consistent with those of an uncorrelated gas of hadrons for temperatures of T<160 MeV and that they can be described by a quasi particle model at best for T>2Tc. To extract this information we use all cumulants of net baryon number and net strangeness fluctuations and their correlations up to the fourth order. In particular we propose observables that serve as indicator for the validity of the hadronic degrees of freedom and show that the partial pressures from different strangeness sectors agree separately with HRG model predictions.

preprint2012arXiv

Baryon number and charge fluctuations from lattice QCD

We calculate electric and baryonic charge fluctuations on the lattice. Results have been obtained with the highly improved staggered quark action (HISQ) and almost physical quark masses on lattices with spacial extent of $N_τ=6,8,12$. Higher order cumulants of the net-charge distributions are increasingly dominated by a universal scaling behavior, which is arising due to a critical point of QCD in the chiral limit. Considering cumulants up to the sixth order, we observe that they generically behave as expected from universal scaling laws, which is quite different from the cumulants calculated within the hadron resonance gas model. Taking ratios of these cumulants, we obtain volume independent results that can be compared to the experimental measurements. Such a comparison will unambiguously relate the QCD transition temperature that has been determined on the lattice with the freeze out temperature of heavy ion collision at LHC and RHIC.

preprint2012arXiv

Fluctuations and Correlations of net baryon number, electric charge, and strangeness: A comparison of lattice QCD results with the hadron resonance gas model

We calculate the quadratic fluctuations of net baryon number, electric charge and strangeness as well as correlations among these conserved charges in (2+1)-flavor lattice QCD at zero chemical potential. Results are obtained using calculations with tree level improved gauge and the highly improved staggered quark (HISQ) actions with almost physical light and strange quark masses at three different values of the lattice cut-off. Our choice of parameters corresponds to a value of 160 MeV for the lightest pseudo scalar Goldstone mass and a physical value of the kaon mass. The three diagonal charge susceptibilities and the correlations among conserved charges have been extrapolated to the continuum limit in the temperature interval 150 MeV <T < 250 MeV. We compare our results with the hadron resonance gas (HRG) model calculations and find agreement with HRG model results only for temperatures T<= 150 MeV. We observe significant deviations in the temperature range 160 MeV < T < 170 MeV and qualitative differences in the behavior of the three conserved charge sectors. At T < 160 MeV quadratic net baryon number fluctuations in QCD agree with HRG model calculations while, the net electric charge fluctuations in QCD are about 10% smaller and net strangeness fluctuations are about 20% larger. These findings are relevant to the discussion of freeze-out conditions in relativistic heavy ion collisions.

preprint2012arXiv

QCD bulk thermodynamics and conserved charge fluctuations with HISQ fermions

After briefly reviewing recent progress by the HotQCD collaboration in studying the 2+1 flavor QCD equation of state, we will focus on results on fluctuations of conserved charges by the BNL-Bielefeld and HotQCD collaborations. Higher order cumulants of the net-charge distributions are increasingly dominated by a universal scaling behavior, which arises due to a critical point of QCD in the chiral limit. Considering cumulants up to the 6^{th} order, we observe that they generically behave as expected from universal scaling laws, which is quite different from cumulants calculated within the hadron resonance gas model. Taking ratios of these cumulants, we obtain volume independent results that can be compared to the experimental measurements. We will argue that the freeze-out chemical potentials and the freeze-out temperature, usually obtained by a HRG model fit to the measured hadronic yields, can also be obtained in a model independent way from ab-initio lattice QCD calculations by utilizing observables related to conserved charge fluctuations. Further, we will show that the freeze-out strangeness and electric charge chemical potentials can be fixed by imposing strangeness neutrality and isospin asymmetry constraints in the lattice QCD calculations, in order to accommodate conditions met in heavy ion collisions. All results have been obtained with the highly improved staggered quark action (HISQ) and almost physical quark masses on lattices with temporal extent of N_τ=6,8,10,12.

preprint2011arXiv

On the universal critical behavior in 3-flavor QCD

We analyze the universal critical behavior at the chiral critical point in QCD with three degenerate quark masses. We confirm that this critical point lies in the universality class of the three dimensional Ising model. The symmetry of the Ising model, which is Z(2), is not directly realized in the QCD Hamiltonian. After making an ansatz for the magnetization- and energy-like operators as linear admixtures of the chiral condensate and the gluonic action, we determine several non-universal mixing and normalization constants. These parameters determine an unambiguous mapping of the critical behavior in QCD to that of the 3d-Ising model. We verify its validity by showing that the thus obtained orderparameter scales in accordance with the magnetic equation of state of the 3d-Ising model.

preprint2010arXiv

Net-baryon number fluctuations in (2+1)-flavor QCD

We present a lattice study of net-baryon number fluctuations in (2+1)-flavor QCD. The results are based on a Taylor expansion of the pressure with respect to the baryon chemical potential. We calculate higher moments of the net-baryon number fluctuations and compare with the corresponding resonance gas results. We find that for temperature below 0.9T_c the fluctuations seem to agree with the hadron resonance gas predictions. Close to T_c, higher moments are increasingly more sensitive to the critical behavior of the QCD phase transition. Furthermore, we estimate the radius of convergence of the Taylor series as well as the curvature of the transition line in the temperature chemical potential plane.

preprint2010arXiv

Universal critical behavior and the transition temperature in (2+1)-flavor QCD

We discuss the universal critical behavior in (2+1)-flavor QCD by analyzing lattice data from improved staggered fermions generated by the HotQCD Collaboration. We present recent results from two different lattice discretizations and various lattice spacings ($N_τ=6,8,12$) at fixed physical strange quark mass ($m_s$) but varying light quark mass ($m_l$). We find that the chiral order-parameter, i.e. the chiral condensate, shows the expected universal scaling that is associated with the critical point in the chiral limit already for light quark masses $m_l/m_s \lsim 0.05$. From an analysis of the disconnected chiral susceptibility we estimate a preliminary value of the QCD transition temperature.

preprint2010arXiv

Universal O(N) scaling and the chiral critical line in (2+1)-flavor QCD with small chemical potentials

We show that for small values of the chemical potential the curvature of the phase transition line can be deduced from an analysis of scaling properties of the chiral condensate and its susceptibilities. We make use of a recent analysis of the magnetic equation of state in (2+1)-flavor QCD where a connection between the QCD parameters and the universal scaling fields could be established. The remaining dependence of the reduced temperature on the chemical potential can be fixed by an analysis of a mixed susceptibility, obtained from a derivative with respect to quark mass and chemical potential. We extract this dependence which describes the curvature of the phase transition line, at two values of the cut-off, $aT=1/4$ and $1/8$. We find that cut-off effects are small for the curvature parameter and determine the transition line in the chiral limit to leading order in the light quark chemical potential. We obtain $T_c(μ_B)/T_c(0) = 1 - 0.00656(66) (μ_B/T)^2 +{\cal O}(μ_B^4)$.

preprint2009arXiv

Development of a GEM-TPC prototype

The use of GEM foils for the amplification stage of a TPC instead of a con- ventional MWPC allows one to bypass the necessity of gating, as the backdrift is suppressed thanks to the asymmetric field configuration. This way, a novel continuously running TPC, which represents one option for the PANDA central tracker, can be realized. A medium sized prototype with a diameter of 300 mm and a length of 600 mm will be tested inside the FOPI spectrometer at GSI using a carbon or lithium beam at intermediate energies (E = 1-3AGeV). This detector test under realistic experimental conditions should allow us to verify the spatial resolution for single tracks and the reconstruction capability for displaced vertexes. A series of physics measurement implying pion beams is scheduled with the FOPI spectrometer together with the GEM-TPC as well.

preprint2009arXiv

Hadronic fluctuations and correlations at nonzero chemical potential

We present a lattice study of fluctuations and correlations among the conserved charges baryon number and strangeness in (2+1)-flavor QCD. The lattice calculations are based on a Taylor expansion of the pressure. Results are presented at zero and nonzero density on lattices with four and six time slices, corresponding to a lattice spacing of $a\approx 0.25$ fm and $a\approx 0.16$ fm, respectively. The quark masses are almost physical, i.e. the light quark mass has been chosen to one tenth of the physical strange quark mass while the strange quark mass has been set to its physical value. We find that all analyzed fluctuations and correlations develop a peak at nonzero density, when they are treated in an expansion up to sixth order. Especially below the critical temperature ($T_c$) fluctuations and correlations increase drastically, whereas above $T_c$ they are rather unaffected.

preprint2008arXiv

Calculating Quark Number Susceptibilities with Domain-Wall Fermions

We present results from calculations of different quark number and hadronic susceptibilities on 2+1-flavor dynamical domain wall ensembles. We find that the iso-spin and electric charge susceptibilities are especially well suited to determine the transition temperature, as these quantities show only small statistical errors. Moreover, the transition values of the coupling obtained from iso-spin and electrical charge susceptibilities are in good agreement with the one obtained from the chiral condensate.

preprint2008arXiv

Non-zero density QCD by the Taylor expansion method: The isentropic equation of state, hadronic fluctuations and more

We discuss the Taylor expansion approach to non-zero baryon chemical potential ($μ_B$) and present results on expanion coefficients of the pressure and energy density up to the 6-th order in $μ_B$. Calculations have been performed with (2+1)-flavor of improved staggerd fermions (p4fat3) and almost physical masses on lattices with temporal extent of four and six time slices. We use the expansion coefficients to construct the isentropic equation of state on lines of constant entropy per baryon number and various different hadronic fluctuations. Furthermore, we estimate the radius of convegence of our expansion, which can be seen as a method to analyze the structure of the QCD phase diagram.

preprint2008arXiv

QCD thermodynamics at zero and non-zero density

We present results on the QCD equation of state, obtained with two different improved dynamical staggered fermion actions and almost physical quark masses. Lattice cut-off effects are discussed in detail as results for three different lattice spacings are available now, i.e. results have been obtained on lattices with temporal extent of $N_τ=4,6$ and 8. Furthermore we discuss the Taylor expansion approach to non-zero baryon chemical potential by means of an expansion of the pressure. We use the expansion coefficients to calculate various fluctuations and correlations among hadronic charges. We find that the correlations reproduce the qualitative behavior of the resonance gas model below $T_c$ and start to agree with the free gas predictions for $T\gsim 1.5T_c$.

preprint2007arXiv

The density of states method at non-zero chemical potential

We study the QCD phase diagram by first principle lattice calculations at so far unreached high densities. For this purpose we employ the density of states method. Unimproved staggered fermions, which describe four quark flavors in the continuum are used in this analysis. Though the method is quite expensive, small lattices show an indication for a triple-point connecting three different phases on the phase diagram.

preprint2004arXiv

Chiral perturbation theory for twisted mass QCD

Quantum Chromodynamics on a lattice with Wilson fermions and a chirally twisted mass term for two degenerate quark flavours is considered in the framework of chiral perturbation theory. The pion masses and decay constants are calculated in next-to-leading order including terms linear in the lattice spacing a. We treat both unquenched and partially quenched QCD. We also discuss the phase structure of twisted mass lattice QCD.