Multicritical Behavior of Two Coupled Ising Models in the Presence of a Random Field
Abstract
A system defined by two coupled Ising models, with a bimodal random field acting in one of them, is investigated. The interactions among variables of each Ising system are infiniteranged, a limit where mean field becomes exact. This model is studied at zero temperature, as well as for finite temperatures, representing physical situations which are appropriate for describing real systems, such as plastic crystals. A very rich critical behavior is found, depending directly on the particular choices of the temperature, couplings, and randomfield strengths. Phase diagrams exhibiting ordered, partiallyordered, and disordered phases are analyzed, showing the sequence of transitions through all these phases, similarly to what occurs in plastic crystals. Due to the wide variety of critical phenomena presented by the model, its usefulness for describing critical behavior in other substances is also expected.
Keywords: Multicritical Phenomena, RandomField Ising Model, Plastic Crystals.
pacs:
05.70.Fh, 05.70.Jk, 64.60.i, 64.60.KwI Introduction
The study of magnetic models has generated considerable progresses in the understanding of magnetic materials, and lately, it has overcome the frontiers of magnetism, being considered in many areas of knowledge. Certainly, the Ising model represents one of the most studied and important models of magnetism and statistical mechanics huang (); reichl (), and it has been employed also to typify a wide variety of physical systems, like lattice gases, binary alloys, and proteins (with a particular interest in the problem of protein folding). Although real magnetic systems should be properly described by means of Heisenberg spins (i.e., threedimensional variables), many materials are characterized by anisotropy fields that make these spins prefer given directions in space, explaining why simple models, characterized by binary variables, became so important for the area of magnetism. Particularly, models defined in terms of Ising variables have shown the ability for exhibiting a wide variety of multicritical behavior by introducing randomness, and/or competing interactions, has attracted the attention of many researchers (see, e.g., Refs. aharony (); mattis (); kaufman (); nogueira98 (); nuno08a (); nuno08b (); salmon1 (); morais12 ()).
Certainly, the simplicity of Ising variables, which are very suitable for both analytical and numerical studies, has led to proposals of important models outside the scope of magnetism, particularly in the area of complex systems. These models have been successful for describing a wide variety of relevant features in such systems, and have raised interest in many fields, like financial markets, optimization problems, biological membranes, and social behavior. In some cases, more than one Ising variable have been used, especially by considering a coupling between them, as proposed within the framework of choice theories fernandez (), or in plastic crystals plastic1 (); brandreview (); folmer (). In the former case, each set of Ising variables represents a group of identical individuals, all of which can make two independent binary choices.
The socalled plastic crystals plastic1 (); brandreview (); folmer (); michel85 (); michel87 (); galam87 (); galam89 (); salinas1 (); salinas2 () appear as states of some compounds considered to be simpler than those of canonical glasses, but still presenting rather nontrivial relaxation and equilibrium properties. Such a plastic phase corresponds to an intermediate stable state, between a hightemperature (disordered) liquid phase, and a lowtemperature (ordered) solid phase and both transitions, namely, liquidplastic and plasticsolid, are first order. In this intermediate phase, the rotational disorder coexists with a translationally ordered state, characterized by the centers of mass of the molecules forming a regular crystalline lattice with the molecules presenting disorder in their orientational degrees of freedom, as shown in Fig. 1. Many materials undergo a liquidplastic phase transition, where the lowertemperature phase presents such a partial orientational order, like the plasticcrystal of Fig. 1. The property of translational invariance makes the plastic crystals much simpler to be studied from both analytical and numerical methods, becoming very useful towards a proper understanding of the glass transition plastic1 (); brandreview (); folmer (). In some plasticcrystal models one introduces a coupling between two Ising models, associating each of these systems respectively, to the translational and rotational degrees of freedom galam87 (); galam89 (); salinas1 (); salinas2 (), as a proposal for explaining satisfactorily thermodynamic properties of the plastic phase.
Accordingly, spin variables and are introduced in such a way to mimic translational and rotational degrees of freedom of each molecule , respectively. The following Hamiltonian is considered galam87 (); galam89 (); salinas1 (); salinas2 (),
(1) 
where represents a sum over distinct pairs of nearestneighbor spins. In the first summation, the Ising variables may characterize two lattices A and B (or occupied and vacant sites). One notices that the rotational variables could be, in principle, continuous variables, although the fact that the minimization of the coupling contribution is achieved for (), or for (), suggests the simpler choice of binary variables () to be appropriate, based on the energy minimization requirement.
In the present model the variables and represent very different characteristics of a molecule. Particularly, the rotational variables are expected to change more freely than the translational ones; for this reason, one introduces a random field acting only on the rotational degrees of freedom. In fact, the whole contribution is known to play a fundamental role for the plastic phase of ionic plastic crystals, like the alkalicyanides KCN, NaCN and RbCN. In spite of its simplicity, the above Hamiltonian is able to capture the most relevant features of the plasticcrystal phase, as well as the associated phase transitions, namely, liquidplastic and plasticsolid ones michel85 (); michel87 (); galam87 (); galam89 (); salinas1 (); salinas2 (); vives ().
A system described by a Hamiltonian slightly different from the one of Eq. (1), in which the whole contribution was replaced by , i.e., with no random field acting on variable separately, was considered in Ref. salinas2 (). In such a work one finds a detailed analysis of the phase diagrams and orderparameter behavior of the corresponding model. However, to our knowledge, previous investigations on the model defined by Eq. (1) have not considered thoroughly the effects of the random field , with a particular attention to the phase diagrams for the case of a randomly distributed bimodal one, ; this represents the main motivation of the present work. In the next section we define the model, determine its freeenergy density, and describe the numerical procedure to be used. In Section III we exhibit typical phase diagrams and analyze the behavior of the corresponding order parameters, for both zero and finite temperatures; the ability of the model to exhibit a rich variety of phase diagrams, characterized by multicritical behavior, is shown. Finally, in Section IV we present our main conclusions.
Ii The Model and FreeEnergy Density
Based on the discussion of the previous section, herein we consider a system composed by two interacting Ising models, described by the Hamiltonian
(2) 
where represent sums over all distinct pairs of spins, a limit for which the meanfield approximation becomes exact. Moreover, and () depict Ising variables, stands for a real parameter, whereas both and are positive coupling constants, which will be restricted herein to the symmetric case, . Although this later condition may seem as a rather artificial simplification of the Hamiltonian in Eq. (1), the application of a random field acting separately on one set of variables, will produce the expected distinct physical behavior associated with and . The random fields will be considered as following a symmetric bimodal probability distribution function,
(3) 
The infiniterange character of the interactions allows one to write the above Hamiltonian in the form
(4) 
from which one may calculate the partition function associated with a particular configuration of the fields ,
(5) 
where and indicates a sum over all spin configurations. One can now make use of the HubbbardStratonovich transformation dotsenkobook (); nishimoribook () to linearize the quadratic terms,
(6) 
where depends on the random fields , as well as on the spin variables, being given by
(7) 
Performing the trace over the spins and defining new variables, related to the respective order parameters,
(8) 
one obtains
(9) 
where
(10)  
Now, one takes the thermodynamic limit (), and uses the saddlepoint method to obtain
(11) 
where the freeenergy density functional results from a quenched average of in Eq. (10), over the bimodal probability distribution of Eq. (3),
(12) 
with
(13) 
The extremization of the freeenergy density above with respect to the parameters and yields the following equations of state,
(15) 
where
(16) 
In the following section we present numerical results for the order parameters and phase diagrams of the model, at both zero and finite temperatures. All phase diagrams are represented by rescaling conveniently the energy parameters of the system, namely, , and . Therefore, for given values of these dimensionless parameters, the equations of state [Eqs.(II) and (15)] are solved numerically for and . In order to avoid metastable states, all solutions obtained for and are substituted in Eq. (12), to check for the minimization of the freeenergy density. The continuous (second order) critical frontiers are found by the set of input values for which the order parameters fall continuously down to zero, whereas the firstorder frontiers were found through Maxwell constructions.
Both ordered ( and ) and partiallyordered ( and ) phases have appeared in our analysis, and will be labeled accordingly. The usual paramagnetic phase (P), given by , always occurs for sufficiently high temperatures. A wide variety of critical points appeared in our analysis (herein we follow the classification due to Griffiths griffiths ()): (i) a tricritical point signals the encounter of a continuous frontier with a firstorder line with no change of slope; (ii) an ordered critical point corresponds to an isolated critical point inside the ordered region, terminating a firstorder line that separates two distinct ordered phases; (ii) a triple point, where three distinct phases coexist, signaling the encounter of three firstorder critical frontiers. In the phase diagrams we shall use distinct symbols and representations for the critical points and frontiers, as described below.

