Charge Susceptibility Research Articles

Magnetic field effect on hadron yield ratios and fluctuations in a hadron resonance gas

This work studies the influence of an external magnetic field on hadron yields and fluctuations in a hadron resonance gas by performing calculations within an updated version of the open-source package. The presence of the magnetic field has a sizable influence on certain hadron yield ratios. Most notably, it leads to enhanced p/π and suppressed n/p ratios, which may serve as a magnetometer in heavy-ion collisions. By attributing the centrality dependence of the p/π ratio in Pb-Pb collisions at 5.02 TeV measured by the ALICE Collaboration entirely to the magnetic field, its maximal strength at freezeout is estimated to be eB≃0.12GeV2≃6.3mπ2 in peripheral collisions. The magnetic field also enhances various conserved charge susceptibilities, which is qualitatively consistent with recent lattice QCD data and is driven in the HRG model by the increase of hadron densities. However, the variances of hadrons do not show any enhancement when normalized by the means, therefore, measurements of second-order fluctuations in heavy-ion collisions appear to offer limited additional information about the magnetic field relative to mean multiplicities. Published by the American Physical Society 2024

Internal consistency of multi-tier GW+EDMFT

The multi-tier GW+EDMFT scheme is an ab-initio method for calculating the electronic structure of correlated materials. While the approach is free from ad-hoc parameters, it requires a selection of appropriate energy windows for describing low-energy and strongly correlated physics. In this study, we test the consistency of the multi-tier description by considering different low-energy windows for a series of cubic SrXO3 (X = V, Cr, Mn) perovskites. Specifically, we compare the 3-orbital t2g model, the 5-orbital t2g + eg model, the 12-orbital t2g + Op model, and (in the case of SrVO3) the 14-orbital t2g + eg + Op model and compare the results to available photoemission and X-ray absorption measurements. The multi-tier method yields consistent results for the t2g and t2g + eg low-energy windows, while the models with Op states produce stronger correlation effects and mostly agree well with experiment, especially in the unoccupied part of the spectrum. We also discuss the consistency between the fermionic and bosonic spectral functions and the physical origin of satellite features, and present momentum-resolved charge susceptibilities.

Thermodynamic Stability at the Two-Particle Level.

We show how the stability conditions for a system of interacting fermions that conventionally involve variations of thermodynamic potentials can be rewritten in terms of one- and two-particle correlators. We illustrate the applicability of this alternative formulation in a multiorbital model of strongly correlated electrons at finite temperatures, inspecting the lowest eigenvalues of the generalized local charge susceptibility in proximity of the phase-separation region. Additionally to the conventional unstable branches, we address unstable solutions possessing a positive, rather than negative, compressibility. Our stability conditions require no derivative of free-energy functions with conceptual and practical advantages for actual calculations and offer a clear-cut criterion for analyzing the thermodynamics of correlated complex systems.

Gauge field fluctuation corrected QED3 effective action by fermionic particle-vortex duality

We present a non-perturbative framework for incorporating gauge field fluctuations into effective actions of QED3 in the infrared using fermionic particle-vortex duality. This approach is demonstrated through the applications to models containing N species of 2-component Dirac fermions in solvable and interpretable electromagnetic backgrounds, focusing on N = 1 or 2. For the N = 1 model, we establish a correspondence between fermion Casimir energy at finite density and the magnetic Euler-Heisenberg Lagrangian, and further evaluate the corrections to their amplitudes. This predicts amplification of charge susceptibility and reduction of magnetic permeability. We additionally provide physical interpretations for each component of our calculation and offer alternative derivations based on energy density measurements in different characteristic lengths. For N = 2, we show that magnetic catalysis is erased in a U(1) × U(1) QED3, indicating no breakdown of chiral symmetry. Reasoning is offered based on the properties of the lowest Landau level wave functions.

Denoising and Extension of Response Functions in the Time Domain.

Response functions of quantum systems, such as electron Green's functions, magnetic, or charge susceptibilities, describe the response of a system to an external perturbation. They are the central objects of interest in field theories and quantum computing and measured directly in experiment. Response functions are intrinsically causal. In equilibrium and steady-state systems, they correspond to a positive spectral function in the frequency domain. Since response functions define an inner product on a Hilbert space and thereby induce a positive definite function, the properties of this function can be used to reduce noise in measured data and, in equilibrium and steady state, to construct positive definite extensions for data known on finite time intervals, which are then guaranteed to correspond to positive spectra.

