Large N approach to Kaon decays and mixing 28 years later: \Delta I=1/2 rule, \hat{B}_K,and \Delta M_K

We review and update our results for $K\rightarrow \pi \pi $ decays and $K^0$ – $\bar{K}^0$ mixing obtained by us in the 1980s within an analytic approximate approach based on the dual representation of QCD as a theory of weakly interacting mesons for large $N$ , where $N$ is the number of colors. In our analytic approach the Standard Model dynamics behind the enhancement of $\hbox {Re}A_0$ and suppression of $\hbox {Re}A_2$ , the so-called $\Delta I=1/2$ rule for $K\rightarrow \pi \pi $ decays, has a simple structure: the usual octet enhancement through the long but slow quark–gluon renormalization group evolution down to the scales $\mathcal{O}(1\, {\hbox { GeV}})$ is continued as a short but fast meson evolution down to zero momentum scales at which the factorization of hadronic matrix elements is at work. The inclusion of lowest-lying vector meson contributions in addition to the pseudoscalar ones and of Wilson coefficients in a momentum scheme improves significantly the matching between quark–gluon and meson evolutions. In particular, the anomalous dimension matrix governing the meson evolution exhibits the structure of the known anomalous dimension matrix in the quark–gluon evolution. While this physical picture did not yet emerge from lattice simulations, the recent results on $\hbox {Re}A_2$ and $\hbox {Re}A_0$ from the RBC-UKQCD collaboration give support for its correctness. In particular, the signs of the two main contractions found numerically by these authors follow uniquely from our analytic approach. Though the current–current operators dominate the $\Delta I=1/2$ rule, working with matching scales $\mathcal{O}(1 \, {\hbox { GeV}})$ we find that the presence of QCD-penguin operator $Q_6$ is required to obtain satisfactory result for $\hbox {Re}A_0$ . At NLO in $1/N$ we obtain $R=\hbox {Re}A_0/\hbox {Re}A_2= 16.0\pm 1.5$ which amounts to an order of magnitude enhancement over the strict large $N$ limit value $\sqrt{2}$ . We also update our results for the parameter $\hat{B}_K$ , finding $\hat{B}_K=0.73\pm 0.02$ . The smallness of $1/N$ corrections to the large $N$ value $\hat{B}_K=3/4$ results within our approach from an approximate cancelation between pseudoscalar and vector meson one-loop contributions. We also summarize the status of $\Delta M_K$ in this approach.     

Description of the Translation-Invariant Splitting Gibbs Measures for the Potts Model on a Cayley Tree

For the \(q\) -state Potts model on a Cayley tree of order \(k\ge 2\) it is well-known that at sufficiently low temperatures there are at least \(q+1\) translation-invariant Gibbs measures which are also tree-indexed Markov chains. Such measures are called translation-invariant splitting Gibbs measures (TISGMs). In this paper we find all TISGMs, and show in particular that at sufficiently low temperatures their number is \(2^{q}-1\) . We prove that there are \([q/2]\) (where \([a]\) is the integer part of \(a\) ) critical temperatures at which the number of TISGMs changes and give the exact number of TISGMs for each intermediate temperature. For the binary tree we give explicit formulae for the critical temperatures and the possible TISGMs. While we show that these measures are never convex combinations of each other, the question which of these measures are extremals in the set of all Gibbs measures will be treated in future work.     

The Distribution of the Area Under a Bessel Excursion and its Moments

A Bessel excursion is a Bessel process that begins at the origin and first returns there at some given time \(T\) . We study the distribution of the area under such an excursion, which recently found application in the context of laser cooling. The area \(A\) scales with the time as \(A \sim T^{3/2}\) , independent of the dimension, \(d\) , but the functional form of the distribution does depend on \(d\) . We demonstrate that for \(d=1\) , the distribution reduces as expected to the distribution for the area under a Brownian excursion, known as the Airy distribution, deriving a new expression for the Airy distribution in the process. We show that the distribution is symmetric in \(d-2\) , with nonanalytic behavior at \(d=2\) . We calculate the first and second moments of the distribution, as well as a particular fractional moment. We also analyze the analytic continuation from \(d<2\) to \(d>2\) . In the limit where \(d\rightarrow 4\) from below, this analytically continued distribution is described by a one-sided Lévy \(\alpha \) -stable distribution with index \(2/3\) and a scale factor proportional to \([(4-d)T]^{3/2}\) .     

