## Brazilian Journal of Physics

##
*Print version* ISSN 0103-9733

*On-line version* ISSN 1678-4448

### Braz. J. Phys. vol.33 no.4 São Paulo Dec. 2003

#### http://dx.doi.org/10.1590/S0103-97332003000400046

**REGULAR ARTICLES**

**On the entropy of the Viana-Bray model **

**J.R.L. de Almeida **

Departamento de Física, Universidade Federal de Pernambuco, 50670-901, Recife-PE, Brazil

**ABSTRACT**

The entropy of the Viana-Bray model at zero temperature and external field is calculated within the solution which takes into account only delta functions for the global order parameter *P*(*h*). It is shown that such solution is unsatisfactory both from the viewpoint of stability analysis and for not reproducing the well known Sherrington-Kirkpatrick result in the large connectivity limit thus pointing out the relevance of considering solutions with continuous part in *P*(*h*) for such model and possibly related models.

**1 Introduction**

Recently a lot of work has been devoted to the study of finite connectivity spin glasses and related models associated to the satisfiability problem (see [1] and references therein). A prototype of these finite connectivity models is the Viana-Bray (VB) model [2] designed to study diluted magnetic systems and applicable to optimization problems [3], specially the 2-SAT problem [4]. One of the most usual analytical technique used to treat these problems is the replica method which for the Sherrington-Kirkpatrick (SK) model [5] yields a negative entropy at low temperatures within the replica symmetric (RS) solution, signaling a problem for the unbroken symmetry ansatz. Eventually, stability study of the RS solution [6] pointed out the necessity to introduce broken replica simmetry solution which led to a physically acceptable solution [7]. Although the VB model has been intensely studied, the stability study of its solutions has been thoroughly analyzed mainly close to the percolation threshold both for the global order parameter with only discrete components [8] or including continuous component [9]. In this note the entropy of the VB model for the discrete component solution at zero temperature and external field is considered and shown to be rather unsatisfactory for it is plainly unstable against longitudinal variations and moreover does not yield SK's result in the large coordination limit. The VB model is described by the Hamiltonian

where s_{i} = ±1 (*i* = 1,2,....,*N*), and the *J _{ij}*'s are infinite-ranged random interactions with probability distribution given by

It thus may describe a highly diluted system with average connectivity *p*. The distribution of the active bonds shall here be taken as a bimodal: *f*(*J _{ij}*) = (d(

*J*- 1) + d(

_{ij}*J*+ 1)) / 2. The variational free energy

_{ij}*f*within the replica method takes the form [8,9]

where *g _{n}*(s

_{a}) is a generalized global order parameter, with 2

^{n}components (n is the number of replicas which is made go to zero at the end of the calculation) and b = 1/

*T*the inverse temperature. It is related to the functional probability distribution of the local fields by a Fourier-like transformation [10]. For our purposes it is more convenient express

*f*in terms of the probability distribution of the local fields

*P*(

*h*) as in [3], valid for replica symmetric solution, which reads

The global order parameters *g _{n}*(s

_{a}) or

*P*(

*h*) obey their equation of motion [3,8,9] which for the latter is

and the relationship between the two order parameters is

Depending on the quantity of interest it may be easier to work with one or the other of these two parameters although the field distribution has a seemingly clear physical appeal. To infer the form of the general replica symmetric solution is seems better to work with g(y) as originally put forward by Katsura [11] for the Bethe lattice and extended for the VB model by de Almeida et al [9]. The simplest solution at very low temperatures is the one assuming that [3]

which we shall take as the one valid close to zero temperature, except for vanishingly small exponential corrections, for as shown in [12] it is found that *P*(*h*) below the spin glass temperature is almost constant. Using the above equation of motion for *P*(*h*) it is straigthforward to show that the parameters in (7) satisfy

and

where *I _{l}*(

*x*) are modified Bessel functions of order

*l*. From equation (4) it is easy to obtain the entropy per spin

*s*= -

*df*/

*dT*, at

*T*= 0, which reads

and as *I _{l}*(

*x*) @ exp(

*x*)/, for large

*x*, the s limiting value as

*p*® ¥ is positive and given by s @ 0.117 . It should not be surprising that the entropy is positive in view of the discrete nature of the solution (7) but in this limit SK's result should be recovered, which is obviously not the case, pointing out the necessity of including the continuous part of

