Equilibrium Instability of Chiral Mesons in External Electromagnetic Field via AdS/CFT
Abstract
We study the equilibrium instability of chiral quarkonia in a plasma in the presence of constant magnetic and electric field and at finite axial chemical potential using AdS/CFT duality. The model in use is a supersymmetric QCD at large ’tHooft coupling and number of colors. We show that the presence of the magnetic field and the axial chemical potential even in the absence of the electric field make the system unstable. In a gapped system, a stable/unstable equilibrium state phase transition is observed and the initial transition amplitude of the equilibrium state to the nonequilibrium state is investigated. We demonstrate that at zero temperature and large magnetic field the instability grows linearly by increasing the quarkonium binding energy. In the constant electric and magnetic field, the system is in a equilibrium state if the Ohm’s law and the chiral magnetic effect cancel their effects. This happens in a subspace of space with constraint equation , where and are called electric and chiral magnetic conductivity, respectively. We analyze the decay rate of a gapless system when this constraint is slightly violated.
a]Seyed Farid Taghavi,
P.O. Box 193955531, Tehran, Iran
\affiliation[b]Department of physics, Kharazmi university,
P.O.Box 3197937551, Tehran, Iran
\emailAdds.f.taghavi@ipm.ir
\emailAddvahedi@ipm.ir
\keywordsQuarkonium mesons, AdS/CFT, Schwinger effect, Chiral magnetic effect, Instability
1 Introduction
The vacuum instability in the presence of strong electric field, called Schwinger effect, was predicted long time ago by Heisenberg and Euler and more formally developed by Schwinger [1, 2]. As an intuitive picture, if the potential between two plates is above a critical value then the (QED) vacuum between two conducting plates sparks similar to a dielectric breakdown in a capacitor. It is well known that the decay rate of the vacuum can be obtained by [3]. Here, the effective Lagrangian is achieved by integrating out the fermionic degrees of freedom of a QED theory. The result is called EulerHeisenberg Lagrangian and its imaginary part for a constant background electric field is [2]
(1) 
where is the fermion mass. The above is related to the transition amplitude of the tunneling though a potential barrier for producing a pairparticle from the vacuum. The required energy for a pair creation is and it should be supplied by separating the particles by a distance greater than , the Compton’s wavelength, via electric field. Using this picture, a critical value of the electric field can be simply estimated to be () where the exponential term in \eqrefSchwingerOrigin becomes significantly large. In the QED theory, the ’tHooft coupling is and therefore .
Recall the similarity of Schwinger effect with capacitor breakdown. The atoms of an insulator are ionized in a strong electric field and the decay rate of a band (and also Mott) insulator groundstate in a Zener breakdown is similar to the relation \eqrefSchwingerOrigin [4]. Beside QED, we expect to observe the same phenomenon in the QCD bound states. Due to the confinement, the perturbation is not applicable in this case and nonperturbative methods such as lattice QCD or AdS/CFT correspondence are needed to study this effect in the hadrons dissociation.
In order to deal with such a nonperturbative phenomenon, we will use AdS/CFT correspondence [5, 6, 7]. According to this correspondence, the theory of super YangMills (SYM) is dual to the type IIB string theory on . We will use a special limit of the duality where the number of colors goes to infinity and ’tHooft coupling is large. In this limit, the strongly coupled SYM theory is dual to the classical type IIB supergravity on . All fields in this case are in the adjoint representation of the gauge group. We add flavors of ”quarks” in the fundamental representation by introducing D7branes in the background metric. The result can be considered as a QCDlike supersymmetric theory where quarks are in a hypermultiplet [8, 9, 10].
A holographic picture of the Schwinger effect in the limit and large ’tHooft coupling is proposed in [11] where the pair production of ”W bosons” in the SYM field theory has been explored. It has been shown that demanding the reality of DiracBornInfeld (DBI) action of the probe D3brane in the background leads to a critical value for electric field, . Note that the value of is similar to our previous intuitive estimation up to a constant. In the context of AdS/CFT, the dissociation of the hadrons in external electric field has been studied in [12, 13]. In this work, the probe D3branes in [11] are replaced by the D7branes to add flavor quarks in the system. In fact, the DBIaction of the D7branes plays the role of the EulerHeisenberg Lagrangian. It is also shown that if we consider the confining force as fermion mass then the nontrivial subleading term in imaginary part of the D7brane action coincides with the supersymmetric Schwinger effect. It is worth mentioning that the holographic picture of Schwinger effect in the different confining backgrounds is also studied in [14].
The meson spectroscopy has been studied in [15] by using the D3/D7 brane holographic model. This model is dual to the mesons in a plasma with deconfined gluons. According to this study, the meson mass can be estimated to be proportional to where stands for the quark mass. At finite temperature, above the critical value , the mesons are melted into the thermal gluonic plasma [16]. The gravity dual picture of the meson melting is related to the topology of the D7branes embedding [26, 17, 18, 19, 20, 21, 22, 23, 24, 25]. Let us mention that the meson melting might happen due to the other sources rather than temperature. Specifically, we are interested in the phase transition in the presence of axial chemical potential in the present work which has been studied in [24]. Using AdS/CFT correspondence, the decay width of the heavy quarkonia in a stronglycoupled thermal plasma has been studied by investigating the imaginary part of the heavy quark potential [27, 28, 29].
The aim of this work is to study whether the magnetic field can dissociate mesons in a chiral supersymmetric plasma. Note that the electric field plays an active role in the meson dissociation by compensating the confining force between quarks in the meson. However, the magnetic field alone does not play such a role in the system. According to the meson spectroscopy, the mesons mass is changed by the influence of electric and magnetic field [33, 16, 30, 31, 32, 26]. Although the meson spectrum is altered in the background magnetic field, it can not dissociate the mesons. In fact, there is a critical value for magnetic field where even at finite temperature there is no melted meson in the system in contrast with electric field [31, 34, 35, 36].
However, we can look at the effect of magnetic field on the plasma differently. Due to the chiral anomaly, there is an anomalous transport coefficient which leads to the induction of electrical current along the magnetic field in a fluid with nonzero axial chemical potential. This effect is called chiral magnetic effect (CME),
(2) 
and has been studied extensively in the literature [37, 38, 39, 40, 41]. Here is the axial chemical potential. One may wonder whether the CME gives rise to the instability similar to the confined mesons in the external electric field. In this article, we will see that the same instability occurs with a different mechanism. Unlike the electric field, the magnetic field is not destructive here as we expected, but, the internal energy enhancement by increasing the axial chemical potential above a critical value will liberate the confined quarks. More precisely, an equilibrium state of mesons in the background magnetic field and axial chemical potential (and temperature) above a critical value decays to a system with nonzero electric current. The critical values of axial chemical potential and temperature can be depicted in stable/unstable equilibrium phase diagram which exhibits similar features compared to the  phase diagram. This has already been studied by using NambuJonaLasinio model with polyakov loop (PNJL model) [42], linear sigma model coupled to quarks and the Polyakov loop () [43] and also lattice QCD simulation [44]. It is worth noting that studying the full time evolution of the equilibrium state to a steady state with nonzero electric current is beyond the scope of this paper. Instead, the decay rate of the initial state will be considered here.
The AdS/CFT setup we will use is developed in [24, 25, 36, 46, 45]. In this setup, the , the symmetry of the supersymmetry in the field theory, is assumed as the axial symmetry and the angular velocity of the rotating D7branes around the D3branes is dual to the axial chemical potential. In [46] this picture is used to investigate the CME where the is broken explicitly by quark mass. The effect of anisotropic plasma on CME [47] and stringy corrections on the holographic picture of CME [48] are also studied in this setup. Here, we investigate the decay rate of the plasma instability by studying the imaginary part of the D7brane action. It is interpreted as the EulerHeisenberg Lagrangian such that the axial chemical potential is encoded in it.
Our paper is organized as follows: In section 2, we review the holographic picture of supersymmetric QCD with flavor using D3/D7 brane construction. In section 3, we study the equilibrium instability for massless quarkonia and also stable/unstable equilibrium phase and its instability for a gapped system. The instability of gapless mesons in the presence of electric and magnetic field and axial chemical potential is investigated in section 4. Finally, in section 5, we summarize the main points of the paper.
2 A Brief Review on the Holographic picture of Supersymmetric QCD
Using D3/D7 branes, a holographic description of a QCDlike field theory has been studied extensively in the past years [10, 49, 50] (for review see [51, 52]). In this section, we briefly review the gravity picture of the supersymmetric field theory at finite temperature and axial chemical potential where the external electric and magnetic field are present
Consider stack of D3branes and stack of the D7branes embedded in the ten dimensional space in the following form,
(3) 
In the limit and for the large ’tHooft coupling, the D3branes are replaced by the background metric and a self dual RR five form field . The background metric is protected against the D7brane backreaction since we are in the probe limit . The dual field theory of this supergravity picture is the super YangMills (SYM) theory in adjoint representation coupled to a hypermultiplet in the fundamental representation.
In the low energy limit, the action of the D7brane for a generic background metric can be written as
(4) 
where the DBI and ChernSimons (CS) actions are given by
{subequations}
{alignat}3
S_DBI &= N_f T_D7 ∫d^8ξ\textdet(g_ab^D7+2πα’F_ab),
S_CS &= N_f T_D7∫P[ΣC^(n)]e^2πα’F_ab.
In the above, is the induced metric on the D7brane and is the tension of the brane. The symbol stands for the pullback of the form field . Moreover, the field strength of the gauge field living on the brane is shown by .
In the absence of the probe D7branes, the rotational symmetry in the transverse direction of the D3branes corresponds to , the symmetry of the SYM. The D7branes break into a rotational symmetry in the plane and plane which is . This isometry corresponds to where is a global symmetry and is the symmetry. For the separated D3 and D7branes, the rotation in plane is also broken and correspondingly the is broken explicitly in the field theory side.
The hypermultiplet is in the fundamental representation of . This hypermultiplet consists of two chiral superfield and with opposite chirally in the notation. The two Weyl fermions of chiral superfields, and are transformed under with and charges and we can combine these two Weyl fermions into a Dirac fermion . In this picture, plays the role of the axial symmetry . It is worth noting that unlike QCD there are some charged mesons under transformation [51].
Here we are interested in a system at finite temperature, finite axial chemical potential and in the presence of external electric and magnetic field. For those reasons, the brane construction by the following design is chosen,