Toward a Mathematical Holographic Principle

In work started in [17] and continued in this paper our objective is to study selectors of multivalued functions which have interesting dynamical properties, such as possessing absolutely continuous invariant measures. We specify the graph of a multivalued function by means of lower and upper boundary maps \(\tau _{1}\) and \(\tau _{2}.\) On these boundary maps we define a position dependent random map \(R_{p}=\{\tau _{1},\tau _{2};p,1-p\},\) which, at each time step, moves the point \(x\) to \(\tau _{1}(x)\) with probability \(p(x)\) and to \(\tau _{2}(x)\) with probability \(1-p(x)\) . Under general conditions, for each choice of \(p\) , \(R_{p}\) possesses an absolutely continuous invariant measure with invariant density \(f_{p}.\) Let \(\varvec{\tau }\) be a selector which has invariant density function \(f.\) One of our objectives is to study conditions under which \(p(x)\) exists such that \(R_{p}\) has \(f\) as its invariant density function. When this is the case, the long term statistical dynamical behavior of a selector can be represented by the long term statistical behavior of a random map on the boundaries of \(G.\) We refer to such a result as a mathematical holographic principle. We present examples and study the relationship between the invariant densities attainable by classes of selectors and the random maps based on the boundaries and show that, under certain conditions, the extreme points of the invariant densities for selectors are achieved by bang-bang random maps, that is, random maps for which \(p(x)\in \{0,1\}.\)      

Quantum Measures on Finite Effect Algebras with the Riesz Decomposition Properties

One kind of generalized measures called quantum measures on finite effect algebras, which fulfil the grade-2 additive sum rule, is considered. One basis of vector space of quantum measures on a finite effect algebra with the Riesz decomposition property (RDP for short) is given. It is proved that any diagonally positive symmetric signed measure \(\lambda \) on the tensor product \(E\otimes E\) can determine a quantum measure \(\mu \) on a finite effect algebra \(E\) with the RDP such that \(\mu (x)=\lambda (x\otimes x)\) for any \(x\in E\) . Furthermore, some conditions for a grade-2 additive measure \(\mu \) on a finite effect algebra \(E\) are provided to guarantee that there exists a unique diagonally positive symmetric signed measure \(\lambda \) on \(E\otimes E\) such that \(\mu (x)=\lambda (x\otimes x)\) for any \(x\in E\) . At last, it is showed that any grade- \(t\) quantum measure on a finite effect algebra \(E\) with the RDP is essentially established by the values on a subset of \(E\) .     

Quantum billiards in multidimensional models with branes

A gravitational $D$ -dimensional model with $l$ scalar fields and several forms is considered. When a cosmological-type diagonal metric is chosen, an electromagnetic composite brane ansatz is adopted and certain restrictions on the branes are imposed; the conformally covariant Wheeler–DeWitt (WDW) equation for the model is studied. Under certain restrictions asymptotic solutions to WDW equation are found in the limit of the formation of the billiard walls which reduce the problem to the so-called quantum billiard on the $(D+ l -2)$ -dimensional Lobachevsky space. Two examples of quantum billiards are considered. The first one deals with $9$ -dimensional quantum billiard for $D = 11$ model with $330$ four-forms which mimic space-like $M2$ - and $M5$ -branes of $D=11$ supergravity. The second one deals with the $9$ -dimensional quantum billiard for $D =10$ gravitational model with one scalar field, $210$ four-forms and $120$ three-forms which mimic space-like $D2$ -, $D4$ -, $FS1$ - and $NS5$ -branes in $D = 10$ $II A$ supergravity. It is shown that in both examples wave functions vanish in the limit of the formation of the billiard walls (i.e. we get a quantum resolution of the singularity for $11D$ model) but magnetic branes could not be neglected in calculations of quantum asymptotic solutions while they are irrelevant for classical oscillating behavior when all $120$ electric branes are present.     

A Relationship Between Fixed Time Wiener Measures and Wiener Measures with Fixed Endpoints

Wiener measures are measures on curves that are derived from two-dimensional Brownian motion. We prove a relationship between two types of Wiener measures: measures on paths with fixed starting point (say the origin \(0\) ) and fixed time duration (say \(1\) ); and measures on paths with fixed endpoints (say \(0\) and \(i\) ). The relationship is that if we take a curve from the first type, weight it by a suitable power of the distance to the endpoint of the curve and then apply the conformal map that takes the endpoint to \(i\) , then we get the curve from the second type.     

