Heavy Quark Potential in a static and strong homogeneous magnetic field

TL;DR

This paper studies how strong magnetic fields influence the potential between heavy quarks in a hot QCD medium, revealing magnetic field thresholds for quarkonium dissociation, with implications for understanding quark-gluon plasma under extreme conditions.

Contribution

It introduces a method to calculate the heavy quark potential in a strong magnetic field, showing significant effects on long-range interactions and quarkonium dissociation thresholds.

Findings

Long-range quarkonium potential is highly affected by magnetic fields.

J/ψ dissociates at magnetic field ~10 m_π^2.

Υ dissociates at magnetic field ~100 m_π^2.

Abstract

We have investigated the properties of quarkonia in a thermal QCD medium in the background of strong magnetic field. For that purpose, we employ the Schwinger proper-time quark propagator in the lowest Landau level to calculate the one-loop gluon self-energy, which in the sequel gives the the effective gluon propagator. As an artifact of strong magnetic field approximation ($eB>>T^2$ and $eB>>m^2$), the Debye mass for massless flavors is found to depend only on the magnetic field which is the dominant scale in comparison to the scales prevalent in the thermal medium. However, for physical quark masses, it depends on both magnetic field and temperature in a low temperature and high magnetic field but the temperature dependence is very meagre and becomes independent of temperature beyond a certain temperature and magnetic field. With the above mentioned ingredients, the potential between heavy quark ($Q$) and anti-quark ($\bar Q$) is obtained in a hot QCD medium in the presence of strong magnetic field by correcting both short and long range components of the potential in real-time formalism. It is found that the long range part of the quarkonium potential is affected much more by magnetic field as compared to the short range part. This observation facilitates us to estimate the magnetic field beyond which the potential will be too weak to bind $Q\bar Q$ together. For example, the $J/\psi$ is dissociated at $eB \sim$ 10 $m_\pi^2$ and $\Upsilon$ is dissociated at $eB \sim$ 100 $m_\pi^2$ whereas its excited states, $\psi^\prime$ and $\Upsilon^\prime$ are dissociated at smaller magnetic field $eB= m_\pi^2$, $13 m_\pi^2$, respectively.

Equations

$$S(x,x^{\prime })=\varphi (x,x^{\prime })\int \frac{d^{4}p}{(2\pi {)}^4}{e}^{ip(x-x^{\prime })}S(p),$$

$$\varphi (x,x^{\prime })=e^{ie{\int }_x^{x^{\prime }}A(\xi )d\xi },$$

$$iS(p)={\int }_0^{\infty }\frac{1}{eB}\frac{ds}{\cos (s)}{e}^{-is\left[m^2-p_{\Vert }^2+\frac{\tan (s)}{s}{p}_{\perp }^2\right]}\left[\left\{\cos (s)+{\gamma }_1{\gamma }_2\sin (s)\right\}(m+\gamma \cdot p_{\Vert })-\frac{\gamma \cdot p_{\perp }}{\cos (s)}\right],$$

$$iS(p)=\sum _n\frac{-i{d}_n(\alpha )D+d_n^{\prime }(\alpha )\bar{D}}{p_L^2+2neB}+i\frac{\gamma \cdot p_{\perp }}{p_{\perp }^2},$$

$$D$$

$$\bar{D}$$

$$d_n(\alpha )$$

$$C_n(2\alpha )$$

$$d_n^{{}^{\prime }}(\alpha )$$

$$p_L^2$$

$$\alpha$$

$$p_{\Vert }^2$$

$$p_{\perp }^2$$

$$i{S}_0(p)=\frac{(1+{\gamma }^0{\gamma }^3{\gamma }^5)({\gamma }^0{p}_0-{\gamma }^3{p}_z+m)}{p_{\parallel }^2-m^2+iϵ}{e}^{-\frac{p_{\perp }^2}{\mid qB\mid }},$$

$$S_{ab}(p)=U_{ac}(p){\left(\begin{array}{cc}{S}_0(p)& 0\\ 0& -S_0^{\ast }(p)\end{array}\right)}_{cd}{U}_{db}(p),$$

$$U(p)=\left(\begin{array}{cc}{\tilde{N}}_2& -{\tilde{N}}_1{e}^{-\beta \mu /2}\\ {\tilde{N}}_1{e}^{-\beta \mu /2}& {\tilde{N}}_2\end{array}\right),$$

$$\widetilde{N_1}(p_0)$$

$$\widetilde{N_2}(p_0)$$

$$n_p^{\pm }(p_0)$$

$$S(p)=\left(\begin{array}{cc}{S}_0(p){\tilde{N}}_2^2+S_0^{\ast }(p){\tilde{N}}_1^2& -S_0(p)\widetilde{N_1}\widetilde{N_2}+S_0^{\ast }(p)\widetilde{N_1}\widetilde{N_2}\\ S_0(p)\widetilde{N_1}\widetilde{N_2}-S_0^{\ast }(p)\widetilde{N_1}\widetilde{N_2}& -S_0(p){\tilde{N}}_1^2-S_0^{\ast }(p){\tilde{N}}_2^2\end{array}\right).$$

$$S_{11}(p)=S_0(p){\tilde{N}}_2^2+S_0^{\ast }(p){\tilde{N}}_1^2,$$