*P*(

*h*), as locally carried out in [13]. Another way to study the validity of the solution (7) is to consider the fluctuations of the variational free energy around it as carried out in [8],[9] and [14], although in this last one only local stability is probed. This has been done thoroughly for

*p*close to the percolation threshold [8] and the results show that the solution is unstable even for just longitudinal variations . However, one may easily show that (7) is an unstable solution for any

*p*, at zero temperature, by considering the l- eigenvalue equations in [8,9] which for longitudinal variations in

*g*(s

_{n}_{a}) is

where

and for a stable solution, at least for longitudinal variations, l > 0 for all possible solutions of (11). For the solution (7), we find the eigenvalues l_{1} = 1/*p*,l_{2} = (1 - *p*)/*p*(1 - *Q*), and l_{3} = 2(1 - *p*)/*p* clearly showing that the solution (7) is unstable for *p* > 1. So it seems that the continuous part of *P*(*h*) must be taken in account for the VB model if well known results are to be recovered and possibly similar models specially those relevant to the K-SAT problem and neural networks.

** Acknowledgments **

The author is grateful to PRONEX/MCT for partial financial support and B. Wemmenhove [15] for discussion concerning this work.

**References**

[1] M. Mezard, F. Ricci-Tersenghi, R. Zecchina, Alternative solutions to diluted p-spin models and XORSAT problems, cond-mat/0207140; M. Mezard, G. Parisi and R. Zecchina, Science **297**, 812 (2002). [ Links ]

[2] L. Viana and A.J. Bray, J. Phys. C: Solid State Phys. **18**, 3037 (1985). [ Links ]

[3] M. Mezard and G. Parisi, Europhys. Lett. 3, 1067 (1987); [ Links ]I. Kanter and H. Sompolinsky, Phys. Rev. Lett. **58**, 164 (1987). [ Links ]

[4] R. Monasson, Phil. Mag. B**77**, 1515 (1998). [ Links ]

[5] D. Sherrington and S. Kirkpatrick, Phys. Rev. Lett. **35**, 1792, (1975). [ Links ]

[6] J.R.L. de Almeida and D.J. Thouless, J. Phys. A: Math. Gen. **11**, 983 (1978). [ Links ]

[7] M. Mezard, G. Parisi and M. Virasoro, *Spin Glass Theory and Beyond, * World Scientific, Singapore, 1987. [ Links ]

[8] P. Mottishaw and C. de Dominicis, J. Phys. A: Math. Gen. **20**, L375 (1987). [ Links ]

[9] J.R.L. de Almeida, P. Mottishaw and C. de Dominicis, J.Phys. A: Math. Gen. **21**, L693 (1988); [ Links ]J.R.L. de Almeida, P. Mottishaw and C. de Dominicis, New Trends in Magnetism.

[10] Y.Y. Goldschmidt and C. de Dominicis, Phys. Rev. B **41**, 2184 (1990). [ Links ]

[11] S. Katsura, Physica A**141**, 556 (1986). [ Links ]

[12] M. Seino and S. Katsura, Prog. Theor. Phys. Suppl. **115**, 237 (1994). [ Links ]

[13] M.J. Oliveira and M.S. de Toledo Piza, Physica A **169**, 317 (1990). [ Links ]

[14] M.J. Oliveira, Physica A **187**, 112 (1992). [ Links ]

[15] B. Wemmenhove and A.C.C. Coolen, Finite Connectivity Attractor Neural Networks, cond-mat/0304282. [ Links ]

Received on 21 August, 2003

**Erratum**

**On the Entropy of the Viana-Bray model,**

J. R. L. de Almeida

Braz. J. Phys. **33**, 892 (2003)

The eigenvalues for the solution of equation (7) should be written as

and

in the longitudinal sector, and

and

in the transversal sector.

These changes have no effect on the final results and conclusions. The details of the analysis leading to these eigenvalues may also be found in the work of Hase et al.[1] and in some of the references quoted by them.

**References**

[1] M. O. Hase, J. R. L. de Almeida, and S. R. Salinas, cond-mat/0411682.

This article has received corrections in agreement with the ERRATUM published in Volume 35 Number 1.