Continuous (second order) critical frontier: continuous line;

Firstorder critical frontier: dotted line;

Tricritical point: located by a black circle;

Ordered critical point: located by a black asterisk;

Triple point: located by an empty triangle.
Iii Phase Diagrams and Behavior of Order Parameters
iii.1 ZeroTemperature Analisis
At , one has to analyze the different spin orderings that minimize the Hamiltonian of Eq. (4). Due to the coupling between the two sets of spins, the minimumenergy configurations will correspond to and antiparallel (), or parallel (). Therefore, in the absence of random fields () one should have (), and (), where . However, when random fields act on the spins, there will be a competition between these fields and the coupling parameter , leading to several phases, as represented in Fig. 2, in the plane versus . One finds three ordered phases for sufficiently low values of and , in addition to P phases for and . All frontiers shown in Fig. 2 are firstorder critical lines.
When one finds ordered phases for all values of , with a vertical straight line at separating the symmetric state (), where , from the antisymmetric one (), characterized by . Two critical frontiers (symmetric under a reflection operation) emerge from the triple point at and , given, respectively, by for , and for . These critical frontiers terminate at and separate the low randomfieldvalued ordered phases from a partiallyordered phase, given by and . As shown in Fig. 2, three triple points appear, each of them signaling the encounter of three firstorder lines, characterized by a coexistence of three phases, defined by distinct values of the magnetizations and , as described below.