Asymptotic Behavior for a Version of Directed Percolation on the Triangular Lattice

We consider a version of directed bond percolation on the triangular lattice such that vertical edges are directed upward with probability $y$ , diagonal edges are directed from lower-left to upper-right or lower-right to upper-left with probability $d$ , and horizontal edges are directed rightward with probabilities $x$ and one in alternate rows. Let $\tau (M,N)$ be the probability that there is at least one connected-directed path of occupied edges from $(0,0)$ to $(M,N)$ . For each $x \in [0,1]$ , $y \in [0,1)$ , $d \in [0,1)$ but $(1-y)(1-d) \ne 1$ and aspect ratio $\alpha =M/N$ fixed for the triangular lattice with diagonal edges from lower-left to upper-right, we show that there is an $\alpha _c = (d-y-dy)/[2(d+y-dy)] + [1-(1-d)^2(1-y)^2x]/[2(d+y-dy)^2]$ such that as $N \rightarrow \infty $ , $\tau (M,N)$ is $1$ , $0$ and $1/2$ for $\alpha > \alpha _c$ , $\alpha < \alpha _c$ and $\alpha =\alpha _c$ , respectively. A corresponding result is obtained for the triangular lattice with diagonal edges from lower-right to upper-left. We also investigate the rate of convergence of $\tau (M,N)$ and the asymptotic behavior of $\tau (M_N^-,N)$ and $\tau (M_N^+ ,N)$ where $M_N^-/N\uparrow \alpha _c$ and $M_N^+/N\downarrow \alpha _c$ as $N\uparrow \infty $ .     

Revisiting radiation patterns in e^+e^- collisions

We propose four simple event-shape variables for semi-inclusive $e^+e^- \rightarrow 4$ -jet events. The observables and cuts are designed to be especially sensitive to subleading aspects of the event structure, and allow to test the reliability of phenomenological QCD models in greater detail. Three of them, $\theta _{14}$ , $\theta ^*$ , and $C_2^{(1/5)}$ , focus on soft emissions off 3-jet topologies with a small opening angle, for which coherence effects beyond the leading QCD dipole pattern are expected to be enhanced. A complementary variable, $M_L^2/M_H^2$ , measures the ratio of the hemisphere masses in 4-jet events with a compressed scale hierarchy (Durham $y_{23}$ – $y_{34}$ ), for which subleading $1\rightarrow 3$ splitting effects are expected to be enhanced. We consider several different parton-shower models, spanning both conventional and dipole/antenna ones, all tuned to the same $e^+e^-$ reference data, and show that a measurement of the proposed observables would allow for additional significant discriminating power between the models.     

Study of Alpha Decay Chains of Superheavy Nuclei and Magic Number Beyond Z = 82 and N = 126

We study various $\alpha $ -decay chains on the basis of the preformed cluster decay model. Our work targets the superheavy elements, which are expected to show extra stability at shell closure. Our computations identify the following combinations of proton and neutron numbers as the most stable nuclei: $Z=112$ , $N=161, 163$ ; $Z=114$ , $N=171, 178, 179$ ; and $Z=124$ , $N=194$ . We also investigate the alternative of heavy cluster emissions in the decay chain of ³⁰¹120, instead of $\alpha $ decay. Our study of cluster radioactivity shows that the half-life for ¹⁰Be decay in ²⁸⁹114 is larger, indicating enhanced stability at $Z=114$ , $N=175$ . Similar calculations concerning the emission of $\ ^{14}{\rm C}$ and $\ ^{34}{\rm Si}$ from ³⁰¹120 find the more stable combinations $Z=114$ , $N=173$ , and $Z=106$ , $N=161$ , respectively. From the same parent, ³⁰¹120, the emission of a $\ ^{49-51}{\rm Ca}$ cluster yielding a $Z=100$ , $N=152$ daughter is the most probable.     

Scaling arguments and gravitating solitons

Scaling arguments are generalized to discuss the existence and the stability of static, spherically symmetric self-gravitating solitons in asymptotically flat and asymptotically anti-de Sitter spacetimes in $D$ dimensions. The formulas obtained from the scaling arguments are applied to the Einstein–Yang–Mills solitons with higher order terms in $D=4$ and $D=5$ dimensions and the Einstein–Yang–Mills–Chern–Simons solitons in $D=5$ dimension. We have shown that the negative cosmological constant, in general, improves the stability of the solitons. In $D=4$ , in fact, the negative cosmological constant is essential for the stability of the solitons. In $D=5$ , however, the higher order terms or the Chern–Simons term play more crucial role in the stability of the solitons; the possibility of having stable solitons exists in an asymptotically flat spacetime.     

