This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Upper bound on the second derivative of the quenched pressure in spin-glass models: weak Griffiths second inequality

Manaka Okuyama1 and Masayuki Ohzeki1,2,3 [email protected]1Graduate School of Information Sciences1Graduate School of Information Sciences Tohoku University Tohoku University Sendai 980-8579 Sendai 980-8579 Japan
2Institute of Innovative Research Japan
2Institute of Innovative Research Tokyo Institute of Technology Tokyo Institute of Technology Yokohama 226-8503 Yokohama 226-8503 Japan
3Sigma-i Co. Japan
3Sigma-i Co. Ltd. Ltd. Tokyo 108-0075 Tokyo 108-0075 Japan Japan
Abstract

The Griffiths first and second inequalities have played an important role in the analysis of ferromagnetic models. In spin-glass models, although the counterpart of the Griffiths first inequality has been obtained, the counterpart of the Griffiths second inequality has not been established. In this study, we generalize the method in the previous work [J. Phys. Soc. Jpn. 76, 074711 (2007)] to the case with multi variables for both symmetric and non-symmetric distributions of the interactions, and derive some correlation inequalities for spin-glass models. Furthermore, by combining the acquired equalities in symmetric distributions, we show that there is a non-trivial positive upper bound on the second derivative of the quenched pressure with respect to the strength of the randomness, which is a weak result of the counterpart of the Griffiths second inequality in spin-glass models for general symmetric distributions.

1 Introduction

The Griffiths inequalities, which characterize that the correlation of ferromagnetic models is always positive, are very important and indispensable in the rigorous analysis of ferromagnetic models. Towards rigorous analyses of finite-dimensional spin-glass models, there are some previous studies [1, 2, 3, 4, 6, 5, 7] which aim to establish the counterpart of the Griffiths inequalities in spin-glass models. These attempts have been partially successful, and, when the distribution function of the interactions is symmetric, the counterpart of the Griffiths first inequality has been obtained in spin-glass models [1, 2]. Besides, it was shown that the ferromagnetic counterpart of the Griffiths second inequality [3, 4] holds on the Nishimori line [8]. However, no counterpart of the Griffiths second inequality has been established for other parameter regions [5]. Rigorous analyses based on the concept of correlation inequalities have not been sufficiently advanced [6, 7].

On the other hand, by focusing only on the distribution of a single interaction among all of the interactions, the previous studies [9, 10] showed that the quenched average of the local energy for Ising models with quenched randomness is always larger than or equal to one in the absence of all the other interactions. Although this is a non-trivial result that holds regardless of any other interaction, an extension to the case with multi variables was not apparent. In this study, we give another derivation of their result [9, 10] and extend it to the case with multi variables. Then, we obtain some correlation inequalities for the Ising models with quenched randomness. Moreover, combining the acquired inequality and our previous work [11] for symmetric distributions, we find that there is a simple upper bound on the second derivative of the quenched pressure with respect to the strength of the randomness. This bound can be regarded as a weak result of the counterpart of the Griffiths second inequality in spin-glass models for general symmetric distributions.

The organization of the paper is as follows. In Sec. II, we define the model and explain the counterpart of the Griffiths inequalities in spin-glass models. In Sec. III, we prove some inequalities for quenched averages in preparation for the next section. Section IV is devoted to obtaining some correlation inequalities for spin glass, which is a systematic extension of the previous study [9, 10]. Furthermore, we provide a non-trivial upper bound on the second derivative of the quenched pressure with respect to the strength of the randomness. Finally, our conclusion is given in Sec. V.

2 Ising model with quenched randomness and counterpart of Griffiths inequalities

Following Ref. [9], we consider a generic form of the Ising model,

H\displaystyle H =\displaystyle= AVλAJAσA,\displaystyle-\sum_{A\subset{V}}\lambda_{A}J_{A}\sigma_{A}, (1)
σA\displaystyle\sigma_{A} \displaystyle\equiv iAσi,\displaystyle\prod_{i\in A}\sigma_{i}, (2)

where VV is the set of sites, the sum over AA is over all the subsets of VV in which interactions exist, and the lattice structure adopts any form. The probability distribution of a random interaction JAJ_{A} is represented as PA(JA)P_{A}(J_{A}). The probability distributions can be generally different from each other, i.e., PA(x)PB(x)P_{A}(x)\neq P_{B}(x), and are also allowed to present no randomness, i.e., PA(JA)=δ(JJA)P_{A}(J_{A})=\delta(J-J_{A}). The parameter λA\lambda_{A} plays a role in controlling the strength of the randomness.

The partition function Z{JA}Z_{\{J_{A}\}} and correlation function σB{JA}\langle\sigma_{B}\rangle_{\{J_{A}\}} for a set of fixed interactions, {JA}\{J_{A}\}, are given by

Z{JA}\displaystyle Z_{\{J_{A}\}} =\displaystyle= \Trexp(βAVλAJAσA),\displaystyle\Tr\exp\left(\beta\sum_{A\subset{V}}\lambda_{A}J_{A}\sigma_{A}\right), (3)
σB{JA}\displaystyle\langle\sigma_{B}\rangle_{\{J_{A}\}} =\displaystyle= \TrσBexp(βAVλAJAσA)Z{JA}.\displaystyle\frac{\Tr\sigma_{B}\exp\left(\beta\sum_{A\subset{V}}\lambda_{A}J_{A}\sigma_{A}\right)}{Z_{\{J_{A}\}}}. (4)

The configurational average over the distribution of the interactions is written as

𝔼[g({JA})]=(AV𝑑JAPA(JA))g({JA}).\displaystyle\mathbb{E}\left[g(\{J_{A}\})\right]=\left(\prod_{A\subset{V}}\int_{-\infty}^{\infty}dJ_{A}P_{A}(J_{A})\right)g(\{J_{A}\}). (5)

For example, the quenched average of the correlation function is obtained as

𝔼[σB{JA}]\displaystyle\mathbb{E}\left[\langle\sigma_{B}\rangle_{\{J_{A}\}}\right] =\displaystyle= (AV𝑑JAPA(JA))\TrσBexp(βAVλAJAσA)Z{JA}.\displaystyle\left(\prod_{A\subset{V}}\int_{-\infty}^{\infty}dJ_{A}P_{A}(J_{A})\right)\frac{\Tr\sigma_{B}\exp\left(\beta\sum_{A\subset{V}}\lambda_{A}J_{A}\sigma_{A}\right)}{Z_{\{J_{A}\}}}. (6)

In addition, the quenched pressure PP is defined as

