Correlations and fluctuations in a magnetized three-flavor PNJL model with and without inverse magnetic catalysis effect

Abstract

The correlations $\chi^{BQ}_{11},\ \chi^{BS}_{11},\ \chi^{QS}_{11}$ and quadratic (quartic) fluctuations $\chi^{B}_{2,4},\ \chi^{Q}_{2,4},\ \chi^{S}_{2,4}$ of baryon number $B$, electric charge $Q$ and strangeness $S$ are investigated in a three-flavor Polyakov loop extended Nambu-Jona-Lasinio model at finite temperature and magnetic field. The inverse magnetic catalysis (IMC) effect is introduced through the magnetic field dependent parameters $G(eB)$ or $T_{0}(eB)$, and we make comparison of the results in the cases with and without IMC effect. Since including IMC effect does not change the strength of phase transition under external magnetic field, it does not lead to qualitative difference in the correlations and fluctuations, but modifies their values. Under vanishing and nonvanishing magnetic field, the correlations and fluctuations increase with temperature, and then show the peak around the pseudocritical temperatures of chiral restoration and deconfinement phase transitions. The peak structure in $\chi^{BQ}_{11}$, $\chi^{B}_{4}$ and $\chi^{Q}_{4}$ are much more apparent than in other correlations and fluctuations. The correlations and fluctuations along the phase transition line under external magnetic field are characterized by the scaled correlations ${\hat{\chi}}_{11}^{BQ}=\frac{\chi_{11}^{BQ}(eB,T_{pc}^{c}(eB))}{\chi_{11}^{BQ}$, ${\hat{\chi}}_{11}^{BS}=\frac{\chi_{11}^{BS}(eB,T_{pc}^{c}(eB))}{\chi_{11}^{BS}$, ${\hat{\chi}}_{11}^{QS}=\frac{\chi_{11}^{QS}(eB,T_{pc}^{c}(eB))}{\chi_{11}^{QS}$ and scaled quadratic (quartic) fluctuations ${\hat{\chi}}_{2,4}^{B}=\frac{\chi_{2,4}^{B}(eB,T_{pc}^{c}(eB))}{\chi_{2,4}^{B}$, ${\hat{\chi}}_{2,4}^{Q}=\frac{\chi_{2,4}^{Q}(eB,T_{pc}^{c}(eB))}{\chi_{2,4}^{Q}$, ${\hat{\chi}}_{2,4}^{S}=\frac{\chi_{2,4}^{S}(eB,T_{pc}^{c}(eB))}{\chi_{2,4}^{S}$ at the pseudocritical temperature $T_{pc}^{c}$ of chiral restoration phase transition. They increase with magnetic fields due to the increase of phase transition strength under magnetic fields. Among them, ${\hat{\chi}}_{11}^{BQ}$ increases fastest, which may serve as the magnetometer of QCD.

I Introduction

Motivated by the strong magnetic field in the core of compact stars and in the initial stage of relativistic heavy ion collisions, the study on Quantum Chromodynamics (QCD) phase structure under external electromagnetic fields has attracted much attention review0 ; review1 ; review2 ; review3 ; review4 ; review5 ; lattice1 ; lattice2 ; lattice2sep ; lattice4 ; lattice5 ; lattice6 ; lattice7 ; lattice9 ; ding2025 ; lattice8 ; fukushima ; mao ; maosep1 ; maosep2 ; maosep3 ; maosep4 ; kamikado ; bf1 ; bf13 ; bf2 ; bf3 ; bf5 ; bf51 ; bf52 ; bf8 ; bf9 ; bf11 ; db1 ; db1sep ; db2 ; db3 ; db5 ; db6 ; pnjl1 ; pnjl1sep ; pnjl2 ; pnjl3 ; pnjl4 ; pqm ; ferr1 ; ferr2 ; mhuang ; meimao1 ; t0effect ; meihuangmao ; t0effectmao ; efield1 ; efield2 ; efield3 ; efield4 .

The Lattice QCD (LQCD) simulations lattice1 ; lattice2 ; lattice2sep ; lattice4 ; lattice5 ; lattice6 ; lattice7 ; lattice9 ; ding2025 observe the inverse magnetic catalysis (IMC) phenomena of $u$ and $d$ quarks, which means the decreasing chiral condensates near the pseudocritical temperature $T^{c}_{pc}$ of chiral restoration phase transition and the decreasing pseudocritical temperature $T^{c}_{pc}$ under external magnetic field. Meanwhile, it is reported that the renormalized Polyakov loop increases with magnetic fields and the transition temperature of deconfinement decreases as the magnetic field grows lattice1 ; lattice2 ; lattice2sep ; lattice4 ; lattice5 ; lattice9 . On analytical side, how to explain the inverse magnetic catalysis phenomena is the open question. Many scenarios are proposed fukushima ; mao ; maosep1 ; maosep2 ; maosep3 ; maosep4 ; kamikado ; bf1 ; lattice9 ; bf13 ; bf2 ; bf3 ; bf5 ; bf51 ; bf52 ; bf8 ; bf9 ; bf11 ; db1 ; db1sep ; db2 ; db3 ; db5 ; db6 ; pnjl1 ; pnjl1sep ; pnjl2 ; pnjl3 ; pnjl4 ; pqm ; ferr1 ; ferr2 ; mhuang ; meimao1 ; t0effect ; meihuangmao ; t0effectmao , such as magnetic inhibition of mesons, sphalerons, gluon screening effect, weakening of strong coupling, and anomalous magnetic moment.

The thermodynamical properties of QCD matter are also influenced by the external magnetic field. Among the thermodynamic quantities, the correlations and fluctuations of the conserved charges are accessible in both theoretical calculations and experimental measurements. They can serve as the useful probes to study the QCD phase transitions, such as to identify the critical end point (CEP) of QCD phase diagram in the temperature-baryon chemical potential plane xwork24 ; xwork25 ; xwork26 ; xwork27 ; xwork28 ; xwork29 ; xwork30 and to serve as the possible magnetometer of QCD (probing the presence of magnetic fields in QCD matter) dingref31 ; addding ; ding2022 ; dingprl2024 ; ding2025 . However, they are much less explored at the chiral restoration phase transition (crossover) with finite temperature, magnetic field and vanishing chemical potential. Except for the LQCD calculations addding ; ding2022 ; dingprl2024 ; ding2025 , the analytical investigations at vanishing chemical potential, finite temperature and magnetic field dingref31 ; dingref33 ; dingref34 ; dingref35 ; dingref36 ; mao2pnjl have been conducted in frame of the hadron resonance gas model, Polyakov loop extended Nambu-Jona-Lasinio (PNJL) model and Polyakov loop extended chiral quark model. It should be mentioned that with three-flavor quarks, the IMC effect is not well considered dingref31 ; dingref33 ; dingref34 ; dingref35 ; dingref36 .

In our current paper, the correlations $\chi^{BQ}_{11},\ \chi^{BS}_{11},\ \chi^{QS}_{11}$ and quadratic (quartic) fluctuations $\chi^{B}_{2,4},\ \chi^{Q}_{2,4},\ \chi^{S}_{2,4}$ of baryon number $B$, electric charge $Q$ and strangeness $S$ are studied in a three-flavor PNJL model at finite temperature and magnetic field. The IMC effect is introduced through the magnetic field dependent coupling between quarks $G(eB)$ bf8 ; bf9 ; mao2pnjl ; geb1 ; su3meson4 ; mao11 and magnetic field dependent interaction between quarks and Polyakov loop $T_{0}(eB)$ t0effectmao ; t0effect ; pnjl3 ; mao2pnjl , respectively. The comparison between the results in the cases with and without IMC effect are made. After the brief introduction, Sec.II presents our three-flavor magnetized PNJL model and the definition of correlations $\chi^{BQ}_{11},\ \chi^{BS}_{11},\ \chi^{QS}_{11}$ and quadratic (quartic) fluctuations $\chi^{B}_{2,4},\ \chi^{Q}_{2,4},\ \chi^{S}_{2,4}$. Section III discusses the numerical results of correlations and fluctuations at finite temperature and magnetic field in the cases without IMC effect and with IMC effect. Finally, we give the summary in Sec.IV.

II theoretical framework

The three-flavor PNJL model under external magnetic field is defined with the Lagrangian density pnjl5 ; pnjl6 ; pnjl7 ; pnjl8 ; pnjl9 ; pnjl10 ; pnjl12 ,

$$\mathcal{L} = \bar{\psi}\left(i\gamma^{\mu}D_{\mu}-\hat{m}_{0}\right)\psi+$$(1)

$$\mathcal{L}_{4} = G\sum_{\alpha=0}^{8}\left[(\bar{\psi}\lambda_{\alpha}\psi)^{2}+($$

$$\mathcal{L}_{6} = -K\left[\text{det}\bar{\psi}(1+\gamma_{5})\psi+\text{det}\bar{$$

$${\cal U}(\Phi,{\bar{\Phi}}) = T^{4}\left[-{b_{2}(t)\over 2}\bar{\Phi}\Phi-{b_{3}\over 6}\left($$

The covariant derivative $D^{\mu}=\partial^{\mu}+iQA^{\mu}-i{\cal A}^{\mu}$ couples quarks to the two external fields, the magnetic field ${\bf B}=\nabla\times{\bf A}$ and the temporal gluon field ${\cal A}^{\mu}=\delta^{\mu}_{0}{\cal A}^{0}$ with ${\cal A}^{0}=g{\cal A}^{0}_{a}\lambda_{a}/2=-i{\cal A}_{4}$ in Euclidean space. The gauge coupling $g$ is combined with the SU(3) gauge field ${\cal A}^{0}_{a}(x)$ to define ${\cal A}^{\mu}(x)$, and $\lambda_{a}$ are the Gell-Mann matrices in color space. We consider magnetic field ${\bf B}=(0,0,B)$ along the $z$-axis by setting $A_{\mu}=(0,0,xB,0)$ in Landau gauge, which couples quarks of electric charge $Q=\text{diag}(Q_{u},Q_{d},Q_{s})=\text{diag}(2/3e,-1/3e,-1/3e)$. $\hat{m}_{0}=\text{diag}(m^{u}_{0},m_{0}^{d},m_{0}^{s})$ is the current quark mass matrix in flavor space. The four-fermion interaction $\mathcal{L}_{4}$ represents the interaction in scalar and pseudo-scalar channels, with Gell-Mann matrices $\lambda_{\alpha},\ \alpha=1,2,...,8$ and $\lambda_{0}=\sqrt{2/3}\mathbf{I}$ in flavor space. The six-fermion interaction or Kobayashi-Maskawa-’t Hooft term $\mathcal{L}_{6}$ is related to the $U_{A}(1)$ anomaly tHooft1 ; tHooft2 ; tHooft3 ; tHooft4 ; tHooft5 . The Polyakov potential ${\cal U}(\Phi,\bar{\Phi})$ describes deconfinement at finite temperature, where $\Phi$ is the trace of the Polyakov loop $\Phi=\left({\text{Tr}}_{c}L\right)/N_{c}$, with $L({\bf x})={\cal P}\text{exp}[i\int^{\beta}_{0}d\tau{\cal A}_{4}({\bf x},\tau)$ and $\beta=1/T$, the coefficient $b_{2}(t)=a_{0}+a_{1}t+a_{2}t^{2}+a_{3}t^{3}$ with $t=T_{0}/T$ is temperature dependent, and the other coefficients $b_{3}$ and $b_{4}$ are constants.

It is useful to convert the six-fermion interaction into an effective four-fermion interaction in the mean field approximation, and the Lagrangian density can be rewritten as 3njlrehberg ; 3pnjlmei

$$\mathcal{L} = \bar{\psi}\left(i\gamma^{\mu}D_{\mu}-\hat{m}_{0}\right)\psi-{\cal U$$

$$+ \sum_{a=0}^{8}\left[K_{a}^{-}\left(\bar{\psi}\lambda^{a}\psi$$

$$+ K_{30}^{-}\left(\bar{\psi}\lambda^{3}\psi\right)\left(\bar{\psi}$$

$$+ K_{03}^{-}\left(\bar{\psi}\lambda^{0}\psi\right)\left(\bar{\psi}$$

$$+ K_{80}^{-}\left(\bar{\psi}\lambda^{8}\psi\right)\left(\bar{\psi}$$

$$+ K_{08}^{-}\left(\bar{\psi}\lambda^{0}\psi\right)\left(\bar{\psi}$$

$$+ K_{83}^{-}\left(\bar{\psi}\lambda^{8}\psi\right)\left(\bar{\psi}$$

$$+ K_{38}^{-}\left(\bar{\psi}\lambda^{3}\psi\right)\left(\bar{\psi}$$

with the effective coupling constants

$$K_{0}^{\pm}=G\pm\frac{1}{3}K\left(\sigma_{u}+\sigma_{d}+\sigma_{$$(3)

$$K_{1}^{\pm}=K_{2}^{\pm}=K_{3}^{\pm}=G\mp\frac{1}{2}K\sigma_{s},$$

$$K_{4}^{\pm}=K_{5}^{\pm}=G\mp\frac{1}{2}K\sigma_{d},$$

$$K_{6}^{\pm}=K_{7}^{\pm}=G\mp\frac{1}{2}K\sigma_{u},$$

$$K_{8}^{\pm}=G\mp\frac{1}{6}K\left(2\sigma_{u}+2\sigma_{d}-\sigma$$

$$K_{03}^{\pm}=K_{30}^{\pm}=\mp\frac{1}{2\sqrt{6}}K\left(\sigma_{u$$

$$K_{08}^{\pm}=K_{80}^{\pm}=\mp\frac{\sqrt{2}}{12}K\left(\sigma_{u$$

$$K_{38}^{\pm}=K_{83}^{\pm}=\pm\frac{1}{2\sqrt{3}}K\left(\sigma_{u$$

and chiral condensates

$$\sigma_{u}=\langle\bar{u}u\rangle,\ \sigma_{d}=\langle\bar{d}d$$(4)

The thermodynamic potential in mean field level contains the mean field part and quark part

$$\Omega = 2G(\sigma_{u}^{2}+\sigma_{d}^{2}+\sigma_{s}^{2})-4K\sigma_{u}$$

$$\Omega_{q} = -\sum_{f=u,d,s}\frac{|Q_{f}B|}{2\pi}\sum_{l}\alpha_{l}\int\frac{$$

$$+ T\ln\left(1+3\Phi e^{-\beta E_{f}^{+}}+3{\bar{\Phi}}e^{-2\beta E$$

$$+ T\ln\left(1+3{\bar{\Phi}}e^{-\beta E_{f}^{-}}+3{\Phi}e^{-2\beta E$$

where $f=u,d,s$ means quark flavors, $l$ Landau levels, $\alpha_{l}=2-\delta_{l0}$ spin factor and $E_{f}^{\pm}=E_{f}\pm\mu_{f}$ contains quark energy $E_{f}=\sqrt{p^{2}_{z}+2l|Q_{f}B|+m_{f}^{2}}$ of longitudinal momentum $p_{z}$ and effective quark masses $m_{u}=m_{0}^{u}-4G\sigma_{u}+2K\sigma_{d}\sigma_{s}$, $m_{d}=m_{0}^{d}-4G\sigma_{d}+2K\sigma_{u}\sigma_{s}$, $m_{s}=m_{0}^{s}-4G\sigma_{s}+2K\sigma_{u}\sigma_{d}$, and quark chemical potential $\mu_{u}=\frac{\mu_{B}}{3}+\frac{2\mu_{Q}}{3}$, $\mu_{d}=\frac{\mu_{B}}{3}-\frac{\mu_{Q}}{3}$, $\mu_{s}=\frac{\mu_{B}}{3}-\frac{\mu_{Q}}{3}+\frac{\mu_{S}}{3}$, with $\mu_{B},\ \mu_{Q}\ ,\mu_{S}$ the chemical potential corresponding to the baryon number $B$, electric charge $Q$ and strangeness $S$, respectively.

The ground state at finite temperature, chemical potential and magnetic field is determined by minimizing the thermodynamic potential

$$\frac{\partial\Omega}{\partial\sigma_{f}} = 0,\ (f=u,d,s),$$

$$\frac{\partial\Omega}{\partial\Phi} = 0,$$(5)

$$\frac{\partial\Omega}{\partial{\bar{\Phi}}} = 0.$$

The thermodynamic potential $\Omega$ is a function of order parameters (chiral condensates $\sigma_{f}$ and Polyakov loop $\Phi,{\bar{\Phi}}$), and hence we obtain five coupled gap equations.

At vanishing chemical potential and finite temperature, the chiral symmetry restoration and deconfinement process are smooth crossover. By considering the derivative of the order parameters (chiral condensates and Polyakov loop) with respect to the temperature, the pseuocritical temperatures $T_{pc}^{c}$ and $T_{pc}^{d}$ of chiral restoration and deconfinement phase transitions are determined by the location of the peak of $\frac{d\sigma_{ud}}{dT}\ \text{with}\ \sigma_{ud}=\frac{\sigma_{u}+\sigma_{d}}$ and $\frac{d\Phi}{dT}$, respectively. The strength of chiral restoration and deconfinement phase transitions is characterized by the peak value of $\frac{d\sigma_{ud}}{dT}$ and $\frac{d\Phi}{dT}$, respectively. Note that at vanishing chemical potential, we have $\Phi={\bar{\Phi}}$.

The fluctuations and correlations of baryon number $B$, electric charge $Q$ and strangeness $S$ can be obtained by taking the derivatives of the thermodynamic potential $\Omega$ with respect to the chemical potentials $\hat{\mu}_{X}=\mu_{X}/T,\ (X=B,\ Q,\ S)$, evaluated at zero chemical potential

$$\chi_{i,j}^{B,Q,S} = -\frac{\partial^{i+j+k}(\Omega/T^{4})}{\partial\hat{\mu}_{B}^{i}$$(6)

In this work, we focus on the correlations $\chi^{BQ}_{11}$, $\chi^{BS}_{11}$, $\chi^{QS}_{11}$, quadratic fluctuations $\chi^{B}_{2},\ \chi^{Q}_{2},\ \chi^{S}_{2}$ and quartic fluctuations $\chi^{B}_{4},\ \chi^{Q}_{4},\ \chi^{S}_{4}$ at finite temperature, magnetic field and vanishing chemical potential.

Figure 1: (first two rows) The chiral condensates $\sigma_{u}/\sigma_{u0},\ \sigma_{d}/\sigma_{d0},\ \sigma_{s}/\sigma_{s0}$ and Polyakov loop $\Phi$ as functions of temperature with vanishing chemical potential and fixed magnetic field $eB/m^{2}_{\pi}=0,\ 10,\ 20$. Here, $\sigma_{u0},\ \sigma_{d0},\ \sigma_{s0}$ means up, down, strange quark chiral condensate in vacuum with vanishing temperature, chemical potential and magnetic field, respectively. (last row) The derivative of chiral condensate and Polyakov loop with respect to the temperature $\frac{d\sigma_{ud}}{dT}$ (left panel), $\frac{d\Phi}{dT}$ (right panel) as functions of temperature at vanishing chemical potential and fixed magnetic field $eB/m^{2}_{\pi}=0,\ 10,\ 20$.

Figure 2: Quadratic fluctuations $\chi^{B}_{2},\ \chi^{Q}_{2},\ \chi^{S}_{2}$ (first row), quartic fluctuations $\chi^{B}_{4},\ \chi^{Q}_{4},\ \chi^{S}_{4}$ (second row) and correlations $\chi^{BQ}_{11},\ -\chi^{BS}_{11},\ \chi^{QS}_{11}$ (last row) as functions of temperature with vanishing chemical potential and fixed magnetic field $eB/m^{2}_{\pi}=0,\ 10,\ 20$.

Figure 3: The scaled correlations ${\hat{\chi}}^{BQ}_{11},\ {\hat{\chi}}^{BS}_{11},\ {\hat{\chi}}^{QS}_{11}$ (left panel), scaled quadratic fluctuations ${\hat{\chi}}^{B}_{2},\ {\hat{\chi}}^{Q}_{2},\ {\hat{\chi}}^{S}_{2}$ (middle panel) and scaled quartic fluctuations ${\hat{\chi}}^{B}_{4},\ {\hat{\chi}}^{Q}_{4},\ {\hat{\chi}}^{S}_{4}$ (right panel) at the pseudocritical temperature of chiral restoration phase transition as functions of magnetic field. The vertical lines are the LQCD results ding2022 ; ding2025 ; dingprl2024 .

Figure 4: (upper panel) Magnetic field dependent parameters $G(eB)$ (red line) and $T_{0}(eB)$ (blue line) fitted from LQCD reported decreasing pseudocritical temperature of chiral restoration phase transition $T_{pc}^{c}(eB)/T_{pc}^{c}(eB=0)$ under external magnetic field lattice1 . (lower panel) Strength of chiral restoration phase transition ($\frac{d\sigma_{ud}}{dT}$ at $T=T_{pc}^{c}$) under external magnetic field with IMC effect (blue and red lines) and without IMC effect (black line).

Figure 5: (First row) The correlation $\chi^{BQ}_{11}$ as a function of temperature with fixed magnetic field $eB=10m^{2}_{\pi}$ (left panel) and $eB=20m^{2}_{\pi}$ (middle panel), and the scaled correlation ${\hat{\chi}}^{BQ}_{11}$ (right panel) at $T=T_{pc}^{c}$ as a function of magnetic field. (Second row) The correlation $-\chi^{BS}_{11}$ as a function of temperature with fixed magnetic field $eB=10m^{2}_{\pi}$ (left panel) and $eB=20m^{2}_{\pi}$ (middle panel), and the scaled correlation ${\hat{\chi}}^{BS}_{11}$ (right panel) at $T=T_{pc}^{c}$ as a function of magnetic field. (Third row) The correlation $\chi^{QS}_{11}$ as a function of temperature with fixed magnetic field $eB=10m^{2}_{\pi}$ (left panel) and $eB=20m^{2}_{\pi}$ (middle panel), and the scaled correlation ${\hat{\chi}}^{QS}_{11}$ (right panel) at $T=T_{pc}^{c}$ as a function of magnetic field. The black vertical lines in the right panels are the LQCD results ding2022 ; ding2025 ; dingprl2024 . For all panels, the blue (red) lines are the results in case of considering IMC effect in $T_{0}(eB)$$(G(eB))$ scheme, and the black lines are the results without IMC effect.

Figure 6: (First row) The quadratic fluctuation $\chi^{B}_{2}$ as a function of temperature with fixed magnetic field $eB=10m^{2}_{\pi}$ (left panel) and $eB=20m^{2}_{\pi}$ (middle panel), and the scaled quadratic fluctuation ${\hat{\chi}}^{B}_{2}$ (right panel) at $T=T_{pc}^{c}$ as a function of magnetic field. (Second row) The quadratic fluctuation $\chi^{Q}_{2}$ as a function of temperature with fixed magnetic field $eB=10m^{2}_{\pi}$ (left panel) and $eB=20m^{2}_{\pi}$ (middle panel), and the scaled quadratic fluctuation ${\hat{\chi}}^{Q}_{2}$ (right panel) at $T=T_{pc}^{c}$ as a function of magnetic field. (Third row) The quadratic fluctuation $\chi^{S}_{2}$ as a function of temperature with fixed magnetic field $eB=10m^{2}_{\pi}$ (left panel) and $eB=20m^{2}_{\pi}$ (middle panel), and the scaled quadratic fluctuation ${\hat{\chi}}^{S}_{2}$ (right panel) at $T=T_{pc}^{c}$ as a function of magnetic field. The black vertical lines in the right panels are the LQCD results ding2022 ; ding2025 . For all panels, the blue (red) lines are the results in case of considering IMC effect in $T_{0}(eB)$$(G(eB))$ scheme, and the black lines are the results without IMC effect.

Figure 7: (First row) The quartic fluctuation $\chi^{B}_{4}$ as a function of temperature with fixed magnetic field $eB=10m^{2}_{\pi}$ (left panel) and $eB=20m^{2}_{\pi}$ (middle panel), and the scaled quadratic fluctuation ${\hat{\chi}}^{B}_{4}$ (right panel) at $T=T_{pc}^{c}$ as a function of magnetic field. (Second row) The quartic fluctuation $\chi^{Q}_{4}$ as a function of temperature with fixed magnetic field $eB=10m^{2}_{\pi}$ (left panel) and $eB=20m^{2}_{\pi}$ (middle panel), and the scaled quadratic fluctuation ${\hat{\chi}}^{Q}_{4}$ (right panel) at $T=T_{pc}^{c}$ as a function of magnetic field. (Third row) The quartic fluctuation $\chi^{S}_{4}$ as a function of temperature with fixed magnetic field $eB=10m^{2}_{\pi}$ (left panel) and $eB=20m^{2}_{\pi}$ (middle panel), and the scaled quadratic fluctuation ${\hat{\chi}}^{S}_{4}$ (right panel) at $T=T_{pc}^{c}$ as a function of magnetic field. Here, the blue (red) lines are the results in case of considering IMC effect in $T_{0}(eB)$$(G(eB))$ scheme, and the black lines are the results without IMC effect.

III Numerical Results

III.1 without IMC effect

Because of the contact interaction in NJL model, the ultraviolet divergence cannot be eliminated through renormalization, and a proper regularization scheme is needed. In this part, we apply the covariant Pauli-Villars regularization mao . By fitting the physical quantities, pion mass $m_{\pi}=138\ \text{MeV}$, pion decay constant $f_{\pi}=93\ \text{MeV}$, kaon mass $m_{K}=495.7\ \text{MeV}$, $\eta^{\prime}$ meson mass $m_{\eta\prime}=957.5\ \text{MeV}$ in vacuum, we fix the current masses of light quarks $m_{0}^{u}=m_{0}^{d}=5.5\ \text{MeV}$, and obtain the parameters $m_{0}^{s}=154.7\ \text{MeV}$, $G\Lambda^{2}=3.627$, $K\Lambda^{5}=92.835$, $\Lambda=1101\ \text{MeV}$ 3pnjlmei . For the Polyakov potential, the parameters are chosen as pnjl6 $a_{0}=6.75$, $a_{1}=-1.95$, $a_{2}=2.625$, $a_{3}=-7.44$, $b_{3}=0.75$, $b_{4}=7.5$, and $T_{0}=270$ MeV.

The typical results of PNJL model with the original parameters determined from the vacuum properties present the catalysis phenomena at finite temperature and vanishing chemical potential. The order parameters, chiral condensates $\sigma_{u}/\sigma_{u0},\ \sigma_{d}/\sigma_{d0},\ \sigma_{s}/\sigma_{s0}$ and Polyakov loop $\Phi$, increase with the magnetic field in the whole temperature region, as shown in Fig.1 first two rows. The last row of Fig.1 depicts the derivative of the order parameters with respect to the temperature $\frac{d\sigma_{ud}}{dT}$ and $\frac{d\Phi}{dT}$ as functions of temperature with fixed magnetic fields. They firstly increase and then decrease with temperature, showing apparent peaks. The location of the peak determines the pseudocritical temperature for chiral restoration and deconfinement phase transitions, respectively, which increase with the magnetic field. The peak value of $\frac{d\sigma_{ud}}{dT}$ and $\frac{d\Phi}{dT}$ determines the strength of chiral restoration and deconfinement phase transitions, respectively, which also increase with the magnetic field. It should be mentioned that the catalysis phenomena of chiral condensates at high temperature and pseudocritical temperature of chiral restoration phase transition under external magnetic field are contrary to LQCD results lattice1 ; lattice2 ; lattice2sep ; lattice4 ; lattice5 ; lattice6 ; lattice7 ; lattice9 ; ding2025 , but the increasing strength of chiral restoration phase transition with increasing magnetic field is consistent with LQCD results lattice1 ; lattice2 ; lattice2sep .

Figure 2 plots quadratic fluctuations $\chi^{B}_{2},\ \chi^{Q}_{2},\ \chi^{S}_{2}$ (first row), quartic fluctuations $\chi^{B}_{4},\ \chi^{Q}_{4},\ \chi^{S}_{4}$ (second row) and correlations $\chi^{BQ}_{11},\ -\chi^{BS}_{11},\ \chi^{QS}_{11}$ (last row) as functions of temperature with vanishing chemical potential and fixed magnetic field $eB/m^{2}_{\pi}=0,\ 10,\ 20$. At vanishing magnetic field, $\chi^{B}_{2},\ \chi^{Q}_{2},\ \chi^{S}_{2}$ firstly increase and then decrease with temperature, with the wide peak around the chiral restoration and deconfinement phase transitions. Turning on the magnetic field, the peak structure becomes more pronounced. In high temperature region, $\chi^{B}_{2}$ and $\chi^{Q}_{2}$ have the higher values as magnetic field grows, but in low temperature region they show non-monotonic behavior with increasing magnetic field. $\chi^{S}_{2}$ increases with magnetic field in the whole temperature region. For quartic fluctuations $\chi^{B}_{4},\ \chi^{Q}_{4},\ \chi^{S}_{4}$, we observe the apparent peak around the chiral restoration and deconfinement phase transitions with and without magnetic field. Their values around the peak are enhanced by the magnetic field. In low temperature region, $\chi^{B}_{4}$ and $\chi^{Q}_{4}$ show non-monotonic dependence on the magnetic field, and in high temperature region they decrease under magnetic fields. $\chi^{S}_{4}$ is not sensitive to the magnetic field in the temperature region away from the phase transitions. The correlations $\chi^{BQ}_{11},\ -\chi^{BS}_{11},\ \chi^{QS}_{11}$ show peak structure around phase transitions, which are enhanced by the magnetic field. In high temperature region, $\chi^{BQ}_{11}$ and $-\chi^{BS}_{11}$ have the higher values as magnetic field grows, but in low temperature region they show non-monotonic behavior with increasing magnetic field. $\chi^{QS}_{11}$ increases with magnetic field in the whole temperature region. The peak structure of fluctuations and correlations is caused by the phase transition (crossover). When turning on the external magnetic field, the strength of phase transition increases (see Fig.1 last row), and this leads to the enhancement of the fluctuations and correlations around the phase transition. It is noticeable that fluctuations and correlations related with strangeness $(\chi^{S}_{2},\ \chi^{S}_{4},\ -\chi^{BS}_{11},\ \chi^{QS}_{11})$ are less sensitive to the magnetic field than others $(\chi^{B}_{2},\ \chi^{Q}_{2},\ \chi^{B}_{4},\ \chi^{Q}_{4},\ \chi^{BQ}_{11})$, and the peak structure of $\chi^{S}_{2},\ \chi^{S}_{4},\ -\chi^{BS}_{11},\ \chi^{QS}_{11}$ is less apparent than $\chi^{B}_{2},\ \chi^{Q}_{2},\ \chi^{B}_{4},\ \chi^{Q}_{4},\ \chi^{BQ}_{11}$, either. The reason is that all up, down and strange quarks make contribution to $\chi^{B}_{2},\ \chi^{Q}_{2},\ \chi^{B}_{4},\ \chi^{Q}_{4},\ \chi^{BQ}_{11}$, but $\chi^{S}_{2},\ \chi^{S}_{4},\ -\chi^{BS}_{11},\ \chi^{QS}_{11}$ are mainly related to strange quarks, which are much more heavier than up and down quarks. The quadratic fluctuations $\chi^{B}_{2},\ \chi^{Q}_{2},\ \chi^{S}_{2}$ and correlation $\chi^{BQ}_{11},\ -\chi^{BS}_{11},\ \chi^{QS}_{11}$ are qualitatively consistent with LQCD simulations addding ; ding2025 . Comparing with the results in Ref dingref35 , where peak structure is observed only in correlation $\chi^{BQ}_{11}$ and quartic fluctuations $\chi^{B}_{4},\ \chi^{Q}_{4},\ \chi^{S}_{4}$, the difference is owing to the different regularization schemes, where the regularization is operated in the vacuum term $(T=\mu_{f}=0)$ of thermodynamical potential in Ref dingref35 , but the regularization is applied in both the vacuum and medium term of thermodynamical potential in our current work.

What is the property of correlations and fluctuations along the phase transition line under external magnetic field? As shown in Fig.3, the scaled correlations ${\hat{\chi}}_{11}^{BQ}=\frac{\chi_{11}^{BQ}(eB,T_{pc}^{c}(eB))}{\chi_{11}^{BQ}$, ${\hat{\chi}}_{11}^{BS}=\frac{\chi_{11}^{BS}(eB,T_{pc}^{c}(eB))}{\chi_{11}^{BS}$, ${\hat{\chi}}_{11}^{QS}=\frac{\chi_{11}^{QS}(eB,T_{pc}^{c}(eB))}{\chi_{11}^{QS}$ and scaled quadratic (quartic) fluctuations ${\hat{\chi}}_{2,4}^{B}=\frac{\chi_{2,4}^{B}(eB,T_{pc}^{c}(eB))}{\chi_{2,4}^{B}$, ${\hat{\chi}}_{2,4}^{Q}=\frac{\chi_{2,4}^{Q}(eB,T_{pc}^{c}(eB))}{\chi_{2,4}^{Q}$, ${\hat{\chi}}_{2,4}^{S}=\frac{\chi_{2,4}^{S}(eB,T_{pc}^{c}(eB))}{\chi_{2,4}^{S}$, at the pseudocritical temperature $T_{pc}^{c}$ of chiral restoration phase transition are plotted as functions of magnetic field. They are all increasing functions of magnetic field, which is related with the increasing strength of phase transition under external magnetic field, shown in Fig.1 last row. Among the scaled correlations and fluctuations, ${\hat{\chi}}^{BQ}_{11}$ increases fastest, which may serve as a magnetometer of QCD dingprl2024 ; ding2025 . These properties are qualitatively consistent with the LQCD results ding2022 ; ding2025 ; dingprl2024 , depicted in the vertical lines in Fig.3 left and middle panels. At the quantitative level, ${\hat{\chi}}^{BS}_{11}$, ${\hat{\chi}}^{QS}_{11}$ and ${\hat{\chi}}^{S}_{2}$ are reasonably good, ${\hat{\chi}}^{BQ}_{11}$ and ${\hat{\chi}}^{B}_{2}$ slightly undershoot the LQCD results, and ${\hat{\chi}}^{Q}_{2}$ overshoots the LQCD result.

III.2 with IMC effect

From LQCD simulations lattice1 ; lattice2 ; lattice2sep ; lattice4 ; lattice5 ; lattice6 ; lattice7 ; lattice9 , the IMC phenomenon of chiral symmetry restoration of $u$ and $d$ quarks can be characterized either by the chiral condensates or by the pseudocritical temperature of chiral restoration phase transition. To simulate the IMC effect in the effect model, one approach is to fit the LQCD results of chiral condensates geb1 ; meson ; geb3 ; geb4 , and another approach is to fit the LQCD results of pseudocritical temperature bf8 ; bf9 ; geb1 ; mao2pnjl ; su3meson4 ; mao11 ; t0effectmao . The two methods give consistent results with each other.

By fitting the LQCD reported decreasing pseudocritical temperature of chiral symmetry restoration $T_{pc}^{c}(eB)/T_{pc}^{c}(eB=0)$ under external magnetic field lattice1 , we introduce the IMC effect in our three-flavor PNJL model through a magnetic field dependent parameter $G(eB)$ and $T_{0}(eB)$, respectively, which represent the influence of external magnetic field to the quark-gluon interaction. On one side, the coupling between quarks plays a significant role in determining the spontaneous breaking and restoration of chiral symmetry. Considering the direct interaction between quarks and external magnetic field, a magnetic field dependent coupling $G(eB)$ bf8 ; bf9 ; mao2pnjl ; geb1 ; su3meson4 ; mao11 is introduced into the PNJL model. On the other side, the interaction between the Polyakov loop and sea quarks may be important for the mechanism of IMC lattice9 . A magnetic field dependent Polyakov loop scale parameter $T_{0}(eB)$ t0effectmao ; t0effect ; pnjl3 ; mao2pnjl is introduced into the PNJL model to mimic the reaction of the gluon sector to the presence of magnetic fields. As plotted in Fig.4$,G(eB)$ and $T_{0}(eB)$ are both monotonic decreasing functions of magnetic field. We have checked that, with our fitted parameter $G(eB)$ or $T_{0}(eB)$, the increase (decrease) of chiral condensates with magnetic fields at the low (high) temperature, the increase of Polyakov loop with magnetic fields in the whole temperature region and the reduction of pseudocritical temperature of deconfinement phase transition under magnetic fields can be realized. Moreover, the inclusion of IMC effect does not change the strength of chiral restoration phase transition under external magnetic field, see Fig.4 lower panel, where $\frac{d\sigma_{ud}}{dT}$ at the pseudocritical temperature of chiral restoration phase transition increases with magnetic fields and are almost coincident in the cases with and without IMC effect. This indicates that including IMC effect will not lead to qualitative difference in the results of correlations and fluctuations.

In the following part, we use the $G(eB)$ and $T_{0}(eB)$ scheme to consider the IMC effect, respectively, and make detailed comparison between the results of correlations and fluctuations with and without IMC effect.

Figure 5 depicts the correlations $\chi^{BQ}_{11}$, $-\chi^{BS}_{11}$, $\chi^{QS}_{11}$ with IMC effect (blue and red lines) and without IMC effect (black lines). At fixed magnetic fields $eB/m^{2}_{\pi}=10,\ 20$, when considering the IMC effect with $G(eB)$ and $T_{0}(eB)$ schemes, the correlations $\chi^{BQ}_{11}$, $-\chi^{BS}_{11}$, $\chi^{QS}_{11}$ increase with temperature and have the peak around the pseudocritical temperatures of chiral restoration and deconfinement phase transitions, which are similar as the case without IMC effect, see the left and middle panels. The value of $\chi^{BQ}_{11}$ with IMC effect are larger (smaller) than those without IMC effect in the low (high) temperature region. The value of $-\chi^{BS}_{11}$ and $\chi^{QS}_{11}$ with IMC effect are larger than those without IMC effect in the whole temperature region. $\chi^{BQ}_{11}$ shows sharper peak than $-\chi^{BS}_{11}$ and $\chi^{QS}_{11}$. In the right panels, we plot the scaled correlations ${\hat{\chi}}_{11}^{BQ}$, ${\hat{\chi}}_{11}^{BS}$, ${\hat{\chi}}_{11}^{QS}$ at the pseudocritical temperature of chiral restoration phase transition as functions of magnetic field with IMC effect (blue and red lines) and without IMC effect (black lines), and LQCD results ding2022 ; ding2025 ; dingprl2024 are depicted in the black vertical lines for references. ${\hat{\chi}}_{11}^{BQ}$, ${\hat{\chi}}_{11}^{BS}$ and ${\hat{\chi}}_{11}^{QS}$ increase with magnetic field. Different methods to consider IMC effect do not lead to qualitative difference in the results, and only cause some quantitative changes. Including the IMC effect with $G(eB)$ and $T_{0}(eB)$ schemes, ${\hat{\chi}}^{BQ}_{11}$ increases faster. ${\hat{\chi}}_{11}^{BS}$ and ${\hat{\chi}}_{11}^{QS}$ increase faster (slower) in the $T_{0}(eB)\ (G(eB))$ scheme.

Figure 6 plots the quadratic fluctuations $\chi^{B}_{2},\ \chi^{Q}_{2},\ \chi^{S}_{2}$ with IMC effect (blue and red lines) and without IMC effect (black lines). At fixed magnetic fields $eB/m^{2}_{\pi}=10,\ 20$, when considering the IMC effect with $G(eB)$ and $T_{0}(eB)$ schemes, $\chi^{B}_{2},\ \chi^{Q}_{2},\ \chi^{S}_{2}$ increase with temperature and have the peak around the pseudocritical temperatures of chiral restoration and deconfinement phase transitions, which are similar as the case without IMC effect, see the left and middle panels. Their value with IMC effect are larger than those without IMC effect in the whole temperature region. In the right panels, we plot the scaled quadratic fluctuations ${\hat{\chi}}_{2}^{B}$, ${\hat{\chi}}_{2}^{Q}$, ${\hat{\chi}}_{2}^{S}$ at the pseudocritical temperature of chiral restoration phase transition as functions of magnetic field with IMC effect (blue and red lines) and without IMC effect (black lines), and LQCD results ding2022 ; ding2025 are depicted in the black vertical lines for references. ${\hat{\chi}}_{2}^{B}$, ${\hat{\chi}}_{2}^{Q}$ and ${\hat{\chi}}_{2}^{S}$ increase with magnetic field. Different methods to consider IMC effect do not lead to qualitative difference in the results, and only cause some quantitative changes. Including the IMC effect with $G(eB)$ and $T_{0}(eB)$ schemes, ${\hat{\chi}}^{B}_{2}$ and ${\hat{\chi}}^{Q}_{2}$ increase faster than without IMC effect. ${\hat{\chi}}_{2}^{S}$ increases faster in the $T_{0}(eB)$ scheme and increases slower in the $G(eB)$ scheme.

Figure 7 shows the quartic fluctuations $\chi^{B}_{4}$, $\chi^{Q}_{4}$, $\chi^{S}_{4}$ with IMC effect (blue and red lines) and without IMC effect (black lines). At fixed magnetic fields $eB/m^{2}_{\pi}=10,\ 20$, when considering the IMC effect with $G(eB)$ and $T_{0}(eB)$ schemes, $\chi^{B}_{4}$, $\chi^{Q}_{4}$ and $\chi^{S}_{4}$ increase with temperature and have the peak around the pseudocritical temperatures of chiral restoration and deconfinement phase transitions, which are similar as the case without IMC effect, see the left and middle panels. Their values with IMC effect are larger (smaller) than those without IMC effect in the low (high) temperature region. $\chi^{B}_{4}$ and $\chi^{Q}_{4}$ show sharper peak than $\chi^{S}_{4}$. In the right panels, we plot the scaled quartic fluctuations ${\hat{\chi}}_{4}^{B}$, ${\hat{\chi}}_{4}^{Q}$, ${\hat{\chi}}_{4}^{S}$ at the pseudocritical temperature of chiral restoration phase transition as functions of magnetic field with IMC effect (blue and red lines) and without IMC effect (black lines). They increase with magnetic field. Different methods to consider IMC effect do not lead to qualitative difference in the results, and only cause some quantitative changes. Including the IMC effect does not alter ${\hat{\chi}}_{4}^{B}$ too much. ${\hat{\chi}}_{4}^{Q}$ increases faster (slower) in the $G(eB)\ (T_{0}(eB))$ scheme. Including the IMC effect with $G(eB)$ and $T_{0}(eB)$ schemes, ${\hat{\chi}}_{4}^{S}$ increases slower than without IMC effect.

IV summary

The correlations $\chi^{BQ}_{11},\ \chi^{BS}_{11},\ \chi^{QS}_{11}$ and fluctuations $\chi^{B}_{2,4},\ \chi^{Q}_{2,4},\ \chi^{S}_{2,4}$ at finite temperature and magnetic field are investigated in frame of a three-flavor PNJL model. The IMC effect is introduced into the PNJL model by two methods, the magnetic field dependent coupling between quarks $G(eB)$, and magnetic field dependent interaction between quarks and Polyakov loop $T_{0}(eB)$. We make comparison of the results in the cases with and without IMC effect.

The PNJL model with IMC effect recovers the decreasing pseudocritical temperature of phase transitions under external magnetic field and the properties of order parameters reported in LQCD simulations. But the increasing strength of chiral restoration phase transition under external magnetic field is obtained in the cases with and without IMC effect. Therefore, the inclusion of IMC effect does not lead to qualitative difference in the results of correlations and fluctuations, and only causes some quantitative changes.

The correlations $\chi^{BQ}_{11},\ \chi^{BS}_{11},\ \chi^{QS}_{11}$ and fluctuations $\chi^{B}_{2,4},\ \chi^{Q}_{2,4},\ \chi^{S}_{2,4}$ show similar properties at finite temperature and magnetic field. With and without magnetic fields, the correlations and fluctuations first increase with temperature, and then show the peak around the pseudocritical temperatures of chiral restoration and deconfinement phase transitions. The peak structure in $\chi^{BQ}_{11}$, $\chi^{B}_{4}$ and $\chi^{Q}_{4}$ is much more apparent than in other correlations and fluctuations. To demonstrate the correlations and fluctuations along the phase transition line under external magnetic field, we calculate the scaled correlations ${\hat{\chi}}_{11}^{BQ}$, ${\hat{\chi}}_{11}^{BS}$, ${\hat{\chi}}_{11}^{QS}$ and scaled fluctuations ${\hat{\chi}}_{2,4}^{B}$, ${\hat{\chi}}_{2,4}^{Q}$, ${\hat{\chi}}_{2,4}^{S}$ at the pseudocritical temperature of chiral restoration phase transition. They increase with magnetic fields, which is caused by the increase of phase transition strength under external magnetic field. Among them, ${\hat{\chi}}_{11}^{BQ}$ increases fastest, which may serve as the magnetometer of QCD. These properties of correlations and quadratic fluctuations are qualitatively consistent with the LQCD results.

Acknowledgement: The work is supported by the NSFC grant 12275204.