Finite Density Research Articles

QCD equation of state and thermodynamic observables from computationally minimal Dyson-Schwinger equations

We study the QCD equation of state and other thermodynamic observables including the isentropic trajectories and the speed of sound. These observables are of eminent importance for the understanding of experimental results in heavy ion collisions and also provide a QCD input for studies of the timeline of heavy-ion-collisions with hydrodynamical simulations. They can be derived from the quark propagator whose gap equation is solved within a minimal approximation to the Dyson-Schwinger equations (DSEs) of QCD at finite temperature and density. This minimal approximation aims at a combination of computational efficiency and simplification of the truncation scheme while maintaining quantitative precision. This minimal DSE scheme is confronted and benchmarked with results for correlation functions and observables from lattice QCD at vanishing density and quantitative functional approaches at finite density. Published by the American Physical Society 2024

Axionlike quasiparticles and topological states of matter: Finite density corrections of the chiral anomaly vertex

We investigate the general structure of the chiral anomaly AVV/AAA and (LLL,RRR) vertices, in the presence of chemical potentials in perturbation theory. The study finds application in anomalous transport, whenever chirally unbalanced matter is present, with propagating external currents that are classically conserved. Examples are topological materials and the chiral magnetic effect in the plasma state of matter of the early universe. We classify the minimal number of form factors of the AVV parametrization, by a complete analysis of the Schouten identities in the presence of a heat bath. We show that the longitudinal (anomaly) sector in the axial-vector channel, for on-shell and off-shell photons, is protected against corrections coming from the insertion of a chemical potential in the fermion loop. When the photons are on-shell, we prove that the transverse sector, in the same channel, is also μ-independent and vanishes. The related effective action is shown to be always described by the exchange of a massless anomaly pole, as in the case of vanishing chemical potentials. The pole is interpreted as an interpolating axionlike quasiparticle generated by the anomaly. In each axial-vector channel, it is predicted to be a correlated fermion/antifermion pseudoscalar (axionlike) quasiparticle appearing in the response function, once the material is subjected to an external chiral perturbation. The cancellation of the μ dependence extends to any chiral current within the Standard Model, including examples such as B (baryon), L (lepton), and B−L. This holds true irrespective of whether these currents exhibit anomalies. Published by the American Physical Society 2024

Robust Effective Ground State in a Nonintegrable Floquet Quantum Circuit.

An external periodic (Floquet) drive is believed to bring any initial state to the featureless infinite temperature state in generic nonintegrable isolated quantum many-body systems in the thermodynamic limit, irrespective of the driving frequency Ω. However, numerical or analytical evidence either proving or disproving this hypothesis is very limited and the issue has remained unsettled. Here, we study the initial state dependence of Floquet heating in a nonintegrable kicked Ising chain of length up to L=30 with an efficient quantum circuit simulator, showing a possible counterexample: the ground state of the effective Floquet Hamiltonian is exceptionally robust against heating, and could stay at finite energy density even after infinitely many Floquet cycles, if the driving period is shorter than a threshold value. This sharp energy localization transition or crossover does not happen for generic excited states. The exceptional robustness of the ground state is interpreted by (i)its isolation in the energy spectrum and (ii)the fact that those states with L-independent ℏΩ energy above the ground state energy of any generic local Hamiltonian, like the approximate Floquet Hamiltonian, are atypical and viewed as a collection of noninteracting quasiparticles. Our finding paves the way for engineering Floquet protocols with finite driving periods realizing long-lived, or possibly even perpetual, Floquet phases by initial state design.

Gravitational chiral anomaly at finite temperature and density

We investigate the gravitational anomaly vertex ⟨TTJ5⟩ (graviton—graviton—axial current) under conditions of finite density and temperature. Through a direct analysis of perturbative contributions, we demonstrate that neither finite temperature nor finite fermion density affects the gravitational chiral anomaly. These results find application in several contexts, from topological materials to the early Universe plasma. They affect the decay of any axion or axionlike particle into gravitational waves, in very dense and hot environments. Published by the American Physical Society 2024