P\displaystyle P =\displaystyle= 𝔼[logZ{JA}].\displaystyle\mathbb{E}\left[\log Z_{\{J_{A}\}}\right]. (7)

Recent studies showed [1, 2], when the probability distribution of a random interaction JBJ_{B} is symmetric, PB(JB)=PB(JB)P_{B}(J_{B})=P_{B}(-J_{B}), the first derivative of the quenched pressure with respect to the strength of the randomness is always positive:

1βPλB\displaystyle\frac{1}{\beta}\frac{\partial P}{\partial\lambda_{B}} =\displaystyle= 𝔼[JBσB{JA}]0,\displaystyle\mathbb{E}\left[J_{B}\langle\sigma_{B}\rangle_{\{J_{A}\}}\right]\geq 0, (8)

which is regarded as the counterpart of the Griffiths first inequality in spin-glass models.

Next, we consider the counterpart of the Griffiths second inequality in spin-glass models. Previous studies [2, 3, 4, 5] have focused on the second derivative of the quenched pressure with respect to the strength of the randomness,

1β22PλBλC\displaystyle\frac{1}{\beta^{2}}\frac{\partial^{2}P}{\partial\lambda_{B}\partial\lambda_{C}} =\displaystyle= 𝔼[JBJC(σBσC{JA}σB{JA}σC{JA})].\displaystyle\mathbb{E}\left[J_{B}J_{C}(\langle\sigma_{B}\sigma_{C}\rangle_{\{J_{A}\}}-\langle\sigma_{B}\rangle_{\{J_{A}\}}\langle\sigma_{C}\rangle_{\{J_{A}\}})\right]. (9)

In ferromagnetic models, the Griffiths second inequality means that the second derivative of the pressure with respect to the ferromagnetic interactions is always positive. Fortunately, on the Nishimori line, a similar relation also holds in spin-glass models and Eq.(9) is always positive [3, 4], which is the ferromagnetic counterpart of the Griffiths second inequality on the Nishimori line in spin-glass models.

In the case of symmetric distributions, however, studies on the counterpart of the Griffiths second inequality have not satisfactorily progressed. When PB(JB)P_{B}(J_{B}) and PC(JC)P_{C}(J_{C}) follow the symmetric Gaussian distributions with the variance ΛB2\Lambda_{B}^{2} and ΛC2\Lambda_{C}^{2}, respectively, by integration by parts, Eq. (9) is deformed to

1β22PλBλC\displaystyle\frac{1}{\beta^{2}}\frac{\partial^{2}P}{\partial\lambda_{B}\partial\lambda_{C}} =\displaystyle= ΛB2βλC𝔼[σB{JA}2].\displaystyle\frac{\Lambda_{B}^{2}}{\beta}\frac{\partial}{\partial\lambda_{C}}\mathbb{E}\left[-\langle\sigma_{B}\rangle_{\{J_{A}\}}^{2}\right]. (10)

This means that, if Eq. (9) is negative for σBσC\sigma_{B}\neq\sigma_{C}, the overlap expectation 𝔼[σB{JA}2]\mathbb{E}\left[\langle\sigma_{B}\rangle_{\{J_{A}\}}^{2}\right] is monotonic non decreasing with the system size [2]. The overlap expectation 𝔼[σB{JA}2]\mathbb{E}\left[\langle\sigma_{B}\rangle_{\{J_{A}\}}^{2}\right] tends to increase with increasing randomness. Then, as the counterpart of the Griffiths second inequality, it is expected that Eq. (9) is always negative for general symmetric distributions; however, an explicit counterexample exists [5] and it was shown that Eq. (9) takes both positive and negative values depending on the details of the model. Thus, the counterpart of the Griffiths second inequality has not been established even for symmetric distributions in spin-glass models and it is an important problem to investigate when Eq. (9) is negative.

On the other hand, it was shown that, for PB(JB)=PB(JB)P_{B}(J_{B})=P_{B}(-J_{B}), the first derivative of the quenched pressure with respect to the strength of the randomness has the following upper bound [9],

1βPλB=𝔼[JBσB{JA}]𝔼[JBtanh(βλBJB)].\displaystyle\frac{1}{\beta}\frac{\partial P}{\partial\lambda_{B}}=\mathbb{E}\left[J_{B}\langle\sigma_{B}\rangle_{\{J_{A}\}}\right]\leq\mathbb{E}\left[J_{B}\tanh(\beta\lambda_{B}J_{B})\right]. (11)

We note that Eq. (11) is independent of any other interaction, which is a non-trivial result. The proof of Eq. (11) was obtained by focusing only on the distribution of a single interaction among all of the interactions; however, it is not clear how to extend it to the case with multi variables. In the following sections, we give another proof of Eq. (11) and extend it to the case with multi variables. Then, we obtain some correlation inequalities in the Ising models with quenched randomness. Furthermore, in Sec. IV, we show that Eq. (9) has a non-trivial positive upper bound.

3 Inequalities for expectations by inequality of arithmetic and geometric means

In this section, we prove three inequalities for expectations, which play an important role in the next section.

We consider the case that the distribution functions of JBJ_{B} and JCJ_{C} satisfy the following relations:

PB(JB)\displaystyle P_{B}(-J_{B}) =\displaystyle= exp(2βNL,BJB)PB(JB),\displaystyle\exp(-2\beta_{\text{NL},B}J_{B})P_{B}(J_{B}), (12)
PC(JC)\displaystyle P_{C}(-J_{C}) =\displaystyle= exp(2βNL,CJC)PC(JC),\displaystyle\exp(-2\beta_{\text{NL},C}J_{C})P_{C}(J_{C}), (13)

where βNL,B\beta_{\text{NL},B} and βNL,C\beta_{\text{NL},C} are allowed to be any real values. For example, in the case of the Gaussian distribution

PB(JB)\displaystyle P_{B}(J_{B}) =\displaystyle= 12πμ2exp((JBJ0)22μ2),\displaystyle\frac{1}{\sqrt{2\pi\mu^{2}}}\exp\left(-\frac{(J_{B}-J_{0})^{2}}{2\mu^{2}}\right), (14)

and the binary distribution

PB(JB)\displaystyle P_{B}(J_{B}) =\displaystyle= pδ(JBJ)+(1p)δ(JB+J),\displaystyle p\delta(J_{B}-J)+(1-p)\delta(J_{B}+J), (15)

βNL,B\beta_{\text{NL},B} is given as follows, respectively,