$$iS_{11}(p)=\Bigg{[}\frac{1}{{p_{\parallel}^{2}-m^{2}+i\epsilon}}+2\pi in_{p}\delta(p_{\parallel}^{2}-m^{2})\Bigg{]}(1+\gamma^{0}\gamma^{3}\gamma^{5})(\gamma^{0}p_{0}-\gamma^{3}p_{z}+m)e^{\frac{-p_{\perp}^{2}}{\mid qB\mid}},$$

$$n_p(p_0)=\frac{1}{e^{\beta \mid p_0\mid }+1}.$$

$${\alpha }_s^{\Vert }(k_3)=\frac{1}{{{\alpha }_s^0({\mu }_0)}^{-1}+\frac{11N_c}{12\pi }\ln (\frac{k_3^2+M_B^2}{{\mu }_0^2})+\frac{1}{3\pi }{\sum }_f\frac{|q_{f}B|}{\sigma }},$$

$${\alpha }_s^0({\mu }_0)=\frac{12\pi }{11N_c\ln (\frac{({\mu }_0^2+M_B^2)}{{\Lambda }_V^2})}.$$

$${\Pi }^{\mu \nu }(k)$$

$$L^{\mu \nu }=8\left[p_{\parallel }^{\mu }{q}_{\parallel }^{\nu }+p_{\parallel }^{\nu }{q}_{\parallel }^{\mu }-g_{\parallel }^{\mu \nu }((p.q{)}_{\parallel }-m_f^2)\right].$$

$${\Pi }^{\mu \nu }(k)=\sum _f{\Pi }_{\parallel }^{\mu \nu }(k_{\parallel })A_f(k_{\perp }),$$

$$A_f(k_{\perp })$$

$${\Pi }_{\parallel }^{\mu \nu }=({\Pi }_{\parallel }^{\mu \nu }{)}_V+({\Pi }_{\parallel }^{\mu \nu }{)}_n+({\Pi }_{\parallel }^{\mu \nu }{)}_{n^2},$$

$$({\Pi }_{\parallel }^{\mu \nu }{)}_V$$

$$({\Pi }_{\parallel }^{\mu \nu }{)}_n$$

$$({\Pi }_{\parallel }^{\mu \nu }{)}_{n^2}$$

$$\frac{1}{x\mp iϵ}=\text{P}\left(\frac{1}{x}\right)\pm i\pi \delta (x),$$

$$\Re{\Pi^{\mu\nu}(k)\mid}_{V}=\Big{(}g_{\parallel}^{\mu\nu}-\frac{k_{\parallel}^{\mu}k_{\parallel}^{\nu}}{k_{\parallel}^{2}}\Big{)}\Pi(k^{2}),$$

$$\displaystyle\Pi(k^{2})=\frac{g^{2}}{4\pi^{2}}\sum_{f}\mid q_{f}B\mid e^{-\frac{k_{\perp}^{2}}{2|q_{f}|B}}\left[\frac{2m_{f}^{2}}{k_{\parallel}^{2}}\left(1-\frac{4m_{f}^{2}}{k_{\parallel}^{2}}\right)^{-1/2}\left\{\ln\frac{1-{\Big{(}1-\frac{4m_{f}^{2}}{k_{\parallel}^{2}}\Big{)}}^{1/2}}{1+{\Big{(}1+\frac{4m_{f}^{2}}{k_{\parallel}^{2}}\Big{)}}^{1/2}}+i\pi\right\}-1\right].$$

Full text

PACS:   12.39.-x,11.10.St,12.38.Mh,12.39.Pn

Keywords: Quantum Chromodynamics, Schwinger proper-time method, Debye mass, strong magnetic field, string tension, dielectric permittivity, Heavy quark potential.

1 Introduction

Lattice gauge theory at very high temperatures and/or baryon densities predicts an interesting window onto the properties of Quantum Chromodynamics (QCD) in guise of a new phase, Quark-gluon Plasma (QGP), which pervaded the early universe, and may be present in the core of neutron stars. To realize this predicted phase, current experimental program of ultra relativistic heavy ion collisions (URHIC) have been designed at different colliders with different center of mass energies, viz. Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory (BNL) at $\sqrt{s}$= 200 GeV per nucleon in Au + Au collisions and Large Hadron Collider (LHC) at European Organization for Nuclear Research (CERN) at $\sqrt{s}$= 2.76 TeV per nucleon in Pb + Pb collisions. Recent analysis suggests that the events of URHIC should be analyzed by incorporating the effect of magnetic field because an intensely strong magnetic field, perpendicular to the reaction plane, is expected to be produced at very early stages of collisions when the event is off-central [1, 2, 3, 4, 5]. Depending on the centrality, the strength of the magnetic field may reach between $m_{\pi}^{2}$ ($\simeq 10^{18}$ Gauss) at RHIC [6] to 10 $m_{\pi}^{2}$ at LHC [7]. At extreme cases it may reach values of 50 $m_{\pi}^{2}$. A very strong magnetic field ($\sim 10^{23}$ Gauss) was also produced in the early universe during the electroweak phase transition due to the gradients in Higgs field [8]. Ultimately such strong magnetic field might significantly affect the production of particles and alter their dynamics at very early stage of the collisions. Since magnetic field induces an anisotropy to the momentum of the affected particles, we might expect it to affect the anisotropic flow of the particles.

