J/ψ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J/\psi $$\end{document} and ηc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta _c$$\end{document} in asymmetric hot magnetized nuclear matter: a unified approach of Chiral SU(3) model and QCD sum rules

We investigate the mass-shift of vector channel J/ψ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J/\psi $$\end{document} and pseudoscalar channel ηc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta _c$$\end{document} in strongly magnetized asymmetric nuclear medium at finite temperature using the conjunction of chiral SU(3) model with the QCD sum rules. The magnetic field dependence of scalar gluon condensate αsπGμνaGaμν\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left\langle \frac{\alpha _{s}}{\pi } G^a_{\mu \nu } {G^a}^{\mu \nu } \right\rangle $$\end{document} as well as the twist-2 gluon condensate αsπGμσaGaνσ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left\langle \frac{\alpha _{s}}{\pi } G^a_{\mu \sigma } {{G^a}_\nu }^{\sigma } \right\rangle $$\end{document} calculated from chiral SU(3) model, are implemented in QCD sum rules to calculate the magnetic field dependence of J/ψ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J/\psi $$\end{document} and ηc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta _c$$\end{document} meson masses. The effects of constant external magnetic field at finite density and temperature of the medium are found to be appreciable in symmetric and asymmetric nuclear matter. The results of the present investigation may be helpful to understand the experimental observables arising from the Compressed Baryonic Matter produced in asymmetric non-central heavy ion collision experiments.