βNL,B\displaystyle\beta_{\text{NL},B} =\displaystyle= J0μ2,\displaystyle\frac{J_{0}}{\mu^{2}}, (16)
βNL,B\displaystyle\beta_{\text{NL},B} =\displaystyle= p1p.\displaystyle\frac{p}{1-p}. (17)

We note that we do not impose any constraint on all the other interactions than JBJ_{B} and JCJ_{C}. In addition, for f({JA})f(\{J_{A}\}), when we focus on the interaction JBJ_{B}, we denote f({JA})f(\{J_{A}\}) as f(JB,{JA}JB)f(J_{B},\{J_{A}\}\setminus J_{B}). Similarly, when we are interested in JBJ_{B} and JCJ_{C}, we represent f({JA})f(\{J_{A}\}) as f(JB,JC,{JA}{JB,JC})f(J_{B},J_{C},\{J_{A}\}\setminus\{J_{B},J_{C}\})

Our first result in this section is as follows.

Lemma 3.1.

We assume that X({JA})X(\{J_{A}\}) satisfies

X({JA})>0,\displaystyle X(\{J_{A}\})>0, (18)
X(JB,{JA}JB)=1X(JB,{JA}JB).\displaystyle X(-J_{B},\{J_{A}\}\setminus J_{B})=\frac{1}{X(J_{B},\{J_{A}\}\setminus J_{B})}. (19)

Then, for any function f({JA})f(\{J_{A}\}) satisfying

f({JA})0,\displaystyle f(\{J_{A}\})\geq 0, (20)
f(JB,{JA}JB)=f(JB,{JA}JB),\displaystyle f(J_{B},\{J_{A}\}\setminus J_{B})=f(-J_{B},\{J_{A}\}\setminus J_{B}), (21)

the following inequality holds

𝔼[f({JA})X({JA})]𝔼[f({JA})exp(βNL,BJB)].\displaystyle\mathbb{E}\left[f(\{J_{A}\})X(\{J_{A}\})\right]\geq\mathbb{E}\left[f(\{J_{A}\})\exp(-\beta_{\text{NL},B}J_{B})\right]. (22)
Proof.

By dividing the integration interval of JBJ_{B} and summing up them, we obtain

𝔼[f({JA})X({JA})]\displaystyle\mathbb{E}\left[f(\{J_{A}\})X(\{J_{A}\})\right] (23)
=\displaystyle= 0𝑑JBPB(JB)𝔼[f({JA})(X(JB,{JA}JB)+X(JB,{JA}JB)e2βNL,BJB)]\displaystyle\int_{0}^{\infty}dJ_{B}P_{B}(J_{B})\mathbb{E}\left[f(\{J_{A}\})\left(X(J_{B},\{J_{A}\}\setminus J_{B})+X(-J_{B},\{J_{A}\}\setminus J_{B})e^{-2\beta_{\text{NL},B}J_{B}}\right)\right]^{\prime}
\displaystyle\geq 0𝑑JBPB(JB)𝔼[2f({JA})X(JB,{JA}JB)X(JB,{JA}JB)e2βNL,BJB]\displaystyle\int_{0}^{\infty}dJ_{B}P_{B}(J_{B})\mathbb{E}\left[2f(\{J_{A}\})\sqrt{X(J_{B},\{J_{A}\}\setminus J_{B})X(-J_{B},\{J_{A}\}\setminus J_{B})e^{-2\beta_{\text{NL},B}J_{B}}}\right]^{\prime}
=\displaystyle= 0𝑑JBPB(JB)𝔼[2f({JA})eβNL,BJB]\displaystyle\int_{0}^{\infty}dJ_{B}P_{B}(J_{B})\mathbb{E}\left[2f(\{J_{A}\})e^{-\beta_{\text{NL},B}J_{B}}\right]^{\prime}
=\displaystyle= 𝔼[f({JA})eβNL,BJB],\displaystyle\mathbb{E}\left[f(\{J_{A}\})e^{-\beta_{\text{NL},B}J_{B}}\right],

where 𝔼[]\mathbb{E}[\cdots]^{\prime} denotes the configurational average over the randomness of the interactions other than JBJ_{B}, and we used the inequality of arithmetic and geometric means. Thus, we prove Lemma 3.1. \Box

Using Lemma 3.1, we give another derivation of the inequality for the local energy[9, 10, 11],

𝔼[f({JA})JBσB{JA}]𝔼[f({JA})JB(tanh(βλBJB)+1eβNL,BJBsinh(2βλBJB))],\displaystyle\mathbb{E}\left[-f(\{J_{A}\})J_{B}\langle\sigma_{B}\rangle_{\{J_{A}\}}\right]\geq\mathbb{E}\left[-f(\{J_{A}\})J_{B}\left(\tanh(\beta\lambda_{B}J_{B})+\frac{1-e^{-\beta_{\text{NL},B}J_{B}}}{\sinh(2\beta\lambda_{B}J_{B})}\right)\right], (24)

where f({JA})f(\{J_{A}\}) satisfies

f({JA})\displaystyle f(\{J_{A}\}) \displaystyle\geq 0,\displaystyle 0, (25)
f(JB,{JA}JB)\displaystyle f(J_{B},\{J_{A}\}\setminus J_{B}) =\displaystyle= f(JB,{JA}JB).\displaystyle f(-J_{B},\{J_{A}\}\setminus J_{B}). (26)

We note that, for βNL,B=0\beta_{\text{NL},B}=0 and f({JA})=1f(\{J_{A}\})=1, this inequality coincides with Eq. (11).

Proof.

For simplicity, we denote Z(JB)Z(J_{B}) as

Z(JB)\displaystyle Z(J_{B}) =\displaystyle= \Trexp(βAVBλAJAσA+βλBJBσB),\displaystyle\Tr\exp\left(\beta\sum_{A\subset{V}\setminus B}\lambda_{A}J_{A}\sigma_{A}+\beta\lambda_{B}J_{B}\sigma_{B}\right), (27)

We note that Z(JB)=Z{JA}Z(J_{B})=Z_{\{J_{A}\}} but Z(JB)Z{JA}Z(-J_{B})\neq Z_{\{J_{A}\}}. Then, putting X({JA})X(\{J_{A}\}) as Z(JB)/Z(JB)=cosh(2βλBJB)σB{JA}sinh(2βλBJB)Z(-J_{B})/Z(J_{B})=\cosh(2\beta\lambda_{B}J_{B})-\langle\sigma_{B}\rangle_{\{J_{A}\}}\sinh(2\beta\lambda_{B}J_{B}) in Lemma. 3.1, we obtain