Anatomy of critical fluctuations in hadronic matter

Critical phenomena in phase transitions of strongly interacting matter, governed by quantum chromodynamics, are inherently encoded in the fluctuations of conserved charges. In this work, we study the net-baryon number density fluctuations, including the lowest-lying nucleon and the baryonic resonance Δ(1232), based on the parity doublet model in the mean-field approximation. We focus on the qualitative features of the second-order susceptibility of the net-baryon number density in dense hadronic matter and how the inclusion of Δ(1232) affects it. We demonstrate that the fluctuations of the individual baryons do not necessarily reflect the total net-baryon number fluctuations at finite density, due to the nontrivial correlations between different particle species. Our results highlight the role of baryonic correlations in the interpretation of data from heavy ion collision experiments. Published by the American Physical Society 2024

Dense QCD2 with matrix product states

We study one-flavor SU(2) and SU(3) lattice QCD in (1 + 1) dimensions at zero temperature and finite density using matrix product states and the density matrix renormalization group. We compute physical observables such as the equation of state, chiral condensate, and quark distribution function as functions of the baryon number density. As a physical implication, we discuss the inhomogeneous phase at nonzero baryon density, where the chiral condensate is inhomogeneous, and baryons form a crystal. We also discuss how the dynamical degrees of freedom change from hadrons to quarks through the formation of quark Fermi seas.

Flexible Bayesian estimation of incubation times.

The incubation period is of paramount importance in infectious disease epidemiology as it informs about the transmission potential of a pathogenic organism and helps to plan public health strategies to keep an epidemic outbreak under control. Estimation of the incubation period distribution from reported exposure times and symptom onset times is challenging as the underlying data is coarse. We develop a new Bayesian methodology using Laplacian-P-splines that provides a semi-parametric estimation of the incubation density based on a Langevinized Gibbs sampler. A finite mixture density smoother informs a set of parametric distributions via moment matching and an information criterion arbitrates between competing candidates. Algorithms underlying our method find a natural nest within the EpiLPS package, which has been extended to cover estimation of incubation times. Various simulation scenarios accounting for different levels of data coarseness are considered with encouraging results. Applications to real data on COVID-19, MERS and Mpox reveal results that are in alignment with what has been obtained in recent studies. The proposed flexible approach is an interesting alternative to classic Bayesian parametric methods for estimation of the incubation distribution.

Matrix product state approximations to quantum states of low energy variance

We show how to efficiently simulate pure quantum states in one dimensional systems that have both finite energy density and vanishingly small energy fluctuations. We do so by studying the performance of a tensor network algorithm that produces matrix product states whose energy variance decreases as the bond dimension increases. Our results imply that variances as small as ∝1/log⁡N can be achieved with polynomial bond dimension. With this, we prove that there exist states with a very narrow support in the bulk of the spectrum that still have moderate entanglement entropy, in contrast with typical eigenstates that display a volume law. Our main technical tool is the Berry-Esseen theorem for spin systems, a strengthening of the central limit theorem for the energy distribution of product states. We also give a simpler proof of that theorem, together with slight improvements in the error scaling, which should be of independent interest.

Polarization dynamics in the paramagnet of a charged quark-gluon plasma

It is commonly understood that the strong magnetic field produced in heavy ion collisions is short-lived. The electric conductivity of the quark-gluon plasma is unable to significantly extend the lifetime of the magnetic field. We propose an alternative scenario to achieve this: with finite baryon density and spin polarization by the initial magnetic field, the quark-gluon plasma behaves as a paramagnet, which may continue to polarize the quark after fading of the initial magnetic field. We confirm this picture by calculations in both quantum electrodynamics and quantum chromodynamics. In the former case, we find a splitting in the damping rates of probe fermions with an opposite spin component along the magnetic field. In the latter case, we find a similar splitting in damping rate of the probe quark in quark-gluon plasma in both high and low density limits. The splitting provides a way of polarizing strange quarks by the quark-gluon plasma paramagnet consisting of light quarks, which effectively extends the lifetime of the magnetic field in heavy ion collisions. Published by the American Physical Society 2024