Charmonium rescattering effects in \psi(2S) \rightarrow \gamma \eta_{c}^{} decay

Charmonium rescattering effects in the M1 transition of $ \psi$ (2S) $ \rightarrow$ $ \gamma$ $ \eta_{c}^{}$ are investigated by modeling a $ \chi_{{cJ}}^{}$ or J/ $ \psi$ rescattering into a $ \eta_{c}^{}$ final state. The absorptive and dispersive part of the transition amplitudes for the rescattering loops of $ \eta$ $ \psi$ ( $ \gamma^{{\ast}}_{}$ ) and $ \gamma$ $ \chi$ ( $ \psi$ ) are separately evaluated. The numerical results show that the contribution from the $ \gamma$ $ \chi$ ( $ \psi$ ) rescattering process is negligible. Compared with the virtual D $ \bar{{D}}$ (D ^*) rescattering processes, the $ \eta$ $ \psi$ ( $ \gamma^{{\ast}}_{}$ ) process may be regarded as the next-leading order of the hadronic loop mechanism, which only offers the partial decay width of ～ 0.045 keV to the $ \psi$ (2S) $ \rightarrow$ $ \gamma$ $ \eta_{c}^{}$ .     

Fluctuation–Dissipation Relation for Systems with Spatially Varying Friction

When a particle diffuses in a medium with spatially dependent friction coefficient \(\alpha (r)\) at constant temperature \(T\) , it drifts toward the low friction end of the system even in the absence of any real physical force \(f\) . This phenomenon, which has been previously studied in the context of non-inertial Brownian dynamics, is termed “spurious drift”, although the drift is real and stems from an inertial effect taking place at the short temporal scales. Here, we study the diffusion of particles in inhomogeneous media within the framework of the inertial Langevin equation. We demonstrate that the quantity which characterizes the dynamics with non-uniform \(\alpha (r)\) is not the displacement of the particle \(\Delta r=r-r^0\) (where \(r^0\) is the initial position), but rather \(\Delta A(r)=A(r)-A(r^0)\) , where \(A(r)\) is the primitive function of \(\alpha (r)\) . We derive expressions relating the mean and variance of \(\Delta A\) to \(f\) , \(T\) , and the duration of the dynamics \(\Delta t\) . For a constant friction coefficient \(\alpha (r)=\alpha \) , these expressions reduce to the well known forms of the force-drift and fluctuation–dissipation relations. We introduce a very accurate method for Langevin dynamics simulations in systems with spatially varying \(\alpha (r)\) , and use the method to validate the newly derived expressions.     

Maximal Distance Travelled by N Vicious Walkers Till Their Survival