and : .

and : .

and : .
Such a rich critical behavior shown for suggests that interesting phase diagrams should occur when the temperature is taken into account. From now on, we investigate the model defined by the Hamiltonian of Eq. (4) for finite temperatures.
iii.2 FiniteTemperature Analysis
As shown above, the zerotemperature phase diagram presents a reflection symmetry with respect to (cf. Fig. 2). The only difference between the two sides of this phase diagram concerns the magnetization solutions characterizing the ordered phases for low randomfield values, where one has (), or (). These results come as a consequence of the symmetry of the Hamiltonian of Eq. (4), which remains unchanged under the operations, , , or , , . Hence, the finitetemperature phase diagrams should present similar symmetries with respect to a change . From now on, for the sake of simplicity, we will restrict ourselves to the case , for which the zerotemperature and lowrandomfield magnetizations present opposite signals, as shown in Fig. 2, i.e., .
In Fig. 3 we exhibit phase diagrams of the model in two particular cases, namely, in the absence of fields [Fig. 3(a)] and for zero coupling [Fig. 3(b)]. These figures provide useful reference data in the numerical procedure to be employed for constructing phase diagrams in more general situations, e.g., in the plane versus , for several values of .
In Fig. 3(a) we present the phase diagram of the model in the plane of dimensionless variables versus , in the absence of random fields (), where one sees the point that corresponds to two noninteracting Ising models, leading to the wellknown meanfield critical temperature of the Ising model []. Also in Fig. 3(a), the ordered solution minimizes the free energy at low temperatures for any ; a secondorder frontier separates this ordered phase from the paramagnetic one that appears for sufficiently high temperatures. For high values of one sees that this critical frontier approaches asymptotically . Since the application of a random field results in a decrease of the critical temperature, when compared with the one of the case aharony (); mattis (); kaufman (), the result of Fig. 3(a) shows that no ordered phase should occur for and .
The phase diagram for is shown in the plane of dimensionless variables versus in Fig. 3(b). The P phase occurs for , whereas for two phases appear, namely, the ordered one (characterized by and , with ), as well as the partiallyordered phase ( and ). Since the two Ising models are uncorrelated for and the random fields act only on the variables, one finds that the critical behavior associated with variables and occur independently: (i) The variables order at , for all values of ; (ii) The critical frontier shown in Fig. 3(b), separating the two lowtemperature phases, is characteristic of an Ising ferromagnet in the presence of a bimodal random field aharony (). The black circle denotes a tricritical point, where the highertemperature continuous frontier meets the lowertemperature firstorder critical line. The type of phase diagram exhibited in Fig. 3(b) will be referred herein as topology I.
The effects of a small interaction [] between the variables and are presented in Fig. 4, where one sees that the topology I [Fig. 3(b)] goes through substantial changes, as shown in Fig. 4(a) (to be called herein as topology II). As expected from the behavior presented in Fig. 3(a), one notices that the border of the P phase (a continuous frontier) is shifted to higher temperatures. However, the most significant difference between topologies I and II consists in the lowtemperature frontier separating the ordered and partiallyordered phases. Particularly, the continuous frontier, as well as the tricritical point shown in Fig. 3(b), give place to an ordered critical point griffiths (), at which the lowtemperature firstorder critical frontier terminates. Such a topology has been found also in some random magnetic systems, like the Ising and BlumeCapel models, subject to random fields and/or dilution kaufman (); salmon1 (); salmon2 (); benyoussef (); carneiro (); kaufmankanner (). In the present model, we verified that topology II holds for any , with the firstorder frontier starting at zero temperature and , which in Fig. 4(a) corresponds to . Such a firstorder line essentially affects the parameter , as will be discussed next.
In Figs. 4(b)–(d) the order parameters and are exhibited versus for conveniently chosen values of , corresponding to distinct physical situations of the phase diagram for . A curious behavior is presented by the magnetization by varying , and more particularly, around the firstorder critical line. For [Fig. 4(c)], one starts at low temperatures essentially to the left of the critical frontier and by increasing one crosses this critical frontier at , very close to the ordered critical point. At this crossing point, presents an abrupt decrease, i.e., a discontinuity, corresponding to a change to the partiallyordered phase; on the other hand, the magnetization remains unaffected when going through this critical frontier. For higher temperatures, becomes very small, but still finite, turning up zero only at the P boundary; in fact, the whole region around the ordered critical point is characterized by a finite small value of . Another unusual effect is presented in Fig. 4(d), for which , i.e., slightly to the right of the firstorder critical frontier: the order parameter is zero for low temperatures, but becomes nonzero by increasing the temperature, as one becomes closer to the critical ordered point. This rather curious phenomenon is directly related to the correlation between the variables and : since for the magnetization is still very close to its maximum value, a small value for is induced, so that both order parameters go to zero together only at the P frontier.
Behind the results presented in Figs. 4(a)–(d) one finds a very interesting feature, namely, the possibility of going continuously from the ordered phase to the partiallyordered phase by circumventing the ordered critical point. This is analogous to what happens in many substances, e.g., water, where one goes continuously (with no latent heat) from the liquid to the gas phase by circumventing a critical end point huang (); reichl ().
In Fig. 5 the freeenergy density of Eq. (12) is analyzed at two different points along the firstorder critical frontier of Fig. 4(a), namely, a lowtemperature one [Fig. 5(a)], and a point at a higher temperature [Fig. 5(b)]. In both cases the free energy presents four minima associated with distinct pairs of solutions . The point at presents , whereas the point at presents . The lowertemperature point represents a coexistence of the two phases shown in the case [cf. Fig. 3(b)], namely, the ordered () and partiallyordered () phases. However, the highertemperature point typifies the phenomenon discussed in Fig. 4, where distinct solutions with coexist, leading to a jump in this order parameter as one crosses the critical frontier, like illustrated in Fig. 4(c) for the point . Although the magnetization presents a very curious behavior in topology II [cf., e.g., Figs. 4(b)–(d)], remains essentially unchanged by the presence of the firstorder critical frontier of Fig. 4(a), as shown also in Fig. 5.
In Fig. 6 we present two other possible phase diagrams, namely, the cases [Fig. 6(a), called herein topology III] and [Fig. 6(b), called herein topology IV]. Whereas topology III represents a special situation that applies only for , exhibiting the richest critical behavior of the present model, topology IV holds for any . In Fig. 6(a) one observes the appearance of several multicritical points, denoted by the black circle (tricritical point), black asterisk (ordered critical point), and empty triangles (triple points): (i) The tricritical point, which signals the encounter of the highertemperature continuous phase transition with the lowertemperature firstorder phase transition, found in the phase diagram [cf. Fig. 3(b)], have curiously disappeared for , and emerged again for ; (ii) The ordered critical point exists for any [as shown in Fig. 4(a)]; (iii) Two triple points, one at a finite temperature, whereas the other one occurs at zero temperature. It should be mentioned that such a zerotemperature triple point corresponds precisely to the one of Fig. 2, at and . The value is very special and will be considered as a threshold for both multicritical behavior and correlations between the two systems. We have observed that for , the critical points shown in Fig. 6(a) disappear, except for the tricritical point that survives for [as shown in Fig. 6(b)]. Changes similar to those occurring herein between topologies II and III, as well as topologies III and IV, were found also in some magnetic systems, like the Ising and BlumeCapel models, subject to random fields and/or dilution kaufman (); salmon1 (); salmon2 (); benyoussef (); carneiro (); kaufmankanner (). Particularly, the splitting of the lowtemperature firstorder critical frontier into two highertemperature firstorder lines that terminate in the ordered and tricritical points, respectively [as exhibited in Fig. 6(a)], is consistent with results found in the BlumeCapel model under a bimodal random magnetic, by varying the intensity of the crystal field kaufmankanner ().
Another important feature of topology III concerns the lack of any type of magnetic order at finite temperatures for , in contrast to the phase diagrams for , for which there is for all [see, e.g., Figs. 3(b) and 4(a)]. This effect shows that represents a threshold value for the coupling between the variables and , so that for the correlations among these variables become significant. As a consequence of these correlations, the fact of no magnetic order on the system () drives the the magnetization of the system to zero as well, for . It is important to notice that the phase diagram of Fig. 2 presents a firstorder critical line for and , at which , whereas in the system both and minimize the Hamiltonian. By analyzing numerically the freeenergy density of Eq. (12) at low temperatures and , we have verified that for any infinitesimal value of destroys such a coexistence of solutions, leading to a minimum free energy at . Consequently, one finds that the lowtemperature region in the interval becomes part of the P phase. Hence, the phase diagram in Fig. 6(a) presents a reentrance phenomena for . In this region, by lowering the temperature gradually, one goes from a P phase to the ordered phase ( ; ), and then back to the P phase. This effect appears frequently in both theoretical and experimental investigations of disordered magnets dotsenkobook (); nishimoribook ().
In Fig. 7 we analyze the behavior of the and for topology III, in the region of multicritical points of the phase diagram for , along three typical thermodynamic paths, as shown in Fig. 7(a). In Fig. 7(b) we exhibit the behavior of and along path (1), where one sees that both parameters go through a jump by crossing the firstorder critical line [], expressing a coexistence of different types of solutions for and at this point. One notices a larger jump in , so that to the right of the ordered critical point one finds a behavior similar to the one verified in topology II, where becomes very small, whereas still presents significant values. Then, by further increasing the temperature, these parameters tend smoothly to zero at the continuous critical frontier separating the ordered and P phases. In Fig. 7(c) we show the magnetizations and along path (2), within the region of the phase diagram where the reentrance phenomenon occurs; along this path, one increases the temperature, going from the P phase to the ordered phase and then to the P phase again. Both parameters are zero for low enough temperatures, jumping to nonzero values at , as one crosses the firstorder critical line. After such jumps, by increasing the temperature, these parameters tend smoothly to zero at the border of the P phase. The behavior shown in Fig. 7(c) confirm the reentrance effect, discussed previously. Finally, in Fig. 7(d) we exhibit the order parameters along thermodynamic path (3), for which the temperature is fixed at , with the field varying in the range . One sees that both magnetizations and display jumps as one crosses each of the two firstorder lines, evidencing a coexistence of different ordered states at the lowertemperature jump, as well as a coexistence of the ordered and P states at the highertemperature jump.
The behavior presented by the order parameters in Figs. 7(b)–(d) shows clearly the fact that represents a threshold value for the coupling between the variables and . In all these cases, one sees that jumps in the magnetization are correlated with corresponding jumps in . These results should be contrasted with those for the cases , as illustrated in Fig. 4(c), where a discontinuity in does not affect the smooth behavior presented by .
The phase diagram shown in Fig. 6(b), which corresponds to topology IV, is valid for any for any . Particularly, the critical point where the lowtemperature firstorder critical frontier touches the zerotemperature axis is kept at , for all , in agreement with Fig. 2. We have verified only quantitative changes in such a phase diagram by increasing the coupling between the variables and . Essentially, the whole continuous critical frontier moves towards higher temperatures, leading to an increase in the values of the critical temperature for , as well as in the temperature associated with the tricritical point, whereas the abscissa of this point remains typically unchanged. Moreover, in what concerns the order parameters, the difference between and decreases, in such a way that for , one obtains . This later effect is illustrated in Fig. 8, where we represent the order parameters and versus temperature, for a sufficiently large value of , namely, , in two typical choices of , close to the tricritical point. In Fig. 8(a) and are analyzed slightly to the left of the tricritical point, exhibiting the usual continuous behavior, whereas in Fig. 8(b) they are considered slightly to the right of the tricritical point, presenting jumps as one crosses the firstorder critical frontier. However, the most important conclusion from Fig. 8 concerns the fact that in both cases one has essentially , showing that the random field applied solely to the system influences the system in a similar way, due to the high value of considered. We have verified that for the two systems become so strongly correlated, such that holds along the whole phase diagram, within our numerical accuracy.
Iv Conclusions
We have analyzed the effects of a coupling between two Ising models, defined in terms of variables and . The model was considered in the limit of infiniterange interactions, where all spins in each system interact by means of an exchange coupling , typical of ferromagnetic interactions. Motivated by a qualitative description of systems like plastic crystals, the variables and would represent rotational and translational degrees of freedom, respectively. Since the rotational degrees of freedom are expected to change more freely than the translational ones, a random field acting only on the variables was considered. For this purpose, a bimodal random field, , with equal probabilities, was defined on the system. The model was investigated through its free energy and its two order parameters, namely, and .
We have shown that such a system presents a very rich critical behavior, depending on the particular choices of and . Particularly, at zero temperature, the phase diagram in the plane versus exhibits ordered, partiallyordered, and disordered phases. This phase diagram is symmetric around , so that for sufficiently low values of one finds ordered phases characterized by () and (). We have verified that plays an important role in the present model, such that at zero temperature one has the disordered phase () for and . Moreover, the partiallyordered phase, where and , occurs for and . In this phase diagram all phase transitions are of the firstorder type, and three triple points were found. In the case of plastic crystals, the sequence of transitions from the disordered to the partiallyordered, and then to the ordered phases, would correspond to the sequence of transitions from the liquid to the plastic crystal, and then to ordered crystal phases.
Due to the symmetry around , the finitetemperature phase diagrams were considered only for , for which the ordered phase was identified by and , whereas the partiallyordered phase by and (equivalent solutions also exist by inverting the signs of these order parameters). Several phase diagrams in the plane versus were studied, by varying gradually . We have found four qualitatively different types of phase diagrams, denominated as topologies I [], II [], III [], and IV []. Such a classification reflects the fact that represents a threshold value for the coupling between the variables and , so that for the correlations among these variables become significant, as verified through the behavior of the order parameters and . From all these cases, only topology IV typifies a wellknown phase diagram, characterized by a tricritical point, where the highertemperature continuous frontier meets the lowertemperature firstorder critical line. This phase diagram is qualitatively similar to the one found for the Ising ferromagnet in the presence of a bimodal random field aharony (), and it does not present the herein physically relevant partiallyordered phase. For , even though the random field is applied only in the system, the correlations lead the system to follow a qualitatively similar behavior.
The phase diagrams referred as topologies I and II exhibit all three phases. In the later case we have found a firstorder critical line terminating at an ordered critical point, leading to the potential physical realization of going continuously from the ordered phase to the partiallyordered phase by circumventing this critical point. In these two topologies, the sequence of transitions from the disordered to the partiallyordered, and then to the ordered phase, represents the physical situation that occurs in plastic crystals. For conveniently chosen thermodynamic paths, i.e., varying temperature and random field appropriately, one may go from the liquid phase (), to a plasticcrystal phase (; ), where the rotational degrees of freedom are found in a disordered state, and then, to an ordered crystal phase (; ).
From the point of view of multicritical behavior, topology III [] corresponds to the richest type of phase diagram, being characterized by several critical lines and multicritical points; one finds its most complex criticality around , signaling a great competition among the different types of orderings. Although the partiallyordered phase does not appear in this particular case, one has also the possibility of circumventing the ordered critical point, such as to reach a region of the phase diagram along which becomes very small, resembling a partiallyordered phase.
Since the infiniterange interactions among variables of each Ising system correspond to a limit where meanfield approach becomes exact, an immediate question concerns whether some of the results obtained above represent an artifact of such limit. Certainly, such a relevant point is directly related with the existence of some of these features in the associated shortrange threedimensional magnetic models. For example, the tricritical point found in topologies III and IV is essentially the same that appears within the meanfield approach of the Ising model in the presence of a bimodal random field. This later model has been extensively investigated on a cubic lattice through different numerical approaches, where the existence of this tricritical point is still very controversial. On the other hand, a firstorder critical frontier terminating at an ordered critical point, and the fact that one can go from one phase to another by circumventing this point, represents a typical physical situation that occurs in real substances. The potential for exhibiting such a relevant feature represents an important advantage of the present model.
Finally, we emphasize that the rich critical behavior presented in the phase diagrams corresponding to topologies II and III suggest the range as appropriate for describing plastic crystals. The potential of exhibiting successive transitions from the ordered to the partiallyordered and then to the disordered phase should be useful for a better understanding of these systems. Furthermore, the characteristic of going continuously from the ordered phase to the partiallyordered phase by circumventing an ordered critical point represents a typical physical situation that occurs in many substances, and opens the possibility for the present model to describe a wider range of materials.
Acknowledgments
The partial financial supports from CNPq, FAPEAMProjetoUniversalAmazonas, and FAPERJ (Brazilian agencies) are acknowledged.
Footnotes
 thanks: Email address:octaviors@gmail.com
 thanks: Corresponding author: Email address: fdnobre@cbpf.br