Non-perturbative intertwining between spin and charge correlations: A ``smoking gun'' single-boson-exchange result

We study the microscopic mechanism controlling the interplay between local charge and local spin fluctuations in correlated electron systems via a thorough investigation of the generalized on-site charge susceptibility of several fundamental many-electron models, such as the Hubbard atom, the Anderson impurity model, and the Hubbard model. By decomposing the numerically determined generalized susceptibility in terms of physically transparent single-boson exchange processes, we unveil the microscopic mechanisms responsible for the breakdown of the self-consistent many-electron perturbation expansion. In particular, we unambiguously identify the origin of the significant suppression of its diagonal entries in (Matsubara) frequency space and the slight increase of the off-diagonal ones which cause the breakdown. The suppression effect on the diagonal elements originates directly from the electronic scattering on local magnetic moments, reflecting their increasingly longer lifetime as well as their enhanced effective coupling with the electrons. Instead, the slight and diffuse enhancement of the off-diagonal terms can be mostly ascribed to multiboson scattering processes. The strong intertwining between spin and charge sectors is partly weakened at the Kondo temperature due to a progressive reduction of the effective spin-fermion coupling of local magnetic fluctuations in the low frequency regime. Our analysis, thus, clarifies the precise mechanism through which the physical information is transferred between different scattering channels of interacting electron problems and highlights the pivotal role played by such an intertwining in the physics of correlated electrons beyond the perturbative regime.

Dynamical evolution of spinodal decomposition in holographic superfluids

We study the nonlinear dynamical evolution of spinodal decomposition in a first-order superfluid phase transition using a simple holographic model in the probe limit. We first confirm the linear stability analysis based on quasinormal modes and verify the existence of a critical length scale related to a gradient instability — negative speed of sound squared — of the superfluid sound mode, which is a consequence of a negative thermodynamic charge susceptibility. We present a comparison between our case and the standard Cahn-Hilliard equation for spinodal instability, in which a critical length scale can be also derived based on a diffusive instability. We then perform several numerical tests which include the nonlinear time evolution directly from an unstable state and fast quenches from a stable to an unstable state in the spinodal region. Our numerical results provide a real time description of spinodal decomposition and phase separation in one and two spatial dimensions. We reveal the existence of four different stages in the dynamical evolution, and characterize their main properties. Finally, we investigate the strength of dynamical heterogeneity using the spatial variance of the local chemical potential and we correlate the latter to other features of the dynamical evolution.

Optical conductivity and damping of plasmons due to electron-electron interaction

We revisit the issue of plasmon damping due to electron-electron interaction. The plasmon linewidth can be related to the imaginary part of the charge susceptibility or, equivalently, to the real part of the optical conductivity, Reσ(q,ω). Approaching the problem first via a standard semiclassical Boltzmann equation, we show that Reσ(q,ω) of a two-dimensional (2D) electron gas scales as q2T2/ω4 for ω≪T, which agrees with the results of Principi [] and Sharma [] but disagrees with that of Mishchenko [], according to which Reσ(q,ω)∝q2T2/ω2. To resolve this disagreement, we rederive Reσ(q,ω) using the original method of Mishchenko for an arbitrary ratio ω/T and show that while the last term is, indeed, present, it is subleading to the q2T2/ω4 term. We give a physical interpretation of both leading and subleading contributions in terms of the shear and bulk viscosities of an electron liquid, respectively. We also calculate Reσ(q,ω) for a three-dimensional electron gas and doped monolayer graphene. We find that, all other parameters being equal, finite temperature has the strongest effect on the plasmon linewidth in graphene, where it scales as T4lnT for ω≪T. Published by the American Physical Society 2024

Effect of finite volume on thermodynamics of quark-hadron matter