𝔼[f({JA})(cosh(2βλBJB)σB{JA}sinh(2βλBJB))]𝔼[f({JA})eβNL,BJB].\displaystyle\mathbb{E}\left[f(\{J_{A}\})\left(\cosh(2\beta\lambda_{B}J_{B})-\langle\sigma_{B}\rangle_{\{J_{A}\}}\sinh(2\beta\lambda_{B}J_{B})\right)\right]\geq\mathbb{E}\left[f(\{J_{A}\})e^{-\beta_{\text{NL},B}J_{B}}\right]. (28)

Furthermore, replacing f({JA})f(\{J_{A}\}) by f({JA})JB/sinh(2βJB)f(\{J_{A}\}){J_{B}}/{\sinh(2\beta J_{B})}, we prove Eq. (24). \Box

Our second result is a extension of Lemma. 3.1 to two-variable case.

Lemma 3.2.

We assume that X({JA})X(\{J_{A}\}) satisfies

X({JA})>0,\displaystyle X(\{J_{A}\})>0, (29)
X(JB,JC,{JA}{JB,JC})=1X(JB,JC,{JA}{JB,JC}).\displaystyle X(-J_{B},-J_{C},\{J_{A}\}\setminus\{J_{B},J_{C}\})=\frac{1}{X(J_{B},J_{C},\{J_{A}\}\setminus\{J_{B},J_{C}\})}. (30)

Then, for any function f({JA})f(\{J_{A}\}) satisfying

f({JA})\displaystyle f(\{J_{A}\}) \displaystyle\geq 0,\displaystyle 0, (31)
f(JB,JC,{JA}{JB,JC}))\displaystyle f(J_{B},J_{C},\{J_{A}\}\setminus\{J_{B},J_{C}\})) =\displaystyle= f(JB,JC,{JA}{JB,JC}))\displaystyle f(-J_{B},J_{C},\{J_{A}\}\setminus\{J_{B},J_{C}\})) (32)
=\displaystyle= f(JB,JC,{JA}{JB,JC}))\displaystyle f(J_{B},-J_{C},\{J_{A}\}\setminus\{J_{B},J_{C}\}))
=\displaystyle= f(JB,JC,{JA}{JB,JC})),\displaystyle f(-J_{B},-J_{C},\{J_{A}\}\setminus\{J_{B},J_{C}\})),

the following inequality holds

𝔼[f({JA})X({JA})]𝔼[f({JA})exp(βNL,BJB)exp(βNL,CJC)].\displaystyle\mathbb{E}\left[f(\{J_{A}\})X(\{J_{A}\})\right]\geq\mathbb{E}\left[f(\{J_{A}\})\exp\left(-\beta_{\text{NL},B}J_{B}\right)\exp\left(-\beta_{\text{NL},C}J_{C}\right)\right]. (33)
Proof.

For simplicity, we denote X(JB,JC,{JA}{JB,JC})X(J_{B},J_{C},\{J_{A}\}\setminus\{J_{B},J_{C}\}) as X(JB,JC)X(J_{B},J_{C}). By dividing the integration interval of JBJ_{B} and JCJ_{C} and summing up them, we find

𝔼[f({JA})X({JA})]\displaystyle\mathbb{E}\left[f(\{J_{A}\})X(\{J_{A}\})\right] (34)
=\displaystyle= 0𝑑JBPB(JB)0𝑑JCPC(JC)\displaystyle\int_{0}^{\infty}dJ_{B}P_{B}(J_{B})\int_{0}^{\infty}dJ_{C}P_{C}(J_{C})
𝔼[f({JA})(X(JB,JC)+X(JB,JC)e2βNL,BJBe2βNL,CJC\displaystyle\mathbb{E}\left[f(\{J_{A}\})\left(X(J_{B},J_{C})+X(-J_{B},-J_{C})e^{-2\beta_{\text{NL},B}J_{B}}e^{-2\beta_{\text{NL},C}J_{C}}\right.\right.
+X(JB,JC)e2βNL,BJB+X(JB,JC)e2βNL,CJC)]′′\displaystyle\left.\left.+X(-J_{B},J_{C})e^{-2\beta_{\text{NL},B}J_{B}}+X(J_{B},-J_{C})e^{-2\beta_{\text{NL},C}J_{C}}\right)\right]^{\prime\prime}
\displaystyle\geq 0𝑑JBPB(JB)0𝑑JCPC(JC)𝔼[4f({JA})eβNL,BJBeβNL,CJC]′′\displaystyle\int_{0}^{\infty}dJ_{B}P_{B}(J_{B})\int_{0}^{\infty}dJ_{C}P_{C}(J_{C})\mathbb{E}\left[4f(\{J_{A}\})e^{-\beta_{\text{NL},B}J_{B}}e^{-\beta_{\text{NL},C}J_{C}}\right]^{\prime\prime}
=\displaystyle= 𝔼[f({JA})eβNL,BJBeβNL,CJC],\displaystyle\mathbb{E}\left[f(\{J_{A}\})e^{-\beta_{\text{NL},B}J_{B}}e^{-\beta_{\text{NL},C}J_{C}}\right],

where 𝔼[]′′\mathbb{E}[\cdots]^{\prime\prime} denotes the configurational average over the randomness of the interactions other than JBJ_{B} and JCJ_{C}, and we used the inequality of arithmetic and geometric means. Therefore, we prove Lemma 3.2.

\Box

Our third result is another extension of Lemma. 3.1 to two-variable case.

Lemma 3.3.

We assume that X({JA})X(\{J_{A}\}) and Y({JA})Y(\{J_{A}\}) satisfy

X({JA})>0,\displaystyle X(\{J_{A}\})>0, (35)
Y({JA})>0,\displaystyle Y(\{J_{A}\})>0, (36)
X(JB,JC,{JA}{JB,JC})\displaystyle X(-J_{B},J_{C},\{J_{A}\}\setminus\{J_{B},J_{C}\}) =\displaystyle= 1X(JB,JC,{JA}{JB,JC}),\displaystyle\frac{1}{X(J_{B},J_{C},\{J_{A}\}\setminus\{J_{B},J_{C}\})}, (37)
Y(JB,JC,{JA}{JB,JC})\displaystyle Y(J_{B},-J_{C},\{J_{A}\}\setminus\{J_{B},J_{C}\}) =\displaystyle= 1Y(JB,JC,{JA}{JB,JC}).\displaystyle\frac{1}{Y(J_{B},J_{C},\{J_{A}\}\setminus\{J_{B},J_{C}\})}. (38)