Naive classical estimates predict that the magnetic field may be very strong typically up to $t_{B}\simeq 0.2$ fm [9]. However, the realistic calculations on the charge transport properties of the plasma (namely, conductivity) may suggest that the magnetic field may remain substantial for significantly longer time [10]. Simultaneously heavy quark and antiquarks pairs also develop into a physical resonances over a formation time $t_{\rm form}\sim 1/E_{\rm bind}$ related to the binding energy of the state, e.g. the charm-anti charm ($c\bar{c}$) pairs form resonances at $t_{c\bar{c}}\sim 0.3$ fm. Thus it becomes reasonable to assume that charmonium production may get significantly influenced by the strong magnetic field. The same argument applies to the bottomonium production. A large number of studies on the in-medium properties of $Q\bar{Q}$ bound states has been carried out using the phenomenological potential models [11, 12, 13], where the effects of the medium are encoded in a temperature dependent potential with non perturbative inputs from the lattice simulations. However, lattice calculations of free energies and other quantities [14] obtained from the correlation functions of Polyakov loop are often taken as input for the potential. Although these quantities have been thought to be related to the color-singlet and color-octet heavy quark potentials at finite temperature, a precise answer is still missing. Recently quarkonia at finite temperature has been studied by taking the advantage of the hierarchies between the non-relativistic scales associated with quarkonia and the thermal energy scales characterizing the system through the effective field theories, viz. NRQCD, pNRQCD etc. [15]. The in-medium modifications of the quarkonium states can be studied from the first principle of QCD by the spectral functions [16] but the reconstruction of the spectral function from the lattice meson correlators turns out to be very difficult. Recently in a new theoretical developments, the heavy quark potential have also been synthesized in strong coupling regime through a novel idea of gauge-gravity duality  [17, 18, 19].

Recently the properties of quarkonia states in a hot medium are explored in perturbative thermal QCD framework by correcting both the perturbative and non perturbative terms of the $Q\bar{Q}$ potential through the dielectric function in real-time formalism [20, 21] in both isotropic as well as anisotropic hot QCD medium, where the anisotropy in the momentum space arose at the very early stages of the collisions due to the different expansion rate in the longitudinal and transverse direction [22]. As mentioned earlier, magnetic field is also produced at the early stages of the collisions, thus it becomes worthwhile to examine the effects of magnetic field on the properties of quarkonia bound states, which is the central theme of our present work. Quantum mechanically both the quarkonium and heavy meson spectra have been analyzed through the solution of non relativistic Schrödinger equation with both harmonic oscillator and Cornell potential with an additional spin-spin interaction term [23, 24]. Moreover lattice studies have also recently explored the possible anisotropies emerging in the static quark-anti quark potentials both at $T=0$ and $T\neq 0$ through the Wilson loop expectation value and Polyakov loop correlator, respectively, in the presence of a magnetic background with respect to the direction of magnetic field [25, 26].

Here we have tried to investigate the effect of strong and homogeneous magnetic field on the properties of quarkonia states. We have first calculated the gluon self energy at finite temperature in a strong magnetic background and then obtain the heavy quark potential by taking the static limit of the effective gluon propagator to see the effects of magnetic field alone on the quarkonium states even in a thermal medium. This is due to the fact that the magnetic field is assumed much stronger than the temperature as well as the mass of the quarks in quark-loop of gluon self-energy ($eB\gg T^{2}$ and $eB\gg m^{2}$), known as “Strong Magnetic Field Approximation (SMFA)”. As a consequence the Debye mass obtained from the static limit of gluon self-energy becomes almost independent of temperature, hence the potential even in the thermal medium depends mainly on the magnetic field because the medium dependence in the potential enters through the Debye mass. Moreover there is another condition for the non relativistic potential approach for heavy quarkonia to be valid is that the mass of the heavy quark (either charm or bottom quark) should be larger than the dominant scale available in the problem ($m_{Q}>>\sqrt{eB}$) (dimensionally) because $\sqrt{eB}$ is the most dominant scale available in the strongly magnetized thermal medium. Thus the above mentioned two conditions constraint the lower and upper limit of the magnetic fields, respectively, for which our work is valid. As a bi-product of this constraint, the magnetic field expected to be produced at RHIC may not satisfy the condition of SMFA. So our work which is valid only in SMFA will be suitable to the LHC events where the magnetic field expected to be produced are well within the limit of the validity of our work.

Our work is arranged in the following way: In subsection 2.1 we will discuss the quark propagator at finite temperature within SMFA. In subsection 2.2 we will calculate the gluon self energy at finite temperature in presence of strong magnetic field. In subsection 2.3 we will compute the screening mass in SMFA by taking the static limit of gluon-self energy. In Section 3, we will obtain the potential from the inverse Fourier transform of the effective propagator in the static limit and explore how the properties of quarkonia could be affected by the presence of strong magnetic field. Finally we will conclude in Section 4.

2 One-loop gluon self-energy and the screening mass in SMFA

The gluon self energy can be affected by the magnetic field in two ways: First, the quark propagator gets affected due to the Landau quantization of the energy levels (known as Landau levels) in the presence of magnetic field. Second, the strong coupling runs with both the magnetic field and temperature. However, in SMFA, it runs exclusively with the magnetic field as we have discussed in the introduction.