The effects of a finite system volume on thermodynamic quantities, such as the pressure, energy density, specific heat, speed of sound, conserved charge susceptibilities, and correlations, in hot and dense strongly interacting matter are studied within the parity-doublet chiral mean field model. Such an investigation is motivated by relativistic heavy-ion collisions, which create a blob of hot QCD matter of a finite volume, consisting of strongly interacting hadrons and potentially deconfined quarks and gluons. The effect of the finite volume of the system is incorporated by introducing lower momentum cutoffs in the momentum integrals appearing in the model, the numerical value of the momentum cutoff being related to the de Broglie wavelength of the given particle species. It is found that some of these quantities show a significant volume dependence, in particular, those sensitive to pion degrees of freedom, and the crossover transition is generally observed to become smoother in finite volume. These findings are relevant for the effective equation of state used in fluid dynamical simulations of heavy-ion collisions and efforts to extract the freeze-out properties of heavy-ion collisions with susceptibilities involving electric charge and strangeness. Published by the American Physical Society 2024

Finite‐Frequency Noise and Dynamical Charge Susceptibility in Single and Double Quantum Dot Systems

AbstractThis study reports on finite‐frequency noise and dynamical charge susceptibility in out‐of‐equilibrium quantum dot systems. Both single and double quantum dots connected to one or two reservoirs of electrons are considered, and these quantities are calculated by using the non‐equilibrium Green function technique. The results are discussed in the light of experimental results, particularly in the low‐frequency limit for which an interpretation in terms of an equivalent RC‐circuit is made. Anti‐symmetrized noise is also studied, defined as the difference between absorption and emission noises, and its relationship with the dynamical charge susceptibility in single quantum dots is established. In double quantum dots, the similarities between the dynamical charge susceptibility, the absorption noise, and the dot occupancy, are highlighted by comparing their respective variations with the bias voltage applied between the two reservoirs, and the detuning energy defined as the difference between the lowest level energies in the two dots.

ACFlow: An open source toolkit for analytic continuation of quantum Monte Carlo data