Then, for any function f({JA})f(\{J_{A}\}) satisfying

f({JA})\displaystyle f(\{J_{A}\}) \displaystyle\geq 0,\displaystyle 0, (39)
f(JB,JC,{JA}{JB,JC}))\displaystyle f(J_{B},J_{C},\{J_{A}\}\setminus\{J_{B},J_{C}\})) =\displaystyle= f(JB,JC,{JA}{JB,JC}))\displaystyle f(-J_{B},J_{C},\{J_{A}\}\setminus\{J_{B},J_{C}\})) (40)
=\displaystyle= f(JB,JC,{JA}{JB,JC}))\displaystyle f(J_{B},-J_{C},\{J_{A}\}\setminus\{J_{B},J_{C}\}))
=\displaystyle= f(JB,JC,{JA}{JB,JC})),\displaystyle f(-J_{B},-J_{C},\{J_{A}\}\setminus\{J_{B},J_{C}\})),

the following inequality holds

𝔼[f({JA})X({JA})Y({JA})]𝔼[f({JA})exp(βNL,BJB)exp(βNL,CJC)].\displaystyle\mathbb{E}\left[f(\{J_{A}\})X(\{J_{A}\})Y(\{J_{A}\})\right]\geq\mathbb{E}\left[f(\{J_{A}\})\exp\left(-\beta_{\text{NL},B}J_{B}\right)\exp\left(-\beta_{\text{NL},C}J_{C}\right)\right]. (41)
Proof.

For simplicity, we denote X(JB,JC,{JA}{JB,JC})X(J_{B},J_{C},\{J_{A}\}\setminus\{J_{B},J_{C}\}) and Y(JB,JC,{JA}{JB,JC})Y(J_{B},J_{C},\{J_{A}\}\setminus\{J_{B},J_{C}\}) as X(JB,JC)X(J_{B},J_{C}) and Y(JB,JC)Y(J_{B},J_{C}). Similarly to the derivation of Lemma 3.2, we find

𝔼[f({JA})X({JA})Y({JA})]\displaystyle\mathbb{E}\left[f(\{J_{A}\})X(\{J_{A}\})Y(\{J_{A}\})\right] (42)
=\displaystyle= 0𝑑JBPB(JB)𝔼[f({JA})(X(JB,JC)Y(JB,JC)+X(JB,JC)Y(JB,JC)e2βNL,BJB)]\displaystyle\int_{0}^{\infty}dJ_{B}P_{B}(J_{B})\mathbb{E}\left[f(\{J_{A}\})\left(X(J_{B},J_{C})Y(J_{B},J_{C})+X(-J_{B},J_{C})Y(-J_{B},J_{C})e^{-2\beta_{\text{NL},B}J_{B}}\right)\right]^{\prime}
\displaystyle\geq 0𝑑JBPB(JB)𝔼[2f({JA})X(JB,JC)X(JB,JC)Y(JB,JC)Y(JB,JC)e2βNL,BJB]\displaystyle\int_{0}^{\infty}dJ_{B}P_{B}(J_{B})\mathbb{E}\left[2f(\{J_{A}\})\sqrt{X(J_{B},J_{C})X(-J_{B},J_{C})Y(J_{B},J_{C})Y(-J_{B},J_{C})e^{-2\beta_{\text{NL},B}J_{B}}}\right]^{\prime}
=\displaystyle= 0𝑑JBPB(JB)𝔼[2f({JA})eβNL,BJBY(JB,JC)Y(JB,JC)]\displaystyle\int_{0}^{\infty}dJ_{B}P_{B}(J_{B})\mathbb{E}\left[2f(\{J_{A}\})e^{-\beta_{\text{NL},B}J_{B}}\sqrt{Y(J_{B},J_{C})Y(-J_{B},J_{C})}\right]^{\prime}
=\displaystyle= 0dJBPB(JB)0dJCPC(JC)𝔼[2f({JA})eβNL,BJBY(JB,JC)Y(JB,JC)\displaystyle\int_{0}^{\infty}dJ_{B}P_{B}(J_{B})\int_{0}^{\infty}dJ_{C}P_{C}(J_{C})\mathbb{E}\left[2f(\{J_{A}\})e^{-\beta_{\text{NL},B}J_{B}}\sqrt{Y(J_{B},J_{C})Y(-J_{B},J_{C})}\right.
+2f({JA})eβNL,BJBe2βNL,CJCY(JB,JC)Y(JB,JC)]′′\displaystyle\left.+2f(\{J_{A}\})e^{-\beta_{\text{NL},B}J_{B}}e^{-2\beta_{\text{NL},C}J_{C}}\sqrt{Y(J_{B},-J_{C})Y(-J_{B},-J_{C})}\right]^{\prime\prime}
\displaystyle\geq 0𝑑JBPB(JB)0𝑑JCPC(JC)\displaystyle\int_{0}^{\infty}dJ_{B}P_{B}(J_{B})\int_{0}^{\infty}dJ_{C}P_{C}(J_{C})
𝔼[4f({JA})eβNL,BJBY(JB,JC)Y(JB,JC)e2βNL,CJCY(JB,JC)Y(JB,JC)]′′\displaystyle\mathbb{E}\left[4f(\{J_{A}\})e^{-\beta_{\text{NL},B}J_{B}}\sqrt{\sqrt{Y(J_{B},J_{C})Y(-J_{B},J_{C})}e^{-2\beta_{\text{NL},C}J_{C}}\sqrt{Y(J_{B},-J_{C})Y(-J_{B},-J_{C})}}\right]^{\prime\prime}
=\displaystyle= 0𝑑JBPB(JB)0𝑑JCPC(JC)𝔼[4f({JA})eβNL,BJBeβNL,CJC]′′\displaystyle\int_{0}^{\infty}dJ_{B}P_{B}(J_{B})\int_{0}^{\infty}dJ_{C}P_{C}(J_{C})\mathbb{E}\left[4f(\{J_{A}\})e^{-\beta_{\text{NL},B}J_{B}}e^{-\beta_{\text{NL},C}J_{C}}\right]^{\prime\prime}
=\displaystyle= 𝔼[f({JA})exp(βNL,BJB)exp(βNL,CJC)],\displaystyle\mathbb{E}\left[f(\{J_{A}\})\exp\left(-\beta_{\text{NL},B}J_{B}\right)\exp\left(-\beta_{\text{NL},C}J_{C}\right)\right],

where we used the inequality of arithmetic and geometric means twice. Thus, we arrive at Lemma 3.3. \Box

4 Some correlation inequalities in spin-glass models