Numerical modelling of sawteeth and sawtooth-free regime

To better understand the sawteeth physics and the sawtooth-free regime associated with the hybrid scenario in tokamak experiments, numerical calculations up to quasi-steady state have been carried out for realistic middle-size tokamak plasma parameters, including the bootstrap current perturbation and basing on both the single- and two-fluid equations with the large aspect ratio approximation. Two types of the sawtooth crash are found in multiple sawteeth simulations: (1) For a low equilibrium bootstrap current fraction, the crash is caused by the internal kink mode, as expected; (2) When the bootstrap current density fraction is larger than 10% in the core region, however, the crash is caused by the non-ideal double kink mode, in contrary to the conventional understanding. In this case, a non-monotonic radial profile of the safety factor q with two q=1 surfaces emerges before the crash, caused by the bootstrap current density and plasma resistivity perturbations, although the original equilibrium has only single q = 1 surface. In both types of sawtooth crashes, the crash time in two-fluid simulations is tens of microseconds, as observed in experiments. Furthermore, for a relatively low ion density and finite bootstrap current density fraction, a transition from the sawtooth to the sawtooth-free regime is found, in which flat q profiles with the q value being about unity in the central region, similar to that observed in hybrid scenario experiments, are maintained by the dynamo effect. To enter into the sawtooth-free regime in two-fluid simulations, a much larger Alfvén velocity than that in single-fluid simulations is required due to the diamagnetic drift.

First order phase transition in the D3-D7 model from the point of view of the fermionic spectral functions

We consider the D3-D7 model and use the spectral function of a probe fermion on D7 to analyze the first order phase transition from the black-hole embedding phase to another black-hole embedding phase in the presence of the finite density and temperature. From the fermionic spectral functions, we study the temperature dependence of the decay rate, and we observe various phenomena that support the first order phase transition including jump in it at the critical temperature that corresponds to the first order phase transition. Published by the American Physical Society 2024

Nonlinear electrodynamics for the vacuum of Dirac materials: Photon magnetic properties and radiation pressures

We investigate the magnetic properties of photons propagating through Dirac materials in a magnetic field, considering both vacuum and medium contributions. Photon propagation properties are obtained through a second-order expansion of nonlinear Euler-Heisenberg electrodynamics at finite density and temperature considering Dirac material parameters (Dirac fine structure constant, band gap, and Fermi velocity). Total magnetization (including electron and photon contributions) and photon-effective magnetic moment are computed. Observables such as photon energy density, radiation pressure, and Poynting vector are obtained by an average of components of the energy-momentum tensor. All quantities are expressed in terms of Lagrangian derivatives. Those related to the vacuum are valid for any value of the external magnetic field, and both the weak- and strong-field limits are recovered. We discuss some ideas of experiments that may contribute to testing in Dirac materials the phenomenology of the strong magnetic field in the quantum electrodynamic (QED) vacuum and how nonlinear corrections on the magnetization, radiation pressure, and birefringence are amplified up to 103 times QED corrections. Published by the American Physical Society 2024

Hysteresis between gas breakdown and plasma discharge

In direct-current (DC) discharge, it is well known that hysteresis is observed between the Townsend (gas breakdown) and glow regimes. Forward and backward voltage sweep is performed using a one-dimensional particle-in-cell Monte Carlo collision (PIC-MCC) model considering a ballast resistor. When increasing the applied voltage after reaching the breakdown voltage (Vb), transition from Townsend to glow discharges is observed. When decreasing the applied voltage from the glow regime, the discharge voltage (Vd) between the anode–cathode gap can be smaller than the breakdown voltage, resulting in a hysteresis, which is consistent with experimental observations. Next, the PIC-MCC model is used to investigate the self-sustaining voltage (Vs) in the presence of finite initial plasma densities between the anode and cathode gap. It is observed that the self-sustaining voltage coincides with the discharge voltage obtained from the backward voltage sweep. In addition, the self-sustaining voltage decreases with increased initial plasma density and saturates above a certain initial plasma density, which indicates a change in plasma resistivity. The decrease in self-sustaining voltage is associated with the electron heat loss at the anode for the low pd (rarefied) regime. In the high pd (collisional) regime, the ion energy loss toward the cathode due to the cathode fall and the inelastic collision loss of electrons in the bulk discharge balance out. Finally, it is demonstrated that the self-sustaining voltage collapses to a singular value, despite the presence of a initial plasma, for microgaps when field emission is dominant, which is also consistent with experimental observations.

