Limit theorems for the cubic mean-field Ising model

We study a mean-field spin model with three- and two-body interactions. The equilibrium measure for large volumes is shown to have three pure states, the phases of the model. They include the two with opposite magnetization and an unpolarized one with zero magnetization, merging at the critical point. We prove that the central limit theorem holds for a suitably rescaled magnetization, while its violation with the typical quartic behavior appears at the critical point.


Introduction
In this paper, we investigate the mean-field Ising spin model with quadratic and cubic interactions.The interest in such a model comes from two large fields of research.The first is condensed matter physics, where the three-body interaction plays a role in the description of the phase separation phenomena of some magnetic alloys [1] lacking spin-flip symmetry.Those physical systems cannot be described by the sole use of a two-body interaction, while a three-body term captures some features of their behavior [2].This fact is well paralleled by the Ginibre theorem about functions of spin configurations that are fully classified by an orthonormal base of k-body interactions [3].Those physical phenomena are well described by statistical mechanics models on regular lattices in finite (d=2,3) dimensions.While some of those models have an exact solution in very special cases [4,5], it is well known that the mean-field approximation provides an analytically viable setting and a fair description of the phase separation.In those cases, the term mean-field approximation is understood in the sense of a special class of probability measure where the Boltzmann-Gibbs variational principle is optimised: instead of minimizing the free energy over all probability measures, one restricts it to product measures on single spins [6,7].
The other field in which the three-body interactions came to play a role is that of the applications to complex systems, in particular those of socio-technical nature [20] where the social network structure with long-range interaction represents a realistic description of the phenomenon and not an approximation of its finite-dimensional version [8][9][10][11].In this case, from a mathematical perspective, the introduction of the three-body interaction entails moving from a graph-theoretical environment of vertices and edges to a richer hypergraph setting where the three-body terms, representing the faces of the hypergraph, are also taken into account.
The presence of the cubic interactions brings technical difficulties in the analysis of the model.In particular, the non-convex energy contribution due to the cubic power prevents the use of the Hubbard-Stratonovich transform, which instead is very efficient in the case of quadratic interactions.More precisely, even if the thermodynamic limit of the free energy can be easily computed by large deviation arguments, the fluctuations of the order parameter cannot be analysed with the classical rigorous methods for a mean-field system with pairwise interaction [12][13][14].In order to overcome this obstacle we need a fine control on the N -asymptotic behavior of the partition function that is obtained by a method similar to that recently introduced in [15].This paper presents a rigorous analysis of the mean-field model with three-and two-body interactions in a zero magnetic field.We show that the infinite-volume properties of the model display new phenomena that are absent in the quadratic mean-field case.In particular, we prove that the equilibria of the system include not only positively and negatively polarized states but also an unpolarized stable state in the presence of a non-zero cubic term that breaks the spin-flip symmetry.Finally, we also study the fluctuation of the magnetization in the entire phase space, specifying the behavior at phase separation and at the critical point.The critical exponent for the magnetization, moreover, takes on a value of zero towards the unpolarized directions of the phase space, and phase transitions can occur in the antiferromagnetic region.
This paper is organised as follows: Section 2 contains the formal definition of the model as well as a statement of the main results.In Section 3, we study the properties of the consistency equation that describes the system in its stationary equilibrium state.These properties provide an analytical description of the system's phase diagram and the magnetization's limiting behavior, as well as the computation of the critical exponents.Finally, Section 4 contains conclusions and perspectives, and the Appendices A and B contain technical and concentration results used throughout the work.