We consider N Brownian particles moving on a line starting from initial positions \(\mathbf{{u}}\equiv \{u_1,u_2,\ldots u_N\}\) such that \(0 . Their motion gets stopped at time \(t_s\) when either two of them collide or when the particle closest to the origin hits the origin for the first time. For \(N=2\) , we study the probability distribution function \(p_1(m|\mathbf{{u}})\) and \(p_2(m|\mathbf{{u}})\) of the maximal distance travelled by the \(1^{\text {st}}\) and \(2^{\text {nd}}\) walker till \(t_s\) . For general N particles with identical diffusion constants \(D\) , we show that the probability distribution \(p_N(m|\mathbf{u})\) of the global maximum \(m_N\) , has a power law tail \(p_i(m|\mathbf{{u}}) \sim {N^2B_N\mathcal {F}_{N}(\mathbf{u})}/{m^{\nu _N}}\) with exponent \(\nu _N =N^2+1\) . We obtain explicit expressions of the function \(\mathcal {F}_{N}(\mathbf{u})\) and of the N dependent amplitude \(B_N\) which we also analyze for large N using techniques from random matrix theory. We verify our analytical results through direct numerical simulations.     

A new point of view on the fifth force in 4D physics

The purpose of the paper is to define and study, in the presence of the electromagnetic potentials, the fifth force induced by the fifth dimension in the $4D$ physics. A $4D$ tensor calculus and the Riemannian horizontal connection on a $5D$ general Kaluza–Klein space, are the main tools in our study. We prove that the fifth force does not violate the basic concepts in $4D$ physics. Our approach is different from what is known so far in literature.     

Reaction pp \rightarrow pp \pi \pi \pi as a background for hadronic decays of the \eta^{{\prime}}_{} meson

Isospin violating hadronic decays of the $ \eta$ and $ \eta{^\prime}$ mesons into 3 $ \pi$ mesons are driven by a term in the QCD Lagrangian proportional to the mass difference of the d and u quarks. The source giving large yield of the mesons for such decay studies are pp interactions close to the respective kinematical thresholds. The most important physics background for $ \eta$ , $ \eta{^\prime}$ $ \rightarrow$ $ \pi$ $ \pi$ $ \pi$ is coming from direct three-pion production reactions. In case of the $ \eta$ meson the background for the decays is relatively low ( $ \approx$ 10% . The purpose of this article is to provide an estimate of the direct pion production background for the $ \eta{^\prime}$ $ \rightarrow$ 3 $ \pi$ decays. Using the inclusive data from the COSY-11 experiment we have extracted the differential cross-section for the pp $ \rightarrow$ pp -multipion production reactions with the invariant mass of the pions equal to the $ \eta{^\prime}$ meson mass and estimated an upper limit for the signal to background ratio for studies of the $ \eta{^\prime}$ $ \rightarrow$ $ \pi^{+}_{}$ $ \pi^{-}_{}$ $ \pi^{0}_{}$ decay.     

Optimal N-Point Configurations on the Sphere: “Magic” Numbers and Smale’s 7th Problem

This paper inquires into the concavity of the map \(N\mapsto v_s(N)\) from the integers \(N\ge 2\) into the minimal average standardized Riesz pair-energies \(v_s(N)\) of \(N\) -point configurations on the sphere \(\mathbb {S}^2\) for various \(s\in \mathbb {R}\) . The standardized Riesz pair-energy of a pair of points on \(\mathbb {S}^2\) a chordal distance \(r\) apart is \(V_s(r)= s^{-1}\left( r^{-s}-1 \right) \) , \(s \ne 0\) , which becomes \(V_0(r) = \ln \frac{1}{r}\) in the limit \(s\rightarrow 0\) . Averaging it over the \(\left( \begin{array}{c} N\\ 2\end{array}\right) \) distinct pairs in a configuration and minimizing over all possible \(N\) -point configurations defines \(v_s(N)\) . It is known that \(N\mapsto v_s(N)\) is strictly increasing for each \(s\in \mathbb {R}\) , and for \(s<2\) also bounded above, thus “overall concave.” It is (easily) proved that \(N\mapsto v_{-2}^{}(N)\) is even locally strictly concave, and that so is the map \(2n\mapsto v_s(2n)\) for \(s<-2\) . By analyzing computer-experimental data of putatively minimal average Riesz pair-energies \(v_s^x(N)\) for \(s\in \{-1,0,1,2,3\}\) and \(N\in \{2,\ldots ,200\}\) , it is found that the map \(N\mapsto {v}_{-1}^x(N)\) is locally strictly concave, while \(N\mapsto {v}_s^x(N)\) is not always locally strictly concave for \(s\in \{0,1,2,3\}\) : concavity defects occur whenever \(N\in {\mathcal {C}}^{x}_+(s)\) (an \(s\) -specific empirical set of integers). It is found that the empirical map \(s\mapsto {\mathcal {C}}^{x}_+(s),\ s\in \{-2,-1,0,1,2,3\}\) , is set-theoretically increasing; moreover, the percentage of odd numbers in \({\mathcal {C}}^{x}_+(s),\ s\in \{0,1,2,3\}\) is found to increase with \(s\) . The integers in \({\mathcal {C}}^{x}_+(0)\) are few and far between, forming a curious sequence of numbers, reminiscent of the “magic numbers” in nuclear physics. It is conjectured that these new “magic numbers” are associated with optimally symmetric optimal-log-energy \(N\) -point configurations on \(\mathbb {S}^2\) . A list of interesting open problems is extracted from the empirical findings, and some rigorous first steps toward their solutions are presented. It is emphasized how concavity can assist in the solution to Smale’s \(7\) th Problem, which asks for an efficient algorithm to find near-optimal \(N\) -point configurations on \(\mathbb {S}^2\) and higher-dimensional spheres.     

Semiclassical strings in supergravity PFT

In this paper, we introduce the bulk viscosity in the formalism of modified gravity theory in which the gravitational action contains a general function \(f(R,T)\) , where \(R\) and \(T\) denote the curvature scalar and the trace of the energy–momentum tensor, respectively, within the framework of a flat Friedmann–Robertson–Walker model. As an equation of state for a prefect fluid, we take \(p=(\gamma -1)\rho \) , where \(0 \le \gamma \le 2\) and a viscous term as a bulk viscosity due to the isotropic model, of the form \(\zeta =\zeta _{0}+\zeta _{1}H\) , where \(\zeta _{0}\) and \(\zeta _{1}\) are constants, and \(H\) is the Hubble parameter. The exact non-singular solutions to the corresponding field equations are obtained with non-viscous and viscous fluids, respectively, by assuming a simplest particular model of the form of \(f(R,T) = R+2f(T)\) , where \(f(T)=\alpha T\) ( \(\alpha \) is a constant). A big-rip singularity is also observed for \(\gamma <0\) at a finite value of cosmic time under certain constraints. We study all possible scenarios with the possible positive and negative ranges of \(\alpha \) to analyze the expansion history of the universe. It is observed that the universe accelerates or exhibits a transition from a decelerated phase to an accelerated phase under certain constraints of \(\zeta _0\) and \(\zeta _1\) . We compare the viscous models with the non-viscous one through the graph plotted between the scale factor and cosmic time and find that the bulk viscosity plays a major role in the expansion of the universe. A similar graph is plotted for the deceleration parameter with non-viscous and viscous fluids and we find a transition from decelerated to accelerated phase with some form of bulk viscosity.     

Mössbauer spectroscopy in the investigation of new mineral–related materials

New materials based on the composition of the mineral schafarzikite, FeSb $_{2}\textit {O}_{4}$ , have been synthesised. $^{57}$ Fe- and $^{121}$ Sb- Mössbauer spectroscopy shows that iron is present as Fe $^{2+}$ and that antimony is present as Sb $^{3+}$ . The presence of Pb $^{2+}$ on the antimony sites in materials of composition FeSb $_{1.5}$ Pb $_{0.5}\textit {O}_{4}$ induces partial oxidation of Fe $^{2+}_{}$ to Fe $^{3+}$ . The quasi-one-dimensional magnetic structure of schafarzikite is retained in FeSb $_{1.5}$ Pb $_{0.5}\textit {O}_{4}$ and gives rise to weakly coupled non-magnetic Fe $^{2+}$ ions coexisting with Fe $^{3+}$ ions in a magnetically ordered state. A similar model can be applied to account for the spectra recorded from the compound Co $_{0.5}$ Fe $_{0.5}$ Sb $_{1.5}$ Pb $_{0.5}\textit {O}_{4}$ .     

P–V criticality of topological black holes in Lovelock–Born–Infeld gravity

To understand the effect of third order Lovelock gravity, $P$ – $V$ criticality of topological AdS black holes in Lovelock–Born–Infeld gravity is investigated. The thermodynamics is further explored with some more extensions and in some more detail than the previous literature. A detailed analysis of the limit case $\beta \rightarrow \infty $ is performed for the seven-dimensional black holes. It is shown that, for the spherical topology, $P$ – $V$ criticality exists for both the uncharged and the charged cases. Our results demonstrate again that the charge is not the indispensable condition of $P$ – $V$ criticality. It may be attributed to the effect of higher derivative terms of the curvature because similar phenomenon was also found for Gauss–Bonnet black holes. For $k=0$ , there would be no $P$ – $V$ criticality. Interesting findings occur in the case $k=-1$ , in which positive solutions of critical points are found for both the uncharged and the charged cases. However, the $P$ – $v$ diagram is quite strange. To check whether these findings are physical, we give the analysis on the non-negative definiteness condition of the entropy. It is shown that, for any nontrivial value of $\alpha $ , the entropy is always positive for any specific volume $v$ . Since no $P$ – $V$ criticality exists for $k=-1$ in Einstein gravity and Gauss–Bonnet gravity, we can relate our findings with the peculiar property of third order Lovelock gravity. The entropy in third order Lovelock gravity consists of extra terms which are absent in the Gauss–Bonnet black holes, which makes the critical points satisfy the constraint of non-negative definiteness condition of the entropy. We also check the Gibbs free energy graph and “swallow tail” behavior can be observed. Moreover, the effect of nonlinear electrodynamics is also included in our research.