We have proved three lemmas for expectations in Sec III. In this section, using the acquired inequalities, we obtain some correlation inequalities in spin-glass models which is an extension of Eq. (24) to the case with multi variables.

4.1 Inequalities for double interactions

Although the previous studies [9, 10] has focused only on a single interaction JBJ_{B} among all of the interactions, Lemma. 3.2 enable us to extend their results to double interactions JBJ_{B} and JCJ_{C} among all of the interactions.

Theorem 4.1.

Under the conditions (12) and (13), the following inequality holds:

𝔼[JBJC(σBσC{JA}+1exp(βNL,BJB)exp(βNL,CJC)sinh(2λBβJB)sinh(2λCβJC))]\displaystyle\mathbb{E}\left[J_{B}J_{C}\left(\langle\sigma_{B}\sigma_{C}\rangle_{\{J_{A}\}}+\frac{1-\exp\left(-\beta_{\text{NL},B}J_{B}\right)\exp\left(-\beta_{\text{NL},C}J_{C}\right)}{\sinh(2\lambda_{B}\beta J_{B})\sinh(2\lambda_{C}\beta J_{C})}\right)\right]
\displaystyle\geq 𝔼[JBJC(tanh(βλBJB)tanh(βλCJC)+σB{JA}tanh(βλBJB)tanh(2λCβJC)+σC{JA}tanh(βλCJC)tanh(2λBβJB))].\displaystyle\mathbb{E}\left[J_{B}J_{C}\left(\tanh(\beta\lambda_{B}J_{B})\tanh(\beta\lambda_{C}J_{C})+\frac{\langle\sigma_{B}\rangle_{\{J_{A}\}}-\tanh(\beta\lambda_{B}J_{B})}{\tanh(2\lambda_{C}\beta J_{C})}+\frac{\langle\sigma_{C}\rangle_{\{J_{A}\}}-\tanh(\beta\lambda_{C}J_{C})}{\tanh(2\lambda_{B}\beta J_{B})}\right)\right].
Proof.

For simplicity, we denote Z(JB,JC)Z(J_{B},J_{C}) as

Z(JB,JC)\displaystyle Z(J_{B},J_{C}) =\displaystyle= \Trexp(βAV{B,C}λAJAσA+βλBJBσB+βλCJCσC).\displaystyle\Tr\exp\left(\beta\sum_{A\subset{V}\setminus\{B,C\}}\lambda_{A}J_{A}\sigma_{A}+\beta\lambda_{B}J_{B}\sigma_{B}+\beta\lambda_{C}J_{C}\sigma_{C}\right). (44)

Setting X({JA})X(\{J_{A}\}) as Z(JB,JC)/Z(JB,JC)Z(-J_{B},-J_{C})/Z(J_{B},J_{C}) in Lemma. 3.2, we obtain

𝔼[f({JA})(cosh(2βJB)cosh(2βJC)+sinh(2βJB)sinh(2βJC)σBσC{JA}\displaystyle\mathbb{E}\left[f(\{J_{A}\})\left(\cosh(2\beta J_{B})\cosh(2\beta J_{C})+\sinh(2\beta J_{B})\sinh(2\beta J_{C})\langle\sigma_{B}\sigma_{C}\rangle_{\{J_{A}\}}\right.\right. (45)
sinh(2βJB)cosh(2βJC)σB{JA}cosh(2βJB)sinh(2βJC)σC{JA})]\displaystyle\left.\left.-\sinh(2\beta J_{B})\cosh(2\beta J_{C})\langle\sigma_{B}\rangle_{\{J_{A}\}}-\cosh(2\beta J_{B})\sinh(2\beta J_{C})\langle\sigma_{C}\rangle_{\{J_{A}\}}\right)\right]
\displaystyle\geq 𝔼[f({JA})eβNL,BJBeβNL,CJC].\displaystyle\mathbb{E}\left[f(\{J_{A}\})e^{-\beta_{\text{NL},B}J_{B}}e^{-\beta_{\text{NL},C}J_{C}}\right].

Then, substituting JBJC/(sinh(2λBβJB)sinh(2λCβJC))J_{B}J_{C}/\left(\sinh(2\lambda_{B}\beta J_{B})\sinh(2\lambda_{C}\beta J_{C})\right) into f({JA})f(\{J_{A}\}), we arrive at

𝔼[JBJC(σBσC{JA}+1eβNL,BJBeβNL,CJCsinh(2λBβJB)sinh(2λCβJC))]\displaystyle\mathbb{E}\left[J_{B}J_{C}\left(\langle\sigma_{B}\sigma_{C}\rangle_{\{J_{A}\}}+\frac{1-e^{-\beta_{\text{NL},B}J_{B}}e^{-\beta_{\text{NL},C}J_{C}}}{\sinh(2\lambda_{B}\beta J_{B})\sinh(2\lambda_{C}\beta J_{C})}\right)\right]
\displaystyle\geq 𝔼[JBJC(1cosh(2βJB)cosh(2βJC)sinh(2βJB)sinh(2βJC)+σB{JA}tanh(2λCβJC)+σC{JA}tanh(2λBβJB))]\displaystyle\mathbb{E}\left[J_{B}J_{C}\left(\frac{1-\cosh(2\beta J_{B})\cosh(2\beta J_{C})}{\sinh(2\beta J_{B})\sinh(2\beta J_{C})}+\frac{\langle\sigma_{B}\rangle_{\{J_{A}\}}}{\tanh(2\lambda_{C}\beta J_{C})}+\frac{\langle\sigma_{C}\rangle_{\{J_{A}\}}}{\tanh(2\lambda_{B}\beta J_{B})}\right)\right]
=\displaystyle= 𝔼[JBJC(tanh(βλBJB)tanh(βλCJC)+σB{JA}tanh(βλBJB)tanh(2λCβJC)+σC{JA}tanh(βλCJC)tanh(2λBβJB))].\displaystyle\mathbb{E}\left[J_{B}J_{C}\left(\tanh(\beta\lambda_{B}J_{B})\tanh(\beta\lambda_{C}J_{C})+\frac{\langle\sigma_{B}\rangle_{\{J_{A}\}}-\tanh(\beta\lambda_{B}J_{B})}{\tanh(2\lambda_{C}\beta J_{C})}+\frac{\langle\sigma_{C}\rangle_{\{J_{A}\}}-\tanh(\beta\lambda_{C}J_{C})}{\tanh(2\lambda_{B}\beta J_{B})}\right)\right].

Thus, we obtain Eq. (4.1). \Box