Definitions and main results
Let us consider N spins σ = (σ i ) i≤N ∈ {−1, +1} N interacting through an Hamiltonian of the form where (K, J, h) ∈ R 3 , K and J tune the interactions among triples and pairs of spins, respectively, while h represents an external field acting on the system.When K = 0, the previous Hamiltonian reduces to the well-know Curie-Weiss case.In this work we will concentrate on the case h = 0 and use the parameter K as a spin-flip symmetry breaking term reducing (1) to an Hamiltonian that can be represented as where m N is the magnetization per particle: The expression (2) highlights the mean-field nature of the model.The Boltzmann-Gibbs probability measure associated to where is the partition function.In equation ( 4), we set the usual inverse temperature β to 1 without loss since it has been reabsorbed in the parameters of the model.Notice that since the Hamiltonian (2) is invariant under the transformation K → −K, and σ i → −σ i for i = 1, ..., N , one can study the model only for K > 0 without loss.Our aim is to obtain a complete characterization of the model's phase diagram, an analysis of the asymptotic distribution of the magnetization in the presence and absence of phase transitions, the fluctuations of the suitably rescaled magnetization (3) w.r.t. the Boltzmann-Gibbs measure (4) at and away from the critical point, and the computation of the critical exponents.
All the above properties are strictly related to the analytical properties of the free energy of the system, which is the starting point of our analysis.Let us define the thermodynamic pressure, i.e., the generating functional as: Notice that p N equals the free energy up to a minus sign.The thermodynamic limit of ( 5) can be easily computed applying Varadhan's integral lemma [12,18], obtaining: Proposition 2.1.Given (K, J) ∈ R 2 the limiting pressure of (5) admits the following variational representation: where is the energy contribution and is the binary entropy contribution.
The critical points of (6) satisfy the consistency equation, A careful analysis shows that, among the solutions of ( 9), the function ϕ(m) in ( 6) can have one or two global maximizers in the interval (−1, 1) for fixed (K, J) (see Figure 1).
In particular, we can divide the parameter space (K, J) ∈ R + × R accordingly to the following: Proposition 2.2 (Phase diagram).For any K > 0, there exists J = γ(K) defined in Proposition 3.3 such that the function m → ϕ(m) has a unique maximum point m * for (K, J) ∈ (R + × R)\γ.Moreover, on the curve γ there are two global maximizers, 0 = m 0 < m 1 and the limit as K → 0 of γ(K) identifies the critical point (K c , J c ) = (0, 1) where the magnetization takes the value m c = 0. Figure 1: Stable solutions of the mean-field equation as a function of K and J.There are three stable phases presented here: the positive polarized phase depicted in red, the unpolarized phase given as the gray plateau, and the negative polarized phase denoted by the blue color.At the critical point, (K, J) = (0, 1), the three phases of the cubic model as well as the two phases of the Curie-Weiss plane (K = 0) coalesce.
In physical terms, the presence of two global maximizers corresponds to the existence of two different thermodynamic equilibrium phases, whereas the curve γ represents the coexistence curve.Let's note that m 0 and m 1 represent a stable paramagnetic state and a positively polarised state, respectively.The paramagnetic state is characterized by the absence of spontaneous magnetic order and the presence of symmetry between the up and down spin, with no preference for either direction.The jump from the paramagnetic state to the polarized state, namely when the magnetization jumps from m 0 to m 1 , represents a first-order phase transition, which is markedly different from the quadratic mean-field model (K = 0) having a secondorder phase transition in J.More precisely if we denotes by m * (K, J) the unique maximizer of ϕ, for any K > 0 there exists J = γ( K) ∈ (−∞, 1) such that 0 = lim This behavior is somehow reminiscent of the Curie-Weiss Potts model analyzed in [19] where for any value of the parameter q a first order phase transition is observed.
Numerical simulations of the phase diagram described in Proposition 2.2 can be seen in Figure 2. In the standard Curie-Weiss model, when J > 0 we know that as soon as h > 0 one obtains a positive magnetization.The reason is that the energy contribution due to h favors only spins aligned with sign(h).On the contrary, in our system, J, K > 0, the energy contribution due to K can be minimized by configurations containing both up and down spin signs.This implies that the entropy contribution can dominate also for small but non-zero K, giving a zero magnetization.
The next theorem contains the law of large numbers and the central limit theorem for the distribution of m N (σ) with respect to the Boltzmann-Gibbs measure.
Finally, we study the behavior of the limiting value of the magnetization near the critical point (K c , J c ) = (0, 1) namely the critical exponents of the model.The average value of the magnetization is given by the LLN in Theorem 2.1 and will be denoted by m * (K, J).The following proposition describes the critical behavior of m * (K, J) when (K, J) → (K c , J c ) from various directions.Proposition 2.3.Let m * (K, J) be the unique maximizer of ϕ(m) defined in Corollary 3.1.Given α ∈ R consider the lines and the function m * (K) ≡ m * (K, J(K)).Then, for K → 0 + , the following holds Remark 2.2.Notice that when α < 0 the critical exponent is 0. The case K = 0 and J → 1 + corresponds to the classical Curie-Weiss model and is well known that 3 Proofs This section contains the proofs of the above results and is organised as follows: In Section 3.1, we prove Proposition 2.2 by studying the properties of the function ϕ(m) appearing in the variational problem (6).Section 3.2 contains the proof of Theorem 2.1 and is based on the asymptotic expansion given in Appendix B. Finally, in Section 3.3, we derive the critical exponents of the model.