In order to have finite temperature in the dual field theory, the AdSSchwarzschild is chosen for the background metric in the following coordinate,
(5) where
(6) In the above, and is the metric of a unit 3sphere. Here, the boundary of the AdS is at and indicates the location of the horizon of the black hole. The Hawking temperature, which is also the temperature of the dual field theory, is given by
(7) 
We mentioned earlier that the rotational symmetry in the plane corresponds to a global axial symmetry in the field theory side and the length of the strings stretched between D3 and D7branes corresponds to the mass of the quarks in the fundamental representations.
Moreover, it is shown in [53] that the rotating D7branes around D3branes with constant angular velocity leads to the complex mass in the dual field theory where is the azimuthal coordinate in the plane in our case. Now it is simple to argue that the angular velocity of the D7brane corresponds to the axial chemical potential in the dual field theory due to the chiral anomaly [46].

We are interested in finding out the response of the system to the background electric and magnetic field. For simplicity, we assume and are parallel when both exist in the background. Thus, we turn on the gauge field on the D7brane in the following form,
(8)
The above brane setup fulfills our desires in the dual field theory.
The embedding of the D7branes can be understood more easily in the coordinate
(9) 
where and are defined in the metric \eqrefmetric
Let us go back to the coordinate mentioned in \eqrefmetric. The D7branes are extended along where . In the static gauge, the D7branes fill the coordinates . Nevertheless, the embedding is not completely determined unless we specify the coordinates which is function of in general. However, translational symmetry in direction and rotational symmetry in the restrict the functionality of the remaining coordinates to the form and . We could choose for the later ansatz. However, this choice leads to instability of the D7brane action while the linear velocity at some point of the rotating brane may exceed the local speed of light. The additional term let the brane to ”twist” in direction and cure this problem [46].
In this section, we focus on the massless quark limit where . Later, we will generalize it to the massive quarks by considering certain assumptions. A straightforward calculation leads to the following form for the DBI action,
(10) 
where prime superscript means derivative with respect to . The terms and are given by
(11)  
(12) 
In the above, we have performed the integration over trivially. Moreover, we have redefined where is the volume of the Minkowski spacetime. Recalling that the 3sphere volume is and , the constant behind the integral is given by
(13) 
The angular velocity of rotating D7brane does not appear in the action \eqrefSimplDBI which means that the effect of the axial chemical potential is not encoded in the DBI action and it appears in the CS action. The selfdual five form is where is the volume form of a 5sphere with radius . This field is found by the following RR four form potential,
(14) 
where is the volume form of a 3sphere with unit radius. Using above and gauge field \eqrefgaugeFields, the CS action can be written as
(15) 
where
(16) 
and we have performed a similar action redefinition .
From \eqrefSimplDBI and \eqrefSimplCS, one can see that the only depends on the derivative of and . Hence, by performing two successive Legendre transformation, we can replace and by two corresponding constants of motion. The constant of motion correspond to is and Legendre transformation trivially cancels the term in the action. It means that there is no constraint on the field and we can freely set it equal to zero in the massless limit. However, the constant of motion related to is not trivial. If we define
(17) 
then the following transformation replaces field by its corresponding constant of motion ,
(18) 
In the next section, we will extensively explore the imaginary part of the brane action. However, we are dealing with a steady state situation here. As a result, in this case the action should be real. The reality condition of the above action demands that the quantity and terms inside the round bracket change their sign at the same . The term changes its sign at
(19) 
while it is clear from \eqrefQ2 that is always a positive quantity. To avoid the imaginary action, the term in the round bracket inside the square root should be equal to zero at . It leads to
(20) 
Using holographic renormalization, it has been shown that the electrical current in the field theory side is proportional to , more precisely,
(21) 
where
(22) 
The and are studied previously in [45] and [46], respectively. Setting () equal to zero, the CS action vanishes and consequently . Unlike , the magnetic conductivity do not depend on the ’tHooft coupling. In fact, CME in the massless quark limit and strong coupling is similar to the calculations in the weak coupling. The term is anomalyinduced transport coefficient and we know that the chiral anomaly is oneloop exact (for more discussions refer to [54] and references therein). However, massive quarks break the chiral symmetry explicitly and we do expect that CME in strongly coupled differs from that in weak coupling. It can be seen that CME in the large and ’tHooft coupling receives stringy corrections for finite quark mass [48].
3 Chiral Equilibrium Instability of Quarkonium Mesons
The response of a QCDlike matter to the external electric and magnetic field at finite temperature and axial chemical potential has been studied in the previous section. We should note that the quarkonia are neutral degrees of freedom and we do not expect any interactions between them and the electromagnetic fields. It means that the electric current exists only in the phase of the melted mesons
Let us first review the physical picture of the instability in the presence of electric field [12]. For a fixed quark mass, consider a system with no electric field initially and turn it on in a short interval such that the final value of the electric field is strong enough to dissociate the quarks inside the mesons. In this case, the system is unstable and we expect that long enough time after turning on the field, the system reaches to the steady current . The full time evolution of the system needs a timedependent calculation in the gravity, however, the instability of equilibrium state just after applying the electric field can be studied by the imaginary part of the D7brane action. Both approaches have been studied in [12].
The electric current is also induced along the magnetic field in a plasma with nonzero axial chemical potential. Note that (as we mentioned in the introduction) there are some mesons in the supersymmetric QCD which are charged. Now, assume the axial chemical potential of the system is below a critical value and suddenly is increased above it. It leads to increasing the internal energy of the plasma and a phase transition happens. In this process, the electrically and chirally charged quarks are liberated from meson confinement, and eventually, CME leads to a steady current in the presence of external magnetic field. The increasing of axial chemical potential may happen by a vacuum tunneling in the gluonic sector [40] where the tunneling leads to a strong enough chiral imbalance in the system such that the axial chemical potential is raised above the critical value.
The time evolution of the induced current is not addressed in this study and we will only focus on the decay rate of the equilibrium state exactly after changing the axial chemical potential by studying the imaginary part of the D7brane action. Moreover, we set the electric field equal to zero to suppress the Schwinger effect in this section.
3.1 Instability of Gapless Mesons
In order to clarify the main idea of the paper, we investigate the (almost) gapless system in the large magnetic field in this subsection. The result for the small magnetic field is presented in the appendix A.1.
In our holographic picture, it is not reasonable to speak about mesons in a system with massless quarks. In fact, at zero temperature and axial chemical potential and in the absence of quark mass scale, the system is conformal and there is no mass spectrum for mesons [15]. Here, we assume the quark mass is extremely small compare to the other physical quantities in the system such as temperature and axial chemical potential. Hence, one can set and the resulting D7brane action is \eqrefLegendreAction. The other assumption is that the magnetic field does not reach to its critical value where the meson melting transition disappears [31]. Above this critical value, the is not the preferred solution energetically and our calculations is not reliable anymore. This point has been studied in [36] for the case of nonzero magnetic field and axial chemical potential.
Exactly after increasing the axial chemical potential suddenly, the electric current should be equal to zero. Regarding this argument (and following [12]), we set in the action \eqrefLegendreAction and study its possible imaginary part. After defining the following dimensionless quantities,
(23) 
the explicit form of the action in this limit is
(24) 
where and . Note that all dimensionful parameters are collected behind the integral.
The field theory as well as the gauge/gravity computations show that the Schwinger effect vanishes for the zero electric field even if the magnetic field is nonzero. But, interestingly, the appearance of the new term inside the square root related to the axial chemical potential is a sign of a new source for the equilibrium instability at the presence of and magnetic field and absence of the electric field.
In the zero temperature, the term inside the square root is negative at its minimum . Therefore, the imaginary part of the action comes from the integration in the range where is the root of the polynomial . The problem becomes more simple in the limit where we can find and analytically. In this limit, we should be careful about the critical value of the magnetic field mentioned at the beginning of this subsection. For example, repeating the [36] procedure, we can find for . Hence, in any case, the magnetic field should be always below an upper bound.
Assuming the polynomial has a finite root at , one can simply ignore and find . Around this point, we consider the expansion for small values of . The coefficients can be achieved order by order in this expansion,
(25) 
We could go further in the expansion, but, comparing to the numerical calculations show that it does not improve the result for .
Until the root is approximately , we have in the range of integration and we can expand for large . Using from above, we find the following expansion for the ,
(26) 
Unlike the holographic pictures of Schwinger effect, the above imaginary part depends on the ’tHooft coupling which might be an evidence for the difference of instability in the weak coupling limit. Moreover, the relation \eqrefBigBImaginaryAnal is an expansion for small . Nevertheless, after substituting physical dimensionful quantities, one finds that it is an expansion with respect to the ’tHooft coupling. It shows that this result cannot be calculated by the perturbative methods
In general, one can find and calculate the integration in \eqrefSimpAction numerically. This has been done and the result is plotted in Fig. 1(a) for and different values of the magnetic field (). The numerical values of normalized to the coefficient are plotted by the blue dots in this figure
\includegraphics[scale=0.80]mu5Eq12.pdf 
\includegraphics[scale=0.8]nonZerTemHighB250.pdf 
(a)  (b) 
Let us now study in the large magnetic field and finite temperature. In this limit, if the negative part of the term inside the square root is in a finite range of then the term can be ignored in \eqrefSimpAction. As we will see, it is the case here. This limit brings the magnetic field out of the square root. If one rescales the integration parameter , the imaginary part of the simplified action turns into the following form,
(27) 
where
(28) 
and is the root of . By inspection, one can find that for the function is monotonically decreasing and in the range , the root is as follows,
(29) 
It confirms the approximation of neglecting the term in \eqrefSimpAction.
Note that goes to infinity in the limit . Let us introduce and demand that it has same infinite properties as . Considering \eqrefFinitTemRoot, it can be simply seen that the condition (high temperature limit) leads to . In other words, large values of shrinks the integration range of to a small interval with length . Therefore, one can set and expand the function as follows,
(30) 
As a result, we choose and fix coefficient such that . The final result is the following,
(31) 
It can be simply checked that the difference between and for () in the integration interval is less than . The integration of approximated function can be performed and the result is
(33) 
The interesting point of the above result is that it is not depend on in the first order of high temperature expansion. The other interesting point is that the above result defers by a factor of in comparison with the first term in the relation \eqrefBigBImaginaryAnal at zero temperature. It shows that raising the temperature in the large magnetic field decreases the instability, meanwhile, the rate of decreasing becomes smaller at higher temperature and finally approaches to \eqrefnonThighB.
In more general case, we can study the equilibrium instability at high magnetic field numerically. The temperature dependence of (normalized to constant and large magnetic field ) for three different values of is shown in Fig. 1(b). It can be seen that the equilibrium instability decreases by increasing the temperature and approach to \eqrefnonThighB. The initial points in the numerical results are also in agreement with the first term \eqrefBigBImaginaryAnal.
3.2 Stable/Unstable Equilibrium State for Gapped Mesons
In this subsection, we extend our study to the massive quarks where the function is not equal to zero anymore. One needs to calculate its functionality by solving the equation of motion. The asymptotic expansion of the solution is where is proportional to the quark mass and is proportional to the quark condensate .
Before exploring the equilibrium instability of chiral mesons, let us mention some remarks on holographic CME [46]. The D7brane embedding (or equivalently in the coordinate \eqrefOldCoord) is achieved by the equation of motion calculated by varying the action \eqrefactionColli. Due to the rotation of the branes, a worldvolume horizon might induce on the D7brane which is given by
(34) 
At the finite temperature and axial chemical potential, there are two possible horizons, background horizon of AdSSchwarzschild and worldvolume horizon. Accordingly, there are three different embeddings for D7branes which depend on the induction of background or worldvolume horizons on them [46, 48, 47]. The CME is nonzero for embeddings that intercept the above locus. Using this picture, one can relate the current to the quark mass as follows. By increasing the quark mass, there is a critical mass (critical embedding) where the D7brane do not cross the locus \eqrefWVHorizon which means CME is zero for (see figure 2 of Ref.[46]). In the field theory picture, for small mass, there are some melted mesons or equivalently quarks that move along the magnetic field via CME. However, for larger mass, there are neutral bound state of quarks in the system, therefore, there is no current.
\includegraphics[scale=0.7]LowTem.pdf  \includegraphics[scale=0.7]midTem.pdf  \includegraphics[scale=0.7]HighTem.pdf 
According to the above picture, we can plot a stable/melted phase diagram for a given quark mass. One should fix the mass and check for what values of and the current becomes nonzero. We have shown three different critical embeddings for same quark mass in Fig. 2. The red dashed curve indicates worldvolume horizon and black curve depicts the background horizon and the blue curve is the critical embedding. All the curves are in \eqrefOldCoord coordinate. Different values of and for same quark mass are shown by red triangles in Fig. 3(a). Here the value of magnetic field is chosen which is much lower than its critical value. For and the critical magnetic field is . One can check that the phase diagram in Fig. 3(a) does not change drastically by choosing and .
Now we would like to find out the decay rate of a system at temperature when the value of suddenly increases. Basically, we could assume that the quarkonia at given , and external magnetic field is in the equilibrium state. It leads to a numerical function for . Using it, we calculate the terms inside the square root in the action and study its value for nonzero and investigate whether it becomes imaginary. Note that the existence of nonzero deforms the D7brane shape, but, if it happens very rapidly then we can assume that the initial shape of the brane does not deviate from . The same logic is used for studying the Schwinger effect in the gapped systems in the presence of electric field [12]. In [12], the initial embedding is the supersymmetric solution for D7brane because the temperate is assumed to be zero initially. In , the supersymmetric solution for D7brane embedding is
(35) 
where . This solution in the coordinate \eqrefOldCoord is simply . Following this logic, we assume at finite temperature and zero chemical potential the solution is approximately what is mentioned in \eqrefsuperSolu. This is a good approximation in the following conditions: ) The magnetic field is much smaller than its critical value, or ) The increase in the axial chemical potential and the magnetic field happens simultaneously. In fact, this approximation becomes exact if all variables turn on at the same time, however, we assume the first option, .
For the massive quarks, the parameters , and in \eqrefQ1, \eqrefQ2 and \eqrefPPdefini change to the following forms,
(36) 
By defining and using redefinitions \eqrefdimlessQuant, the action can be rewritten as
(37) 
where .
We would like to find a stable region in the , and space. First assume which means . In \eqrefMoreGenAction01, the term in the square bracket should be positive in the range to avoid the negative values inside the square root. It happens if we choose or equivalently,
(38) 
In the region , however, there is always a window that the action \eqrefMoreGenAction01 becomes imaginary. In this region, has always a positive value, but the term inside the square bracket goes to at and at which means it changes its sign somewhere in the middle. The relation \eqrefstableRegion shows the stable/unstable region for the given and . Note that in this relation is absent. It is compatible with the previous remark that the magnetic field has not active role in the quarkonia instability. However, we should mention that the approximation \eqrefsuperSolu leads to the magnetic field independence in the stable/unstable regions. We expect the value of magnetic field (well below its critical value) changes mildly the stable/unstable region in a more accurate study. There is a critical value for the temperature where above that, the system becomes unstable for any values of axial chemical potential. Similarly a critical value for is obtained,
(39) 
\includegraphics[scale=0.73]criticalRegion02.pdf  \includegraphics[scale=0.75]MasGapPlot.pdf 
(a)  (b) 
The stable region \eqrefstableRegion is the region where the supersymmetric solution \eqrefsuperSolu for a fixed mass does not intercept the worldvolume horizon \eqrefWVHorizon for different values of and . In other words, one can find \eqrefstableRegion by using the supersymmetric solution and the method we have found the stable/melted meson phases. This is the shaded region in Fig. 3(a). The difference between stable region (shaded region) and stable meson phase (the region below the red triangles) is that the later is achieved in the equilibrium state where and are assumed to be fixed, while the former is calculated for the case is suddenly increases from zero. In the gravity side, it can be understood by the difference between two solutions of in two cases.
Now we will study the decay of the equilibrium state when increases suddenly such that the system goes into the unstable region in  space. In the limit and , it is possible to find an analytical result. The imaginary part of the action \eqrefMoreGenAction01 in this limit is
(40) 
which can be written in terms of the field theory quantities as,
(41) 
Interestingly, only the term in the square bracket is added compare to the gapless system (see \eqrefBigBImaginaryAnal). Recall that adding mass in the Schwinger effect \eqrefSchwingerOrigin leads to a summation of exponentials which is interpreted as a sum over instantons.
The other feature appears if we write in terms of the quarkonium mass and its binding energy . In this case, the ’tHooft coupling does not appear explicitly in the equation. The mass of the lightest meson in the model under consideration is given by [15]
According to the above, the meson mass is much smaller than its quark constituents. Consequently, the binding energy is very large in this model and is given by . Using these facts, we have
(42) 
This is true for large number of colors and large ’tHooft coupling of SYM theory and what we could ask is whether it is true for QCD. Unlike the mesons in our model, the binding energy is not large in QCD. As an example the meson binding energy is around which is one order smaller than its mass and constituent quarks. In any case, according to the above result, the system is more unstable for mesons with larger binding energies.
In more general cases, a numerical analysis is needed to find . It can be found in Fig. 3(b) for (red), (black) and (blue). The quark mass is chosen and . It is seen that the is nonzero when the . The red curve corresponds to the analytical result for instability \eqrefGappedAnal01. There is a good agreement between the analytical and the numerical results for .
4 Melted Meson Equilibrium Point and its Instability
In this section, we again study a gapless system, but this time, in the presence of electric and magnetic field at the same time. One can imagine a case that both Ohm’s law and CME produce current in opposite directions such that the net current is equal to zero. In this case, there is a subspace in space in which the system is always in equilibrium
Choosing , the D7brane action is written as
(43) 
where
(44) 
Let us define and where and are roots of and , respectively. Then the imaginary part of the action comes from integration of in the range . When the so the system is stable. This happens in a special point where and have common roots. Such a point exists if the equation is satisfied for at least one point in space. Using \eqrefConductis, this equation reduces to
(45) 
where . Here, we will try to find when we slightly violate the constraint .
In general, numerical calculations are needed to find the . However, we can simplify it in some special limits and find analytical results. In the limit , the electrical conductivity is
(46) 
Now the equation \eqrefOHMCME can be written as
(47) 
in the large magnetic field where
(48) 
Regarding this equation, mimics the axial chemical potential in large magnetic field and the equation \eqrefStablePoint reduces to .
Let us calculate for a large magnetic field and zero temperature in a small interval with radius around stable region . In this limit, we have
(49) 
where and and . The term inside the square root can be factorized in the following form
(50) 
Inside the interval of integration we have . It indicates that the leading term inside the first bracket in the above relation does not depend on for finite and . However, the leading term in the second bracket is in the order of . We approximate the first bracket by evaluating it around . One can check that the value of the first bracket in this approximation is equal . The integration of the second bracket inside the square root is nothing but the area of a semicircle with radios . The final result is the following
(51) 
According to the above result, the leading term of the instability is quadratically increasing when we move away from the stable subspace in large magnetic field and zero temperature.
In Fig. 4, the numerical calculations are compared with analytical results for . The stable point occurs at . The red dashed curve in this figure is analytical result \eqrefPerturbFromStable and blue dots are obtained by numerical calculations. It can be seen that two results have good agreement around the stable point . Note that at the leading term inside the first square bracket is not large compare to anymore and the analytical result is not valid. The has been already calculated in \eqrefBigBImaginaryAnal for the case and it shows correct value in comparison with numerical solutions at this point.
As a final remark, let us mention that by setting in \eqrefSimpAction the diverges (for details please refer to [13]). Interestingly, the existence of in our study leads to disappearance of this divergence even at the zero temperature.
5 Summary
Using AdS/CFT, equilibrium instability of chiral quarkonia in a plasma in the presence of external electric and magnetic field and at finite axial chemical potential has been studied. We have found that despite the system is stable when there exist magnetic field only, in the case that it is accompanied by axial chemical potential the system may become unstable. This instability has been investigated via imaginary part of the D7branes action.
The stable/unstable phase transition and the instability of mesons (see Fig. 3) have been investigated in the following picture. Consider a stable system of mesons in a plasma at finite temperature such that its axial chemical potential is much smaller than a critical value (e.g. initially). Also, a nonzero external magnetic field exists in the background. Recall that quarkonia in the system are charged with respect to but they are not electrically charged. For that reason, they can not couple to magnetic field and no current can be induced by CME. However, by assuming an external mechanism (a vacuum tunneling for example) a sudden change in the axial chemical potential occurs such that goes beyond the critical value . Then the system becomes unstable and constituent quarks of the quarkonium are liberated. In this case, the new degrees of freedom are electrically and axially charged and CME produces the current. It means that the initial equilibrium decays to another state which is ultimately a steady state with constant current. For quarkonia with mass and binding energy the at zero temperature and large magnetic field has been presented in \eqrefbindingMesonInst. Also the numerical study for more general cases have been depicted in Fig. 3(b).
Finally, the massless mesons have been investigated in a system that two Ohm’s law and CME cancel their effects. These are two different mechanisms in producing current. Basically, it happens for specific choices of , and (and ) such that . Any deviation from this subspace makes the equilibrium state unstable. We have studied when the constraint is slightly violated.
We would like to thank M. AliAkbari, H. Ebrahim, S. Fayazbakhsh, A. E. Mosaffa, H. Arfaei, A. Davodi, N. Abbasi, Y. Ayazi, D. Allahbakhshi for useful discussions. We are also grateful to D. Kaviani for helpful and fruitful comments.
Appendix A More About Instability of Gapless Systems
a.1 Small Magnetic Field
Recall the action \eqrefSimpAction. At zero temperature, the term inside the square root is which is negative in the range . In the limit , we can easily ignore the term and immediately find the root . Now obtaining the imaginary part of the action \eqrefSimpAction is straightforward,
(52) 
where is the beta special function.
\includegraphics[scale=0.8]nonZerTemLowB.pdf 
For the finite temperature case, we do as follows. By sending the magnetic field to zero and defining the imaginary part of the action \eqrefSimpAction can be written as
(53) 
where
(54) 
Similar to the case at zero temperature, the root of can be found perturbatively as the following,