We mention that we had found the following inequality for βNL,B=βNL,C=0\beta_{\text{NL},B}=\beta_{\text{NL},C}=0 in Ref. [11],

𝔼[JBJCσBσC{JA}]0.\displaystyle\mathbb{E}\left[J_{B}J_{C}\langle\sigma_{B}\sigma_{C}\rangle_{\{J_{A}\}}\right]\geq 0. (47)

Numerical calculation suggests that the right-hand side in Eq. (4.1) does not have a definite sign for βNL,B=βNL,C=0\beta_{\text{NL},B}=\beta_{\text{NL},C}=0. Thus, it is considered that Eq. (4.1) is independent of Eq. (47).

Interestingly, an inequality of the same form as in Eq. (4.1) holds for 𝔼[JBJCσB{JA}σC{JA}]\mathbb{E}\left[J_{B}J_{C}\langle\sigma_{B}\rangle_{\{J_{A}\}}\langle\sigma_{C}\rangle_{\{J_{A}\}}\right] as well.

Theorem 4.2.

Under the conditions (12) and (13), the following inequality holds:

𝔼[JBJC(σB{JA}σC{JA}+1exp(βNL,BJB)exp(βNL,CJC)sinh(2λBβJB)sinh(2λCβJC))]\displaystyle\mathbb{E}\left[J_{B}J_{C}\left(\langle\sigma_{B}\rangle_{\{J_{A}\}}\langle\sigma_{C}\rangle_{\{J_{A}\}}+\frac{1-\exp\left(-\beta_{\text{NL},B}J_{B}\right)\exp\left(-\beta_{\text{NL},C}J_{C}\right)}{\sinh(2\lambda_{B}\beta J_{B})\sinh(2\lambda_{C}\beta J_{C})}\right)\right]
\displaystyle\geq 𝔼[JBJC(tanh(βλBJB)tanh(βλCJC)+σB{JA}tanh(βλBJB)tanh(2λCβJC)+σC{JA}tanh(βλCJC)tanh(2λBβJB))].\displaystyle\mathbb{E}\left[J_{B}J_{C}\left(\tanh(\beta\lambda_{B}J_{B})\tanh(\beta\lambda_{C}J_{C})+\frac{\langle\sigma_{B}\rangle_{\{J_{A}\}}-\tanh(\beta\lambda_{B}J_{B})}{\tanh(2\lambda_{C}\beta J_{C})}+\frac{\langle\sigma_{C}\rangle_{\{J_{A}\}}-\tanh(\beta\lambda_{C}J_{C})}{\tanh(2\lambda_{B}\beta J_{B})}\right)\right].
Proof.

The proof is very similar to Theorem 4.1. We denote Z(JB)Z(J_{B}) and Z(JC)Z(J_{C}) as

Z(JB)\displaystyle Z(J_{B}) =\displaystyle= \Trexp(βAVBλAJAσA+βλBJBσB),\displaystyle\Tr\exp\left(\beta\sum_{A\subset{V}\setminus B}\lambda_{A}J_{A}\sigma_{A}+\beta\lambda_{B}J_{B}\sigma_{B}\right), (49)
Z(JC)\displaystyle Z(J_{C}) =\displaystyle= \Trexp(βAVCλAJAσA+βλCJCσC).\displaystyle\Tr\exp\left(\beta\sum_{A\subset{V}\setminus C}\lambda_{A}J_{A}\sigma_{A}+\beta\lambda_{C}J_{C}\sigma_{C}\right). (50)

Putting X({JA})X(\{J_{A}\}), Y({JA})Y(\{J_{A}\}) and f({JA})f(\{J_{A}\}) as Z(JB)/Z(JB)Z(-J_{B})/Z(J_{B}), Z(JC)/Z(JC)Z(-J_{C})/Z(J_{C}) and JBJC/(sinh(2λBβJB)sinh(2λCβJC))J_{B}J_{C}/\left(\sinh(2\lambda_{B}\beta J_{B})\sinh(2\lambda_{C}\beta J_{C})\right) in Lemma. 3.3, we obtain Eq. (4.2). \Box

4.2 Upper bound on second derivative of quenched pressure

We have extended the previous studies [9, 10] to two-variable cases. However, at first glance, Eqs. (4.1) and (4.2) take complex forms and it is not clear what these inequalities mean. Here, we show that Eq. (4.2) enable us to obtain a simple upper bound on the second derivative of the quenched pressure with respect to the strength of the randomness.

For βNL,B=βNL,C=0\beta_{\text{NL},B}=\beta_{\text{NL},C}=0, we have already shown the following inequality [11],

𝔼[JBJCtanh(βλBJB)tanh(βλCJC)]\displaystyle\mathbb{E}\left[J_{B}J_{C}\tanh(\beta\lambda_{B}J_{B})\tanh(\beta\lambda_{C}J_{C})\right] \displaystyle\geq 𝔼[JBJCσBσC{JA}].\displaystyle\mathbb{E}\left[J_{B}J_{C}\langle\sigma_{B}\sigma_{C}\rangle_{\{J_{A}\}}\right]. (51)

Then, combining Eqs. (4.2) and (51), we arrive at the following result.

Corollary 4.3.

For symmetric distribution βNL,B=βNL,C=0\beta_{\text{NL},B}=\beta_{\text{NL},C}=0 and σBσC\sigma_{B}\neq\sigma_{C}, the second derivative of the quenched pressure with respect to the strength of the randomness has a non-trivial upper bound,

1β22λBλC𝔼[logZ]\displaystyle\frac{1}{\beta^{2}}\frac{\partial^{2}}{\partial\lambda_{B}\partial\lambda_{C}}\mathbb{E}\left[\log Z\right] =\displaystyle= 𝔼[JBJC(σBσC{JA}σB{JA}σC{JA})]\displaystyle\mathbb{E}\left[J_{B}J_{C}(\langle\sigma_{B}\sigma_{C}\rangle_{\{J_{A}\}}-\langle\sigma_{B}\rangle_{\{J_{A}\}}\langle\sigma_{C}\rangle_{\{J_{A}\}})\right]
\displaystyle\leq 𝔼[JBJC(tanh(βλBJB)σB{JA}tanh(2βλCJC)+tanh(βλCJC)σC{JA}tanh(2βλBJB))].\displaystyle\mathbb{E}\left[J_{B}J_{C}\left(\frac{\tanh(\beta\lambda_{B}J_{B})-\langle\sigma_{B}\rangle_{\{J_{A}\}}}{\tanh(2\beta\lambda_{C}J_{C})}+\frac{\tanh(\beta\lambda_{C}J_{C})-\langle\sigma_{C}\rangle_{\{J_{A}\}}}{\tanh(2\beta\lambda_{B}J_{B})}\right)\right].
Remark 4.4.