References
 K. Huang, Statistical Mechanics, second edition (John Wiley and Sons, New York, 1987).
 L. E. Reichl, A Modern Course in Statistical Physics, second edition (John Wiley and Sons, New York, 1998).
 A. Aharony, Phys. Rev. B 18, 3318 (1978).
 D. C. Mattis, Phys. Rev. Lett. 55, 3009 (1985).
 M. Kaufman, P. E. Kluzinger and A. Khurana, Phys. Rev. B 34, 4766 (1986).
 E. Nogueira Jr., F. D. Nobre, F. A. da Costa, and S. Coutinho, Phys. Rev. E 57, 5079 (1998); Erratum, Phys. Rev. E 60, 2429 (1999).
 N. Crokidakis and F. D. Nobre, J. Phys. Condens. Matter 20, 145211 (2008).
 N. Crokidakis and F. D. Nobre, Phys. Rev. E 77, 041124 (2008).
 O. R. Salmon, N. Crokidakis and F. D. Nobre, J. Phys.: Condens. Matter 21, 056005 (2009).
 C. V. Morais, S. G. Magalhães, and F. D. Nobre, J. Stat. Mech. P01013 (2012).
 A. Fernandez del Río, E. Korutcheva and J. de la Rubia, Complexity, 17, 31 (2012).
 J. N. Sherwood, Ed., The Plastically Crystalline State (John Wiley and Sons, New York, 1979).
 R. Brand, P. Lunkenheimer, and A. Loidl, J. Chem. Phys. 116, 10386 (2002).
 J. C. W. Folmer, R. L. Withers, T. R. Welberry and J. D. Martin, Phys. Rev. B 77, 144205 (2008).
 K. H. Michel, and J. M. Rowe, Phys. Rev. B 32, 5818 (1985).
 K. H. Michel, and J. Naudts, J. Chem. Phys. 67, 547 (1987).
 S. Galam, Phys. Lett. A 122, 271 (1987).
 S. Galam and M. Gabay, Europhys. Lett. 8, 167 (1989).
 S. Galam, V. B. Henriques and S. R. Salinas, Phys. Rev. B 42, 6720 (1990).
 S. Galam, S. R. Salinas and Y. Shapir, Phys. Rev. B 51, 2864 (1995).
 E. Vives and A. Planes, Phys. Rev. B 43, 13335 (1991).
 V. Dotsenko, Introduction to the Replica Theory of Disordered Statistical Systems, Cambridge University Press, Cambridge, 2001.
 H. Nishimori, Statistical Physics of Spin Glasses and Information Processing, Oxford University Press, Oxford 2001.
 R. B. Griffiths, Phys. Rev. B 12, 345 (1975).
 O. R. Salmon and J. Rojas, J. Phys. A: Math. Theor. 43, 125003 (2010).
 A. Benyoussef, T. Biaz, M. Saber and M. Touzani, J. Phys. C 20, 5349 (1987).
 C. E. I. Carneiro, V. B. Henriques and S. R. Salinas, J. Phys. A: Math. Gen. 23, 3383 (1990).
 M. Kaufman and M. Kanner, Phys. Rev. B 42, 2378 (1990).