-
Dynamic critical behavior of the chiral phase transition from the real-time functional renormalization group
Authors:
Johannes V. Roth,
Yunxin Ye,
Sören Schlichting,
Lorenz von Smekal
Abstract:
In the chiral limit the complicated many-body dynamics around the second-order chiral phase transition of two-flavor QCD can be understood by appealing to universality. We present a novel formulation of the real-time functional renormalization group that describes the stochastic hydrodynamic equations of motion for systems in the same dynamic universality class, which corresponds to Model G in the…
▽ More
In the chiral limit the complicated many-body dynamics around the second-order chiral phase transition of two-flavor QCD can be understood by appealing to universality. We present a novel formulation of the real-time functional renormalization group that describes the stochastic hydrodynamic equations of motion for systems in the same dynamic universality class, which corresponds to Model G in the Halperin-Hohenberg classification. Our approach preserves all relevant symmetries of such systems with reversible mode couplings. We show that the calculations indeed produce the non-trivial value $z=d/2$ for the dynamic critical exponent, where $d$ is the number of spatial dimensions. From the momentum and temperature dependence of the diffusion coefficient of the conserved charge densities, we extract the dimensionless universal scaling function.
△ Less
Submitted 7 March, 2024;
originally announced March 2024.
-
Medium modification of pion production in low energy Au+Au collisions
Authors:
Christian Kummer,
Kai Gallmeister,
Lorenz von Smekal
Abstract:
There is a major mismatch between the charged pion yields in Au+Au collisions at low energies calculated by various transport models and the experimental measured values from the Hades collaboration. In this work, reasonable improvements on the equation of state, in-medium modification of cross sections, and the influence of the nuclear potential for Delta resonances will be investigated in the fr…
▽ More
There is a major mismatch between the charged pion yields in Au+Au collisions at low energies calculated by various transport models and the experimental measured values from the Hades collaboration. In this work, reasonable improvements on the equation of state, in-medium modification of cross sections, and the influence of the nuclear potential for Delta resonances will be investigated in the framework of the GiBUU transport model. As a result, we demonstrate that theoretical calculations can indeed describe the charged pion yields measured by Hades for Au+Au collisions rather well, but that a mismatch then remains between calculations and data for the yields of neutral pions extracted from dileptons within the same experiment.
△ Less
Submitted 10 January, 2024; v1 submitted 16 September, 2023;
originally announced September 2023.
-
Medium induced mixing, spatial modulations and critical modes in QCD
Authors:
Maximilian Haensch,
Fabian Rennecke,
Lorenz von Smekal
Abstract:
The mixing between the chiral condensate and the density in hot and dense QCD matter is familiar. We show that the mixing relevant for the ground state is considerably more extensive, and in particular also involves gluonic degrees of freedom. As a result, the Hessian of the QCD effective action is non-Hermitian, but retains a symmetry under combined charge- and complex conjugation. This can lead…
▽ More
The mixing between the chiral condensate and the density in hot and dense QCD matter is familiar. We show that the mixing relevant for the ground state is considerably more extensive, and in particular also involves gluonic degrees of freedom. As a result, the Hessian of the QCD effective action is non-Hermitian, but retains a symmetry under combined charge- and complex conjugation. This can lead to complex-conjugate pairs of eigenvalues of this Hessian, signaling regimes with spatially modulated correlations. Furthermore, based on the analytic structure of the quark determinant at a chiral critical point, we demonstrate that the corresponding massless critical mode is composed of the chiral condensate, the density and the Polyakov loops. Due to an avoided crossing, the critical mode turns out to be disconnected from the chiral condensate in vacuum. We present general arguments for all these features and illustrate them through explicit model calculations.
△ Less
Submitted 5 August, 2024; v1 submitted 30 August, 2023;
originally announced August 2023.
-
QCD Anderson transition with overlap valence quarks on a twisted-mass sea
Authors:
Robin Kehr,
Dominik Smith,
Lorenz von Smekal
Abstract:
In this work we probe the QCD Anderson transition by studying spectral distributions of the massless overlap operator on gauge configurations created by the twisted mass at finite temperature collaboration (tmfT) with 2+1+1 flavors of dynamical quarks and the Iwasaki gauge action. We assess finite-size and discretization effects by considering two different lattice spacings and several physical vo…
▽ More
In this work we probe the QCD Anderson transition by studying spectral distributions of the massless overlap operator on gauge configurations created by the twisted mass at finite temperature collaboration (tmfT) with 2+1+1 flavors of dynamical quarks and the Iwasaki gauge action. We assess finite-size and discretization effects by considering two different lattice spacings and several physical volumes, and mimic the approach to the continuum limit through stereographic projection. Fitting the inflection points of the participation ratios of the overlap Dirac eigenmodes, we obtain estimates of the temperature dependence of the mobility edge, below which quark modes are localized. We observe that it is well-described by a quadratic polynomial and systematically vanishes at temperatures below the pseudocritical one of the chiral transition. In fact, our best estimates within errors overlap with that of the chiral phase transition temperature of QCD in the chiral limit.
△ Less
Submitted 29 April, 2024; v1 submitted 26 April, 2023;
originally announced April 2023.
-
Critical dynamics in a real-time formulation of the functional renormalization group
Authors:
Johannes V. Roth,
Lorenz von Smekal
Abstract:
We present first calculations of critical spectral functions of the relaxational Models A, B, and C in the Halperin-Hohenberg classification using a real-time formulation of the functional renormalization group (FRG). We revisit the prediction by Son and Stephanov that the linear coupling of a conserved density to the non-conserved order parameter of Model A gives rise to critical Model-B dynamics…
▽ More
We present first calculations of critical spectral functions of the relaxational Models A, B, and C in the Halperin-Hohenberg classification using a real-time formulation of the functional renormalization group (FRG). We revisit the prediction by Son and Stephanov that the linear coupling of a conserved density to the non-conserved order parameter of Model A gives rise to critical Model-B dynamics. We formulate both 1-loop and 2-loop self-consistent expansion schemes in the 1PI vertex functions as truncations of the effective average action suitable for real-time applications, and analyze in detail how the different critical dynamics are properly incorporated in the framework of the FRG on the closed-time path. We present results for the corresponding critical spectral functions, extract the dynamic critical exponents for Models A, B, and C, in two and three spatial dimensions, respectively, and compare the resulting values with recent results from the literature.
△ Less
Submitted 21 March, 2023;
originally announced March 2023.
-
Dilepton Decay of Low-mass $ρ$ Mesons
Authors:
K. Gallmeister,
U. Mosel,
L. von Smekal
Abstract:
The HADES collaboration has extracted dilepton mass spectra for $π^-$ induced reactions on the proton from a comparison of data taken on C and CH$_2$ targets. The spectra were interpreted in terms of different versions of vector meson dominance. Here we present results obtained from the theory and generator GiBUU. We first check the subtraction procedure used and then discuss the obtained mass spe…
▽ More
The HADES collaboration has extracted dilepton mass spectra for $π^-$ induced reactions on the proton from a comparison of data taken on C and CH$_2$ targets. The spectra were interpreted in terms of different versions of vector meson dominance. Here we present results obtained from the theory and generator GiBUU. We first check the subtraction procedure used and then discuss the obtained mass spectra for the proton target. We point out that any conclusions on the version of VMD requires the knowledge of the $ρ$ spectral function in the interesting mass region.
△ Less
Submitted 10 October, 2022;
originally announced October 2022.
-
Quarks and Triality in a Finite Volume
Authors:
Milad Ghanbarpour,
Lorenz von Smekal
Abstract:
In order to understand the puzzle of the free energy of an individual quark in QCD, we explicitly construct ensembles with quark numbers $N_V\neq 0\!\mod 3$, corresponding to non-zero triality in a finite subvolume $V$ on the lattice. We first illustrate the basic idea in an effective Polyakov-loop theory for the heavy-dense limit of QCD, and then extend the construction to full Lattice QCD, where…
▽ More
In order to understand the puzzle of the free energy of an individual quark in QCD, we explicitly construct ensembles with quark numbers $N_V\neq 0\!\mod 3$, corresponding to non-zero triality in a finite subvolume $V$ on the lattice. We first illustrate the basic idea in an effective Polyakov-loop theory for the heavy-dense limit of QCD, and then extend the construction to full Lattice QCD, where the electric center flux through the surface of $V$ has to be fixed at all times to account for Gauss's law. This requires introducing discrete Fourier transforms over closed center-vortex sheets around the spatial volume $V$ between all subsequent time slices, and generalizes the construction of 't Hooft's electric fluxes in the purge gauge theory. We derive this same result from a dualization of the Wilson fermion action, and from the transfer matrix formulation with a local $\mathbb Z_3$-Gauss law to restrict the dynamics to sectors with the required center charge in $V$.
△ Less
Submitted 23 June, 2022;
originally announced June 2022.
-
Real-time methods for spectral functions
Authors:
Johannes V. Roth,
Dominik Schweitzer,
Leon J. Sieke,
Lorenz von Smekal
Abstract:
In this paper we develop and compare different real-time methods to calculate spectral functions. These are classical-statistical simulations, the Gaussian state approximation (GSA), and the functional renormalization group (FRG) formulated on the Keldysh closed-time path. Our test-bed system is the quartic anharmonic oscillator, a single self-interacting bosonic degree of freedom, coupled to an e…
▽ More
In this paper we develop and compare different real-time methods to calculate spectral functions. These are classical-statistical simulations, the Gaussian state approximation (GSA), and the functional renormalization group (FRG) formulated on the Keldysh closed-time path. Our test-bed system is the quartic anharmonic oscillator, a single self-interacting bosonic degree of freedom, coupled to an external heat bath providing dissipation analogous to the Caldeira-Leggett model. As our benchmark we use the spectral function from exact diagonalization with constant Ohmic damping. To extend the GSA for the open system, we solve the corresponding Heisenberg-Langevin equations in the Gaussian approximation. For the real-time FRG, we introduce a novel general prescription to construct causal regulators based on introducing scale-dependent fictitious heat baths. Our results explicitly demonstrate how the discrete transition lines of the quantum system gradually build up the broad continuous structures in the classical spectral function as temperature increases. At sufficiently high temperatures, classical, GSA and exact-diagonalization results all coincide. The real-time FRG is able to reproduce the effective thermal mass, but overestimates broadening and only qualitatively describes higher excitations, at the present order of our combined vertex and loop expansion. As temperature is lowered, the GSA follows the ensemble average of the exact solution better than the classical spectral function. In the low-temperature strong-coupling regime, the qualitative features of the exact result are best captured by our real-time FRG calculation, with quantitative improvements to be expected at higher truncation orders.
△ Less
Submitted 31 May, 2022; v1 submitted 23 December, 2021;
originally announced December 2021.
-
Critical dynamics of relativistic diffusion
Authors:
Dominik Schweitzer,
Sören Schlichting,
Lorenz von Smekal
Abstract:
We study the dynamics of self-interacting scalar fields with $Z_2$ symmetry governed by a relativistic Israel-Stuart type diffusion equation in the vicinity of a critical point. We calculate spectral functions of the order parameter in mean-field approximation as well as using first-principles classical-statistical lattice simulations in real-time. We observe that the spectral functions are well-d…
▽ More
We study the dynamics of self-interacting scalar fields with $Z_2$ symmetry governed by a relativistic Israel-Stuart type diffusion equation in the vicinity of a critical point. We calculate spectral functions of the order parameter in mean-field approximation as well as using first-principles classical-statistical lattice simulations in real-time. We observe that the spectral functions are well-described by single Breit-Wigner shapes. Away from criticality, the dispersion matches the expectations from the mean-field approach. At the critical point, the spectral functions largely keep their Breit-Wigner shape, albeit with non-trivial power-law dispersion relations. We extract the characteristic time-scales as well as the dynamic critical exponent $z$, verifying the existence of a dynamic scaling regime. In addition, we derive the universal scaling functions implied by the Breit-Wigner shape with critical power-law dispersion and show that they match the data. Considering equations of motion for a system coupled to a heat bath as well as an isolated system, we perform this study for two different dynamic universality classes, both in two and three spatial dimensions.
△ Less
Submitted 4 September, 2022; v1 submitted 4 October, 2021;
originally announced October 2021.
-
Effects of chiral symmetry restoration on meson and dilepton production in relativistic heavy-ion collisions
Authors:
A. B. Larionov,
L. von Smekal
Abstract:
We include effects of chiral symmetry and its restoration in the kinetic equations for baryon propagation and explore the consequences for $η$, $π^0$, $ρ$ and dilepton production in heavy-ion collisions at 1-2A GeV. Numerical calculations are performed using the GiBUU microscopic transport model supplemented by the parity-doublet model for the mean fields of the nucleon and the $N^*(1535)$ resonan…
▽ More
We include effects of chiral symmetry and its restoration in the kinetic equations for baryon propagation and explore the consequences for $η$, $π^0$, $ρ$ and dilepton production in heavy-ion collisions at 1-2A GeV. Numerical calculations are performed using the GiBUU microscopic transport model supplemented by the parity-doublet model for the mean fields of the nucleon and the $N^*(1535)$ resonance. In this chiral model, a strong dropping of the Dirac mass of the $N^*(1535)$ in the high-density stage of a collision leads to a considerable enhancement in the production of this resonance as compared to the standard (non-linear) Walecka model. As the system expands, the Dirac masses of these abundant soft $N^*(1535)$ resonances gradually increase and ultimately cross the $N η$ decay threshold. As a result, an enhanced low-energy $η$ production is observed in the calculations with chiral mean fields. Comparing with TAPS data on $η$ and $π^0$ production we find that the chiral model improves the agreement for the $m_t$-spectra of $η$'s at small $m_t$ in heavy colliding systems. A similar enhancement is also observed in the soft $ρ$ production due to chiral symmetry and its partial restoration. The resulting dilepton yields at low and intermediate invariant masses are slightly enhanced due to these chiral effects which further improve the agreement between GiBUU transport simulations and HADES data for C+C at 1A GeV.
△ Less
Submitted 24 March, 2022; v1 submitted 8 September, 2021;
originally announced September 2021.
-
Self-consistent $O(4)$ model spectral functions from analytically continued FRG flows
Authors:
Christopher Jung,
Jan-Hendrik Otto,
Ralf-Arno Tripolt,
Lorenz von Smekal
Abstract:
In this paper we explore practicable ways for self-consistent calculations of spectral functions from analytically continued functional renormalization group (aFRG) flow equations. As a particularly straightforward one we propose to include parametrizations of self-energies based on explicit analytic one-loop expressions. To exemplify this scheme we calculate the spectral functions of pion and sig…
▽ More
In this paper we explore practicable ways for self-consistent calculations of spectral functions from analytically continued functional renormalization group (aFRG) flow equations. As a particularly straightforward one we propose to include parametrizations of self-energies based on explicit analytic one-loop expressions. To exemplify this scheme we calculate the spectral functions of pion and sigma meson of the $O(4)$ model at vanishing temperature in the broken phase. Comparing the results with those from previous aFRG calculations, we explicitly demonstrate how self-consistency at all momenta fixes the tight relation between particle masses and decay thresholds. In addition, the two-point functions from our new semi-analytic FRG scheme have the desired domain of holomorphy built in and can readily be studied in the entire cut-complex frequency plane, on physical as well as other Riemann sheets. This is very illustrative and allows, for example, to trace the flow of the resonance pole of the sigma meson across an unphysical sheet. In order to assess the limitations due to the underlying one-loop structure, we also introduce a fully self-consistent numerical scheme based on spectral representations with scale-dependent spectral functions. The most notable improvement of this numerically involved calculation is that it describes the three-particle resonance decay of an off-shell pion, $π^* \to σπ\to3π$. Apart from this further conceptual improvement, overall agreement with the results from the considerably simpler semi-analytic one-loop scheme is very encouraging, however. The latter can therefore provide a sound and practicable basis for self-consistent calculations of spectral functions in more realistic effective theories for warm and dense matter.
△ Less
Submitted 20 October, 2021; v1 submitted 22 July, 2021;
originally announced July 2021.
-
Vector and Axial-Vector Mesons in Nuclear Matter
Authors:
Ralf-Arno Tripolt,
Christopher Jung,
Lorenz von Smekal,
Jochen Wambach
Abstract:
As a first step towards a realistic phenomenological description of vector and axial-vector mesons in nuclear matter, we calculate the spectral functions of the $ρ$ and the $a_1$ meson in a chiral baryon-meson model as a low-energy effective realization of QCD, taking into account the effects of fluctuations from scalar mesons, nucleons, and vector mesons within the Functional Renormalization Grou…
▽ More
As a first step towards a realistic phenomenological description of vector and axial-vector mesons in nuclear matter, we calculate the spectral functions of the $ρ$ and the $a_1$ meson in a chiral baryon-meson model as a low-energy effective realization of QCD, taking into account the effects of fluctuations from scalar mesons, nucleons, and vector mesons within the Functional Renormalization Group (FRG) approach. The phase diagram of the effective hadronic theory exhibits a nuclear liquid-gas phase transition as well as a chiral phase transition at a higher baryon-chemical potential. The in-medium $ρ$ and $a_1$ spectral functions are calculated by using the previously introduced analytically-continued FRG (aFRG) method. Our results show strong modifications of the spectral functions in particular near the critical endpoints of both phase transitions which may well be of relevance for electromagnetic rates in heavy-ion collisions or neutrino emissivities in neutron-star merger events.
△ Less
Submitted 20 August, 2021; v1 submitted 3 May, 2021;
originally announced May 2021.
-
Static magnetic susceptibility in finite-density SU(2) lattice gauge theory
Authors:
P. V. Buividovich,
D. Smith,
L. von Smekal
Abstract:
We study static magnetic susceptibility $χ(T, μ)$ in $SU(2)$ lattice gauge theory with $N_f = 2$ light flavours of dynamical fermions at finite chemical potential $μ$. Using linear response theory we find that $SU(2)$ gauge theory exhibits paramagnetic behavior in both the high-temperature deconfined regime and the low-temperature confining regime. Paramagnetic response becomes stronger at higher…
▽ More
We study static magnetic susceptibility $χ(T, μ)$ in $SU(2)$ lattice gauge theory with $N_f = 2$ light flavours of dynamical fermions at finite chemical potential $μ$. Using linear response theory we find that $SU(2)$ gauge theory exhibits paramagnetic behavior in both the high-temperature deconfined regime and the low-temperature confining regime. Paramagnetic response becomes stronger at higher temperatures and larger values of the chemical potential. For our range of temperatures $0.727 \leq T/T_c \leq 2.67$, the first coefficient of the expansion of $χ(T, μ)$ in even powers of $μ/T$ around $μ=0$ is close to that of free quarks and lies in the range $(2 \ldots 5) \cdot 10^{-3}$. The strongest paramagnetic response is found in the diquark condensation phase at $μ> m_π/2$.
△ Less
Submitted 20 October, 2021; v1 submitted 20 April, 2021;
originally announced April 2021.
-
Numerical Study of the Chiral Separation Effect in Two-Color QCD at Finite Density
Authors:
P. V. Buividovich,
D. Smith,
L. von Smekal
Abstract:
We study the Chiral Separation Effect (CSE) in finite-density SU(2) lattice gauge theory with dynamical quarks. We find that the CSE is well described by the free quark result in the high-temperature quark-gluon plasma phase. As one enters the confinement regime with broken chiral symmetry at chemical potential smaller than half of the pion mass, the CSE response is gradually suppressed towards lo…
▽ More
We study the Chiral Separation Effect (CSE) in finite-density SU(2) lattice gauge theory with dynamical quarks. We find that the CSE is well described by the free quark result in the high-temperature quark-gluon plasma phase. As one enters the confinement regime with broken chiral symmetry at chemical potential smaller than half of the pion mass, the CSE response is gradually suppressed towards low temperatures in comparison to the free quark result. This suppression can be approximately described by assuming that the CSE current is proportional to the charge density, rather than the chemical potential, as suggested in the literature (ArXiv:1712.01256, Phys.Rev.D97(2018)085020). We also provide an upper bound on the contribution of disconnected fermionic diagrams to the CSE, which is consistent with zero within our statistical errors and small compared to that of the connected diagrams. Our results are obtained mainly in the QCD-like regime of SU(2) gauge theory at low densities, and hence should be at least qualitatively applicable to QCD as well.
△ Less
Submitted 16 August, 2021; v1 submitted 9 December, 2020;
originally announced December 2020.
-
Dilepton production in microscopic transport theory with in-medium $ρ$-meson spectral function
Authors:
A. B. Larionov,
U. Mosel,
L. von Smekal
Abstract:
We use the microscopic GiBUU transport model to calculate dilepton ($e^+e^-$) production in heavy-ion collisions at SIS18 energies focusing on the effect of collisional broadening of the $ρ$-meson. The collisional width of the $ρ$-meson at finite temperature and baryon density in nuclear matter is calculated on the basis of the collision integral of the GiBUU model. A systematic comparison with HA…
▽ More
We use the microscopic GiBUU transport model to calculate dilepton ($e^+e^-$) production in heavy-ion collisions at SIS18 energies focusing on the effect of collisional broadening of the $ρ$-meson. The collisional width of the $ρ$-meson at finite temperature and baryon density in nuclear matter is calculated on the basis of the collision integral of the GiBUU model. A systematic comparison with HADES data on dilepton production in heavy-ion collisions is performed. The collisional broadening of the $ρ$ improves the agreement between theory and experiment for the dilepton invariant-mass distributions near the $ρ$ pole mass and for the excess radiation in Au+Au at $1.23 A$ GeV. We furthermore show that some remaining underprediction of the experimental dilepton spectra in C+C at $1 A$ GeV and Au+Au at $1.23 A$ GeV at intermediate invariant masses $0.2-0.4$ GeV can be accounted for by adjusting the $pn$ bremsstrahlung cross section in a way to agree with the inclusive dilepton spectrum from $dp$ collisions at $1.25 A$ GeV.
△ Less
Submitted 4 March, 2021; v1 submitted 24 September, 2020;
originally announced September 2020.
-
Electric conductivity in finite-density SU(2) lattice gauge theory with dynamical fermions
Authors:
P. V. Buividovich,
D. Smith,
L. von Smekal
Abstract:
We study the dependence of the electric conductivity on chemical potential in finite-density $SU(2)$ gauge theory with $N_f = 2$ flavours of rooted staggered sea quarks, in combination with Wilson-Dirac and Domain Wall valence quarks. The pion mass is reasonably small with $m_π/m_ρ \approx 0.4$. We concentrate in particular on the vicinity of the chiral crossover, where we find the low-frequency e…
▽ More
We study the dependence of the electric conductivity on chemical potential in finite-density $SU(2)$ gauge theory with $N_f = 2$ flavours of rooted staggered sea quarks, in combination with Wilson-Dirac and Domain Wall valence quarks. The pion mass is reasonably small with $m_π/m_ρ \approx 0.4$. We concentrate in particular on the vicinity of the chiral crossover, where we find the low-frequency electric conductivity to be most sensitive to small changes in fermion density. Working in the low-density QCD-like regime with spontaneously broken chiral symmetry, we obtain an estimate of the first nontrivial coefficient $c(T)$ of the expansion of conductivity $σ(T,μ) = σ(T,0) \left(1 + c(T) (μ/T)^2 + O(μ^4)\right)$ in powers of $μ$, which has rather weak temperature dependence and takes its maximal value $c(T) \approx 0.10 \pm 0.07$ around the critical temperature. At larger densities and lower temperatures, the conductivity quickly grows towards the diquark condensation phase, and also becomes closer to the free quark result. As a by-product of our study we confirm the conclusions of previous studies with heavier pion that for $SU(2)$ gauge theory the ratio of crossover temperature to pion mass $T_c/m_π \approx 0.4$ at $μ=0$ is significantly smaller than in real QCD.
△ Less
Submitted 17 November, 2020; v1 submitted 10 July, 2020;
originally announced July 2020.
-
Spectral functions and dynamic critical behavior of relativistic $Z_2$ theories
Authors:
Dominik Schweitzer,
Sören Schlichting,
Lorenz von Smekal
Abstract:
We investigate the dynamic critical behaviour of a relativistic scalar field theory with $Z_2$ symmetry by calculating spectral functions of the order parameter at zero and non-vanishing momenta from first-principles classical-statistical lattice simulations in real-time. We find that at temperatures above the critical point $(T > T_c)$, the spectral functions are well described by relativistic qu…
▽ More
We investigate the dynamic critical behaviour of a relativistic scalar field theory with $Z_2$ symmetry by calculating spectral functions of the order parameter at zero and non-vanishing momenta from first-principles classical-statistical lattice simulations in real-time. We find that at temperatures above the critical point $(T > T_c)$, the spectral functions are well described by relativistic quasi-particle peaks. Close to the transition temperature $(T \sim T_c)$, we observe strong infrared contributions building up. In the ordered phase at low temperatures $(T < T_c)$, in addition to the quasi-particle peak, we observe a soft mode with a dispersion relation indicative of collective excitations. Investigating the spectral functions close to $T_c$, we demonstrate that the behavior in the vicinity of the critical point is controlled by dynamic scaling functions and the dynamic critical exponent $z$, which we determine from our simulations. By considering the equations of motion for a closed system and a system coupled to a heat bath, we extract the dynamic critical behavior for two different dynamic universality classes (Models A & C) in two and three spatial dimensions.
△ Less
Submitted 25 November, 2020; v1 submitted 7 July, 2020;
originally announced July 2020.
-
A density of states approach to the hexagonal Hubbard model at finite density
Authors:
Michael Körner,
Kurt Langfeld,
Dominik Smith,
Lorenz von Smekal
Abstract:
We apply the Linear Logarithmic Relaxation (LLR) method, which generalizes the Wang-Landau algorithm to quantum systems with continuous degrees of freedom, to the fermionic Hubbard model with repulsive interactions on the honeycomb lattice. We compute the generalized density of states of the average Hubbard field and divise two reconstruction schemes to extract physical observables from this resul…
▽ More
We apply the Linear Logarithmic Relaxation (LLR) method, which generalizes the Wang-Landau algorithm to quantum systems with continuous degrees of freedom, to the fermionic Hubbard model with repulsive interactions on the honeycomb lattice. We compute the generalized density of states of the average Hubbard field and divise two reconstruction schemes to extract physical observables from this result. By computing the particle density as a function of chemical potential we assess the utility of LLR in dealing with the sign problem of this model, which arises away from half filling. We show that the relative advantage over brute-force reweighting grows as the interaction strength is increased and discuss possible future improvements.
△ Less
Submitted 14 August, 2020; v1 submitted 8 June, 2020;
originally announced June 2020.
-
Fermionic excitations at finite temperature and density
Authors:
Ralf-Arno Tripolt,
Dirk H. Rischke,
Lorenz von Smekal,
Jochen Wambach
Abstract:
We study fermionic excitations in a hot and dense strongly interacting medium consisting of quarks and (pseudo-)scalar mesons. In particular, we use the two-flavor quark-meson model in combination with the Functional Renormalization Group approach, which allows to take into account the effects from thermal and quantum fluctuations. The resulting fermionic excitation spectrum is investigated by cal…
▽ More
We study fermionic excitations in a hot and dense strongly interacting medium consisting of quarks and (pseudo-)scalar mesons. In particular, we use the two-flavor quark-meson model in combination with the Functional Renormalization Group approach, which allows to take into account the effects from thermal and quantum fluctuations. The resulting fermionic excitation spectrum is investigated by calculating the quark spectral function at finite temperature, quark chemical potential, and spatial momentum. This involves an analytic continuation from imaginary to real energies by extending the previously introduced analytically-continued FRG (aFRG) method to the present case. We identify three different collective excitations in the medium: the ordinary thermal quark, the plasmino mode, and an ultra-soft "phonino" mode. The dispersion relations of these modes are extracted from the quark spectral function. When compared to corresponding results from an FRG-improved one-loop calculation remarkable agreement is found.
△ Less
Submitted 14 May, 2020; v1 submitted 26 March, 2020;
originally announced March 2020.
-
Quantum phase transitions on the hexagonal lattice
Authors:
Dominik Smith,
Pavel Buividovich,
Michael Körner,
Maksim Ulybyshev,
Lorenz von Smekal
Abstract:
Hubbard-type models on the hexagonal lattice are of great interest, as they provide realistic descriptions of graphene and other related materials. Hybrid Monte Carlo simulations offer a first-principles approach to study their phase structure. Here, we review the present status of our work in this direction.
Hubbard-type models on the hexagonal lattice are of great interest, as they provide realistic descriptions of graphene and other related materials. Hybrid Monte Carlo simulations offer a first-principles approach to study their phase structure. Here, we review the present status of our work in this direction.
△ Less
Submitted 20 January, 2021; v1 submitted 28 December, 2019;
originally announced December 2019.
-
Continuum Goldstone spectrum of two-color QCD at finite density with staggered quarks
Authors:
Jonas Wilhelm,
Lukas Holicki,
Dominik Smith,
Björn Wellegehausen,
Lorenz von Smekal
Abstract:
We carry out lattice simulations of two-color QCD and spectroscopy at finite density with two flavors of rooted-staggered quarks and a diquark source term. As in a previous four-flavor study, for small values of the inverse gauge coupling we observe a Goldstone spectrum which reflects the symmetry-breaking pattern of a Gaussian symplectic chiral random-matrix ensemble (GSE) with Dyson index…
▽ More
We carry out lattice simulations of two-color QCD and spectroscopy at finite density with two flavors of rooted-staggered quarks and a diquark source term. As in a previous four-flavor study, for small values of the inverse gauge coupling we observe a Goldstone spectrum which reflects the symmetry-breaking pattern of a Gaussian symplectic chiral random-matrix ensemble (GSE) with Dyson index $β_D=4$, which corresponds to any-color QCD with adjoint quarks in the continuum instead of QC$_2$D wih fundamental quarks. We show that this unphysical behavior occurs only inside of the bulk phase of $SU(2)$ gauge theory, where the density of $Z_2$ monopoles is high. Using an improved gauge action and a somewhat larger inverse coupling to suppress these monopoles, we demonstrate that the continuum Goldstone spectrum of two-color QCD, corresponding to a Gaussian orthogonal ensemble (GOE) with Dyson index $β_D=1$, is recovered also with rooted-staggered quarks once simulations are performed away from the bulk phase. We further demonstrate how this change of random-matrix ensemble is reflected in the distribution of eigenvalues of the Dirac operator. By computing the unfolded level spacings inside and outside of the bulk phase, we demonstrate that, starting with the low-lying eigenmodes which determine the infrared physics, the distribution of eigenmodes continuously changes from the GSE to the GOE one as monopoles are suppressed.
△ Less
Submitted 10 October, 2019;
originally announced October 2019.
-
Fluctuating vector mesons in analytically continued FRG flow equations
Authors:
Christopher Jung,
Lorenz von Smekal
Abstract:
In this work we study contributions due to vector and axial-vector meson fluctuations to their in-medium spectral functions in an effective low-energy theory inspired by the gauged linear sigma model. In particular, we show how to describe these fluctuations in the effective theory by massive (axial-)vector fields in agreement with the known structure of analogous single-particle or resonance cont…
▽ More
In this work we study contributions due to vector and axial-vector meson fluctuations to their in-medium spectral functions in an effective low-energy theory inspired by the gauged linear sigma model. In particular, we show how to describe these fluctuations in the effective theory by massive (axial-)vector fields in agreement with the known structure of analogous single-particle or resonance contributions to the corresponding conserved currents. The vector and axial-vector meson spectral functions are then computed by numerically solving the analytically continued functional renormalization group (aFRG) flow equations for their retarded two-point functions at finite temperature and density in the effective theory. We identify the new contributions that arise due to the (axial-)vector meson fluctuations, and assess their influence on possible signatures of a QCD critical endpoint and the restoration of chiral symmetry in thermal dilepton spectra.
△ Less
Submitted 30 September, 2019;
originally announced September 2019.
-
Spectral functions and critical dynamics of the $O(4)$ model from classical-statistical lattice simulations
Authors:
Sören Schlichting,
Dominik Smith,
Lorenz von Smekal
Abstract:
We calculate spectral functions of the relativistic $O(4)$ model from real-time lattice simulations in classical-statistical field theory. While in the low and high temperature phase of the model, the spectral functions of longitudinal $(σ)$ and transverse $(π)$ modes are well described by relativistic quasi-particle peaks, we find a highly non-trivial behavior of the spectral functions in the cro…
▽ More
We calculate spectral functions of the relativistic $O(4)$ model from real-time lattice simulations in classical-statistical field theory. While in the low and high temperature phase of the model, the spectral functions of longitudinal $(σ)$ and transverse $(π)$ modes are well described by relativistic quasi-particle peaks, we find a highly non-trivial behavior of the spectral functions in the cross over region, where additional structures appear. Similarly, we observe a significant broadening of the quasi-particle peaks, when the amount explicit $O(4)$ symmetry breaking is reduced. We further demonstrate that in the vicinity of the $O(4)$ critical point, the spectral functions develop an infrared power law associated with the critical dynamics, and comment on the extraction of the dynamical critical exponent $z$ from our simulations.
△ Less
Submitted 23 July, 2019;
originally announced August 2019.
-
Numerical evidence of conformal phase transition in graphene with long-range interactions
Authors:
Pavel Buividovich,
Dominik Smith,
Maksim Ulybyshev,
Lorenz von Smekal
Abstract:
Using state of the art Hybrid-Monte-Carlo (HMC) simulations we carry out an unbiased study of the competition between spin-density wave (SDW) and charge-density wave (CDW) order in suspended graphene. We determine that the realistic inter-electron potential of graphene must be scaled up by a factor of roughly 1.6 to induce a semimetal-SDW phase transition and find no evidence for CDW order. A stud…
▽ More
Using state of the art Hybrid-Monte-Carlo (HMC) simulations we carry out an unbiased study of the competition between spin-density wave (SDW) and charge-density wave (CDW) order in suspended graphene. We determine that the realistic inter-electron potential of graphene must be scaled up by a factor of roughly 1.6 to induce a semimetal-SDW phase transition and find no evidence for CDW order. A study of critical properties suggests that the universality class of the three-dimensional chiral Heisenberg Gross-Neveu model with two fermion flavors, predicted by renormalization group studies and strong-coupling expansion, is unlikely to apply to this transition. We propose that our results instead favor an interpretation in terms of a conformal phase transition. In addition, we describe a variant of the HMC algorithm which uses exact fermionic forces during molecular dynamics trajectories and avoids the use of pseudofermions. Compared to standard HMC this allows for a substantial increase of the integrator stepsize while achieving comparable Metropolis acceptance rates and leads to a sizable performance improvement.
△ Less
Submitted 17 May, 2019; v1 submitted 16 December, 2018;
originally announced December 2018.
-
The Anderson transition in QCD with $N_f=2+1+1$ twisted mass quarks: overlap analysis
Authors:
Lukas Holicki,
Ernst-Michael Ilgenfritz,
Lorenz von Smekal
Abstract:
Chiral Random Matrix Theory has proven to describe the spectral properties of low temperature QCD very well. However, at temperatures above the chiral symmetry restoring transition it can not provide a global description. The level-spacing distribution in the lower part of the spectrum of the Dirac operator is Poisson-like. There the eigenmodes are localized in space-time and separated from the re…
▽ More
Chiral Random Matrix Theory has proven to describe the spectral properties of low temperature QCD very well. However, at temperatures above the chiral symmetry restoring transition it can not provide a global description. The level-spacing distribution in the lower part of the spectrum of the Dirac operator is Poisson-like. There the eigenmodes are localized in space-time and separated from the rest of the spectrum by a so-called mobility edge. In analogy to Anderson localization in condensed-matter systems with random disorder this has been called the QCD-Anderson transition. Here, we study the localization features of the low-lying eigenmodes of the massless overlap operator on configurations generated with $N_f=2+1+1$ twisted mass Wilson sea quarks and present results concerning the temperature dependence of the mobility edge and the mechanism of the quark-mode localization. We have used various methods to fix the spectral position of the delocalization transition and verify that the mobility edge extrapolates to zero at a temperature within the chiral transition region.
△ Less
Submitted 2 October, 2018;
originally announced October 2018.
-
Fermionic spectral functions with the Functional Renormalization Group
Authors:
Ralf-Arno Tripolt,
Johannes Weyrich,
Lorenz von Smekal,
Jochen Wambach
Abstract:
We present first results on the calculation of fermionic spectral functions from analytically continued flow equations within the Functional Renormalization Group approach. Our method is based on the same analytic continuation from imaginary to real frequencies that was developed and used previously for bosonic spectral functions. In order to demonstrate the applicability of the method also for fe…
▽ More
We present first results on the calculation of fermionic spectral functions from analytically continued flow equations within the Functional Renormalization Group approach. Our method is based on the same analytic continuation from imaginary to real frequencies that was developed and used previously for bosonic spectral functions. In order to demonstrate the applicability of the method also for fermionic correlations we apply it here to the real-time quark propagator in the quark-meson model and calculate the corresponding quark spectral functions in the vacuum.
△ Less
Submitted 1 November, 2018; v1 submitted 31 July, 2018;
originally announced July 2018.
-
Hybrid-Monte-Carlo study of competing order in the extended fermionic Hubbard model on the hexagonal lattice
Authors:
Pavel Buividovich,
Dominik Smith,
Maksim Ulybyshev,
Lorenz von Smekal
Abstract:
Using first-principle Hybrid-Monte-Carlo (HMC) simulations, we carry out an unbiased study of the competition between spin-density wave (SDW) and charge-density wave (CDW) order in the extended Hubbard model on the two dimensional hexagonal lattice at half filling. We determine the phase diagram in the space of on-site and nearest-neighbor couplings $U$ and $V$ in the region $V<U/3$, which can be…
▽ More
Using first-principle Hybrid-Monte-Carlo (HMC) simulations, we carry out an unbiased study of the competition between spin-density wave (SDW) and charge-density wave (CDW) order in the extended Hubbard model on the two dimensional hexagonal lattice at half filling. We determine the phase diagram in the space of on-site and nearest-neighbor couplings $U$ and $V$ in the region $V<U/3$, which can be simulated without a fermion sign problem, and find that a transition from semimetal to a SDW phase occurs at sufficiently large $U$ for basically all $V$. Tracing the corresponding phase boundary from $V=0$ to the $V=U/3$ line, we find evidence for critical scaling in the Gross-Neveu universality class for the entire boundary. With rather high confidence we rule out the existence of the CDW ordered phase anywhere in the range of parameters considered. We also discuss several improvements of the HMC algorithm which are crucial to reach these conclusions, in particular the improved fermion action with exact sublattice symmetry and the complexification of the Hubbard-Stratonovich field to ensure the ergodicity of the algorithm.
△ Less
Submitted 6 December, 2018; v1 submitted 18 July, 2018;
originally announced July 2018.
-
In-medium spectral functions and dilepton rates with the Functional Renormalization Group
Authors:
Ralf-Arno Tripolt,
Christopher Jung,
Naoto Tanji,
Lorenz von Smekal,
Jochen Wambach
Abstract:
We present recent results on in-medium spectral functions of vector and axial-vector mesons, the electromagnetic (EM) spectral function and dilepton rates using the Functional Renormalization Group (FRG) approach. Our method is based on an analytic continuation procedure that allows us to calculate real-time quantities like spectral functions at finite temperature and chemical potential. As an eff…
▽ More
We present recent results on in-medium spectral functions of vector and axial-vector mesons, the electromagnetic (EM) spectral function and dilepton rates using the Functional Renormalization Group (FRG) approach. Our method is based on an analytic continuation procedure that allows us to calculate real-time quantities like spectral functions at finite temperature and chemical potential. As an effective low-energy model for Quantum Chromodynamics (QCD) we use an extended linear-sigma model including quarks where (axial-)vector mesons as well as the photon are introduced as gauge bosons. In particular, it is shown how the $ρ$ and the $a_1$ spectral function become degenerate at high temperatures or chemical potentials due to the restoration of chiral symmetry. Preliminary results for the EM spectral function and the dilepton production rate are discussed with a focus on the possibility to identify signatures of the chiral crossover and the critical endpoint (CEP) in the QCD phase diagram.
△ Less
Submitted 13 July, 2018;
originally announced July 2018.
-
Numerical analytic continuation of Euclidean data
Authors:
Ralf-Arno Tripolt,
Philipp Gubler,
Maksim Ulybyshev,
Lorenz von Smekal
Abstract:
In this work we present a direct comparison of three different numerical analytic continuation methods: the Maximum Entropy Method, the Backus-Gilbert method and the Schlessinger point or Resonances Via Padé method. First, we perform a benchmark test based on a model spectral function and study the regime of applicability of these methods depending on the number of input points and their statistic…
▽ More
In this work we present a direct comparison of three different numerical analytic continuation methods: the Maximum Entropy Method, the Backus-Gilbert method and the Schlessinger point or Resonances Via Padé method. First, we perform a benchmark test based on a model spectral function and study the regime of applicability of these methods depending on the number of input points and their statistical error. We then apply these methods to more realistic examples, namely to numerical data on Euclidean propagators obtained from a Functional Renormalization Group calculation, to data from a lattice Quantum Chromodynamics simulation and to data obtained from a tight-binding model for graphene in order to extract the electrical conductivity.
△ Less
Submitted 19 November, 2018; v1 submitted 31 January, 2018;
originally announced January 2018.
-
Spectral Functions from the Functional Renormalization Group
Authors:
Jochen Wambach,
Christopher Jung,
Fabian Rennecke,
Ralf-Arno Tripolt,
Lorenz von Smekal
Abstract:
We present results for in-medium spectral functions obtained within the Functional Renormalization Group framework. The analytic continuation from imaginary to real time is performed in a well-defined way on the level of the flow equations. Based on this recently developed method, results for the sigma and the pion spectral function for the quark-meson model are shown at finite temperature, finite…
▽ More
We present results for in-medium spectral functions obtained within the Functional Renormalization Group framework. The analytic continuation from imaginary to real time is performed in a well-defined way on the level of the flow equations. Based on this recently developed method, results for the sigma and the pion spectral function for the quark-meson model are shown at finite temperature, finite quark-chemical potential and finite spatial momentum. It is shown how these spectral function become degenreate at high temperatures due to the restoration of chiral symmetry. In addition, results for vector- and axial-vector meson spectral functions are shown using a gauged linear sigma model with quarks. The degeneration of the $ρ$ and the $a_1$ spectral function as well as the behavior of their pole masses is discussed.
△ Less
Submitted 6 December, 2017;
originally announced December 2017.
-
The role of frequency dependence in dynamical gap generation in graphene
Authors:
M. E. Carrington,
C. S. Fischer,
L. von Smekal,
M. H. Thoma
Abstract:
We study the frequency dependencies of the fermion and photon dressing functions in dynamical gap generation in graphene. We use a low energy effective QED-like description, but within this approximation, we include all frequency dependent effects including retardation. We obtain the critical coupling by calculating the gap using a non-perturbative Dyson-Schwinger approach. Compared to the results…
▽ More
We study the frequency dependencies of the fermion and photon dressing functions in dynamical gap generation in graphene. We use a low energy effective QED-like description, but within this approximation, we include all frequency dependent effects including retardation. We obtain the critical coupling by calculating the gap using a non-perturbative Dyson-Schwinger approach. Compared to the results of our previous calculation [1] which used a Lindhard screening approximation instead of including a self-consistently calculated dynamical screening function, the critical coupling is substantially reduced.
△ Less
Submitted 8 March, 2018; v1 submitted 3 November, 2017;
originally announced November 2017.
-
Spectrum of QCD at Finite Isospin Density
Authors:
Philipp Scior,
Lorenz von Smekal,
Dominik Smith
Abstract:
We study the phase diagram of QCD at finite isospin density using two flavors of staggered quarks. We investigate the low temperature region of the phase diagram where we find a pion condensation phase at high chemical potential. We started a basic analysis of the spectrum at finite isospin density. In particular, we measured pion, rho and nucleon masses inside and outside of the pion condensation…
▽ More
We study the phase diagram of QCD at finite isospin density using two flavors of staggered quarks. We investigate the low temperature region of the phase diagram where we find a pion condensation phase at high chemical potential. We started a basic analysis of the spectrum at finite isospin density. In particular, we measured pion, rho and nucleon masses inside and outside of the pion condensation phase. In agreement with previous studies in two-color QCD at finite baryon density we find that the Polyakov loop does not depend on the density in the staggered formulation.
△ Less
Submitted 20 October, 2017; v1 submitted 17 October, 2017;
originally announced October 2017.
-
The low-temperature behavior of the quark-meson model
Authors:
Ralf-Arno Tripolt,
Bernd-Jochen Schaefer,
Lorenz von Smekal,
Jochen Wambach
Abstract:
We revisit the phase diagram of strong-interaction matter for the two-flavor quark-meson model using the Functional Renormalization Group. In contrast to standard mean-field calculations, an unusual phase structure is encountered at low temperatures and large quark chemical potentials. In particular, we identify a regime where the pressure decreases with increasing temperature and discuss possible…
▽ More
We revisit the phase diagram of strong-interaction matter for the two-flavor quark-meson model using the Functional Renormalization Group. In contrast to standard mean-field calculations, an unusual phase structure is encountered at low temperatures and large quark chemical potentials. In particular, we identify a regime where the pressure decreases with increasing temperature and discuss possible reasons for this unphysical behavior.
△ Less
Submitted 19 February, 2018; v1 submitted 18 September, 2017;
originally announced September 2017.
-
A Hybrid-Monte-Carlo study of monolayer graphene with partially screened Coulomb interactions at finite spin density
Authors:
Michael Körner,
Dominik Smith,
Pavel Buividovich,
Maksim Ulybyshev,
Lorenz von Smekal
Abstract:
We report on Hybrid-Monte-Carlo simulations at finite spin density of the $π$-band electrons in monolayer graphene with realistic inter-electron interactions. Unlike simulations at finite charge-carrier density, these are not affected by a fermion-sign problem. Our results are in qualitative agreement with an interaction-induced warping of the Fermi contours, and a reduction of the bandwidth as ob…
▽ More
We report on Hybrid-Monte-Carlo simulations at finite spin density of the $π$-band electrons in monolayer graphene with realistic inter-electron interactions. Unlike simulations at finite charge-carrier density, these are not affected by a fermion-sign problem. Our results are in qualitative agreement with an interaction-induced warping of the Fermi contours, and a reduction of the bandwidth as observed in angle resolved photoemission spectroscopy experiments on charge-doped graphene systems. Furthermore, we find evidence that the neck-disrupting Lifshitz transition, which occurs when the Fermi level traverses the van Hove singularity (VHS), becomes a true quantum phase transition due to interactions. This is in-line with an instability of the VHS towards the formation of electronic ordered phases, which has been predicted by a variety of different theoretical approaches.
△ Less
Submitted 7 November, 2017; v1 submitted 12 April, 2017;
originally announced April 2017.
-
Inter-electron interactions and the RKKY potential between H adatoms in graphene
Authors:
Pavel Buividovich,
Dominik Smith,
Maksim Ulybyshev,
Lorenz von Smekal
Abstract:
We use first-principles Quantum Monte-Carlo simulations to study the Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction between hydrogen adatoms attached to a graphene sheet. We find that the pairwise RKKY interactions at distances of a few lattice spacings are strongly affected by inter-electron interactions, in particular, the potential barrier between widely separated adatoms and the dimer config…
▽ More
We use first-principles Quantum Monte-Carlo simulations to study the Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction between hydrogen adatoms attached to a graphene sheet. We find that the pairwise RKKY interactions at distances of a few lattice spacings are strongly affected by inter-electron interactions, in particular, the potential barrier between widely separated adatoms and the dimer configuration becomes wider and thus harder to penetrate. We also point out that anti-ferrromagnetic and charge density wave orderings have very different effects on the RKKY interaction. Finally, we analyze the stability of several regular adatom superlattices with respect to small displacements of a single adatom, distinguishing the cases of adatoms which populate either both or only one sublattice of the graphene lattice.
△ Less
Submitted 22 September, 2017; v1 submitted 16 March, 2017;
originally announced March 2017.
-
What we can learn from two-dimensional QCD-like theories at finite density
Authors:
Björn H. Wellegehausen,
Lorenz von Smekal
Abstract:
We study generic properties of strongly interacting matter at finite density as relevant to heavy-ion collisions at moderate beam energies or the physics of neutron stars and their mergers. Because of the fermion-sign problem in lattice QCD, here we simulate QCD-like theories without this problem at finite density. These theories (two-color QCD, G2-QCD, or adjoint QCD) typically contain bosonic ba…
▽ More
We study generic properties of strongly interacting matter at finite density as relevant to heavy-ion collisions at moderate beam energies or the physics of neutron stars and their mergers. Because of the fermion-sign problem in lattice QCD, here we simulate QCD-like theories without this problem at finite density. These theories (two-color QCD, G2-QCD, or adjoint QCD) typically contain bosonic baryons, for example diquarks, or other more exotic states of matter. It is therefore important to understand the effects of such bosonic matter and disentangle them from fermionic baryons where they exist to draw conclusions for QCD. Simulations of these theories, for instance G2-QCD, reveal an interesting and rich phase diagram at zero temperature. Many open questions arise, partly due to the lack of high precision or large volume/continuum data. This is the reason why we study two-dimensional QCD-like theories. In this contribution we shall discuss differences between QCD-like theories at baryon chemical and isospin chemical potential. Furthermore we present simulation results on the phase diagram and spectroscopy at finite density for G2- and two-color-QCD and compare it to free lattice fermions.
△ Less
Submitted 1 February, 2017;
originally announced February 2017.
-
Two-colour QCD at finite density with two flavours of staggered quarks
Authors:
Lukas Holicki,
Jonas Wilhelm,
Dominik Smith,
Björn Wellegehausen,
Lorenz von Smekal
Abstract:
In this contribution we revisit simulations of two-color QCD with rooted staggered quarks at finite density, where baryon-number spontaneously breaks and a diquark condensate forms. We thereby pay special attention to simulating outside the lattice-artifact bulk phase, in which $Z_2$ monopoles condense, and investigate some of the consequences of this, e.g. on the chiral and the diquark condensate…
▽ More
In this contribution we revisit simulations of two-color QCD with rooted staggered quarks at finite density, where baryon-number spontaneously breaks and a diquark condensate forms. We thereby pay special attention to simulating outside the lattice-artifact bulk phase, in which $Z_2$ monopoles condense, and investigate some of the consequences of this, e.g. on the chiral and the diquark condensate which were known to be well described by chiral effective field theory. Not surprisingly, on finer lattices outside the bulk phase the quark condensate now requires additive renormalization before it can be compared with effective field theory predictions. The subtraction must necessarily depend on the chemical potential, however. The diquark condensate is not affected by this problem and remains in good agreement with these predictions. We also compare staggered with Wilson quarks to demonstrate that the two fermion discretizations yield qualitatively different results well below half-filling already. We close with prelimiary results for the Goldstone spectrum to demonstrate that the continuum pattern is recovered also with staggered quarks outside the bulk phase.
△ Less
Submitted 17 January, 2017;
originally announced January 2017.
-
Non-Local effective SU(2) Polyakov-loop models from inverse Monte-Carlo methods
Authors:
Bardiya Bahrampour,
Björn Wellegehausen,
Lorenz von Smekal
Abstract:
The strong-coupling expansion of the lattice gauge action leads to Polyakov-loop models that effectively describe gluodynamics at low temperatures, and together with the hopping expansion of the fermion determinant provides insight into the QCD phase diagram at finite density and low temperatures, although for rather heavy quarks. At higher temperatures the strong-coupling expansion breaks down an…
▽ More
The strong-coupling expansion of the lattice gauge action leads to Polyakov-loop models that effectively describe gluodynamics at low temperatures, and together with the hopping expansion of the fermion determinant provides insight into the QCD phase diagram at finite density and low temperatures, although for rather heavy quarks. At higher temperatures the strong-coupling expansion breaks down and it is expected that the interactions between Polyakov loops become non-local. Here, we therefore test how well pure SU(2) gluodynamics can be mapped onto different non-local Polyakov models with inverse Monte-Carlo methods. We take into account Polyakov loops in higher representations and gradually add interaction terms at larger distances. We are particularly interested in extrapolating the range of non-local terms in sufficiently large volumes and higher representations. We study the characteristic fall-off in strength of the non-local couplings with the interaction distance, and its dependence on the gauge coupling in order to compare our results to existing proposals for non-local effective actions.
△ Less
Submitted 1 December, 2016;
originally announced December 2016.
-
Competing order in the fermionic Hubbard model on the hexagonal graphene lattice
Authors:
Pavel Buividovich,
Dominik Smith,
Maksim Ulybyshev,
Lorenz von Smekal
Abstract:
We study the phase diagram of the fermionic Hubbard model on the hexagonal lattice in the space of on-site and nearest neighbor couplings with Hybrid-Monte-Carlo simulations. With pure on-site repulsion this allows to determine the critical coupling strength for spin-density wave formation with the standard approach of introducing a small mass term, explicitly breaking the sublattice symmetry. The…
▽ More
We study the phase diagram of the fermionic Hubbard model on the hexagonal lattice in the space of on-site and nearest neighbor couplings with Hybrid-Monte-Carlo simulations. With pure on-site repulsion this allows to determine the critical coupling strength for spin-density wave formation with the standard approach of introducing a small mass term, explicitly breaking the sublattice symmetry. The analogous mass term for charge-density wave formation above a critical nearest-neighbor repulsion, on the other hand, would introduce a fermion sign problem. The competition between the two and the phase diagram in the space of the two coouplings can however be studied in simulations without explicit sublattice symmetry breaking. Our results compare qualitatively well with the Hartree-Fock phase diagram. We furthermore demonstrate how spin-symmetry breaking by the Euclidean time discretization can be avoided also, when using an improved fermion action based on an exponetial transfer matrix with exact sublattice symmetry.
△ Less
Submitted 31 October, 2016;
originally announced October 2016.
-
In-Medium Spectral Functions of Vector- and Axial-Vector Mesons from the Functional Renormalization Group
Authors:
Christopher Jung,
Fabian Rennecke,
Ralf-Arno Tripolt,
Lorenz von Smekal,
Jochen Wambach
Abstract:
In this work we present first results on vector and axial-vector meson spectral functions as obtained by applying the non-perturbative functional renormalization group approach to an effective low-energy theory motivated by the gauged linear sigma model. By using a recently proposed analytic continuation method, we study the in-medium behavior of the spectral functions of the $ρ$ and $a_1$ mesons…
▽ More
In this work we present first results on vector and axial-vector meson spectral functions as obtained by applying the non-perturbative functional renormalization group approach to an effective low-energy theory motivated by the gauged linear sigma model. By using a recently proposed analytic continuation method, we study the in-medium behavior of the spectral functions of the $ρ$ and $a_1$ mesons in different regimes of the phase diagram. In particular, we demonstrate explicitly how these spectral functions degenerate at high temperatures as well as at large chemical potentials, as a consequence of the restoration of chiral symmetry. In addition, we also compute the momentum dependence of the $ρ$ and $a_1$ spectral functions and discuss the various time-like and space-like processes that can occur.
△ Less
Submitted 27 October, 2016;
originally announced October 2016.
-
Dynamical gap generation in graphene with frequency dependent renormalization effects
Authors:
M. E. Carrington,
C. S. Fischer,
L. von Smekal,
M. H. Thoma
Abstract:
We study the frequency dependencies in the renormalization of the fermion Greens function for the $π$-band electrons in graphene and their influence on the dynamical gap generation at sufficiently strong interaction. Adopting the effective QED-like description for the low-energy excitations within the Dirac-cone region we self consistently solve the fermion Dyson-Schwinger equation in various appr…
▽ More
We study the frequency dependencies in the renormalization of the fermion Greens function for the $π$-band electrons in graphene and their influence on the dynamical gap generation at sufficiently strong interaction. Adopting the effective QED-like description for the low-energy excitations within the Dirac-cone region we self consistently solve the fermion Dyson-Schwinger equation in various approximations for the photon propagator and the vertex function with special emphasis on frequency dependent Lindhard screening and retardation effects.
△ Less
Submitted 5 May, 2016;
originally announced May 2016.
-
Spectral functions and in-medium properties of hadrons
Authors:
Ralf-Arno Tripolt,
Lorenz von Smekal,
Jochen Wambach
Abstract:
The in-medium modifications of hadron properties such as masses and decay widths have been a major focus of the scientific work of Gerry Brown and the insights gained by him and his collaborators made them major drivers of this field for several decades. Their prediction of experimental signals in di-lepton pair production in relativistic heavy-ion collisions were instrumental in initiating large…
▽ More
The in-medium modifications of hadron properties such as masses and decay widths have been a major focus of the scientific work of Gerry Brown and the insights gained by him and his collaborators made them major drivers of this field for several decades. Their prediction of experimental signals in di-lepton pair production in relativistic heavy-ion collisions were instrumental in initiating large experimental campaigns which continue until today. In this chapter we review recent results which elucidate the relation of hadronic spectral properties at finite temperature and density to the restoration of spontaneously broken chiral symmetry.
△ Less
Submitted 3 May, 2016;
originally announced May 2016.
-
Effective Polyakov loop models for QCD-like theories at finite chemical potential
Authors:
Philipp Scior,
Lorenz von Smekal
Abstract:
We study genuine finite density effects in QCD-like theories with three-dimensional Polyakov-loop effective theories for heavy quarks. These are derived from the full QCD-like theories by combined strong-coupling and hopping expansions. In particular, we investigate the cold and dense regimes of phase diagrams where we expect to find Bose-Einstein-condensation of diquark baryons or a fermionic fir…
▽ More
We study genuine finite density effects in QCD-like theories with three-dimensional Polyakov-loop effective theories for heavy quarks. These are derived from the full QCD-like theories by combined strong-coupling and hopping expansions. In particular, we investigate the cold and dense regimes of phase diagrams where we expect to find Bose-Einstein-condensation of diquark baryons or a fermionic first-order liquid-gas transition, depending on the gauge group of the theory. In two-color QCD, for example, we observe evidence of a continuous zero-temperature transition to finite diquark density when the quark chemical potential $μ$ reaches half the diquark mass, i.e. without binding energy. In G$_2$-QCD we observe, in addition to this "Silver Blaze" onset of diquark density, a second transition in the density towards an exponential increase by roughly $3μ/T$ corresponding to a finite density of G$_2$-nucleons.
△ Less
Submitted 15 February, 2016;
originally announced February 2016.
-
Baryonic Matter Onset in Two-Color QCD with Heavy Quarks
Authors:
Philipp Scior,
Lorenz von Smekal
Abstract:
We study the cold and dense regime in the phase diagram of two-color QCD with heavy quarks within a three-dimensional effective theory for Polyakov loops. This theory is derived from two-color QCD in a combined strong-coupling and hopping expansion. In particular, we study the onset of diquark density as the finite-density transition of the bosonic baryons in the two-color world. In contrast to pr…
▽ More
We study the cold and dense regime in the phase diagram of two-color QCD with heavy quarks within a three-dimensional effective theory for Polyakov loops. This theory is derived from two-color QCD in a combined strong-coupling and hopping expansion. In particular, we study the onset of diquark density as the finite-density transition of the bosonic baryons in the two-color world. In contrast to previous studies of heavy dense QCD, our zero-temperature extrapolations are consistent with a continuous transition without binding energy. They thus provide evidence that the effective theory for heavy quarks is capable of describing the characteristic differences between diquark condensation in two-color QCD and the liquid-gas transition of nuclear matter in QCD.
△ Less
Submitted 3 August, 2015;
originally announced August 2015.
-
On Dyson-Schwinger studies of Yang-Mills theory and the four-gluon vertex
Authors:
Markus Q. Huber,
Anton K. Cyrol,
Lorenz von Smekal
Abstract:
We review the status of calculations of Yang-Mills Green functions from Dyson-Schwinger equations. The role of truncations is discussed and results for the four-gluon vertex are presented.
We review the status of calculations of Yang-Mills Green functions from Dyson-Schwinger equations. The role of truncations is discussed and results for the four-gluon vertex are presented.
△ Less
Submitted 12 June, 2015;
originally announced June 2015.
-
Chiral Mirror-Baryon-Meson Model and Nuclear Matter beyond Mean-Field
Authors:
Johannes Weyrich,
Nils Strodthoff,
Lorenz von Smekal
Abstract:
We consider a chiral baryon-meson model for nucleons and their parity partners in mirror assignment interacting with pions, sigma and omega mesons to describe the liquid-gas transition of nuclear matter together with chiral symmetry restoration in the high density phase. Within the mean-field approximation the model is known to provide a phenomenologically successful description of the nuclear-mat…
▽ More
We consider a chiral baryon-meson model for nucleons and their parity partners in mirror assignment interacting with pions, sigma and omega mesons to describe the liquid-gas transition of nuclear matter together with chiral symmetry restoration in the high density phase. Within the mean-field approximation the model is known to provide a phenomenologically successful description of the nuclear-matter transition. Here, we go beyond this approximation and include mesonic fluctuations by means of the functional renormalization group. While these fluctuations do not lead to major qualitative changes in the phase diagram of the model, beyond mean-field, one is no-longer free to adjust the parameters so as to reproduce the binding energy per nucleon, the nuclear saturation density, and the nucleon sigma term all at the same time. However, the prediction of a clear first-order chiral transition at low temperatures inside the high baryon-density phase appears to be robust.
△ Less
Submitted 10 July, 2015; v1 submitted 10 April, 2015;
originally announced April 2015.
-
Glueballs from the Bethe-Salpeter equation
Authors:
Helios Sanchis-Alepuz,
Christian S. Fischer,
Christian Kellermann,
Lorenz von Smekal
Abstract:
We formulate a framework to determine the mass of glueball states of Landau gauge Yang-Mills theory in the continuum. To this end we derive a Bethe-Salpeter equation for two gluon bound states including the effects of Faddeev-Popov ghosts. We construct a suitable approximation scheme such that the interactions in the bound state equation match a corresponding successful approximation of the Dyson-…
▽ More
We formulate a framework to determine the mass of glueball states of Landau gauge Yang-Mills theory in the continuum. To this end we derive a Bethe-Salpeter equation for two gluon bound states including the effects of Faddeev-Popov ghosts. We construct a suitable approximation scheme such that the interactions in the bound state equation match a corresponding successful approximation of the Dyson-Schwinger equations for the Landau gauge ghost and gluon propagators. Based upon a recently obtained solution for the propagators in the complex momentum plane we obtain results for the mass of the scalar and pseudoscalar glueballs. In the scalar channel we find a mass value in agreement with lattice gauge theory.
△ Less
Submitted 20 March, 2015;
originally announced March 2015.
-
Test of the Atiyah-Singer Index Theorem for Fullerene with a Superconducting Microwave Resonator
Authors:
B. Dietz,
T. Klaus,
M. Miski-Oglu,
A. Richter,
M. Bischoff,
L. von Smekal,
J. Wambach
Abstract:
Experiments have been performed using a spherical superconducting microwave resonator that simulates the geometric structure of the C60 fullerene molecule. The objective was to study with very high resolution the exceptional spectral properties emerging from the symmetries of the icosahedral structure of the carbon lattice. In particular, the number of zero modes has been determined to test the pr…
▽ More
Experiments have been performed using a spherical superconducting microwave resonator that simulates the geometric structure of the C60 fullerene molecule. The objective was to study with very high resolution the exceptional spectral properties emerging from the symmetries of the icosahedral structure of the carbon lattice. In particular, the number of zero modes has been determined to test the predictions of the Atiyah-Singer index theorem, which relates it to the topology of the curved carbon lattice. This is, to the best of our knowledge, the first experimental verification of the index theorem.
△ Less
Submitted 13 March, 2015;
originally announced March 2015.
-
Lattice simulations of $G_2$-QCD at finite density
Authors:
Bjoern H. Wellegehausen,
Lorenz von Smekal
Abstract:
$G_2$-QCD, in which the exceptional Lie group $G_2$ replaces the $SU(3)…
▽ More
$G_2$-QCD, in which the exceptional Lie group $G_2$ replaces the $SU(3)$ gauge group of QCD, does not suffer from a fermion sign problem. It can therefore be simulated also at comparatively low temperatures and high densities on the lattice, which hence allows to map out the phase diagram of this QCD-like theory. We briefly review some of our previous results from such finite density simulations to then present further evidence for a first-order transition to what might be called $G_2$-nuclear matter. In order to isolate diquark condensation effects, we introduce simulations with Majorana fermions and diquark sources. This allows to disentangle states in the spectrum that are connected by charge conjugation. We discuss chiral symmetry in the presence of diquark sources and present first results from our ongoing large-scale simulations.
△ Less
Submitted 27 January, 2015;
originally announced January 2015.
-
Effective SU(2) Polyakov Loop Theories with Heavy Quarks on the Lattice
Authors:
Philipp Scior,
David Scheffler,
Dominik Smith,
Lorenz von Smekal
Abstract:
We compare SU(2) Polyakov loop models with different effective actions with data from full two-color QCD simulations around and above the critical temperature. We then apply the effective theories at finite temperature and density to extract quantities like Polyakov loop correlators, effective Polyakov loop potentials and baryon density.
We compare SU(2) Polyakov loop models with different effective actions with data from full two-color QCD simulations around and above the critical temperature. We then apply the effective theories at finite temperature and density to extract quantities like Polyakov loop correlators, effective Polyakov loop potentials and baryon density.
△ Less
Submitted 22 December, 2014;
originally announced December 2014.