We note that, from Eq. (24), the right-hand side in Eq. (4.3) is always positive. Thus, this bound can be considered as a weak result of the counterpart of the Griffiths second inequality for general symmetric distributions.

In order to investigate the tightness of Eq. (4.3), for example, we consider a closed chain of six spins with one added interaction between 22 and 55,

H\displaystyle H =\displaystyle= λ1,2J1,2σ1σ2λ2,3J2,3σ2σ3λ3,4J3,4σ3σ4λ4,5J4,5σ4σ5\displaystyle-\lambda_{1,2}J_{1,2}\sigma_{1}\sigma_{2}-\lambda_{2,3}J_{2,3}\sigma_{2}\sigma_{3}-\lambda_{3,4}J_{3,4}\sigma_{3}\sigma_{4}-\lambda_{4,5}J_{4,5}\sigma_{4}\sigma_{5} (53)
λ5,6J5,6σ5σ6λ6,1J6,1σ6σ1λ2,5J2,5σ2σ5,\displaystyle-\lambda_{5,6}J_{5,6}\sigma_{5}\sigma_{6}-\lambda_{6,1}J_{6,1}\sigma_{6}\sigma_{1}-\lambda_{2,5}J_{2,5}\sigma_{2}\sigma_{5},

where all the interactions follows independently a symmetric binary distribution

P(J)=12δ(J1)+12δ(J+1),\displaystyle P(J)=\frac{1}{2}\delta(J-1)+\frac{1}{2}\delta(J+1), (54)

and λ1,2=λ2,3=λ3,4=λ4,5=λ5,6=λ2,5=λ\lambda_{1,2}=\lambda_{2,3}=\lambda_{3,4}=\lambda_{4,5}=\lambda_{5,6}=\lambda_{2,5}=\lambda. This model had been investigated in the previous study [5] and, for B={1,2}B=\{1,2\}, C={2,3}C=\{2,3\} and β=1\beta=1, it was shown that the second derivative of the quenched pressure takes a negative value for λ<0.695\lambda<0.695 and a positive value for λ>0.695\lambda>0.695. For this model (53), we numerically calculate both sides in Eq. (4.3) for B={1,2}B=\{1,2\}, C={2,3}C=\{2,3\} and β=1\beta=1 from λ=0.001\lambda=0.001 to λ=3\lambda=3 in Fig. 1. The numerical calculation shows that the acquired inequality is strict in the regions with small randomness (λ<<1\lambda<<1) but give weak evaluation in the regions with strong randomness.

Refer to caption
Figure 1: Randomness dependence of both sides in Eq. (4.3) for B={1,2}B=\{1,2\}, C={2,3}C=\{2,3\} and β=1\beta=1 for the closed chain of six spins with one added interaction (53). The horizontal axis denotes the strength of the randomness, λ\lambda. The dotted line and the thick line show the left-hand side and the right-hand side in Eq. (4.3), respectively.

5 Conclusions

We have obtained several correlation inequalities for the Ising models with quenched randomness. Our main inequalities (4.1) and (4.2) are extensions of previous studies [9, 10] (24) to two-variable case. Our method can be easily extended to the case with more than three variables. Then, it is possible to obtain an infinite number of correlation inequalities in principle and the problem is how to find a meaningful one.

Furthermore, using the obtained inequalities for general symmetric distributions, we have given the positive upper bound (4.3) on the second derivative of the quenched pressure with respect to the strength of the randomness (9). Numerical calculation shows that our bound (4.3) is strict in the regions with small randomness. Equation (9) does not always take a negative value [5] and, thus, the counterpart of the Griffiths second inequality has not been established in spin glass models. Our bound (4.3) is slight progress on this issue and can be regarded as a weak result of the counterpart of the Griffiths second inequality in spin glass models for general symmetric distributions, which is the first non-trivial upper bound on Eq. (9).

It is an important problem to improve our bound (4.3) to a tighter one. One direction for future research is to consider model-dependent properties such as the shape of the lattice and the interaction. Besides, it may also be useful to consider the effects of other interactions. We have only focused on two interactions JBJ_{B} and JCJ_{C}, and have not imposed any constraint on all the other interactions. Incorporating information other than the two interactions may make our bound (4.3) tighter.

{acknowledgment}

The present work was financially supported by JSPS KAKENHI Grant No. 18H03303, 19H01095, 19K23418, and the JST-CREST (No.JPMJCR1402) for Japan Science and Technology Agency.

References

  • [1] P. Contucci and S. Graffi, Monotonicity and Thermodynamic Limit for Short Range Disordered Models, J. Stat. Phys. 115, 581 (2004).
  • [2] P. Contucci and J. Lebowitz, Correlation Inequalities for Spin Glasses, Ann. Henri Poincare 8, 1461(2007).
  • [3] S. Morita, H. Nishimori, and P. Contucci, Griffiths inequalities for the Gaussian spin glass, J. Phys. A 37, L203 (2004).
  • [4] H. Kitatani, Griffiths Inequalities for Ising Spin Glasses on the Nishimori Line, J. Phys. Soc. Jpn. 78, 044714 (2009).
  • [5] P. Contucci, F. Unguendoli, and C. Vernia, Lack of monotonicity in spin glass correlation functions, J. Phys. A: Math. Theor. 41, 385001 (2008).
  • [6] P. Contucci and S. Starr, Thermodynamic Limit for Spin Glasses. Beyond the Annealed Bound, J. Stat. Phys. 135, 1159 (2006).
  • [7] P. Contucci and C. Giardina, Perspective on Spin Glasses, (Cambridge University, Cambridge, 2013).
  • [8] H. Nishimori, Internal Energy, Specific Heat and Correlation Function of the Bond-Random Ising Model, Prog. Theor. Phys. 66, 1169 (1981).
  • [9] H. Kitatani, H. Nishimori, and A. Aoki, Inequalities for the Local Energy of Random Ising Models, J. Phys. Soc. Jpn. 76, 074711 (2007).
  • [10] M. Okuyama and M. Ohzeki, Inequality for local energy of Ising models with quenched randomness and its application, J. Phys. Soc. Jpn. 89, 064704 (2020).
  • [11] M. Okuyama and M. Ohzeki, Some inequalities for correlation functions of Ising models with quenched randomness, arXiv:2004.05832 (2020).