Hawking-Page transition in holographic QCD at finite density

We study the confinement/deconfinement transition of QCD matter for a quark-gluon plasma at finite density using AdS/QCD duality. In order to determine the critical temperature and its dependence on the quark chemical potential, the semi-classical Hawking-Page approach is considered for a charged AdS black hole. The result obtained is consistent with the QCD phase diagram, with the critical temperature decreasing as the chemical potential increases. Using the soft wall holographic model, the value of quark density in which the phase transition occurs at zero temperature is estimated. For higher densities the QCD matter is deconfined, independent of the temperature.

Coupling Nano and Atomic Electric Field Confinement for Robust Alkaline Oxygen Evolution

AbstractThe alkaline oxygen evolution reaction (OER) is a promising avenue for producing clean fuels and storing intermittent energy. However, challenges such as excessive OH− consumption and strong adsorption of oxygen‐containing intermediates hinder the development of alkaline OER. In this study, we propose a cooperative strategy by leveraging both nano‐scale and atomically local electric fields for alkaline OER, demonstrated through the synthesis of Mn single atom doped CoP nanoneedles (Mn SA‐CoP NNs). Finite element method simulations and density functional theory calculations predict that the nano‐scale local electric field enriches OH− around the catalyst surface, while the atomically local electric field improves *O desorption. Experimental validation using in situ attenuated total reflection infrared and Raman spectroscopy confirms the effectiveness of the nano‐scale and atomically electric fields. Mn SA‐CoP NNs exhibit an ultra‐low overpotential of 189 mV at 10 mA cm−2 and stable operation over 100 hours at ~100 mA cm−2 during alkaline OER. This innovative strategy provides new insights for enhancing catalyst performance in energy conversion reactions.

Coupling Nano and Atomic Electric Field Confinement for Robust Alkaline Oxygen Evolution.

The alkaline oxygen evolution reaction (OER) is a promising avenue for producing clean fuels and storing intermittent energy. However, challenges such as excessive OH- consumption and strong adsorption of oxygen-containing intermediates hinder the development of alkaline OER. In this study, we propose a cooperative strategy by leveraging both nano-scale and atomically local electric fields for alkaline OER, demonstrated through the synthesis of Mn single atom doped CoP nanoneedles (Mn SA-CoP NNs). Finite element method simulations and density functional theory calculations predict that the nano-scale local electric field enriches OH- around the catalyst surface, while the atomically local electric field improves *O desorption. Experimental validation using in situ attenuated total reflection infrared and Raman spectroscopy confirms the effectiveness of the nano-scale and atomically electric fields. Mn SA-CoP NNs exhibit an ultra-low overpotential of 189 mV at 10 mA cm-2 and stable operation over 100 hours at ~100 mA cm-2 during alkaline OER. This innovative strategy provides new insights for enhancing catalyst performance in energy conversion reactions.

Probing molecular spectral functions and unconventional pairing using Raman spectroscopy

An impurity interacting with an ultracold Fermi gas can form either a polaron state or a dressed molecular state, the molaron, in which the impurity forms a bound state with one gas particle. This molaron state features rich physics, including a negative effective mass around unitarity and a first-order transition to the polaron state. However, these features have remained so far experimentally inaccessible. In this work we show theoretically how the molaron state can be directly prepared experimentally even in its excited states using Raman spectroscopy techniques. Initializing the system in the ultrastrong coupling limit, where the binding energy of the molaron is much larger than the Fermi energy, our protocol maps out the momentum-dependent spectral function of the molecule. Using a diagrammatic approach we furthermore show that the molecular spectral function serves as a direct precursor of the elusive Fulde-Ferell-Larkin-Ovchinnikov phase, which is realized for a finite density of fermionic impurity particles. Our results pave the way to a systematic understanding of how composite particles form in quantum many-body environments and provide a basis to develop new schemes for the observation of exotic phases of quantum many-body systems. Published by the American Physical Society 2024