Proof of Proposition 2.2
The complete proof of Proposition 2.2 follows from Propositions 3.1, 3.2, 3.3 and 3.4 below.
Let us start studying in detail the variational principle (6) and observe that the function ϕ(m) satisfies Therefore the variational pressure ϕ(m) attains it maximum in at least one point m = m(K, J) ∈ (−1, 1), which satisfy Indeed, from (20) where g(m, K) := arctanh(m) − Km 2 and set J c = 1.Then: a. for J < Ψ(K), there exists a unique solution, m 0 = 0, and it is the maximum Proof.Let us start by noticing that m = 0 is always a solution of (21).Moreover, Now, let's rewrite (21) as The solutions of (21) are the intersections between the line mJ and the function g(m, K).Therefore the function Ψ(K) in ( 22) is a benchmark to study the number of solutions of ϕ ′ (m) = 0 when J varies.Indeed by definition, Ψ(K) represents the smallest value of J in order to have a positive solution for (23).Let us start collecting some properties of the function g(m, K).By definition we have that This implies that, Since the function m → 2m (1 − m 2 ) 2 is strictly increasing on [0, 1), then g ′′ (m, K) = 0 has only one solution, namely g(m, K) has only one inflection point.Moreover, observe that, as m → 1 − , g(m, K) → +∞.a.If J < Ψ(K) then it's clear that (21) has a unique solution m 0 = 0 which is a maximum point since in this case ϕ ′′ (0) < 0 .
b.If Ψ(K) < J < J c , continuity of g and the fact that for m → 1 − , g(m, K) → +∞, imply that (21) has three solutions, m 0 , m 1 and m 3 , where m 1 and m 3 are positive.It's also easy to check using the properties of the function g(m, K) that m 0 and m 1 are local maxima while m 3 is a local minima.
c.If J = Ψ(K), then there is only one intersection point m 4 between the line mJ and the function g(m, K).Standard reasoning allows to conclude that m 4 is an inflection point for ϕ.
d. Finally suppose that J ≥ J c .The fact that g ′ (0, K) = 1 and g ′′ (0, K) = −2K < 0 for K > 0, means that the line mJ starts above the function g.Now, since g has at most one inflection point and g(m, K) → +∞ as m → 1 − , one can conclude that there exist a unique positive solution m 2 ∈ (0, 1) of ϕ ′ (m) = 0.
The solutions made mention in Proposition 3.1 are displayed in Figure 3.
In the next proposition we obtain the differentiability of the solution(s) of the consistency equation (21) with respect to the parameters J and K. (a) m 1 is continuous in its domain namely Moreover, for any i ∈ {0, 1, 2} it holds that Remark 3.1.Notice that (b) implies that there are no degenerate maximum points of ϕ(m) for K > 0. Therefore the only degenerate maximum is obtained for (K, J) = (K c , J c ) = (0, 1), that is the critical point of a Curie-Weiss model, here the magnetization takes the value m c = 0.
Proof.(a) Let's start with m 1 and take (K, J) in its domain, namely . Observe from ( 21) that, Hence, m 1 is the unique maximum point of φ(m), then by the Berge's maximum theorem A.1 (see [21,22]), m 1 is continuous for (K, J) ∈ D. To prove the smoothness of m 1 on the interior of its domain it's enough to show that ϕ ′′ (m 1 ) < 0 and then apply the implicit function theorem A.2 (see [22,23]).Let G(m) := ϕ ′′ (m) then, and hence, We want to prove that G(m 1 ) < 0 if Ψ(K) < J < J c .Clearly since m 1 is a local maximizer it's enough to show that G(m 1 ) ̸ = 0. Recall that m 1 is the biggest positive solution of ϕ ′ (m) = 0. It's easy to check that G(m) = 0 has at most two solutions.Assume by contradiction that G(m 1 ) = 0 if Ψ(K) < J < J c , then G(m) < 0 or G(m) > 0 in a left neighbourhood of m 1 .
• Suppose that G(m) < 0 in a left neighbourhood of m 1 then G(m) cannot be always negative, otherwise ϕ ′ (m) is decreasing and, since ϕ ′ (0) = 0 then ϕ ′ (m) = 0 can not have more than one solution.This contradicts point b) of Proposition 3.1.Therefore there exist an interval where G(m) > 0 but keeping in mind the properties of G in (27) and the fact that G is continuous, this implies that there are at least three solutions for G(m) = 0, but this is impossible because we already observed that G(m) = 0 has at most two solutions.
Therefore by the implicit function theorem A.2, since ϕ ′′ (m) ̸ = 0 on the interior of the domains of m 1 and m 2 , we can conclude that m 1 and m 2 are C ∞ .
(b) We already proved that for any i ∈ {0, 1, 2}, ϕ ′′ (m i ) < 0 for suitable K, J.For the second part a direct computation shows that: and similarly, Using the fact that m i , i = {0, 1, 2} are the stationary points of ϕ(•), we have that and similarly for ∂m i ∂J one obtains and this concludes the proof.Now we study which of the stationary points described by Proposition 3.1 are global maximizers of ϕ(m) and show the existence of a phase transition.These stationary points are: m 0 , m 1 , and m 2 .Let us start by recalling the result of Proposition 3.1:

the only the global maximum point of ϕ(m)
To identify the coexistence of two global maximum points of ϕ(m) when Ψ(K) < J < J c , consider the following function: Notice that ∆(K, J) can be extended by continuity at J = Ψ(K) and J = J c .In the above equation we use ϕ(•, K, J) to emphasis the dependence of ϕ on the parameters.
) and m 2 (K, 1) is the only global maximum for ϕ(m, K, J).Now, by continuity of ϕ(m) and m 1 , we have that J → ∆(K, J) is a continuous function, and then the existence of the wall J = γ(K) follows from the application of the intermediate value theorem.For the uniqueness part we observe that J → ∆(K, J) is strictly increasing.Indeed from Proposition 3.2 we know that ϕ(m 1 ), m 1 are smooth functions and for J ∈ (Ψ(K), J c ).
Note that on the curve γ there are two global maximum points of ϕ(m).Let us define Therefore by Proposition 3.2 one can conclude that m * (K, J) is continuous on its domain (R + × R)\γ and it is C ∞ on (R + × R) \ γ.Moreover the following holds: ) and at least C 1 for K = 0.In particular, Proof.i.We begin by showing that γ(K) where ∆ is defined by equation (32) for Ψ(K) ≤ J < J c and K > 0. Furthermore, observe that ∆ is C ∞ in its domain by the smoothness of ϕ and m 1 .Recall from the proof of Proposition 3.3 that hence, by the implicit function theorem A.
From equations ( 28) and ( 29), we have that, A consequence of this property is that also when J < 0 (antiferromagnetic case) and very large there is always going to be phase transition between a polarized and unpolarized state.
ii.Now we prove that the extended function γ ∈ C 1 (R + ).Recall that γ(K) ∈ [Ψ(K), J c ] and observe that lim K→K + c Ψ(K) = J c then lim which implies that γ is continuous at K c .Now we have that which implies that γ ′ (K c ) = − 2 3 m c = 0 by the application of mean value theorem.