The purpose of analytic continuation is to establish a real frequency spectral representation of single-particle or two-particle correlation function (such as Green's function, self-energy function, spin and charge susceptibilities) from noisy data generated in finite temperature quantum Monte Carlo simulations. It requires numerical solutions of a family of Fredholm integral equations of the first kind, which is indeed a challenging task. In this paper, an open source toolkit (dubbed ACFlow) for analytic continuation of quantum Monte Carlo data is presented. We first give a short introduction to the analytic continuation problem. Next, three popular analytic continuation algorithms, including the maximum entropy method, the stochastic analytic continuation, and the stochastic optimization method, as implemented in this toolkit are reviewed. And then we elaborate on the major features, implementation details, basic usage, inputs, and outputs of this toolkit. Finally, four representative examples, including analytic continuations of Matsubara self-energy function, Matsubara and imaginary time Green's functions, and current-current correlation function, are shown to demonstrate the usefulness and flexibility of the ACFlow toolkit. Program summaryProgram Title: ACFlowCPC Library link to program files:https://doi.org/10.17632/th6w74gwjm.1Developer's repository link:https://github.com/huangli712/ACFlowLicensing provisions: GNU General Public License Version 3Programming language: JuliaNature of problem: Most of the quantum Monte Carlo methods work on imaginary axis. In order to extract physical observables and compare them with the experimental results, analytic continuation must be done in the post-processing stage to convert the quantum Monte Carlo simulated data from imaginary axis to real axis.Solution method: Three well-established analytic continuation methods, including the maximum entropy method, the stochastic analytic continuation (both A. W. Sandvik's and K. S. D. Beach's algorithms), and the stochastic optimization method, have been implemented in the ACFlow toolkit.Additional comments including restrictions and unusual features: The ACFlow toolkit is written in pure Julia language. It is highly optimized and parallelized. It can be executed interactively in a Jupyter notebook environment.

Infinite-layer nickelates as Ni-eg Hund’s metals

The recent and exciting discovery of superconductivity in the hole-doped infinite-layer nickelate Nd1−δSrδNiO2 draws strong attention to correlated quantum materials. From a theoretical view point, this class of unconventional superconducting materials provides an opportunity to unveil a physics hidden in correlated quantum materials. Here we study the temperature and doping dependence of the local spectrum as well as the charge, spin and orbital susceptibilities from first principles. By using ab initio LQSGW+DMFT methodology, we show that onsite Hund’s coupling in Ni-d orbitals gives rise to multiple signatures of Hund’s metallic phase in Ni-eg orbitals. The proposed picture of the nickelates as an eg (two orbital) Hund’s metal differs from the picture of the Fe-based superconductors as a five orbital Hund’s metal as well as the picture of the cuprates as doped charge transfer insulators. Our finding uncover a new class of the Hund’s metals and has potential implications for the broad range of correlated two orbital systems away from half-filling.

Higher-form symmetries, anomalous magnetohydrodynamics, and holography

In U(1)U(1) Abelian gauge theory coupled to fermions, the non-conservation of the axial current due to the chiral anomaly is given by a dynamical operator F_{\mu\nu} \tilde{F}^{\mu\nu}FμνF̃μν constructed from the field-strength tensor. We attempt to describe this physics in a universal manner by casting this operator in terms of the 2-form current for the 1-form symmetry associated with magnetic flux conservation. We construct a holographic dual with this symmetry breaking pattern and study some aspects of finite temperature anomalous magnetohydrodynamics. We explicitly calculate the charge susceptibility and the axial charge relaxation rate as a function of temperature and magnetic field and compare to recent lattice results. At small magnetic fields we find agreement with elementary hydrodynamics weakly coupled to an electrodynamic sector, but we find deviations at larger fields.

Pairing susceptibility of the two-dimensional Hubbard model in the thermodynamic limit

We compute the diagrammatic expansion of the particle-particle susceptibility via algorithmic Matsubara integration and compute the correlated pairing susceptibility in the thermodynamic limit of the 2D Hubbard Model. We study the static susceptibility and its dependence on the pair momentum $\mathbf{q}$ for a range of temperature, interaction strength, and chemical potential. We show that $d_{x^2-y^2}$-wave pairing is expected in the model in the $U/t\to 0^+ $ limit from direct perturbation theory. From this, we identify key second and third-order diagrams that support pairing processes and note that the diagrams responsible are not a part of charge or spin susceptibility expansions. We find two key components for pairing at momenta $(0,0)$ and $(\pi,\pi)$ that can be well fit as separate bosonic modes. We extract amplitudes and correlation length scales where we find a predominantly local $(\pi,\pi)$ pairing and non-local $\mathbf{q}=(0,0)$ pairs and present the relative weights of these modes for variation in temperature, doping, and interaction strength.

Two-dimensional extended Hubbard model at half-filling

We consider the extended Hubbard model on a two-dimensional square lattice at half-filling. The model is investigated using the strong coupling diagram technique. We sum infinite series of ladder diagrams allowing for full-scale charge and spin fluctuations and the actual short-range antiferromagnetic order for nonzero temperatures. In agreement with earlier results, we find the first-order phase transition in the charge subsystem occurring at v = v c ≳ U/4 with v and U the intersite and on-site Coulomb repulsion constants. The transition reveals itself in an abrupt sign change of a sharp maximum in the zero-frequency charge susceptibility at the corner of the Brillouin. States arising at the transition have alternating deviations of electron occupations from the mean value on neighboring sites. Due to fluctuations, these alternating occupation deviations have short-range order. For the considered parameters, such behavior is found for U ≲ 5t with t the hopping constant. For the insulating case U ≳ 6t, in which the transition is not observed, we find a continuous growth of the Mott gap with v. The evolution of the electron density of states with increasing v is also considered.

Evidence of charge susceptibility and multiple f–c hybridization configurations with the La doping in CeGe: a DFT + DMFT study

Kondo coupling has been extensively investigated in several Ce-based systems. However, the search for materials showing the interplay between the Kondo effect, spin–orbit interaction, and crystal-field effect along with the presence of local charge susceptibility; remains a challenge for the condensed matter community. Actually, in Ce-based systems, the strong coupling of the conduction electrons to the local magnetic moments usually hides these properties. Here, we present a detailed investigation of Ce0.6La0.4Ge through a combined density functional theory and dynamic mean-field theory study. Our investigations give evidence of the significant charge susceptibility and the multiple different f–c hybridization configurations. The weakening of the magnetization owing to the dilution of the Ce-site is the main cause for the appearance of such properties, which is believed to occur due to the presence of the relevant local moment and f–c hybridization over the competition with the on-site Coulomb interaction.

Dynamical stability from quasi normal modes in 2nd, 1st and 0th order holographic superfluid phase transitions

We study a simple extension of the original Hartnoll, Herzog and Horowitz (HHH) holographic superfluid model with two nonlinear scalar self-interaction terms λ|ψ|4 and τ|ψ|6 in the probe limit. Depending on the value of λ and τ, this setup allows us to realize a large spectrum of holographic phase transitions which are 2nd, 1st and 0th order as well as the “cave of wind” phase transition. We speculate the landscape pictures and explore the near equilibrium dynamics of the lowest quasinormal modes (QNMs) across the whole phase diagram at both zero and finite wave-vector. We find that the zero wave-vector results of QNMs correctly present the stability of the system under homogeneous perturbations and perfectly agree with the landscape analysis of homogeneous configurations in canonical ensemble. The zero wave-vector results also show that a 0th order phase transition cannot occur since it always corresponds to a global instability of the whole system. The finite wave-vector results show that under inhomogeneous perturbations, the unstable region is larger than that under only homogeneous perturbations, and the new boundary of instability match with the turning point of condensate curve in grand canonical ensemble, indicating a new explanation from the subsystem point of view. The additional unstable section also perfectly match the section with negative value of charge susceptibility.

Electron-Hole Asymmetry of Quantum Collective Excitations in High-Tc Copper Oxides

We carry out a systematic study of collective spin- and charge excitations for the canonical single-band Hubbard, $t$-$J$-$U$, and $t$-$J$ models of high-temperature copper-oxide superconductors, both on electron- and hole-doped side of the phase diagram. Recently developed variational wave function approach, combined with the expansion in inverse number of fermionic flavors, is employed. All three models exhibit a substantial electron-hole asymmetry of magnetic excitations, with a robust paramagnon emerging for hole-doping, in agreement with available resonant inelastic $x$-ray scattering data for the cuprates. The $t$-$J$ model yields additional high-energy peak in the magnetic spectrum that is not unambiguously identified in spectroscopy. For all considered Hamiltonians, the dynamical charge susceptibility contains a coherent mode for both hole- and electron doping, with overall bandwidth renormalization controlled by the on-site Coulomb repulsion. Away from the strong-coupling limit, the antiferromagnetic ordering tendency is more pronounced on electron-doped side of the phase diagram.

Constraints on hadron resonance gas interactions via first-principles Lattice QCD susceptibilities

We investigate extensions of the Hadron Resonance Gas (HRG) Model beyond the ideal case by incorporating both attractive and repulsive interactions into the model. When considering additional states exceeding those measured with high confidence by the Particle Data Group, attractive corrections to the overall pressure in the HRG model are imposed. On the other hand, we also apply excluded-volume corrections, which ensure there is no overlap of baryons by turning on repulsive (anti)baryon-(anti)baryon interactions. We emphasize the complementary nature of these two extensions and identify combinations of conserved charge susceptibilities that allow us to constrain them separately. In particular, we find interesting ratios of susceptibilities that are sensitive to one correction and not the other. This allows us to constrain the excluded volume and particle spectrum effects separately. Analysis of the available lattice results suggests the presence of both the extra states in the baryonstrangeness sector and the repulsive baryonic interaction, with indications that hyperons have a smaller repulsive core than non-strange baryons. We note that these results are interesting for heavy-ion-collision systems at both the LHC and RHIC.

Electrical conductivity of the quark-gluon plasma from the low energy limit of photon and dilepton spectra

Fluid dynamic considerations are used to determine the electric current spectral density in the regime of small energies and momenta. The spectral density in this regime is parameterized by the electric conductivity, the charge susceptibility, and the relaxation time for the electric current, which is needed for relativistic causality. Experimentally, the spectral function can be accessed through the production rates of photons and dileptons in the expanding quark-gluon plasma. We use fluid dynamic simulations of high energy nuclear collisions, together with the transport limit of the spectral density, to obtain photon and dielectron spectra for different values of the conductivity and relaxation times. The yields of photon and dileptons produced in the plasma are compared to the background from decays of short-lived hadrons. We discuss how experiments can constrain the electrical conductivity and associated relaxation time of the quark-gluon plasma.