Schwinger first obtained the fermionic propagator in coordinate space by the proper-time formalism [27], then Tsai has obtained the same in momentum space and used it to calculate the vacuum polarization in magnetic field [28]. The vacuum polarization tensor has also been obtained in a gauge invariant manner for both strong and weak magnetic field limit [29, 30, 31]. We shall now going to extend these calculations to QCD to calculate the gluon self-energy, which will in turn help to study the properties of quarkonia quantum field theoretically in the presence of strong magnetic field.

2.1 Fermionic propagator in presence of magnetic field

2.1.1 Vacuum in a static and homogeneous magnetic

field:

For the sake of simplicity, we assume the magnetic field to be constant and homogeneous. We also assume the magnetic field to be along $z$-direction and of magnitude $B$. Such a magnetic field can be obtained from a vector potential $A_{\mu}=(0,0,Bx,0)$. The choice of vector potential is not unique as the same magnetic field can also be obtained from a symmetric potential given by $A_{\mu}=(0,\frac{-By}{2},\frac{Bx}{2},0)$. Using the proper-time method formulated by Schwinger, the fermion propagator in such a magnetic field can be written in the coordinate space as [27]

[TABLE]

where the the phase factor, $\phi(x,x^{\prime})$ is given by

[TABLE]

which becomes unity for a closed fermion loop with two fermion lines, i.e, $\phi(x,x^{\prime})=1$ [32].

However, the same propagator was first calculated by [28, 32] in the momentum space as

[TABLE]

The propagator (3) in the momentum space can also be expressed in a more convenient way using the associated Laguerre polynomials ($L_{n}$)

[TABLE]

where the following quantities are defined as [32]

[TABLE]

The order of the Laguerre polynomial also corresponds to the number of energy eigenvalues in a magnetic field, known as Landau levels. In SMFA, the particles occupy the lowest Landau level (LLL) ($n=0$) only, thus, in SMFA, the fermion propagator in eq.(4) reduces to the following form

[TABLE]

where $m$ and $q$ are the mass and electric charge of the fermion, respectively.

2.1.2 Heat Bath in a strong homogeneous magnetic field

In thermal medium, the system possesses additional thermal scales, viz. T, $gT$ etc., which are well separated in weak coupling regime ($T>gT>\cdot\cdot$), in addition to the quark masses. So at finite temperature, strong magnetic field approximation implies that both conditions $eB>>T^{2}$ and $eB>>m^{2}$ are to be satisfied. To switch on the temperature in the vacuum propagator (5) in real-time formalism, the matrix propagator is diagonalized by the matric $U$ as

[TABLE]

where the matrix U(p) is given by

[TABLE]

with

[TABLE]

where $\beta$ is the inverse of temperature and $\mu$ is the chemical potential. In the present work we are working for baryonless medium ($\mu=0$), i.e. $n^{+}_{p}=n^{-}_{p}=n_{p}$. Plugging eq.(7) in eq.(6), we get the fermion propagator as

[TABLE]

For calculating the gluon self energy in an equilibrium medium, we need only the “11”-component of the matrix propagator expressed in eq.(11)

[TABLE]

[TABLE]

where the distribution function is given by

[TABLE]

The above description for fermionic propagator can be easily generalized to quarks of $f$-th flavor with which we are going to calculate the gluon self-energy.

2.2 Gluon self-energy in a hot QCD medium in presence of strong

magnetic field

As mentioned earlier, since we are working within SMFA so we may consider the strong coupling to run with the magnetic field only. For this purpose we closely follow the running coupling in [33], where the coupling runs with the momentum parallel and perpendicular to the magnetic field separately. In our case of magnetic field ($\vec{B}=B\hat{z}$), we will use the coupling dependent on the longitudinal component only because the energy of Landau levels for quarks in SMFA depend only on the longitudinal component of momentum. In fact, the coupling dependent on the transverse momentum does not depend on magnetic field at all. So in our calculation, the relevant coupling is given by [33]

[TABLE]

where

[TABLE]

In the above eq.(13), $M_{B}$ is taken $\sim~{}1$ GeV as an infrared mass and the string tension, $\sigma=0.18{\rm{GeV}}^{2}$.

For system in equilibrium, we need only the “11”-component of the gluon self-energy matrix calculated in real time formalism, which is given by

[TABLE]

where the factor 1/2 arises due to trace in color-space and the trace due to $\gamma$ matrices is given by

[TABLE]

Separating the momentum integration into longitudinal ($\parallel$) and transverse ($\perp$) components with respect to the magnetic field, the gluon self-energy can be factorized into $\parallel$ and $\perp$ components of momentum integration

[TABLE]

where the transverse component is given by

[TABLE]

It may be noted that in LLL approximation, the dependence of self-energy on the magnetic field is fully encapsulated in the transverse component whereas the longitudinal part carries no dependence on the magnetic field. We will now calculate the longitudinal component of the self-energy by decomposing eq.(16) into vacuum and thermal parts:

[TABLE]

where $(\Pi^{\mu\nu}_{\parallel})_{V}$ is the vacuum part, $(\Pi^{\mu\nu}_{\parallel})_{n}$ and $(\Pi^{\mu\nu}_{\parallel})_{n^{2}}$ are the thermal contributions due to single and double distribution functions, respectively. They are explicitly given by

[TABLE]

We will now calculate the vacuum term for the gluon self-energy.