Proof of Theorem 2.1
In this section we provide the details of the proof for Theorem 2.1 following closely the argument in [15].
then ϕ(m) has a unique global maximizer m * with ϕ ′′ (m * ) < 0. It's easy to check that ϕ(m) satisfies the hypothesis of Lemma B.1, therefore (64) gives concentration inequality for m N in a suitable neighbourhood of m * under the probability measure (4).More precisely, for any α ∈ (0, 1  6 ] and N large enough one has where Therefore the convergence in distribution (10) follows from (43) by standard approximation arguments.
To obtain the central limit for m N , it is enough to compute the limit of the moment generating function of the rescaled random variable N 1 2 (m N − m * ).For a fixed t ∈ R, the moment generating function of N 1 2 (m N − m * ) can be expressed as which implies (11).
2. Let's recall that by Proposition 2.2 there exist two global maximizers m i of ϕ(m) for i ∈ {0, 1} on γ.Moreover by point b) of Proposition 3.2 we know that ϕ ′′ (m i ) < 0 for i ∈ {0, 1}.Now, following the same argument as before, formula (73) in Lemma B.2 gives the concentration inequality for m N within a suitable neighbourhood of m i with respect to the Gibbs measure (4).Therefore the convergence in distribution ( 13) and ( 12) follows the asymptotic expansions of the (restricted) partition function in Lemma B.2.
To obtain the local central limit theorem for m N around the global maximizers m i , we will show that the moment generating function of 1] is such that m i is the unique maximizer of ϕ(m) on its interior.The moment generating function of Following the asymptotic expansion of the partition function in (74) (see Lemma B.2), the fraction on the right side of equation (52) reduces to Now, taking Taylor's expansion of ϕ N (m i,N (t)) at m i up to the second order, one can repeat the same arguments as in the unique maximum case, obtaining This completes the proof of ( 14).
3. Notice that the critical point (K c , J c ) = (0, 1) is a degenerate maximum point for ϕ(m) in the sense that ϕ ′′ (m * (K, J)) (K,J)=(0,1) = 0.This does not allow the use of the asymptotic expansions in Lemma B.1.However, one can simply notice that the Hamiltonian H N of the model at the critical point (K c , J c ) = (0, 1) coincides at any N ∈ N with the Hamiltonian function of the standard Curie-Weiss model at the critical temperature J = 1 and zero external field.Therefore (15) and ( 16) are a well known results and their proof can be found in [13].

Proof of Proposition 2.3
Proof.Let us start with the case α ≥ 0. This implies from equation (17) that J(K) ≥ J c = 1 and then m * (K) ≡ m 2 (K, J(K)) where m 2 is the only positive solution of the consistency equation (21).
Clearly m * (K) → 0 as K → 0 + , hence by Taylor's expansion we have that Hence From the above equation, neglecting higher order corrections we have Now, if α > 0 then Otherwise if α = 0, then Let's turn on the case α < 0. From Proposition 3.4 we know that γ(K) is at least C 1 at K = 0. Since lim K→0 + γ ′ (K) = 0 we know that if J(K) < γ(K) for K small enough, then m * (K) ≡ m 0 (K, J) = 0.