Introduction
The investigation of the properties of hadrons in hot and dense matter under the effect of strong external magnetic field is a challenging area of research in non-perturbative QCD regime and of very importance from non-central Heavy-Ion Collisions (HICs) point of view which aim at understanding strong interaction physics. Due to highly asymmetric nature (different numbers of neutrons and protons) of heavy-ions, it is very imperative to include the isospin asymmetry while studying the properties of mesons in nuclear matter at finite temperature. In RHIC and LHC experiments, strongly interacting matter is produced at high temperature and very low baryon density. In addition, in non-central HICs, huge magnetic fields are believed to be produced and are estimated to be of the order of eB ∼ 15m 2 π (5 × 10 19 gauss) at Large a e-mail: rajesh.sism@gmail.com b e-mails: iitd.arvind@gmail.com; kumara@nitj.ac.in Hadron Collider (LHC) CERN [eB ∼ 2m 2 π (6.5 × 10 18 gauss) at Relativistic Heavy Ion Collider (RHIC) BNL] [1][2][3] and it decrease rapidly after the collision as the residuals move away from the collision zone. But the interaction with remnant matter leads to the generation of induced circular currents, according to Lenz's rule, which further produce the induced magnetic field along the direction of external magnetic field. These currents affect the electric conductivity of the medium, which slows down the relaxation time, hence decay rate of external magnetic field [4][5][6][7][8][9][10]. The above effect is known as chiral magnetic effect [2,[11][12][13]. The external magnetic field modifies the QCD vacuum properties such as gluon condensates [14,15] and quark condensates [16][17][18][19][20][21] and due to this effect the interaction of heavy quarkonia [9,10,22,23] with strong magnetic field leads to the modification of measurable observables [7,10,24]. In Ref. [25], it is shown that the electrical conductivity of the medium is large even at temperature near T c . The time evolution of magnetic field is a debatable topic and to have a proper estimate of relaxation time, further investigation in the solutions of magnetohydrodynamics equations are required. With the construction of CBM experimental Facility for Antiproton and Ion Research (FAIR) and Nuclotron-based Ion Collider Facility (NICA) at Dubna, Russia, we expect significant progress to understand the in-medium modifications of hadrons at high baryonic density and moderate temperature.
Strong magnetic fields also exist in the astrophysical compact objects namely magnetars. Magnetars are the neutron stars having very strong magnetic field [26,27]. It is suggested that the soft gamma-ray repeaters (SGRs) and anomalous X-ray pulsars (AXPs) reside in the category of magnetars [28,29]. The magnetic field at the interior and exterior of magnetars are believed to be as large as ∼ 10 18 gauss and ∼ 10 14 gauss [30], respectively. The huge magnetic field in the centre is suggested to be due to the dense quark matter magnetic core [31,32] and has been studied substantially in literature [33][34][35][36][37]. Also, note that the magnetic field pro-duced in the heavy-ion collisions is much above the field produced in magnetars and is comparable to QCD scale and hence, will play an important role in understanding the QCD phase diagram. Such strong magnetic field leads to many interesting physical phenomenon, for example, chiral separation effect [11,[38][39][40], heavy quark transport phenomenon [41], magnetic catalysis of chiral symmetry breaking [11], inverse magnetic catalysis [11], medium modification of heavy flavoured mesons [9,10,23], magnetic inhibition at finite temperature and density [42].
In the present paper, we will study the properties of ground state S-wave charmonium, J/ψ and η c , at finite density and temperature under strong external magnetic field using QCD sum rules and chiral SU (3) model. QCD sum rules are used to evaluate hadronic properties, such as mass in term of condensates to get valuable information about the structure of QCD regime [72,75,76]. The mass modifications of heavy quarkonia, notably charmonium, and bottomonium have been investigated using QCD sum rules [10,22,58,[77][78][79][80][81]. This model remarkably predicted the mass of η c in mixed states [82]. Also, within QCD sum rules, the in-medium mass and decay width of open charm mesons have been calculated [83]. The chiral SU (3) model has been used vastly to study the nuclear matter [58], kaons and antikaons in nuclear and hyperonic matter [84], finite nuclei [55] and in-medium properties of vector mesons [84,85]. This model is also generalized to SU (4) to investigate the in-medium properties of D mesons [56][57][58][59]. To study the in-medium mass of J/ψ and η c mesons, we consider the contributions of the scalar gluon condensates, α s π G a μν G a μν and twist-2 tensorial gluon oper-ator, α s π G a μσ G a ν σ up to dimension four [76]. The scalar gluon condensate as well as the twist-2 gluon operator in the asymmetric nuclear medium under the effect of finite magnetic field are calculated from chiral SU (3) model [55] as will be discussed in detail in the next section. In this article, we use mean field approximation, which is a nonperturbative relativistic approach to solve approximately the nuclear many-body problem. In this approximation, all the meson fields are considered as classical fields, hence only the vector and scalar fields contribute to nucleon-meson interaction Lagrangian term, as expectation value of the other mesons is zero [56]. The quantum and thermal fluctuations near chiral phase transition (beyond mean field approximation) has also been studied in literature by using the Functioal Renoramalization Group (FRG) approach [86,87] and Polyakov extended quark meson (PQM) model [88,89]. It was observed that the critical point moves towards the higher chemical potential and lower temperature region [90,91].
The phase diagram become more stabilized after the inclusion of fluctuations [87]. Interest in the in-medium properties of charmonium was ignited when Matsui and Satz proposed that the decrease in the yield of J/ψ state in HICs due to color screening effect should be considered as a probe of the production of Quark Gluon Plasma (QGP) [92]. Imperative results in favor of J/ψ suppression were observed at CERN SPS and RHIC experiment [93][94][95]. The statistical recombination of primordially produced charm quark pairs may lead to the increase in the yield of J/ψ mesons and this picture is more important at LHC energies [96][97][98]. Charmonium is a bound state of charm and an anti-charm quark. The in-medium properties of charmonia are modified in the nuclear matter through gluon condensates [58,99]. The gluon condensates modify feebly, hence the obtained mass-shift of lowest charmonium states, J/ψ and η c is very small in nuclear medium [77]. In Ref.
[99], the in-medium mass of charmonium states have been studied under linear density approximation in the nuclear medium using QCD second order stark effect. A significant mass shift is observed for the excited charmonium states (ψ(3686) and ψ(3770)), but a small mass shift for the ground state J/ψ(3097). In Ref. [78], mass shift of J/ψ and η c is studied using QCD sum rules through the medium modifications of gluon condensates in linear density approximation. These studies are applicable up to nuclear saturation density only. In Refs. [58,59], gluon condensates evaluated in nonlinear chiral SU (3) model, in terms of dilaton field, were used as input in QCD second order stark effect and in QCD sum rules, respectively, to evaluate the masses of charmonium at finite density and temperature. In Refs. [100][101][102], the effect of temperature on the J/ψ and η c was also studied in deconfined phase and observed that the heavy quark bound states can survive in the deconfined plasma. In addition to this, the in-medium properties of excited charmonium states, χ c0 and χ c1 are also studied in the literature [78,103]. Recently, in Ref. [10], authors have investigated the masses of lowest charmonium states J/ψ and η c in strong magnetic field using QCD sum rules. In Ref. [22], magnetically induced mixing between J/ψ and η c have been calculated. Also, chiral SU (3) model along with QCD second order stark effect was extended in Ref. [104] to study the effect of external magnetic field on the masses of J/ψ, ψ(3686) and ψ(3770) mesons at finite density and zero temperature. D, B, and ρ mesons are also studied in the presence of strong magnetic field [9,23,105,106]. Due to the attractive potential of J/ψ and η c , it is also predicted that these mesons can form a bound state with nucleons [107,108].
We organize the paper as follows. In Sect. 2, we will present chiral SU (3) model to describe the hot asymmetric nuclear matter in the presence of an external magnetic field. The scalar density as well as the number density of the charged proton has contributions from the Landau energy levels, whereas, the uncharged neutrons do not have contributions from Landau energy levels. In Sect. 3, we describe the QCD sum rules to calculate the in-medium mass-shift of charmonium. We discuss the results of the present investigation in Sect. 4, and finally in Sect. 5 summary of this work will be given.

The hadronic chiral SU L (3) × SU R (3) model
We use an effective field theoretical approach to describe hadron-hadron interactions, based on the broken scale invariance [55,56,84] and non-linear realization of chiral symmetry [109][110][111] in the presence of external magnetic field at finite density and temperature. In this model, the hadronhadron interactions are expressed in terms of the scalar fields σ , ζ , δ, χ and vector fields ω and ρ. The δ and ρ fields are introduced to incorporate the effect of isospin asymmetry in hadronic medium. Also, the dilaton field χ is proposed to express the scale symmetry breaking, leading to a nonvanishing trace of the energy-momentum tensor [55]. The Lagrangian density of this model under mean-field approximation is given as Individually, and In Eq. (1), L kin represents the kinetic energy term, L N M is the nucleon-meson interaction term, where S and V represents the spin-0 and spin-1 mesons, respectively. Here, the effective mass of nucleons is given as In above, g σ i , g ζ i and g δi represent the coupling strengths of nucleons (i = p, n) with σ , ζ and δ fields respectively and τ 3 is the third component of isospin. The term L vec of Eq. (1) generates the mass of vector mesons through the interactions with scalar mesons and contains the quartic self-interaction terms, L 0 describes the spontaneous chiral symmetry breaking, and L S B describes the explicit chiral symmetry breaking.
In grand canonical ensemble, the partition function for nuclear system is given by [112] where β = 1/T andĤ, μ i andN i are the Hamiltonian density operator, chemical potential and number density operator, respectively. The thermodynamical potential, at given temperature T reads By inserting mean field Hamiltonian density [112] in terms of Lagrangian density, the thermodynamic potential , per unit volume, V in zero magnetic field can be expressed as where the sum runs over neutron and proton, γ i is the spin degeneracy factor for nucleons, 3 ρ are the effective single particle energy of nucleons and effective nucleon chemical potential, respectively. In addition, vacuum potential energy, V vac is subtracted in order to get zero vacuum energy.
In the presence of magnetic field, the Lagrangian density given by Eq. (1) modifies to where In above ψ i is a wave function corresponds to the i th nucleon and the second term represents the tensorial interaction with the electromagnetic field tensor, F μν . Also, k i and μ N are the anomalous magnetic moment of i th nucleon and nuclear magneton, given as μ N = e 2m N , respectively, where m N is the vacuum mass of the nucleon. We choose the magnetic field to be uniform and along the Z -axis, and hence the vector potential A μ = (0, 0, Bx, 0). In the following, we discuss the interaction of charged and uncharged particles with external magnetic field.

Charged proton in magnetic field
In the presence of uniform magnetic field, due to the charged nature of proton, Lorentz force comes into picture. Hence, the transverse momenta of proton with an electric charge q p are confined to discrete Landau levels, ν, with, k 2 ⊥ = 2ν|q p |B, where ν ≥ 0 is an integral quantum number [113]. Thus, the volume integral converts into line integral where k is the momenta along the direction of magnetic field. Also, the summation represents a sum over the discrete orbital angular momentum, n of proton in the perpendicular plane. The orbital quantum number is related to Here, the quantum number s is +1 for spin up and −1 for spin down protons. In addition to this, the effective single particle proton energy also gets quantized [114] and is given bỹ Under these changes, the first term of Eq. (9), for proton will be

Uncharged neutron in magnetic field
For uncharged neutron, there is no Landau quantization in the presence of external magnetic field, hence 3 remains unchanged [113]. The first term of thermodynamic potential, , per unit volume given by Eq. (9), for neutron, will be whereẼ n s is the effective single particle energy of neutron in the presence of magnetic field and is given bỹ The net thermodynamical potential in the presence of external magnetic field can be written as Now, by minimizing the thermodynamical potential /V of the nuclear system, the coupled equations of motion of the non-strange meson field σ , the strange scalar meson field ζ , the scalar iso-vector meson field δ, the vector meson field ω, the vector-isovector meson field ρ, and the dilaton field χ , are determined and given as and respectively. In above, m π , m K and f π , f K denote the mass and decay constant of π , K mesons, respectively and the other parameters k 0 , k 2 and k 4 are fitted to reproduce the vacuum mass of σ , ζ and χ meson and the remaining constants k 1 is fixed to produce the effective nucleon mass at saturation density around 0.65m N and k 3 is the constraint by η and η masses. Furthermore, ρ s i and ρ v i represent the scalar and vector/number densities of ith nucleon respectively.
In the presence of the magnetic field, the vector density as well as the scalar density of proton can be extracted from Eqs. (21) to (18) and given as [33,115] and respectively, wherem p is the effective mass under the effect of magnetic field, defined as Similarly, for the uncharged neutron, the number and scalar densities are given by [33,115] and respectively. In above, f p k,ν,s ,f p k,ν,s , f n k,s andf n k,s represent the finite temperature distribution functions for particles and antiparticles for proton and neutron, and given as If we neglect the effect of magnetic field in Eq. (10), the vector and scalar densities of the nucleons are modified and given as and respectively. As discussed earlier, the effect of isospin asymmetry is introduced by incorporating δ and ρ fields having dependence on the isospin asymmetry parameter, η = (ρ v n − ρ v p )/2ρ N , which is the measure of the abundance of one isospin component over other.
In the present investigation, gluon condensates will be used to calculate the effective mass of charmonia within QCD sum rules. The scalar gluon condensate α s π G a μν G a μν and the twist-2 tensorial gluon operator, α s π G a μσ G a ν σ , given by Eqs. (39) and (40) respectively in the following are derived in terms of medium modified σ , ζ , δ and dilaton field χ within chiral SU (3) model as explained below. Within this model, gluon condensates are obtained from the energy-momentum tensor in terms of χ field [58]. We start with energy-momentum tensor where L χ represents the Lagrangian density term, which is written in the model to incorporate scale breaking property of QCD and is given by To obtain the trace of above energy-momentum tensor, we multiply it by g μν , The energy-momentum tensor, T μν and its trace, T μ μ in massless QCD can be written as [116,117] and respectively. In Eq. (36), G 2 is twist-2 gluon operator.
By comparing Eq. (35) with (37), scalar gluon condensate, G 0 can be written as If we take the effect of finite quark mass term [58], the modified scalar gluon condensate, G 0 is expressed as The comparison of energy-momentum tensor in chiral model with QCD [58], gives the expression for twist-2 gluon operator and is given by The value of d may be taken from QCD beta function, β QC D at the one loop level [55], with N c colors and N f flavors, Here, the first term in parentheses arises from the selfinteraction of the gluons and the second term is proportional to N f , is the contribution from quark pairs. This equation suggests the value d = 6/33 for three flavors and three colors.
Here, q = (ω, q), represents four momentum vector and symbol c stands for the pseudoscalar (P) and vector (V ) mesons. Each current is defined as j P =ψ q γ 5 ψ q and j V μ = ψ q γ μ ψ q with ψ q being the quark operator.
We can render the current-current correlation in the region of positive and large Q 2 = q 2 − ω 2 through an OPE. Hence, left-hand side of above equation can be written as [118] (q) = n W n (q) Ô n .
In above, theÔ n are operators of dimension n, and W n are the perturbative Wilson coefficients. By expressing correlation function as polarization func-tion˜ J , which reduces to vacuum polarization function in the limit ρ N → 0 [77], the nth moment can be expressed as at a fixed Q 2 0 = 4m 2 c ξ . Using operator product expansion, the moment M J n can be written as [76] where A J n (ξ ), a J n (ξ ), b J n (ξ ) and c J n (ξ ) are the Wilson coefficients [72,76] and ξ is the normalization scale. The common factor A J n results from the bare loop diagram. The coefficients a J n take into account perturbative radiative corrections, while the coefficients b J n are associated with the scalar gluon condensate through As described earlier, the contribution of the scalar gluon condensate is incorporated through the σ , ζ , δ and dilaton field χ within the chiral SU (3) model. Using Eq. (39), the above equation can be rewritten in terms of these fields as The coefficients c J n are associated with the value of φ c , which gives the contribution from twist-2 gluon operator, φ c appearing in Eq. (45) is defined by where twist-2 gluon operator α s π G a μσ G a ν σ is given by Eq. (40). Explicitly in terms of σ , ζ , δ and dilaton field χ , φ c will be, The m c and α s parameters are the running charm quark mass and running coupling constant and are ξ dependent [72], and are given by where, m c ≡ m c ( p 2 = −m 2 c ) = 1.26 GeV [75], and with, α s 4m 2 c 0.3 and Q 2 0 = 4m 2 c ξ [72]. The in-medium mass of the J/ψ and η c charmonium states in terms of M J n can be written as In the next section, we shall discuss the results for the J/ψ and η c meson mass modifications in asymmetric nuclear matter in the presence of an external magnetic field, B, at finite temperature.

Results and discussions
The effects of external magnetic field in hot and dense nuclear matter on the properties of J/ψ and η c in the present work are calculated through the in-medium behaviour of scalar gluon condensate α s π G a μν G a μν and the twist-2 tensorial gluon operator, α s π G a μσ G a ν σ . We have divided our discussion in three sections. In Section A, we will discuss the effect of external magnetic field on in-medium behaviour of scalar fields σ , ζ , δ and χ , through which medium dependence of gluon condensates will be evaluated. The in-medium effects on gluon condensates are discussed in Section B. Finally, in Section C, medium modification of J/ψ and η c meson will be discussed.
4.1 Scalar fields σ , ζ , δ, and χ in hot magnetized nuclear medium By solving coupled system of non-linear equations (18-(23) for σ , ζ , δ, ω, ρ and χ , temperature and density dependence of these fields in the presence of magnetic field is calculated. These equations have the dependence on vector density, ρ v i   and scalar density, ρ s i of the nucleons, which further have magnetic field dependency as can be seen from Eqs. (24) to (28). Various parameters used in the present work are listed in Table 1.
In Fig. 1, we show the variation of scalar fields σ, ζ, δ, and χ as a function of temperature T at magnetic field B = 0, 5m 2 π and 7m 2 π (1m 2 π = 2.818 × 10 18 gauss = 5.48 × 10 −2 GeV 2 ), and nucleon density, ρ N = 0. As can be seen from Fig. 1, at ρ N = 0, the effects of magnetic field are more appreciable above certain high value of temperature. For given temperature, the magnitude of scalar fields σ and ζ and the dilaton fields χ are observed to decrease with increase in magnetic field. Further, in Figs. 2 and 3, we plot these scalar fields individually as a function of temperature T at different magnetic field strength, eB and for finite nucleon density. We show the results at nucleon densities ρ N = ρ 0 (left column) and 4ρ 0 (right column) and isospin asymmetry, η = 0, 0.3 and 0.5. Here, e and ρ 0 are the electrostatic unit of charge and nuclear saturation density respectively.
In Fig. 2, we observe that at zero and non-zero values of magnetic field, for finite value of nucleon density ρ N , the magnitude of scalar fields σ and ζ increases with increase in temperature and then start decreasing after certain value. For example, in symmetric nuclear matter, at ρ N = ρ 0 and B = 0, the magnitudes of both σ and ζ fields increase upto temperature T = 145 MeV and then start decreasing with further increase in temperature. For finite value of magnetic field B, the magnitude of both scalar fields decreases as compared to B = 0, but follows the same functional dependence with respect to temperature. The medium modified values of σ and ζ fields in symmetric nuclear matter are listed in Table 2. The decrease in the magnitude of σ and ζ fields as a function of magnetic field strength will result in decrease in the mass of nucleons (see Eq. 6). Thus, finite magnetic field strength at high temperature will further support the restoration of chiral symmetry and may result in decrease in critical temperature. Moreover, with the rise of temperature, the value of thermal distribution function given in Eqs. (29) and (30) decreases, which results in decrease in the scalar density given by Eqs. (25) and (28). In addition, at μ * i = 0, i.e. for finite density, there are contributions from higher momenta. Thus, the variation of σ and ζ field with respect to temperature reflects the competing effects between the thermal distribution functions and the contributions from higher momenta states [84].
Comparing the behaviour of scalar fields at finite asymmetry of the medium with the symmetric situation, we observe that isospin asymmetric effect on scalar field σ are more appreciable than ζ field. This is related to the fact that the σ field contains non zero isospin quark content, while ζ field quark content is independent of isospin. Therefore, former is more sensitive to non-zero isospin asymmetry parameter, η, than the latter (ζ has little dependence as it is solved along with other scalar fields). For given density and temperature, the decrease in magnitude of scalar fields is observed to be small at η = 0.3 as compared to η = 0, as we increase the external magnetic field strength from zero to finite value. However, as we move to η = 0.5, below T ∼ = 80 MeV, the magnitude of scalar field σ increases with increase in the magnetic field strength, whereas, above this temperature, the trend is opposite. The reason behind this crossover behaviour is the fact that in pure neutron matter, η = 0.5, at a particular value of temperature, neutron scalar density, ρ s n decrease more slowly as compared to proton scalar density, ρ s p . The net effect of ρ s n and ρ s p results in the above behavior of σ and ζ fields. In the present investigation, the variation of scalar fields and nucleon mass at finite density, temperature and external magnetic field are in accordance with the calculations done under non-linear Walecka model [26]. In this article, authors have studied the properties of stellar matter at finite density, tempearature and uniform magnetic field.
The scalar isovector field δ having contribution due to isospin asymmetry of the medium is plotted in Fig. 3. It is observed that at η = 0, the δ field is zero for magnetic field B = 0, but non-zero for finite magnetic field strength. This is because the δ field is determined by the difference between the scalar densities of neutron and proton. For zero magnetic field, the value of δ field is zero as scalar density of proton is equal to scalar density of neutron (see Eqs. 31, 32). On the other hand, in the presence of magnetic field, due to charged nature of proton, Landau quantization occurs, this leads to the inequality ρ s p = ρ s n (see Eqs. 24, 28), hence non-zero value of the δ is obtained. In the same figure, the dilaton field χ is also plotted, which is introduced in the chiral SU (3) model to preserve broken scale invariance property of QCD at tree level [119]. It is observed that the magnitude of dilaton field χ shows behaviour similar to σ and ζ fields. It first increases up to certain temperature T , then decreases with further increase in temperature. It is observed from the figure that the dilaton field χ varies least as compared to other mesonic fields. Therefore, frozen glueball limit is assumed in many research papers on chiral SU (3) model [9,55]. In the lower density region, the effects of magnetic field on χ

(d)
field are very less, whereas at high density these become significant.

In-medium gluon condensate and twist-2 gluon operator
In Fig. 4, we show the variations of scalar gluon condensate G 0 (Eq. 39) as well as the twist-2 gluon operator G 2 (Eq. 40) with temperature for different values of magnetic field, nucleon density and isospin asymmetry parameter. The scalar gluon condensate, G 0 consists of two terms. The first term has the dependence on the fourth power of χ field, whereas the second term contributes the finite quark mass effects from the energy-momentum tensor, which have the σ and ζ field dependence. On the other hand, the twist-2 gluon operator, G 2 has contribution of three terms, which further have the dependency on scalar fields. Although, the condensates G 0 and G 2 depends upon all scalar fields (see Eqs. 39,40), but the observed behaviour resembles χ field due to its fourth power dependence. For better understanding of condensate's behaviour, the different terms of G 0 and G 2 are plotted in Fig.  5 considering finite quark mass term contribution to scalar gluon condensates, in symmetric nuclear medium, at temperature T = 0 (100) MeV, the values of G 0 , in units of 10 −2 GeV 4 , are observed to be 1.842 (1.871), 1.809 (1.838) and 1.764 (1.800) at eB = 0, 5m 2 π and 7m 2 π , respectively. However, in massless quark limit (Eq. 38) these values changes to 2.177 (2.210), 2.139 (2.172) and 2.089 (2.130). In addition, the tensorial gluon operator, G 2 varies very less in high density as compared with low density as a function of temperature. This behaviour reflects the competing effects between the second and third term of Eq. (40) plotted in Fig. 5. The crossover behaviour of condensates in highly asymmetric nuclear medium reflects the χ field, which is solved along with the σ , ζ and δ fields. For example, in symmetric nuclear matter, at ρ N = 4ρ 0 and temperature T = 0 (100) MeV, the values of χ field are observed to be 402.2 (403.5), 400.4 (402) and 398.1 (400) MeV for eB = 0, 5m 2 π and 7m 2 π respectively and for η = 0.5, these values changes to 403.2 (404.1), 403.8 (403.9) and 404. 3 (403.3). In vacuum, the twist-2 gluon operator vanishes due to the Lorentz symmetry. But in the current investigation, the presence of finite density, temperature and magnetic field breaks the Lorentz symmetry, hence nonvanishing value of G 2 is obtained [77,117]. The observed behaviour of G 0 and G 2 as a function of temperature at zero magnetic field and nucleon density is consistent with the condensates behaviour observed in Ref. [117].

In-medium masses of J/ψ and η c
In this section, we calculate the in-medium mass shift of J/ψ and η c mesons, using medium modified scalar and tensorial gluon condensate, by applying QCD sum rules. In this application, we use the moment in the range 5 ≤ n ≤ 9 and 7 ≤ n ≤ 10 for J/ψ and η c mesons respectively to get the mass in stable region. We have used ξ = 1 and corresponding running charm quark mass, m c = 1228.52 MeV as well as parameter, α s = 0.2636. In Eqs. (47) and (49), we obtain the value of φ b and φ c using medium modified gluon condensates within the chiral SU (3) model and are used to calculate moment M J n , and hence the in-medium charmonium mass m * C , under QCD sum rules. The obtained values of φ b and φ c are tabulated in the Table 3. At T = 0 and ρ N = ρ 0 , the above values of φ b can be compared with the    Table 3 In the above table, we tabulate the values of φ b and φ c under the effect of magnetic field, temperature (in MeV) and density eB (m 2 π ) φ b in unit of 10 −3 φ c in unit of 10 −5  G 0 Total G 0 without finite quark mass term G 0 with only finite quark mass term G 2 Total G 2 First term G 2 Second term G 2 Third term values 1.7 × 10 −3 and 1.6 × 10 −3 at vacuum and nuclear saturation density respectively in Ref. [77]. We plot in-medium mass of J/ψ and η c in Figs. 6, 7, and 8 for isospin asymmetry η = 0, 0.3 and 0.5 respectively. We show the results for J/ψ and η c at T = 0,100 MeV and ρ N = ρ 0 ,4ρ 0 and for magnetic field eB = 0, 5m 2 π and 7m 2 π . In each subplot the results are compared with vacuum situation, i.e. ρ N = B = T = 0. In Table 4, we summarize the values of mass-shift of J/ψ and η c at eB = 5m 2 π and 7m 2 π from the values at B = 0. From Fig. 6 we observe that, for given density and temperature, the finite magnetic field causes more decrease in masses of J/ψ and η c mesons. However, we observe less decrease in mass of J/ψ and η c , as we move from T = 0 to T = 100 MeV. The negative mass-shift due to magnetic field is more in high density regime. For example, at ρ N = ρ 0 , η = 0 and T = 0, the mass-shift of J/ψ and η c is − 3.39 (− 3.68) and − 4.63 (− 5.07) MeV for eB = 5m 2 π (7m 2 π ) with respect to vacuum mass respectively, whereas at 4ρ 0 , the mass-shift modifies to − 11.63 (− 14.68) and − 11.97 (− 12.96) MeV . This is due to the fact that at high density the values of gluon condensates decrease more with the increase of magnetic field, which further modify φ b Fig. 6 The in-medium mass of J/ψ and η c meson plotted as a function of n for symmetric nuclear matter at different values of temperature, magnetic field and density Fig. 7 The in-medium mass of J/ψ and η c meson plotted as a function of n for asymmetric nuclear matter, η = 0.3 at different values of temperature, magnetic field and density 3180 3200 3220 3240 3260 m Fig. 8 The in-medium mass of J/ψ and η c meson plotted as a function of n for asymmetric nuclear matter, η = 0.5 at different values of temperature, magnetic field and density 3180 3200 3220 3240 3260 m through σ , ζ , δ, and χ fields resulting in the in-medium modifications of gluon condensates and masses of J/ψ and η c . This interaction also accounts for scalar and number density of nucleons under the effect of magnetic field. For proton, since it has non-zero charge, we have contributions from the Landau levels hence it has strong influence, whereas for neutrons (zero charge) there are no Landau level contributions. Mainly, the in-medium modifications of J/ψ and η c owes to the change in scalar gluon condensate and twist-2 gluon operator, which is obtained from the medium modification of scalar fields σ , ζ , δ, and χ . These fields in strongly magnetized hot and dense isospin asymmetric (introduced by the δ and ρ fields) nuclear matter is obtained by solving the coupled equations of motions (Eqs. [18][19][20][21][22][23]. We have observed that the effective mass of J/ψ and η c decrease with the increase in magnetic field except in the case of low-temperature regime of high isospin asymmetry, where increase in effective mass is observed due to the crossover behaviour of scalar fields and scalar densities of protons and neutrons. The mass modifications of these mesons are observed to be appreciable in high density and symmetric nuclear matter. It is also observed that the masses of vector meson, J/ψ and pseudoscalar meson, η c increase as we move towards high temperature. The observed mass-shift may be used to study the decay of charmonia [121] in the presence of magnetic field. The present work will shed light on non-perturbative effects of magnetic field on the properties of mesons and hadrons, which may help to understand the experimental observables such as decay width [120,121] and cross-section [122], in asymmetric noncentral heavy ion collision experiments.