2.2.1 Vacuum contribution ($T=0$, $eB\neq 0$)

The vacuum term in strong magnetic field can be calculated easily as it is similar to the calculation of self energy in vacuum without magnetic field except the fact that the dimension of the momentum integration is now reduced from $4$ to $2$. This dimensional reduction in fact removes the divergences usually encountered in 4-dimension, thus we do not need any regularization any more. Using the identity

[TABLE]

the real part of the vacuum term in the gluon-self energy has been calculated as

[TABLE]

where the form factor, $\Pi(k^{2})$ is given by

[TABLE]

Therefore the “00”-component ($\mu=\nu=0$) of the real part of vacuum term of the gluon self-energy (using the metric $g^{\mu\nu}_{\parallel}=diag(1,0,0,-1)$ is given by

[TABLE]

In the limit of massless quarks $(m_{f}=0)$, the gluon self-energy due to vacuum term in the static limit ($k_{0}=0$, $\vec{k}\rightarrow 0$) is given by the scale available to the magnetic field only in SMFA

[TABLE]

For the physical quark masses $(m_{f}\neq 0$), the vacuum term in the static limit ($k_{0}=0$, $\vec{k}\rightarrow 0$) vanishes

[TABLE]

2.2.2 Medium contribution

The (thermal) medium contribution to the gluon self-energy contains two terms: the first one (22) involves single distribution function and the second one (23) involves the product of two distribution functions. We will first consider the medium contribution due to the single distribution function only. Using the property of Dirac delta function, the gluon self-energy in eq.(22) is reduced to

[TABLE]

Taking $\mu=\nu=0$, the real part of “00”-component of $(\Pi^{\mu\nu}_{\parallel})_{n}$ becomes

[TABLE]

where the “00” component of $L^{\mu\nu}$ is

[TABLE]

and the other notations are

[TABLE]

After performing the $p_{0}$ integration we get from eq.(30)

[TABLE]

where we have defined

[TABLE]

and

[TABLE]

In the limit of massless quarks ($m_{f}=0$), the gluon self-energy in eq.(32) gets simplified into

[TABLE]

Using eq.(18) and multiplying the transverse component, $A(k_{\perp})$ from eq.(19), the contribution to the real part of self-energy from the component having single distribution function becomes

[TABLE]

which, in the static limit ($k_{0}=0$, $\vec{k}\rightarrow 0$) becomes

[TABLE]

However, for the physical quark masses $(m_{f}\neq 0$), the self-energy in eq.(30) reduces to, by putting $k_{0}=0$

[TABLE]

where the integrand, $I_{n}$, is given by

[TABLE]

and the distribution functions are given by

[TABLE]

Further taking the $k_{z}\rightarrow 0$ limit, the integrand, $I_{n}$ is simplified into

[TABLE]

Thus for the physical quark masses $(m_{f}\neq 0$), the contribution to the gluon self-energy having single distribution function in the static limit reduces to

[TABLE]

Finally the medium contribution to the gluon self-energy involving the product of two distribution functions given in eq.(23) does not contribute to the real-part of the gluon self-energy, i.e.

[TABLE]

We have thus so far evaluated the vacuum as well as medium contribution to one-loop gluon self energy, therefore we add them up to obtain the real-part of one-loop gluon self-energy in static limit for massless quarks

[TABLE]

and for the physical quark masses ($m_{f}\neq 0$)

[TABLE]

2.3 Debye screening mass in strong magnetic field:

The Debye screening manifests in the collective oscillation of the medium via the dispersion relation and is obtained by the static limit of the longitudinal part (“00” component) of gluon self-energy, i.e.

[TABLE]

Therefore, eq.(40) gives the very simple form for the square of the Debye mass for massless quarks, which is already derived in [34, 35]

[TABLE]

It shows that $m_{D}^{2}$ depends strongly on the magnetic field and is independent of temperature, thus the collective behavior of the medium gets strongly affected by the presence of strong magnetic field. However, for physical quark masses, the Debye mass is given by from eq.(41)

[TABLE]

which depends on both magnetic field and temperature. However, $m_{D}^{2}$ depends strongly on the magnetic field and the dependence on temperature is very weak and the screening mass becomes temperature-independent beyond a certain temperature.

Now, for the SMFA to be valid, we have to be careful in choosing the range of temperature and magnetic field. For example, for temperatures up to 300 MeV, the starting value of $eB$ has to be much higher than $0.09~{}{\rm{GeV}}^{2}$. Here we have taken the starting magnetic field to be $eB=10~{}m_{\pi}^{2}\sim 0.2~{}GeV^{2}$. However, for the upper bound on the magnetic field, the constraint comes from the heavy quark mass ($m_{Q}\gg\sqrt{eB}$) as discussed in the introduction. So, we have taken the highest magnetic field for charmonium states to be $eB=25~{}m_{\pi}^{2}$ which gives us $\sqrt{eB}\sim 0.7~{}GeV$. Thus, to see the variation of the Debye masses with the strong magnetic field, we have numerically calculated $m_{D}^{2}$ as a function of $eB$ (in units of $m_{\pi}^{2}$) for the temperature range T=200 - 300 MeV in Fig.1-a and noticed that $m_{D}^{2}$ is almost linearly increasing with $eB$ for smaller temperature. For higher temperatures, $m_{D}^{2}$ deviates slightly from the linearly increasing trend.

As we understood earlier in SMFA, the strongly magnetized thermal medium with massless quarks possesses only one scale available related to the magnetic field ($eB$) so by the dimensional arguments the square of the Debye mass is linear in $eB$ whereas for the the medium with physical quark masses, even in SMFA there is a weak competition between the dominant scale, $eB$ and much weaker scales, mass (m) and temperature (T) (rather their ratio, $m/T$) in the form of Boltzmann damping factor ($\exp(-m/T)$) as in eq.(44). This is seen in Fig.1-b, where a comparison of Debye masses with and without incorporating the quarks masses is made.

To see the temperature dependence of the Debye mass explicitly we have plotted $m_{D}$ with the temperature directly with increasing values of $eB=10m_{\pi}^{2}$, $15m_{\pi}^{2}$ and $25m_{\pi}^{2}$ in Fig.2. We can see how weakly the screening mass depends on temperature. It increases very slightly with temperature and beyond a point, the screening mass is practically a constant with magnetic field. The effect of temperature is slightly more pronounced for high magnetic field and low temperature.

Very recently the effects of a magnetic background on color-screening phenomena in QGP is also explored through the estimation of both the magnetic and electric screening masses by measuring the Polyakov loop correlators on the lattice for various temperatures [36]. There, it is found that the magnetic field induces an increase of both the magnetic and the electric screening masses and, to some extent, also the appearance of an anisotropy in Polyakov loop correlators. Both screening masses are found to increase linearly with the magnetic field and the influence of the magnetic field on both the masses is enhanced at lower temperatures and is asymptotically diminished in the higher temperature. Thus our aforesaid results on the Debye mass qualitatively agree with their findings for the electric screening mass, which is of interest to us for the screening of heavy quark potential. However, their lattice estimates for the electric screening masses are approximately larger by an order of magnitude than our results. This large difference may be attributed due to the nonperturbative effects, which is beyond the scope of this work.

3 Heavy quark potential in a hot QCD medium:

The derivation of potential between a heavy quark $Q$ and its anti-quark ($\bar{Q}$) either from EFT (pNRQCD) or from first principle QCD may not be plausible because the hierarchy of non relativistic scales and thermal scales assumed in weak coupling EFT calculations may not be satisfied and the adequate quality of the data is not available in the present lattice correlator studies, respectively, so one may use the potential model to circumvent the problem.

Since the mass of the heavy quark ($m_{Q}$) is very large, so the requirement: $m_{Q}\gg\sqrt{eB}\gg\Lambda_{QCD}$ and $T\ll m_{Q}$ is satisfied for the description of the interactions between a pair of heavy quark and anti-quark at finite temperature in the presence of magnetic field in terms of quantum mechanical potential. Thus we can obtain the medium-modification to the vacuum potential by correcting both its short and long-distance part with a dielectric function $\epsilon(\textbf{k})$ as

[TABLE]

where we have subtracted a $r$-independent term (to renormalize the heavy quark free energy) which is the perturbative free energy of quarkonium at infinite separation. The dielectric function is related to the “00”-component of effective gluon propagator in static limit as

[TABLE]

and $V(\textbf{k})$ is the Fourier transform (FT) of the Cornell potential. To obtain the FT of the potential, we regulate both terms with the same screening scale. However, different scales for the Coulomb and linear pieces were also employed in [37] to include non-perturbative effects in the free energy beyond the deconfinement temperature through a dimension-two gluon condensate.

At present, we regulate both terms by multiplying with an exponential damping factor and is switched off after the FT is evaluated. This has been implemented by assuming $r$- as distribution ($r\rightarrow$ $r\exp(-\gamma r))$. The FT of the linear part - $\sigma r\exp{(-\gamma r)}$ is

[TABLE]

After putting $\gamma=0$, we obtain the FT of the linear term $\sigma r$ as,

[TABLE]

The FT of the Coulomb piece is straightforward and is given by

[TABLE]

thus the FT of the full Cornell potential becomes

[TABLE]

The “00”-component of effective gluon propagator in static limit has been obtained with the help of “00”-component of one-loop gluon self energy. We have already calculated the “00” component of one-loop gluon self energy in presence of strong magnetic field at finite temperature in eq.(41), hence the “00”-component of effective gluon propagator in static limit is given by

[TABLE]

Therefore the real part of the static potential can be obtained by substituting the dielectric permittivity $\epsilon(\textbf{k})$ from eq.(46) and the Fourier transformation from eq.(50) into the definition of the potential (45)

[TABLE]

where the Coulombic and string term of the potential are given by (with the dimensionless quantity $\hat{r}=rm_{D}$)

[TABLE]

respectively. It is thus evident that the medium dependence in the potential enters through the Debye mass, which in turn depends on both temperature and magnetic field for physical quark masses and depends only on magnetic field for massless quarks. This gives a characteristic dependence of the potential on both temperature and magnetic field. The $r$-independent terms in the potential insures $V(r,T)$ to reduce to the Cornell potential in $T\rightarrow 0$ limit [38]. However, such terms could also arise naturally from the basic computations of real time static potential in hot QCD [39] and from the real and imaginary time correlators in a thermal QCD medium [40]. These terms in the potential are needed in computing the masses of the quarkonium states and to compare the results with the lattice studies. It is equally important while comparing our effective potential with the free energy in lattice studies.

Since we are exploring the effect of medium on the potential between $Q$ and $\bar{Q}$ in strong magnetic field approximation so we probe it by varying the strength of magnetic field ($eB$) from $10~{}m_{\pi}^{2}$ to $25~{}m_{\pi}^{2}$ (in Figure 3-a) at a temperature T=150 MeV. It is found that as the strength of the magnetic field increases the potential becomes stronger. To see the competition between the magnetic field and temperature, we have plotted the potential in Figure 3-b in a hotter medium (T=300 MeV). As we have seen earlier in Figure 1-a that the (square) Debye screening mass increases very little with temperature, here also the potential changes little as compared to Figure 3-a. This small dependence on temperature stems from SMFA as we have observed in the Debye screening mass.

Usually potential model studies are limited to the medium-modification of the perturbative part of the potential only where it is assumed that the string-tension vanishes abruptly at the deconfinement point. Since the phase transition in QCD for physical quark masses is found to be a crossover [41], so the string tension may not vanish at the deconfinement temperature. This issue, usually overlooked in the literature where only a screened Coulomb potential was assumed above $T_{c}$ and the linear term was neglected, is certainly worth for an investigation. To see the effect of the linear term on the potential, in addition to the Coulomb term, we have plotted the potential (in Fig.4-a) with ($\sigma\neq 0$) and without string term ($\sigma=0$) in a magnetic field $eB=10~{}m_{\pi}^{2}$. As we know already in vacuum (T=0), the inclusion of the linear term makes the potential in short-distance interaction less attractive and in long-distance interaction the linear term makes the potential more repulsive, compared to the Coulomb term alone. However, the medium modification causes the linear term attractive and overall the medium modifications to both Coulomb and string term makes the potential more attractive (seen in Figure 4-a) compared to the vacuum potential.

To see the effect of the scale (Debye mass) at which the screening takes place on both the linear and Coulombic term we have plotted the potential at a larger magnetic field, $eB=25~{}m_{\pi}^{2}$ in Fig. 4-b, where we found that the increase of the scale (screening mass) makes the linear term less attractive, compared to the lower scale ($eB$). To understand the observations in Figure 4, we have probed the range of interactions, viz. short-range ($r=0.2$ fm), intermediate ($r=0.5$ fm) and long-range ($r=1$ fm) interactions of $Q\bar{Q}$ potential as a function of magnetic field ($eB$) in figure 5 and found that only the long range interaction ($r=1$ fm) has been affected noticeably. Overall observation is that as the strength of the magnetic field increases the long range QCD force becomes more and more short range, thus implying that magnetic field facilitates early dissolution of $Q\bar{Q}$ states.

It is important to mention here that we have not observed any anisotropy in our potential with respect the direction of magnetic field, which is expected as a common sense because the magnetic field breaks the translational invariance of space. In our perturbative framework too, at the starting point the quark propagator in strong magnetic field approximation gets factorized into functions involving the parallel and perpendicular components with respect to the direction of magnetic field ($\vec{B}=B\hat{z}$) and so obviously the gluon self-energy is decoupled into parallel and perpendicular components. But while thermalizing the quark propagator and gluon self-energy at finite temperature in strong magnetic field through the distribution function, there is no scope to introduce the momentum anisotropy in distribution function because the quarks dispersion relation is restricted to the LLL only ($E_{0}(p_{z})=\sqrt{p_{z}^{2}+m_{f}^{2}}$), i.e. only the longitudinal component of momentum is present, hence no anisotropy arises between transverse and longitudinal components. One of us had recently derived an anisotropic heavy quark potential [20, 21] in perturbative thermal QCD, where the anisotropy in the potential in the coordinate space had arisen from the manifested momentum anisotropy with respect to the direction of anisotropy in the distribution function.

However, recently a novel magnetic field-induced anisotropic behavior was first observed in the heavy quark potential in the longitudinal-traverse plane with respect to the direction of $\vec{B}$ [25]. Later the the set-up is extended to measure Polyakov loop correlators on the lattice to extract the potential for both zero and finite temperature in place of earlier Wilson loop expectation values used at zero temperature, for different orientations with respect to $\vec{B}$ [26]. The reason of anisotropy arises due to the averaging of both the Wilson loop expectation value and Polyakov loop correlators differently for different orientations with respect to the magnetic field.

3.1 Dissociation of Heavy Quarkonia in magnetic field

In this section, we shall discuss the dissociation of charmonium and bottomonium states due to an external strong magnetic field in a hot QCD medium. The concept of dissociation temperature becomes irrelevant here because the scale at which the collective oscillations develop depends only on the magnetic field albeit we are considering a hot QCD medium because in strong magnetic field approximation ($eB\gg T^{2}$), the scale at which the collective oscillation sets in is associated with magnetic field only because $eB$ is the most dominant scale in strong magnetic field limit (if the partons are assumed massless), not the thermally generated scales. This in turn makes the potential to depend only on the magnetic field through the dependence of Debye mass on the magnetic field. Thus, it makes sense here to discuss the dissociation of quarkonium states due to the magnetic field only as far as SMFA is valid.

As we know that in the presence of medium, the potential between a heavy quark ($Q$) and its anti-quark ($\bar{Q}$) will be screened, as a result if the screening is strong enough, the potential becomes too weak to form the resonance. Thus we can argue that the quarkonium states will be dissolved in a medium if the Debye screening radius, $r_{D}$ (=$\frac{1}{m_{D}}$) in a given medium is smaller than the bound state radius of a particular resonance state then the medium inhibits the formation of the particular resonance and $Q$ and $\bar{Q}$ will be dissolved into the medium. Since the screening mass in strong magnetic field increases with the magnetic field therefore the (critical) magnetic field at which the $Q\bar{Q}$ potential becomes too feeble to hold $Q\bar{Q}$ together becomes smaller for the excited states. We can thus estimate the lower limit of critical magnetic field for various charmonium and bottomonium states by the criteria: $\sqrt{\langle{r^{i}}^{2}\rangle}={r_{D}}(B_{d}^{i})$, i.e., for magnetic field larger than ${B_{d}}^{i}$, the $i$-th quarkonium states cease to exist. For example, $J/\psi$ will be dissociated at $eB=14m_{\pi}^{2}$ and its excited state, $\psi^{\prime}$ is dissociated at smaller magnetic field $eB=m_{\pi}^{2}$ whereas $\Upsilon$ will be dissociated at $eB=130m_{\pi}^{2}$ and $\Upsilon^{\prime}$ is dissociated at smaller magnetic field $eB=13m_{\pi}^{2}$.

To understand the in-medium properties of the quarkonium states quantitatively, one need to solve the Schrödinger equation with the medium-modified potential, $V(r;B,T)$. There are some numerical methods to solve the Schrödinger equation either in partial differential form (time-dependent) or eigen value form (time-independent) by the finite difference time domain method (FDTD) or matrix method, respectively. In the later method, the stationary Schrödinger equation can be solved in a matrix form through a discrete basis, instead of the continuous real-space position basis spanned by the states $|\overrightarrow{x}\rangle$. Here the confining potential V is subdivided into N discrete wells with potentials $V_{1},V_{2},...,V_{N+2}$ such that for $i^{\rm{th}}$ boundary potential, $V=V_{i}$ for $x_{i-1}<x<x_{i};~{}i=2,3,...,(N+1)$. Therefore for the existence of a bound state, there must be exponentially decaying wave function in the region $x>x_{N+1}$ as $x\rightarrow\infty$ and has the form:

[TABLE]

where, $P_{{}_{E}}=\frac{1}{2}(A_{N+2}-B_{N+2})$, $Q_{{}_{E}}=\frac{1}{2}(A_{N+2}+B_{N+2})$ and, $\gamma_{{}_{N+2}}=\sqrt{2\mu(V_{N+2}-E)}$. The eigenvalues can be obtained by identifying the zeros of $Q_{E}$. Using this method, we have found that $J/\psi$ and $\Upsilon$ is dissociated at $eB=5m_{\pi}^{2}$ and $eB=50m_{\pi}^{2}$, respectively.

Though the dissociation magnetic fields, obtained from the two different methods apparently look different, its easy to see that qualitatively they are similar. Using both methods, we found that the dissociation magnetic field for $\Upsilon$ is roughly an order of magnitude greater than the dissociation magnetic field for $J/\psi$. Even though their absolute value obtained from the two different methods differ, they lie in the same ball park, which is $\sim$ 10 $m_{\pi}^{2}$ for $J/\psi$ and $\sim$ 100 $m_{\pi}^{2}$ for $\Upsilon$.

4 Conclusions

In this article, we have explored the effects of strong and homogeneous magnetic field on the properties of quarkonium states. For that purpose we have derived the potential between a heavy quark and its anti-quark by the medium corrections to both Coulomb and linear term of $Q\bar{Q}$ potential at T=0, unlike the medium correction to the Coulomb term alone. Although the medium considered is thermal but due to strong magnetic field approximation, all other scales present in the thermal medium becomes irrelevant as the scale related to magnetic field dominates over other. This is exactly what happens in the collective oscillation of the medium in the form of Debye mass. In fact, the Debye mass becomes completely independent of temperature for massless quarks and depends very weakly on temperature for massive quarks. However, beyond a certain temperature, the dependence is so weak that it is almost insignificant. As a result the heavy quark potential mainly depends on the magnetic field with a very feeble dependence on the temperature. This is expected as the effect of the medium on the potential enters through the Debye mass. In particular the long distance part of the potential gets significantly affected, whereas the short distance part is mildly affected.

We have then studied the dissociation of quarkonium states in a medium. Since the potential in SMFA depends mainly only on the magnetic field thus we have discussed the dissociation of quarkonium states due to the magnetic field only. We have estimated the critical value of magnetic field beyond which the resonance does not form in two methods. The first one gives a lower limit of critical magnetic field for both charmonium and bottomonium states at which the Debye screening radius becomes smaller than the bound state radius of a particular resonance state. The other one comes from the consideration of the binding energies of a specific state obtained from the energy eigenvalues of the Schrödinger equation. In brief, $J/\psi$ is dissociated at $eB\sim$ 10 $m_{\pi}^{2}$ and $\Upsilon$ is dissociated at $eB\sim$ 100 $m_{\pi}^{2}$.

5 Acknowledgments

We are thankful to Aritra Bandyopadhyay for a fruitful discussion during this work. Bhaswar is thankful to the Ministry of Human Resource Development, Government of India for the financial assistance.