Conclusion and perspectives
In this work, we have studied how the three-body interaction, which provides a spinflip symmetry-breaking parameter, induces phase transitions with novel properties in the mean-field setting.In particular, we derived all the critical exponents and the limiting distribution of a suitably rescaled magnetization in the entire phase space.The presence of a stable paramagnetic phase and the fact that, also in the antiferromagnetic regime, the model presents phase transitions and phase coexistence are interesting for applications in socio-technical environments [20] and possibly in other fields [24,25].
A possible research development will be to extend the results of the present work to multi-populated models [20,[29][30][31][32][33][34][35][36][37].In these models, the invariance of the Hamiltonian with respect to permutations among sites is replaced by a weaker one that takes into account the existence of different species of spins.This setting is particularly useful in social science applications [20,32,35,36].Moreover, as mentioned in the introduction, the mean-field approximation involved in the study of some finitedimensional lattices provides a natural emergence of the multi-populated models.It is well known, for instance, that a system on a simple cubic lattice [26,27] with ferromagnetic and antiferromagnetic couplings has a factorized equilibrium measure that corresponds to a two-populated mean-field model.Similarly, it has been shown in [7] that on a regular square lattice, a system with cubic interaction has a product state equilibrium described by a two-populated mean-field model, while on a regular triangular lattice [6], by a three-populated mean-field model.
We also mention that in the case of quadratic interaction, Stein's method provides stronger results (Berry-Esseen type bounds) on the rate at which the convergence to the normal distribution takes place (see [16,17]).The extension of the above method to our model and more generally to higher order interaction is an interesting open problem.We plan to develop those research directions in the future.

A Technical results
This section of the appendix presents some useful technical results applied in the work.We begin by stating the Berge's maximum theorem in the following Proposition.
(a) The following function is continuous: (b) Suppose that for all x, y ∈ R n the function v → f (v, x) achieves its maximum on [−1, c(y)] in a unique point.Then also the following function is continuous: The following proposition partially states Dini's implicit function theorem.Then we provide two simple corollaries that are used in the paper.

B Concentration results and asymptotic expansions
In this section of the appendix, we state concentration properties of the magnetization and asymptotic expansions of the partition function for a large class of Ising mean-field models and give proofs using the same methods and arguments recently introduced in [15].Consider a mean-field spin model with energy density f N , namely where m N = 1 N i≤N σ i is the magnetization density.We assume that (f N ) is a sequence of continuous functions f N : [−1, 1] N → R converging uniformly to f .We assume also that f N has bounded derivatives up to order 4 converging uniformly to f ′ , f ′′ , f ′′′ , f ′′′′ .We denote the law of the magnetization under the Gibbs measure by We begin by applying the Laplace approximation of an integral over a shrinking interval B N,α via the Riemann approximation of the sum in equation ( 68) with an integral and the binomial coefficient can be approximated by the Stirling's approximation method.Notice that by the Riemann approximation (see Appendix Lemma A.2 and B.7 of [15]) of the sum, we have that and the binomial coefficient in (68) can be approximated as It follows from equations ( 69) and ( 70) and the Laplace approximation (see Appendix Lemma A.3 of [15]) of an integral over a shrinking interval B N,α that: e N ϕ N (x) (1 + O(N −1 ))dx Therefore, This completes the proof of Lemma B.1.
Lemma B.2. Suppose ϕ(x) has S ∈ N global maximizers x i such that ϕ ′′ (x i ) < 0. For i ≤ S, let A i ⊂ [−1, 1] be an interval such that x i ∈ int(A i ) is the unique maximizer of ϕ on cl(A i ).Then for N large enough ϕ N has a unique global maximizer x i,N → x i on A i with ϕ ′′ N (x i,N ) < 0 and for α ∈ 0, 1  6 , one has 2 ) (73 where B N,α,S = i≤S B N,α (x i,N ), moreover the restricted partition function on A i can be expanded as, and the unrestricted partition function can be expanded as, Note that, here, int(A i ) and cl(A i ) denote the interior and closure of A i , respectively.
Proof.The fact that for N large enough ϕ N has a unique maximizer x i,N → x i with ϕ ′′ (x i,N ) < 0 can be proved applying to the function ϕ N restricted to cl(A i ) and using the same argument of Lemma B.1.
Clearly, for N large enough, B N,α (x i,N ) ⊂ A i and 2 ) (76) following a step-by-step argument used to prove equation (64).Now, for i ≤ S and N large enough, one has that A i \B N,α (x i,N ) = A i \B N,α,S and then µ N (m N (σ) ∈ B c N,α (x i,N ) m N (σ) 2 ). (77) This completes the proof of equation ( 73) following from equation (77).
The proof for the asymptotic expansion of the partition function when there are multiple global maximizers of ϕ follows exactly the same argument for the case with unique global maximizer.Note that for fixed i ≤ S and N large, m N (σ) concentrates around x i ∈ A i as it was shown in equation (76).Hence, exp N f N (x) .
(78) Now, following the exact computation and argument in Lemma B.1, we have that the restricted partition function for each of the global maximizers x i can be expanded as e N ϕ N (x i,N ) (x 2 i,N − 1)ϕ ′′ (x i,N ) Assuming that m N (σ) concentrates around S global maximizers x i,N for i ≤ S then, equation (75) follows from (79).Hence, we have

Figure 2 :
Figure 2: Phase diagram of the model with coexistence curve γ and the critical point (K c , J c ) in the (K, J) plane.
Therefore, there exists ϵ > 0 such that ϕ(m) is strictly increasing on [−1, −1 + ϵ] and strictly decreasing on [1 − ϵ, 1].This implies that, the local maximizers of ϕ(m) does not include −1 and +1.Notice also that, since K > 0, if m > 0 then ϕ( m) > ϕ(− m) therefore the supremum of ϕ(m) cannot be reached at negative values.A complete classification of the critical points of ϕ(m) is contained in the following proposition: Proposition 3.1.(Classification of critical points) For all K > 0 and J ∈ R, the solutions to equation (21) can be described as follow:Define the function ) has three solutions, i.e., m 0 , m 1 > m 3 > 0. Furthermore, m 0 , m 1 are local maximum points while m 3 is a local minimum point of ϕ(m), c. for J = Ψ(K), there exist two solutions, m 0 and m 1 > 0.Where m 0 is the maximum point of ϕ(m) and m 1 is an inflection point.d.If J ≥ J c , there exists a unique positive solution m 2 which is the only maximum point of ϕ(m) in equation (6).

Figure 3 :
Figure 3: The points of intersection between the blue curve g(m, K) as defined in (23) and red curve f (m) = Jm.The solution of the equation (21) are the points of intersection between g(m, K) and f (m).