Surface Magnetization in Antiferromagnets: Classification, Example Materials, and Relation to Magnetoelectric Responses

We use symmetry analysis and density-functional theory to determine and characterize surface terminations that have a finite equilibrium magnetization density in antiferromagnetic materials. A nonzero magnetic dipole moment per unit area or “surface magnetization” can arise on particular surfaces of many antiferromagnets due to the bulk magnetic symmetries. Such surface magnetization underlies intriguing physical phenomena like interfacial magnetic coupling and can be used as a readout method of antiferromagnetic domains. However, a universal description of antiferromagnetic surface magnetization is lacking. We first introduce a classification system based on whether the surface magnetization is either sensitive or robust to roughness and on whether the magnetic dipoles at surface of interest are compensated or uncompensated when the bulk magnetic order is retained at the surface. We show that roughness-sensitive categories can be identified by a simple extension of a previously established group-theory formalism for identifying roughness-robust surface magnetization. We then map the group-theory method of identifying surface magnetization to a novel description in terms of bulk magnetic multipoles, which are already established as symmetry indicators for bulk magnetoelectric responses at both linear and higher orders. We use density-functional calculations to illustrate that nominally compensated surfaces in magnetoelectric Cr2O3 and centrosymmetric altermagnetic FeF2 develop a finite magnetization density at the surface, in agreement with our predictions based on both group theory and the ordering of the bulk multipoles. Our analysis provides a comprehensive basis for understanding the surface magnetic properties and their intimate correspondence to bulk magnetoelectric effects in antiferromagnets and has important implications for technologically relevant phenomena such as exchange-bias coupling. Published by the American Physical Society 2024

Fisher and Bayes-Fisher information measures for finite mixture distributions

. In this work, we consider Fisher information and Bayes-Fisher information measures for the mixing parameter vector of a finite mixture density function and develop some associated results for this model. We provide several interesting connections between these measures and some known informational measures, such as chi-square divergence, Shannon entropy, Kullback-Leibler, Jeffreys, and Jensen-Shannon divergences. Finally, to demonstrate the usefulness of the Bayes-Fisher information measure, we apply it to a real example in image processing and present some numerical results. Our findings in this regard show that this information measure is an effective criteria for quantifying the similarity between two images.

Depolarization Field-Induced Photovoltaic Effect in Graphene/α-In2Se3/Graphene Heterostructures.

The ferroelectric photovoltaic effect (FPVE) enables alternate pathways for energy conversion that are not allowed in centrosymmetric materials. Understanding the dominant mechanism of the FPVE at the ultrathin limit is important for defining the ultimate efficiency. In contrast to the wide band gap conventional thin-film ferroelectrics, 2D α-In2Se3 has an ideal band gap of 1.3 eV and enables the fabrication of ultrathin and stable heterostructures, providing the perfect platform to explore FPVE in the nanoscale limit. Here, we study the ferroelectric layer thickness-dependent FPVE in vertical few-layer graphene/α-In2Se3/graphene heterostructures. We find that the short-circuit photocurrent is antiparallel to the ferroelectric polarization and increases exponentially with decreasing thickness. We show that the observed behavior is predicted by the depolarization field model, originating from the unscreened bound charges due to the finite density of states in semimetal few-layer graphene. As a result, the heterostructures show enhancement of the power conversion efficiency, reaching 2.56 × 10-3% under 100 W/cm2 in 18 nm thick α-In2Se3, approximately 275 times more than the 50 nm thick α-In2Se3. These results demonstrate the importance of the depolarization field at the nanoscale and define design principles for the potential of harnessing FPVE at reduced dimension.