Next Article in Journal
Artificial Intelligence (AI) Trust Framework and Maturity Model: Applying an Entropy Lens to Improve Security, Privacy, and Ethical AI
Next Article in Special Issue
Heat-Bath and Metropolis Dynamics in Ising-like Models on Directed Regular Random Graphs
Previous Article in Journal
Shannon Entropy of Ramsey Graphs with up to Six Vertices
Previous Article in Special Issue
Relevant Analytic Spontaneous Magnetization Relation for the Face-Centered-Cubic Ising Lattice
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characteristics of an Ising-like Model with Ferromagnetic and Antiferromagnetic Interactions

by
Boris Kryzhanovsky
1,
Vladislav Egorov
1 and
Leonid Litinskii
2,*
1
Scientific Research Institute for System Analysis RAS, Moscow 117218, Russia
2
Independent Researcher, Haifa 2641418, Israel
*
Author to whom correspondence should be addressed.
Entropy 2023, 25(10), 1428; https://doi.org/10.3390/e25101428
Submission received: 4 September 2023 / Revised: 27 September 2023 / Accepted: 28 September 2023 / Published: 9 October 2023
(This article belongs to the Special Issue Ising Model: Recent Developments and Exotic Applications II)

Abstract

:
In the framework of mean field approximation, we consider a spin system consisting of two interacting sub-ensembles. The intra-ensemble interactions are ferromagnetic, while the inter-ensemble interactions are antiferromagnetic. We define the effective number of the nearest neighbors and show that if the two sub-ensembles have the same effective number of the nearest neighbors, the classical form of critical exponents ( α = 0 , β = 1 / 2 , γ = γ = 1 , δ = 3 ) gives way to the non-classical form ( α = 0 , β = 3 / 2 , γ = γ = 0 , δ = 1 ), and the scaling function changes simultaneously. We demonstrate that this system allows for two second-order phase transitions and two first-order phase transitions. We observe that an external magnetic field does not destroy the phase transitions but only shifts their critical points, allowing for control of the system’s parameters. We discuss the regime when the magnetization as a function of the magnetic field develops a low-magnetization plateau and show that the height of this plateau abruptly rises to the value of one when the magnetic field reaches a critical value. Our analytical results are supported by a Monte Carlo simulation of a three-dimensional layered model.

1. Introduction

Magnetic properties of multilayer superlattices have been subject to intensive theoretical and experimental study. Special interest is given to multilayer structures consisting of thin layers of various ferromagnetic and antiferromagnetic materials [1,2,3]. An example of such structures is designs with alternating ferromagnetic and non-ferromagnetic layers. In these structures, the thickness of the non-ferromagnetic layer can be chosen in such a way that the long-range exchange interaction between the ferromagnetic layers is antiferromagnetic [4]. This sort of structure has a giant magnetoresistance. Ferromagnetic layered structures with antiferromagnetic inter-layer interactions are theoretically investigated in [5,6] in the framework of the Heisenberg model. These studies reveal that this kind of system can have many distinct phases: ferromagnetic, antiferromagnetic, paramagnetic, and a whirling phase. Papers [7,8] consider the effect of weak interaction between two square Ising lattices on the phase transitions in this system. It is shown that if a structure has alternating layers of different ferromagnetic materials, there might exist a compensation temperature, which is a temperature below the critical value, at which the full magnetization of the lattice is zero. The conditions for the appearance of the compensation effect and the critical behavior of such systems near the compensation temperature are investigated in [9,10,11,12,13,14,15]. In recent years, great consideration has been given to the study of spin systems with competing ferromagnetic and antiferromagnetic interactions, or with interactions between different groups of spins [16,17,18,19]. This sort of system usually has a complex landscape of free energy containing many local minima separated from the global minimum by a deep potential barrier.
In the works [20,21,22,23,24], the authors investigated the spin-1 Ising models, where a spin has three possible states. They showed that in these models the critical temperature shifts with increasing crystal field and at the certain values of the field the transition order changes from the second to the first order. In the present paper, we demonstrate that such critical behavior is possible in the two-state spin systems.
This paper considers the mean field approximation for the simplest Ising-like model having not only the global minimum, but also a local minimum of energy. Spins are divided into two subgroups (sub-ensembles). The intra-group interaction of spins is ferromagnetic, while spins from different sub-ensembles can interact in both a ferromagnetic and antiferromagnetic manner. With antiferromagnetic inter-group interaction, the ground state is the ferromagnetic ordering within sub-ensembles and antiferromagnetic ordering among spins of different sub-ensembles. Yet, there is also a local energy minimum corresponding to a complete ferromagnetic order.
To compare the mean-field model characteristics with real systems, we carried out a Monte Carlo simulation of a 3D layered model with ferromagnetic intra-layer interactions and antiferromagnetic inter-group interactions. We took into account only the interactions between the nearest neighbors. This model corresponds to realistic ferromagnetic structures with giant magnetoresistance. Although it is expected that the results of the mean-field approach will differ from the results for real systems, it is well known that this approach approximates the 3D systems quite well. We found that the simulation results agree qualitatively with the predictions of our model. The agreement is expected to be better for lattices with larger dimensions.
The paper is organized as follows. State equations for the mean-field model are given in Section 2. Section 3 investigates the critical properties of the system and derives the values of the critical exponents. The effect of an external magnetic field on the temperature dependencies of magnetization is considered in Section 4. Section 5 compares the results of the Monte Carlo simulation with the analytical results obtained in the previous chapters. The discussion and our conclusions are given in Section 6. To conclude, we would like to mention that in our paper [25], for the first time, we used a similar technique.

2. Basic Expressions

2.1. Equations of State and Critical Temperature

Let us consider a system of N spins with a connection matrix J ^ . The system is divided into two groups I and II holding p 1 N and p 2 N spins, respectively ( N = p 1 N + p 2 N , p 1 + p 2 = 1 ). The inter-spin connection in group I is defined by a quantity J 11 , in group II by J 22 , and cross-connections between the spins of group I and group II are defined by a quantity J 12 :
J 11 = q 11 N ,   J 12 = q 12 N ,   J 22 = q 22 N .
For definiteness, in what follows, we assume that spin interaction within each group is of ferromagnetic nature, the interaction between spins from different groups is antiferromagnetic, and the direction of the external magnetic field H is positive:
q 11 0 ,   q 22 0 ,   q 12 0 ,   H 0 .
The case q 12 > 0 in the absence of the external field was investigated earlier in [26,27].
The Hamiltonian of this system can be written as
E N = ( J 11 i , j = 1 N p 1 S i S j + J 22 i , j = 1 N p 2 S i S j + J 12 i = 1 N p 1 j = 1 N p 2 S i S j ) H ( i = 1 N p 1 S i + j = 1 N p 2 S j ) ,
where S and S are the values of the spins from groups I and II, respectively. The sums in (2) are taken over all the spins of the corresponding groups. Introducing partial magnetizations of the two groups, m 1 and m 2 :
m 1 = 1 N p 1 i = 1 N p 1 S i ,   m 2 = 1 N p 2 j = 1 N p 2 S j ,
We can represent the energy of the system per one spin in the mean-field approximation, E = E N / N , as follows:
E = 1 2 ( q 11 p 1 2 m 1 2 + 2 q 12 p 1 p 2 m 1 m 2 + q 22 p 2 2 m 2 2 ) H ( p 1 m 1 + p 2 m 2 ) ,
and the full magnetization of the system becomes: M = p 1 m 1 + p 2 m 2 .
The statistical sum of this spin system is determined as:
Z = n 1 = 0 N p 1 n 2 = 0 N p 2 ( N p 1 n 1 ) ( N p 2 n 2 ) e N K E ,
where K = 1 / T is the inverse temperature, n 1 and n 2 denote, respectively, the number of downward spins (−1) in groups I and II; n 1 and n 2 are related to partial magnetizations by the following expressions:
n 1 = N p 1 1 + m 1 2 , n 2 = N p 2 1 + m 2 2 .
If the total number of spins N is large, then summation in Formula (4) can be replaced by integration. Using Stirling’s formula to estimate binomial coefficients, we obtain
Z = 1 1 1 1 e N F ( m 1 , m 2 ) d m 1 d m 2 ,
where
F ( m 1 , m 2 ) = p 1 S 1 ( m 1 ) + p 2 S 2 ( m 2 ) + K E ( m 1 , m 2 ) ,
and
S i ( m i ) = 1 + m i 2 ln 1 + m i 2 + 1 m i 2 ln 1 m i 2 ,   i = 1 , 2 .
Let us estimate this integral using the saddle point method. The equations for the saddle point ( F / m 1 = 0 , F / m 2 = 0 ) are
1 2 K ln 1 + m 1 1 m 1 = q 11 p 1 m 1 q 12 p 2 m 2 + H 1 2 K ln 1 + m 2 1 m 2 = q 22 p 2 m 2 q 12 p 1 m 1 + H . ,
These equations define the equilibrium values of the partial magnetizations m 1 and m 2  at a given temperature K . Substituting these values into expression (5), we obtain the free energy per spin F .
Let us determine the critical temperature of a phase transition, K c . Assuming H = 0 and m 1 0 , ( K K c ), in the case of q 12 = 0 from (6), we obtain an obvious result: a system split into two non-interacting subsystems has two critical points:
K c ( 1 ) = 1 q 11 p 1 ,   K c ( 2 ) = 1 q 22 p 2
Expression (7) fully agrees with the mean-field theory [8] because in the notations of (1) and (2), the quantities q 11 p 1 and q 22 p 2 are the numbers of the nearest neighbors in groups I and II.
In the general case of q 12 0 , the system has one critical point, which we can find from (6):
K c = 2 q 11 p 1 + q 22 p 2 + ( q 11 p 1 q 22 p 2 ) 2 + 4 p 1 p 2 q 12 2

2.2. Energy Minima and Critical Value H c

To proceed further, it is necessary to investigate the presence of energy minima to which the system can converge when the temperature approaches zero. It is easy to see that the energy (3) can be at its minimum only in configurations where m 1 = ± 1 and m 2 = ± 1 . The stability of such states (the presence of the minimum) is only possible when two conditions are simultaneously met:
H 1 m 1 > 0   and   H 2 m 2 > 0
where H 1 and H 2 are the local fields acting on the spins of groups I and II, respectively:
H 1 = q 11 p 1 m 1 + q 12 p 2 m 2 + H ,   H 2 = q 22 p 2 m 2 + q 12 p 1 m 1 + H .
When there is no external magnetic field ( H = 1 ), the unconditionally stable configurations are S + = S ( m 1 = + 1 , m 2 = 1 ) and S + = S ( m 1 = 1 ,   m 2 = + 1 ) , because spins of both groups are directed along the local fields acting on them: H 1 m 1 = q 11 p 1 + | q 12 | p 2 > 0 and H 2 m 2 = q 22 p 2 + | q 12 | p 1 > 0 . The stability conditions (8) for configurations S + + = S ( m 1 = + 1 , m 2 = + 1 ) and S = S ( m 1 = 1 , m 2 = 1 ) are H 1 m 1 = q 11 p 1 | q 12 | p 2 > 0 and H 2 m 2 = q 22 p 2 | q 12 | p 1 > 0 . Simple calculations show that these configurations correspond to energy minima only when two conditions are met simultaneously:
q 11 q 22 > q 12 2   and   | q 12 | q 11 + | q 12 | < p 1 < q 22 q 22 + | q 12 | .
It follows from (3) that if H = 0 , the configurations S + and S + correspond to a doubly degenerate global energy minimum, while the configurations S + + and S are local minima if (10) is met.
When H 0 , the picture becomes more complicated—the number of minima can vary from one to four. This situation is illustrated in Figure 1, which demonstrates how configuration energies change with the growing strength of the field H . We are interested in the energy values in the configurations S + + ( m 1 = 1 , m 2 = 1 ), S + ( m 1 = 1 , m 2 = 1 ) and S + ( m 1 = 1 , m 2 = 1 ):
E + + = 1 2 ( q 11 p 1 2 + q 22 p 2 2 ) + | q 12 | p 1 p 2 H ,
E + = 1 2 ( q 11 p 1 2 + q 22 p 2 2 ) | q 12 | p 1 p 2 H ( p 1 p 2 ) ,
E + = 1 2 ( q 11 p 1 2 + q 22 p 2 2 ) | q 12 | p 1 p 2 H ( p 2 p 1 ) .
The energy E in the configuration S ( m 1 = 1 , m 2 = 1 ) is of no interest, because it is a minimum only when H is very small.
When condition (10) is not met and the field H is small, the system has two energy minima corresponding to the antiferromagnetic configurations S + and S + (see Figure 1). When the field reaches the value H + + = max ( | q 12 | p 1 q 22 p 2 , | q 12 | p 2 q 11 p 1 ) , the ferromagnetic configuration S + + also becomes the energy minimum. The minima achieved in the configurations S + and S + disappear, when the field H is larger than the values H + = q 22 p 2 + | q 12 | p 1 and H + = q 11 p 1 + | q 12 | p 2 , correspondingly. The only one minimum remains in the configuration S + + . If the condition (10) is met, for all values of the field H there is the minimum in the configuration S + + and when H < min ( q 22 p 2 | q 12 | p 1 , q 11 p 1 | q 12 | p 2 ) there is the fourth minimum in the configuration S .
Let us introduce the critical value of the magnetic field as follows:
H c = | q 12 | p max ,   p max = max ( p 1 , p 2 ) .
Let p 1 > p 2 . Then, from comparing E + + , E + and E + , we can conclude:
  • when H < H c , the antiferromagnetic configuration m 1 = 1 , m 2 = 1 is the ground state of the system, and the ferromagnetic configuration m 1 = 1 , m 2 = 1 is its local minimum;
  • when H > H c , the ground state is the ferromagnetic state m 1 = 1 , m 2 = 1 , and the antiferromagnetic configuration m 1 = 1 , m 2 = 1 is its local minimum.
If p 1 < p 2 , we should swap indices 1 and 2 in these expressions. The aforesaid is summarized in Table 1.

2.3. Balanced System and “Critical” Value q 12 = q 12

Further analysis reveals that the thermodynamic characteristics of balanced and unbalanced systems are considerably different. We call a system balanced if it has the effective number of spin “nearest neighbors” in the first sub-ensemble is equal to that in the second sub-ensemble. This equality can be expressed as:
q 11 p 1 + q 12 p 2 = q 22 p 2 + q 12 p 1 .
Indeed, a spin from the first sub-ensemble has q 11 p 1 neighbors from this same sub-ensemble and q 12 p 2 neighbors from the other sub-ensemble. The same is true for a spin from the second sub-ensemble. We use the term “effective” number of neighbors because the neighbors with antiferromagnetic interaction make a negative contribution to this number.
The balance (12) is possible only when a certain relationship between the interaction parameters holds. For a symmetric system with p 1 = p 2 and q 11 = q 22 , it follows from (12) that the system is balanced for any q 12 . In a more general case, the balance condition (12) can be rewritten as q 12 = q 12 * , where q 12 * < 0 is a certain “critical” value of antiferromagnetic interaction:
q 12 * = q 11 p 1 q 22 p 2 p 1 p 2 ,   q 12 * < 0 .
Below we will use notations q 12 = q 12 * or q 12 q 12 * with q 12 * defined in (13) to indicate whether a system is balanced or unbalanced.
Note that a system can be balanced only within a specific range of parameters p 1 , 2 and q 11 , q 22 . Indeed, it follows from (13) that the condition q 12 * < 0 can be written as:
q 12 * < 0 i f { q 22 q 11 + q 22 p 1 < 1 2 , q 11 q 22 1 2 < p 1 q 22 q 11 + q 22 , q 11 q 22 .
Relationship (14) is illustrated as a phase diagram in Figure 2.

3. Critical Behavior

Let us introduce the relative temperature t as follows
t = T T c T c = K c K K .
We will consider the behavior of various physical quantities near the critical temperature, paying special attention to the difference between critical exponents in the general case ( q 12 q 12 * ) and critical exponents in the case of the “critical” antiferromagnetic interaction ( q 12 = q 12 * ).

3.1. Spontaneous Magnetization and Critical Exponent β

  • The case q 12 q 12 * . Let us consider the equation of state (6) when H = 0 and m 1 , 2 0 . Let us expand it in small parameters m 1 , 2 with an accuracy up to the terms of the order of m 1 , 2 3 and t . Then, for partial magnetizations near the critical temperature ( K > K c ), we obtain:
    m 1 2 = 3 t D 1 ,   m 2 2 = 3 t D 2 ,
    where
    D 1 = p 2 ( 1 K c q 22 p 2 ) [ 2 K c ( q 11 p 1 + q 22 p 2 ) ] p 1 ( 1 K c q 11 p 1 ) 2 + p 2 ( 1 K c q 22 p 2 ) 2 , D 2 = p 1 ( 1 K c q 11 p 1 ) [ 2 K c ( q 11 p 1 + q 22 p 2 ] p 1 ( 1 K c q 11 p 1 ) 2 + p 2 ( 1 K c q 22 p 2 ) 2 .
The full spontaneous magnetization M 0 is defined as:
M 0 = p 1 m 1 + p 2 m 2 = ± 3 t ( p 1 D 1 p 2 D 2 ) .
Here, we take into account the fact that with H = 0 and q 12 < 0 , the inequality m 1 m 2 < 0 holds.
2.
The case q 12 = q 12 * . Here, p 1 D 1 = p 2 D 2 and the parenthesized expression from (16) become zero. It means that in the expressions resulting from the expansion of the equation of state (6), we must retain the terms of the order of t 2 . Then, the spontaneous magnetization near the critical temperature is described as
M 0 = ± ( t ) 3 2 3 p 1 p 2 ( p 1 p 2 ) K c | q 12 * | ( 1 3 p 1 p 2 ) 3 / 2 .
Thus, in an unbalanced system ( q 12 q 12 * ), the critical exponent β = 1 / 2 , which agrees with the conventional mean-field theory. In a balanced system ( q 12 = q 12 * ), the critical exponent takes a “non-classical” value β = 3 / 2 . Here, by the classic value, we mean β = 1 / 2 for the case where the interaction is ferromagnetic only.

3.2. Jump in Heat Capacity ( H = 0 ) and Critical Exponent α

At the critical point, the heat capacity C experiences a finite jump. Indeed, when t > 0 , the heat capacity C = 0 , and when t 0 , the quantity C = K d E / d K can be easily found using (3) and (23). Then, for the heat capacity jump at the critical point, we have:
Δ C = lim t 0 C = 3 2 K c p 1 p 2 ( q 11 p 1 q 22 p 2 ) 2 + 4 p 1 p 2 q 12 2 p 1 ( 1 K c q 11 p 1 ) 2 + p 2 ( 1 K c q 22 p 2 ) 2 .
When q 12 = q 12 * , expression (18) takes the form:
Δ C = 3 2 K c p 1 p 2 p 1 3 + p 2 3 .
Since the heat capacity tends to a finite value at the critical point, the classical definition of the critical exponent α loses sense. Therefore, we will use its alternative definition [28]:
F + ( K ) F ( K ) t 2 α   when   t 0 ,
where F + ( K ) and F ( K ) are the free energies for K < K c and K > K c , respectively. To compute the exponent α from (19), the functions F + ( K ) and F ( K ) must be analytically extended to the K -axis beyond their domain. Using the expression for free energy (5) and taking into account the equation of state (6) and expressions for spontaneous magnetization (23), we obtain:
F + ( K ) F ( K ) = 3 t 2 4 [ 2 K c ( q 11 p 1 2 D 1 2 q 12 p 1 p 2 D 1 D 2 + q 22 p 2 2 D 2 ) + p 1 D 1 2 + p 2 D 2 2 ] .
Comparing (20) with (19), we obtain α = 0 , which agrees with the classical mean-field model. Note that this result is independent of the magnitude of the antiferromagnetic interaction, i.e., it holds both for q 12 = q 12 * and q 12 q 12 * .
Typical curves of spontaneous partial magnetization m 1 , 2 ( K ) and energy variance σ 2 ( K ) = C ( K ) / K as functions of temperature without an external field are shown in Figure 3.

3.3. Susceptibility χ ( H = 0 ) and Critical Exponents γ and γ

Let us consider the behavior of the susceptibility of the system near the critical point, given that H = 0 . We define the full and partial susceptibilities as follows:
χ = M ( K , H ) H = p 1 χ 1 + p 2 χ 2 ,   χ 1 = m 1 ( K , H ) H ,   χ 2 = m 2 ( K , H ) H .
Differentiating the equations of state (6) with respect to H and solving the resulting equations for χ 1 and χ 2 in view of (15), we obtain the following expressions:
{ χ 1 = t + [ 1 K c p 2 ( q 22 q 12 ) ] t ( q 11 p 1 q 22 p 2 ) 2 + 4 p 1 p 1 q 12 2 i f t > 0 χ 2 = t + [ 1 K c p 1 ( q 11 q 12 ) ] t ( q 11 p 1 q 22 p 2 ) 2 + 4 p 1 p 1 q 12 2 ,
and
{ χ 1 = K c [ 1 K c p 2 ( q 22 c ) ] t ( 1 3 D 2 ) t [ ( 1 K c q 11 p 1 ) ( 1 3 D 1 ) + ( 1 K c q 22 p 2 ) ( 1 3 D 2 ) ] i f t < 0 χ 2 = K c [ 1 K c p 1 ( q 11 c ) ] t ( 1 3 D 1 ) t [ ( 1 K c q 11 p 1 ) ( 1 3 D 1 ) + ( 1 K c q 22 p 2 ) ( 1 3 D 2 ) ] ,
  • When q 12 q 12 * , the terms in the square brackets in the numerators of (21) and (22) are not zero, and the quantity t in the numerators can be neglected. In this event, the susceptibility of the system χ = p 1 χ 1 + p 2 χ 2 takes the classical form:
    χ = 1 K c p 1 p 2 ( q 11 + q 22 2 q 12 ) | t | ( q 11 p 1 q 22 p 2 ) 2 + 4 p 1 p 1 q 12 2 i f t > 0 o r t < 0 ,
    holding true for both t > 0 and t < 0 .
  • When q 12 = q 12 * , we have [ 1 K c p 2 ( q 22 q 12 * ) ] = [ 1 K c p 1 ( q 11 q 12 * ) ] = 0 . Then, χ 1 and χ 2 in (21) and (22) become constant values, and full susceptibility χ = p 1 χ 1 + p 2 χ 2 takes the same “non-classical” form for both t > 0 and t < 0 :
    χ = 1 | q 12 * | .
    From (23) and (24), we can derive critical exponents:
    { γ = γ = 1 i f q 12 q 12 * γ = γ = 0 i f q 12 = q 12 *
It should be noted that in a balanced system (12) (at q 12 = q 12 * ), partial susceptibilities χ 1 and χ 2 experience a finite jump at the critical point; yet K -dependence of the full susceptibility χ remains continuous.

3.4. Scaling Hypothesis and Critical Exponent δ

According to the scaling hypothesis, the field H is a homogeneous function of variables M 1 / β and t near the critical point. Let us consider how the critical exponent δ and the form of the scaling function rely on the interaction parameters.
  • The case q 12 q 12 * (unbalanced system). Expanding the equations of state (6) in terms of the small parameters m 1 , m 2 , t and extracting full magnetization M = p 1 m 1 + p 2 m 2 , we obtain:
    K c H = M 3 R 1 + t M R 2 ,
    where
    R 1 = ( 1 K c q 11 p 1 ) 2 p 1 + ( 1 K c q 22 p 2 ) 2 p 2 3 p 1 p 2 [ 1 K c p 1 p 2 ( q 11 + q 22 2 q 12 ) ] 2 ,   R 2 = K c ( q 11 p 1 q 22 p 2 ) 2 + 4 p 1 p 2 q 12 2 1 K c p 1 p 2 ( q 11 + q 22 2 q 12 ) .
    It is easy to see that the scaling hypothesis in this case is confirmed, because if q 12 q 12 * , it follows from (16) that β = 1 / 2 . Indeed, expression (25) can be represented in the classical form K c H = M | M | δ 1 h s ( t | M | 1 / β ) with the critical exponent δ = 3 and the scaling function h s ( x ) of the form:
    h s ( x ) = R 1 + R 2 x ,   where   x = t / M 2 .
  • The case q 12 = q 12 * (balanced system (12)). Performing similar calculations with regard to the relations 1 K c p 1 ( q 11 + | q 12 * | ) = 1 K c p 2 ( q 22 + | q 12 * | ) = 0 , which take place when q 12 = q 12 * , we obtain:
    K c H = { K c | q 12 * | M , t > 0 K c | q 12 * | M | t | 3 / 2 R 3 , t < 0 ,   where   R 3 = 3 p 1 p 2 | p 1 p 2 | ( 1 3 p 1 p 2 ) 3 / 2 .
    In this event, (17) gives β = 3 / 2 . Correspondingly, expression (26) can be rewritten in the classical form K c H = M | M | δ 1 h s ( t | M | 1 / β ) , with the critical exponent δ = 1 and the scaling function h s ( x ) of the form:
    h s ( x ) = { K c | q 12 * | , t > 0 K c | q 12 * | | x | 3 / 2 R 3 , t < 0 ,   where   x = t / M .
It is easy to see that with q 12 = q 12 * , the Rushbrooke relation, α + 2 β + γ 2 , the Widom relation, γ β ( δ 1 ) , and the Griffiths relation, β ( δ + 1 ) 2 α , hold as strict inequalities. At the same time, the monograph [28] claims that, being a consequence of the scaling hypothesis, these relations should turn into strict equalities. Thus, the validity of the scaling hypothesis remains in question when q 12 = q 12 * .
Summing up the above, we should note that the critical exponents of the model under consideration correspond to the parameters of the classical mean-field model provided that q 12 q 12 * , i.e., when the effective numbers of neighbors in different sub-lattices are not equal (see Table 2). However, in a balanced system, when condition (12) is met (i.e., q 12 = q 12 * ), the critical exponents β , γ and δ take “non-classical” values. Then, the relation α + 2 β + γ = 2 , which is a consequence of the scaling hypothesis, is violated.
At the end of this section let us show how functions of certain physical quantities change when q 12 = q 12 * . Figure 4a presents the temperature dependence of susceptibility in the absence of an external field χ ( K , 0 ) . We can see that when q 12 q 12 * , the susceptibility diverges at the critical point, just as in the case of the classical mean-field model. When q 12 = q 12 * , the susceptibility has a finite critical value defined by expression (24). At the “critical” values of antiferromagnetic interactions and critical temperature, the field-dependences of partial magnetizations m 1 ( K c , H ) and m 2 ( K c , H ) also change (Figure 4b). If q 12 q 12 * and the field H is small, the system is in the antiferromagnetic state ( m 1 m 2 < 0 ), and it makes a cross-over to the ferromagnetic state ( m 1 m 2 > 0 ) when the field increases. If q 12 = q 12 * , the system is in the ferromagnetic state at any value of H , and its partial magnetizations are equal ( m 1 = m 2 ).

4. Dependences of Physical Quantities on Temperature and External Field

In this section, we investigate how the physical parameters of the system depend on temperature in the presence of an external field. We find that these dependences change dramatically when the antiferromagnetic interaction parameter q 12 takes a “critical” value q 12 = q 12 * . In the first subsection, we deal with the situation when q 12 q 12 * . In the second subsection, we consider a symmetric configuration q 11 = q 22 ,   p 1 = p 2 = 1 / 2 , which is a particular case of a system with “critical” antiferromagnetic interaction.

4.1. Unbalanced System ( q 12 q 12 * )

As follows from the equations of state (6), at H H c the right sides of both equations of state are necessarily positive and, therefore, when m 1 , 2 > 0 the system is in the ferromagnetic phase at any temperature. That is why the change in phase is only possible when the external field is weak, H < H c . It follows from expression (9) for the local fields that for small K the local fields H 1 , 2 > 0 and, therefore, when m 1 , 2 > 0 , both sub-groups are in the ferromagnetic phase. If H < H c , when m 1 m 2 < 0 , there is a temperature K p ( H ) , at which point the transition into the antiferromagnetic phase occurs. If q 12 q 12 * , there are four fundamentally different behaviors of partial magnetizations as functions of temperature, m 1 , 2 ( K ) , as determined by the characteristics of the model and by the strength of the external field H (Figure 5 and Figure 6).
Before we proceed to further analysis, let us consider the behavior of partial magnetizations at high temperatures, i.e., at K 0 . Differentiating Equation (6) with respect to K and retaining only the quantities of the first order of smallness, we obtain the following expressions for the derivatives m ˙ 1 and m ˙ 2 at K 0 :
m ˙ 1 > m ˙ 2 i f p 1 ( q 11 + | q 12 | ) > p 2 ( q 22 + | q 12 | ) , m ˙ 2 > m ˙ 1 i f p 1 ( q 11 + | q 12 | ) < p 2 ( q 22 + | q 12 | ) .
This implies that with growing K , the partial magnetization m 1 increases faster than m 2 if the effective number of neighbors in the lattice I is greater than that in the lattice II, and vice versa. For definiteness, in this subsection, we will assume that p 1 ( q 11 + | q 12 | ) > p 2 ( q 22 + | q 12 | ) . If the opposite is true, all considerations will hold up to a permutation of the indices 1 and 2.
1.
The case H < H c . Let us first consider the case H < H c , when the ground state is antiferromagnetic ( m 1 m 2 = 1 at K ).
Let p 1 ( q 11 + | q 12 | ) > p 2 ( q 22 + | q 12 | ) . Then, in accordance with (27), for K < < 1 , the magnetization m 1 grows with K faster than m 2 . In this case, the magnetization m 1 will keep growing as K increases, while the magnetization m 2 will initially grow and then will drop to negative values, passing through a zero value at a particular temperature K = K p . At further growth of K , there are two distinct scenarios for the curve development in the K > K p region:
  • If p 1 > p 2 , the partial magnetization m 1 will monotonously increase, and the partial magnetization m 2 will decrease, tending to 1 and −1, correspondingly (Figure 5a). There is no phase transition in this case.
  • If p 1 < p 2 , the system’s ground state is the configuration with m 1 = 1 , m 2 = 1 . At a certain temperature, a first-order phase transition will occur, as a result of which the magnetization m 1 will become negative, and the magnetization m 2 will be positive (Figure 5b). This transition will be accompanied by an abrupt change in the measurable parameters: magnetization M ( K ) , internal energy U ( K ) , and heat capacity C ( K ) .
Note that the full magnetization of the system M = p 1 m 1 + p 2 m 2 always remains positive, and in the limit K , we have M = | p 1 p 2 | independently of the relationship between p 1 and p 2 .
The temperature K p of the phase change, i.e., the temperature at which m 2 passes zero, can be easily found: at m 2 = 0 , the second part of Equation (6) reduces to | q 12 | p 1 m 1 = H . Substituting this relation in the first part of Equation (6), we obtain:
K p = | q 12 | 2 H ( q 11 + | q 12 | ) ln p 1 | q 12 | + H p 1 | q 12 | H ,   m 2 ( K = K p ) = 0 ,   m 1 ( K = K p ) = H | q 12 | p 1 .
When the external field is weak ( H p 1 | q 12 | ), the temperature of the transition into the antiferromagnetic phase approaches the quantity
K p ( H 0 ) = 1 p 1 ( q 11 + | q 12 | )
Let us compare (29) with the temperature K c . It is easy to see that inequality K p ( H 0 ) K c is identical to condition (27), i.e., when the field is weak, the temperature K p always tends to a value less or equal to K c . The two temperatures can be equal, K p ( H 0 ) = K c , only if q 12 = q 12 * .
Moreover, the expression (28) is correct given that p 1 | q 12 | > H . It means that m 2  can pass zero only if the external field is relatively weak, or if the transition into the antiferromagnetic phase is an abrupt process. Otherwise, m 2 cannot pass zero, and m 1 , 2 ( K ) 0 for any K .
2.
If H > H c , the behavior of the temperature dependences is quite predictable, without any peculiarities (see Figure 6). The partial magnetizations are positive at any K . In particular, if the field H is close enough to the critical value H c , the magnetization m 2 decreases with the growth of K in a certain temperature range, but it never becomes negative (Figure 6a). When H H c , the partial magnetizations m 1 and m 2 are increasing functions of K over the whole temperature range. When H > H c , the ground state of the system is the ferromagnetic state: M = m 1 = m 2 = 1 , K .
Summarizing the results of Section 4.1, we can highlight the following properties of an unbalanced ( q 12 q 12 * ) system: Firstly, the presence of magnetic field suppresses a second-order phase transition in the system. Secondly, if 0 < H < H c , a first-order phase transition can occur for a certain relation between parameters (27) and p 1 , 2 , as illustrated in Figure 5b.

4.2. Balanced System ( q 12 = q 12 * )

If the condition of equal effective numbers of neighbors (12) is met, there is a solution to the equations of state (6): m 1 = m 2 . This solution does not always correspond to the free energy minimum, and the behavior of partial magnetizations is so diverse that the analysis of the general case ( q 11 q 22 , p 1 p 2 ) is not within the scope of this paper.
Here, we consider the simplest (symmetrical) case, when
q 11 = q 22 = q   and   p 1 = p 2 = 1 2 .
Note that even in this simplest case, the behavior of partial magnetizations can be fairly diverse (Figure 7). Although the below formulae are true only for symmetric systems (30), general patterns hold true for all systems with q 12 = q 12 * .
We recall that in the symmetric model, the equality q 12 = q 12 * is true for all negative q 12 . The critical temperature in this case is determined as:
K c = 2 q + | q 12 | ,
and equations of state (6) take the form:
1 K ln 1 + m 1 1 m 1 = q m 1 | q 12 | m 2 + 2 H , 1 K ln 1 + m 2 1 m 2 = q m 2 | q 12 | m 1 + 2 H .
As we see, these equations have the same coefficients and admit a solution m 1 = m 2 . Indeed, analysis of (31) shows that over a certain temperature interval 0 K K S , partial magnetizations m 1 = m 2 = m are defined by the equation:
1 K ln 1 + m 1 m = m ( q | q 12 | ) + 2 H
Magnetization m is positive because the external magnetic field H is positive. If the field is smaller than a certain value H s max , at a temperature K = K s , the curve m = m ( K ) experiences a “soft” splitting into two diverging partial magnetizations m 1 ( K ) and m 2 ( K ) (Figure 7a).
  • Soft splitting and merging. Let us determine the “soft” splitting point K S . Quantities m 1 = m 2 = m 3 at this point are derived from (32) at K = K S . Let us consider the behavior of magnetizations m 1 , 2 in a small vicinity of K S , introducing the small deviation
    t s = T T s T s = K s K K ,   | t s | 0 .
In this case, we seek partial magnetizations in the form m 1 = m 2 = m S δ when K < K s and in the form m 1 , 2 = m S ± δ 1 , 2 when K > K s , where δ 0 and δ 1 , 2 0 . Let us substitute m 1 , 2 of this form into (31) and (32) and carry out expansion in small parameters δ , δ 1 , 2 retaining terms up to δ 1 , 2 3 . Equating the terms of the same order of smallness, we obtain the expression for the splitting point:
K s = K c 1 m s 2 ,
as well as the expressions for small deviations from m S :
δ = t s κ , δ 1 , 2 = t s κ 1 ± κ 2 t s ,
where
κ = ( q | q 12 | ) m s + 2 H 2 | q 12 | κ 1 = 3 ( 1 m s 2 ) ( | q 12 | q m s 2 2 H m s | q 12 | 3 q m s 2 ) κ 2 = ( 1 + 3 m s 2 ) κ 1 3 ( 1 m s 2 ) 2 6 m s ( 1 m s 2 )
2.
Splitting point and “bubble” formation. The quantity m S is a solution to Equation (32) at K = K S . This equation allows for a simple graphical solution if we rewrite it in the form:
H = R ( m s ) 1 2 m s ( | q 12 | q ) + 1 4 ( | q 12 | + q ) ( 1 m s 2 ) ln 1 + m s 1 m s
One example of such a solution for | q 12 | < q , | q 12 | = q and | q 12 | > q is given in Figure 8. The function R ( m s ) has one maximum at m s = m max , where m max  is the solution to the equation R ( m ) / m = 0 :
m max ln 1 + m max 1 m max = 2 | q 12 | q + | q 12 | .
Correspondingly, Equation (36) can have a solution only when H H max , where
H max = ( | q 12 | q m max 2 ) 2 m max .
If | q 12 | = q , the right wing of the curve R ( m s ) is on the abscissa axis. In the case | q 12 | < q , the curve is lower than the abscissa. As one can see, if | q 12 | q , Equation (36) may have two solutions. The first one ( m s m max ) corresponds to the splitting point of the curves m 1 ( K ) and m 2 ( K ) , the second solution ( m s m max ) corresponds to the merging point, and the formation of a “bubble” is shown in Figure 7b,c.
If | q 12 | > q , the right wing of the curve R ( m s ) behaves as R = ( | q 12 | q ) / 2 and is above the abscissa. Therefore, for 0 < H ( | q 12 | q ) / 2 , Equation (36) has only one solution ( m s m max ), which is responsible for the splitting, whereas for ( | q 12 | q ) / 2 < H H max , there are two solutions m s m max and m s m max , i.e., a “bubble” can form (as in Figure 7b,c).
Let us determine the asymptotic values of m max and H max . In the case of weak antiferromagnetic interaction ( | q 12 | < < q ), it follows from (37) and (38) that
m max | q 12 q | 1 / 2 ,   H max | q 12 | 3 / 2 2 q ,
and in the limit case | q 12 | > > q , we have:
m max 0.8336 ,   H max = | q 12 | 2 m max 1.2 H c .
As we see, a soft splitting can be observed when H is varied within a wide enough range. Note that the width of this range grows with | q 12 | .
If soft splitting is not feasible and the antiferromagnetic state ( H < H c ) is the ground state, then passing through the critical point is accompanied by a jump of partial magnetizations (Figure 7b). If the condition of a “bubble” formation is met, a jump in magnetization occurs after the merging point at H < H c (Figure 7c).
Generally, the temperature behavior of magnetizations is very diverse (Figure 7), depending strongly on the relative magnitudes of the ratio of | q 12 | / q and H . This aspect of the model will be investigated in more detail in future papers.
3.
Field-controlled phase transition at K = K s . A second-order phase transition occurs at the point of soft splitting/merging and is accompanied by a jump in heat capacity. The magnitude of the heat capacity jump, ΔC under the splitting can be easily calculated by differentiating expression (3) for the energy, with the account of expression (35):
Δ C = 3 ( q + | q 12 | ) ( | q 12 | q m s 2 2 H m s ) 2 4 | q 12 | ( | q 12 | 3 q m s 2 ) .
At the point of merging of the curves m 1 ( K ) and m 2 ( K ) , Δ C is determined by the same expression, but with a different value of m s .
It follows from what has been said that K s determined by (33) is nothing else but a magnetic-field-dependent critical point: K s = K c ( H ) . In other words, in a symmetrical system, the external field does not suppress a phase transition but only shifts it towards larger temperatures: K c ( H 0 ) > K c ( H = 0 ) K c . The range of the critical point values can be quite large. As follows from (33), given strong antiferromagnetic interaction ( | q 12 | > > q ) the phase transition at the soft-splitting point can be greatly shifted: with the field changing from H = 0 to H = H max ~ 1.2 H c ( m s changing from m s = 0 to m s = m max ~ 0.83 ), K s changes from K s = K c to K s 3.3 K c . The variations at the merging point (if any) cover a yet wider range, where the case of m s 1 , (i.e., Ks →∞) can take place. For instance, Figure 7f demonstrates the case when, with H 0.999 H c , a “bubble” forms, and the system has two critical points: one at the splitting point at K s 1.4 K c and the second at the merging point at K s 8.3 K c . We have deliberately chosen such interaction parameters that reveal the richness of the system: besides two second-order phase transitions, it has two first-order phase transitions (the curve m 2 ( K ) undergoes two abrupt changes).
4.
Critical exponents of a phase transition at the point of soft splitting. Let us consider the critical behavior of the involved physical quantities near the temperature K s .
We start by pointing out that the exponent α s = 0 , since, in accordance with Equation (41), at the splitting point the heat capacity experiences a finite jump. We will define other critical exponents as follows:
m s M ( K ) ( t s ) β s   with   t s 0 ,
m s M ( K ) t s β s   with   t s 0 + ,
χ ( K ) ( t s ) γ s   with   t s 0 ,
χ ( K ) t s γ s   with   t s 0 + ,
M ( K s , H ) m s ( H H s ) 1 / δ s   with   H H s 0 + ,
m s M ( K s , H ) ( H s H ) 1 / δ s   with   H H s 0 ,
where H s is the magnitude of the external field at which the phase transition happens.
It follows from Equation (34) that in the vicinity of the splitting point, the full magnetisation is determined by the following equations:
M = m s + κ 2 t s ,   t s < 0 , M = m s κ t s ,   t s > 0 .
It follows from (42) that the critical exponents β s = β s = 1 . For H 0 ( m s 0 ), we have κ 2 0 , and hence we should take into account the terms of the order of t s 3 / 2 in the expansion of the full magnetization M over t s . This agrees with the results obtained earlier in Section 3.1. We also point out that at large values of | q 12 | , there is such a value of the magnetic field H < H max beyond which κ 2 becomes negative, so that the magnetization increases when the temperature rises above the splitting point.
We differentiate the equations of states (31) and (32) over H and, making use of Equations (33) and (34), we obtain the following for the susceptibilities at the critical points:
χ = 1 | q 12 | , t s 0 + χ = [ | q 12 | + 2 κ 1 m s 2 ( q + | q 12 | ) ( 1 m s 2 + 2 m s ) ( 1 m s 2 ) ( 3 m s 2 + 1 ) κ 1 ] 1 , t s 0
It follows from (43) that the susceptibility has a finite jump in the critical point. Its magnitude depends on m s , and hence on the magnitude of the field H . In the absence of the external field, the jump disappears. From Equation (43), we find the critical exponents γ s = γ s = 0 .
We will use Equation (32) to estimate the exponent δ s . Let us look at how the magnetization changes in response to a small increase in the external field. We introduce the notations Δ M = M m s and Δ H = H H s , and accounting only for the first-order terms in Δ M , we find at the temperature K s :
Δ H = | q 12 | Δ M .
Equation (44) is valid also for negative values of Δ H and Δ M ; therefore, we obtain for the critical exponents: δ s = δ s = 1 .
We see that in this case, the relationships α + 2 β + γ = 2 , γ = β ( δ 1 ) and β ( δ + 1 ) = 2 α , which are the consequences of the scaling hypothesis, are satisfied in the form of equalities. Except for β , all other critical exponents at the soft splitting point are equal to their values in the absence of the external field.
5.
Restrictions in the soft splitting. Above we analyzed the processes of soft splitting and merging with the assumption that these regimes could be reached. However, sometimes this is not possible.
Firstly, the field H should be such that the condition K s < K p is met, i.e., the curves m 1 , 2 ( K ) must split before one of them passes zero.
Secondly, the field cannot be too big, H < H max , otherwise the soft splitting will be replaced by an abrupt jump-like change in partial magnetizations m 1 , 2 ( K ) .
Thirdly, the splitting at point K s is possible if δ 1 and δ 2 take real values. It follows from (35) that δ 1 , δ 2 at t s < 0 , if
( | q 12 | q m s 2 2 H m s ) ( | q 12 | 3 q m s 2 ) > 0 .
Expression (45) is a necessary condition for the soft splitting point to exist. The quantities δ 1 and δ 2 must be real numbers at t s > 0 for K s to be the point of merging. Correspondingly, the necessary condition for the point of merging to exist is as follows:
( | q 12 | q m s 2 2 H m s ) ( | q 12 | 3 q m s 2 ) < 0 .
Examination of (46) in light of (36)–(40) shows that merging of the curves m 1 ( K ) and m 2 ( K ) (i.e., the formation of a “bubble”) usually takes place only when the magnitude of H is very close to H c .

5. Computer Simulation of the Layered Model

Let us compare the results of the mean-field model with the results of the computer simulation of lattices with a finite interaction radius. We want to make sure that in the case of “critical” antiferromagnetic interaction q 12 = q 12 * , which takes place if the condition (12) is satisfied, the critical parameters indeed take non-classical values presented in Table 2. To this end, let us consider a three-dimensional cubic lattice consisting of alternating two-dimensional layers of spins (see Figure 9). The intra-layer spin interaction is described by the interaction constants J 11 > 0  and J 22 > 0  for even and odd layers, respectively. The interlayer interaction is antiferromagnetic, and it is described by the constant J 12 < 0 . The interaction only between the nearest neighbors is taken into account.
This model is described by the Hamiltonian
E N = ( J 11 i , j S i S j + J 22 i , j S i S j + J 12 i , j S i S j ) H ( i = 1 N 1 S i + j = 1 N 2 S j ) ,
where i , j stands for a set of nearest neighbors. The relationship between the interaction parameters of the layered model and the parameters of the mean-field model is defined by the expression:
q 11 = 8 J 11 ,   q 22 = 8 J 22 ,   q 12 = 4 J 12 .
Condition (12), under which q 12 = q 12 * , is set for the simplest case of p 1 = p 2 , q 11 = q 22 .
The Metropolis algorithm is used to calculate the thermodynamic parameters for the lattices of the size N = L × L × L , where L is the dimension of the lattice varying from 6 to 64. The linear dimension L is always an even number to secure the equality p 1 = p 2 . The case of p 1 p 2 was also modeled preliminarily, but the results of the modeling showed that in this case, it is impossible to achieve the equality of the moduli of partial magnetizations, | m 1 | = | m 2 | , over a wide temperature range even if condition (12), the effective neighbor numbers being equal, is met. A possible cause for this is the effect of boundaries. Therefore, to make it possible to study “critical” antiferromagnetic interactions occurring in the balanced model ( q 12 = q 12 * ), the modeling is carried out for the case of p 1 = p 2 only.
1.
Figure 10 shows the temperature dependences of partial magnetizations of the layered model when q 11 = q 22 . The general behavior of these dependences is similar to the dependences obtained in the mean-field model for a balanced system (c). Depending on the magnitude of the external field, we can observe both “soft” (Figure 10b) and “hard” (Figure 10c) splitting. However, we failed to observe magnetization “bubbles” in this model.
When the lattice dimension is big enough and the external field H is close to, but less than, the critical magnitude H c , the system does not proceed into the antiferromagnetic phase by cooling; instead, it stays in the local energy minimum (Figure 10d). We observe the same situation for a fully connected lattice (the analog of the mean-field model) if we use the Metropolis algorithm to compute the thermodynamic parameters.
To avoid misunderstanding, we should add that “hard” splitting accompanied by an abrupt change in magnetization (shown in Figure 10c) was observed only with relatively small lattice dimensions ( L 24 ). A small lattice size makes it possible to overcome the energy barrier between the local and global minima. When dimensions are large ( L = 32 , 64 ) and K , the system always transits into the ferromagnetic state with the magnetization M = 1 , which is a local minimum at H < H c .
2.
Let us consider the temperature dependence of the susceptibility χ = χ ( K ) in the layered model in the absence of the external field: H = 0 . If q 11 = q 22 , i.e., q 12 = q 12 * , the form of the curve χ ( K ) is independent of the lattice dimension, and the susceptibility takes a finite value at the critical point (Figure 11a,b). With large | q 12 | , the curve χ ( K ) does not have a maximum (Figure 11a). It means that the critical exponents γ = γ = 0 , which agrees with the numbers in Table 2.
If q 11 q 22 , the susceptibility peak, as expected, is observed at the critical temperature (Figure 11c,d). The height of the peak increases with the lattice dimensions L , which means that the susceptibility diverges at the critical point. This finding agrees with the classical critical exponents γ = γ = 1 .
3.
To obtain the value of the critical exponent δ, we need to evaluate the critical temperature K c . The temperature dependences of the Binder cumulants [29] for lattices of different dimensions were built for this purpose. The asymptotic value of K c for the case of L is determined as a point of intersection of these dependences. We define the Binder cumulants for partial and full magnetization in the following way:
g 1 = 1 m 1 4 3 m 1 2 2 ,   g 2 = 1 m 2 4 3 m 2 2 2 ,   g = 1 M 4 3 M 2 2 .
When q 11 = q 22 , in the absence of the external field the full magnetization M is zero. For this reason, we had to use partial magnetization cumulants to determine the critical temperature (Figure 12a). The results obtained with the help of the cumulants g 1 и g 2 gave approximately the same value of K c . This allowed us to determine the exponent δ = 1 to within 10−3 (Table 3).
If q 11 q 22 , the cumulants g 1 and g 2 give different estimates of K c . At the same time, the cumulants g for the full magnetization do not intersect at one point (Figure 12b). For this reason, we evaluated K c roughly as the middle of the interval where the cumulant g experiences a jump. As can be seen from Table 3, such inaccuracy in the determination of K c results in large deviations of the exponent δ .
To evaluate the exponent δ , we built relationships between ln M and ln H (Figure 13).
If q 11 = q 22 , the relationships between ln M and ln H for lattices of different dimensions L merge into one line already at a sufficiently small H (see Figure 13a). The slope of this relationship was used to evaluate δ . The exponent δ was found to be independent of the ratio q / | q 12 | , being roughly equal to one (see Table 3). Its deviations from one are about 10−3, which is comparable to the error of our evaluation. Thus, the value of the exponent δ in this case fully agrees with its value δ = 1 in the mean-field model for q 12 = q 12 * .
If q 11 q 22 , the curves ln M as functions of ln H become linear at certain H , whose magnitude decreases with growing L (see Figure 13b). The parameter δ measured by the curve slope differs considerably from the expected value δ = 3 in all the investigated configurations (Table 3), and depends on the parameters of the model.
4.
Let us consider how the critical exponents change with the increase in the external field in a balanced layered model. We estimated the splitting temperature K s by plotting the temperature dependencies of the following cumulants:
g s = 1 ( m 1 m 2 ) 4 3 ( m 1 m 2 ) 2 2 .
The temperature K s was determined as the point of intersection of the cumulants for lattices with various dimensions (Figure 14). We point out that the cumulants g s are usually used for determining the critical temperature of an antiferromagnetic phase transition [30]. The magnetization at the splitting point, m s , was calculated as the magnitude of the local maximum of magnetization for the lattice of the maximum size ( L = 32 ).
To determine the critical exponents β s and β s , we plotted the dependencies of ln ( m s M ) on ln | t s | for the lattices with various dimensions (Figure 15). These curves become linear when t s < 0 , and their slope allows us to estimate the exponent β s (see Table 4). For t s > 0 , the dependencies of ln ( m s M ) on ln t s are nonlinear for the whole range of t s , making it impossible to estimate the value of the exponent β s .
The temperature dependencies of the susceptibility χ near the splitting point are shown in Figure 16. As one can see, when the external field is present, the susceptibility curves depend on the lattice size. The accuracy of our simulations is insufficient to unambiguously attribute this observation to either a divergence of the susceptibility at the temperature K s , or to its finite jump.
In order to determine the critical exponents δ s and δ s , we plotted the dependencies of ln | Δ M | on ln | Δ H | for the lattices with the dimensions ranging from L = 12 to L = 32 (Figure 17). Independently of the sign of Δ H , these curves have a linear asymptotic, the slope of which allowed us to extract the critical exponents (Table 4).
The results of these simulations are shown in Table 4. The splitting temperature K s differs substantially from the one obtained in the framework of the mean-field model (MFM). The corresponding numbers are presented in column 2 of Table 4. However, when in Equation (33) we use the values of K c and m s obtained as a result of this modelling, we obtain K s rather accurately (compare three numbers in the column 5 of Table 4). The obtained critical exponent of the layered model β s takes values in a range of 0.5–0.6, which differs significantly from its value from the mean-field model ( β s = 1 ). The value of the exponent β s decreases with the increase in the external field. However, similarly to the mean-field model, the value of β s exceeds considerably the values of the exponents of partial magnetization, β m 1 , 2 0.3 , measured in the absence of the external field (in the mean-field model, β m 1 , 2 = 1 / 2 ). For small values of the external field, the values of the exponents δ s and δ s are approximately equal to 1. With the increase in the external field, the exponent δ s increases, while δ s decreases.

6. Results and Discussions

We have studied the spin system consisting of two antiferromagnetically interacting sub-lattices. We obtained our results using the mean-field model, but we are not able to assert that they hold for the general case. However, the Monte Carlo simulation for the 3D layered model with a finite interaction radius confirms our conclusions. Our results are as follows:
  • In the absence of an external field, the system has a second-order phase transition. If the relationship between the system’s parameters is such that the effective numbers of neighbors in the sub-lattices are equal (12), then the critical exponents of the transition differ significantly from the usual values of the mean-field model (see Table 2), and the susceptibility has a finite value at the critical point. It should be mentioned that the concept of the universality class is still valid because the critical exponents for each sub-lattice remain the same, and the unusual critical behavior takes place only for total values of magnetization and susceptibility.
    • These conclusions are validated by computer simulations of a three-dimensional layered lattice. In particular, the balanced system ( q 12 = q 12 * ) is shown to have the following critical exponents: γ = γ = 0 and δ = 1 . The dependences m1 = m1(K), m2 = m2(K) and M = M ( K ) (Figure 10) generated in the computer simulation reproduce the theoretical curves (Figure 6 and Figure 7).
    • If the condition of equality of the effective number of neighbors (12) is met, the external field does not destroy the phase transition, but only shifts the transition temperature. If the external field is weak enough, H H max , where H max is defined by (38), the shift of the critical point does not change the qualitative characteristics of the system, and the transition remains a second-order phase transition. In Section 4.2, we show that it is possible to vary the critical temperature over a sufficiently wide range (from K c to 3.3 K c ) by varying H within the interval 0 H H max . Under these conditions, the values of the critical exponents α , γ and δ do not change in the mean-field model, and the exponent β becomes twice as large as its value for the partial magnetizations. In addition, for a finite external field, the system exhibits a jump in the susceptibility at the phase transition point.
The results of the numerical modeling on a three-dimensional layered lattice demonstrate that the critical exponents β and δ do change when the external field increases (see Table 4). Nevertheless, the results of the numerical modeling coincide with the predictions of the mean-field model in two aspects: (1) the exponent β for the full magnetization considerably exceeds its magnitude for the partial magnetization, and (2) the presence of an external field results in a jump of the susceptibility.
If H > H max , the shift of the critical point is accompanied by qualitative changes: at the critical point, a first-order phase transition replaces the second-order phase transition, which is accompanied by jumps in the magnetization and the internal energy (Figure 7d).
Since this system has a lot of free parameters, the overall picture of possible regimes is very diverse. Given a certain relation between the parameters of the model and the external field, there may be two first-order phase transitions and two second-order phase transitions (Figure 7f). The general behavior of the temperature dependences of the involved physical quantities, as well as determining the conditions under which various regimes occur, is a complicated problem that deserves to be considered in a separate paper.
2.
The behavior of the system depends on the relation between the magnitude of the magnetic field and the critical value H c (11). When H < H c , the system’s ground state is an antiferromagnetic state with m 1 = m 2 , | m 1 | = | m 2 | = 1 and the magnetization M = | p 1 p 2 | . When H > H c , the ground state is a ferromagnetic state with M = m 1 = m 2 = 1 (Table 1). Correspondingly, adiabatic cooling of the system ( K ) transfers it into the antiferromagnetic state at H < H c and into the ferromagnetic state at H > H c . The analysis given in Section 4 shows that (except for some special cases) the field dependence of the magnetization M = M ( H ) at a fixed temperature follows the following pattern:
  • If K is smaller than or of the order of K c , then the magnetization M grows monotonously with H , and becomes equal to one when the system reaches the ferromagnetic state at H H c ;
  • If K is several times greater than K c , the picture is different: for 0 < H < H c there is a plateau of a height M = | p 1 p 2 | , which obtains transferred to the ferromagnetic limit M = 1 at H = H c by means of a jump.
The plateau on the curve M = M ( H ) is shown in Figure 18, where we vary the quantity | p 1 p 2 | while keeping the other parameters constant. It is seen that the increase in | p 1 p 2 | results in the growth of the height of the plateau, as well as in the increase in the field magnitude at which the jump to the state M = 1 takes place. Note that the curves in Figure 18 resemble the well-known dependence M = M ( H ) demonstrated by many materials with giant magnetoresistance [31,32,33].

Author Contributions

Conceptualization, L.L. and B.K.; methodology, B.K., V.E. and L.L.; software, V.E.; validation, V.E. and L.L.; formal analysis, B.K. and L.L.; investigation, V.E.; writing—original draft preparation, L.L. and B.K.; writing—review and editing, B.K. and V.E.; funding acquisition, B.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by State Program of Scientific Research Institute for System Analysis, Russian Academy of Sciences, project no. FNEF-2022-0003.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are openly available in FigShare at https://doi.org/10.6084/m9.figshare.24218805.v1.

Acknowledgments

The authors are grateful to Marina Litinskaya for help when preparing this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ilkovič, V. Magnetic Properties of Ising-Type Ferromagnetic Films with a Sandwich Structure. Phys. St. Sol. 1998, 207, 131–137. [Google Scholar] [CrossRef]
  2. Kuzniak-Glanowska, E.; Konieczny, P.; Pełka, R.; Muzioł, T.M.; Kozieł, M.; Podgajny, R. Engineering of the XY Magnetic Layered System with Adeninium Cations: Monocrystalline Angle-Resolved Studies of Nonlinear Magnetic Susceptibility. Inorg. Chem. 2021, 60, 10186–10198. [Google Scholar] [CrossRef] [PubMed]
  3. Grünberg, P. Layered magnetic structures in research and application. Acta Mater. 2000, 48, 239–251. [Google Scholar] [CrossRef]
  4. Taborelli, M.; Allenspach, R.; Boffa, G.; Landolt, M. Magnetic coupling of surface adlayers: Gd on Fe (100). Phys. Rev. Lett. 1986, 56, 2869. [Google Scholar] [CrossRef]
  5. Camley, R.E.; Tilley, D.R. Phase transitions in magnetic superlattices. Phys. Rev. B 1988, 37, 3413. [Google Scholar] [CrossRef]
  6. Camley, R.E. Properties of magnetic superlattices with antiferromagnetic interfacial coupling: Magnetization, susceptibility, and compensation points. Phys. Rev. B 1989, 39, 12316. [Google Scholar] [CrossRef]
  7. Lipowski, A. Critical temperature in the two-layered Ising model. Phys. A 1998, 250, 373–383. [Google Scholar] [CrossRef]
  8. Horiguchi, T.; Lipowski, A.; Tsushima, N. Spin-32 ising model and two-layer Ising model. Phys. A 1996, 224, 626–638. [Google Scholar] [CrossRef]
  9. Diaz, I.J.L.; da Silva Branco, N. Monte Carlo simulations of an Ising bilayer with non-equivalent planes. Phys. A Stat. Mech. Appl. 2017, 468, 158–170. [Google Scholar] [CrossRef]
  10. Diaz, I.; Branco, N. Monte Carlo study of an anisotropic Ising multilayer with antiferromagnetic interlayer couplings. Phys. A 2018, 490, 904–917. [Google Scholar] [CrossRef]
  11. Gharaibeh, M.; Obeidat, A.; Qaseer, M.K.; Badarneh, M. Compensation and critical behavior of Ising mixed spin (1-1/2-1) three layers system of cubic structure. Phys. A 2020, 550, 124147. [Google Scholar] [CrossRef]
  12. Wang, W.; Xue, F.-L.; Wang, M.-Z. Compensation behavior and magnetic properties of a ferrimagnetic mixed-spin (1/2, 1) Ising double layer superlattice. Phys. B 2017, 515, 104–111. [Google Scholar] [CrossRef]
  13. Kaneyoshi, T.; Jascur, M. Compensation temperatures of ferrimagnetic bilayer systems. J. Magn. Magn. Mats. 1993, 118, 17–27. [Google Scholar] [CrossRef]
  14. Szałowski, K.; Balcerzak, T. Normal and inverse magnetocaloric effect in magnetic multilayers with antiferromagnetic interlayer coupling. J. Phys. Cond. Mat. 2014, 26, 386003. [Google Scholar] [CrossRef]
  15. Balcerzak, T.; Szałowski, K. Ferrimagnetism in the Heisenberg–Ising bilayer with magnetically non-equivalent planes. Phys. A 2014, 395, 183–192. [Google Scholar] [CrossRef]
  16. Campa, A.; Gori, G.; Hovhannisyan, V.; Ruffo, S.; Trombettoni, A. Ising chains with competing interactions in the presence of long-range couplings. J. Phys. A 2019, 52, 344002. [Google Scholar] [CrossRef]
  17. Giraldo-Barreto, J.; Restrepo, J. Tax evasion study in a society realized as a diluted Ising model with competing interactions. Phys. A 2021, 582, 126264. [Google Scholar] [CrossRef]
  18. Jin, H.-K.; Pizzi, A.; Knolle, J. Prethermal nematic order and staircase heating in a driven frustrated Ising magnet with dipolar interactions. Phys. Rev. B 2022, 106, 144312. [Google Scholar] [CrossRef]
  19. Stariolo, D.A.; Cugliandolo, L.F. Barriers, trapping times, and overlaps between local minima in the dynamics of the disordered Ising p-spin model. Phys. Rev. E 2020, 102, 022126. [Google Scholar] [CrossRef]
  20. Zhang, Q.; Wei, G.-Z.; Liang, Y.-Q. Phase diagrams and tricritical behavior in spin-1 Ising model with biaxial crystal-field on honeycomb lattice. J. Magn. Magn. Mat. 2002, 253, 45–50. [Google Scholar] [CrossRef]
  21. Ez-Zahraouy, H.; Kassou-Ou-Ali, A. Phase diagrams of the spin-1 Blume-Capel film with an alternating crystal field. Phys. Rev. B 2004, 69, 064415. [Google Scholar] [CrossRef]
  22. Zivieri, R. Critical behavior of the classical spin-1 Ising model for magnetic systems. AIP Adv. 2022, 12, 035326. [Google Scholar] [CrossRef]
  23. Kaneyoshi, T. Contribution to the theory of spin-1 Ising models. J. Phys. Soc. Jpn. 1987, 56, 933–941. [Google Scholar] [CrossRef]
  24. Kaneyoshi, T. The phase transition of the spin-one Ising model with a random crystal field. J. Phys. C Solid State Phys. 1988, 21, L679–L682. [Google Scholar] [CrossRef]
  25. Kryzhanovsky, B.; Litinskii, L.; Egorov, V. Analytical Expressions for Ising Models on High Dimensional Lattices. Entropy 2021, 23, 1665. [Google Scholar] [CrossRef]
  26. Hansen, P.L.; Lemmich, J.; Ipsen, J.H.; Mouritsen, O.G. Two coupled Ising planes: Phase diagram and interplanar force. J. Stat. Phys. 1993, 73, 723–749. [Google Scholar] [CrossRef]
  27. Ferrenberg, A.M.; Landau, D.P. Monte Carlo study of phase transitions in ferromagnetic bilayers. J. Appl. Phys. 1991, 70, 6215–6217. [Google Scholar] [CrossRef]
  28. Baxter, R. Exactly Solved Models in Statistical Mechanics; Dover Books on Physics; Dover Publications: Mineola, NY, USA, 2007. [Google Scholar]
  29. Binder, K. Finite size scaling analysis of Ising model block distribution functions. Z. Phys. B 1981, 43, 119–140. [Google Scholar] [CrossRef]
  30. Godoy, M.; Wagner, F. Nonequilibrium antiferromagnetic mixed-spin Ising model. Phys. Rev. E 2002, 66, 036131. [Google Scholar] [CrossRef]
  31. Drovosekov, A.B.; Kholin, D.I.; Kreinies, N.M. Magnetic Properties of Layered Ferrimagnetic Structures Based on Gd and Transition 3d Metals. J. Exp. Theor. Phys. 2020, 131, 149–159. [Google Scholar] [CrossRef]
  32. Telford, E.J.; Dismukes, A.H.; Lee, K.; Cheng, M.; Wieteska, A.; Bartholomew, A.K.; Chen, Y.; Xu, X.; Pasupathy, A.N.; Zhu, X.; et al. Layered antiferromagnetism induces large negative magnetoresistance in the van der Waals semiconductor CrSBr. Adv. Mats. 2020, 32, 2003240. [Google Scholar] [CrossRef] [PubMed]
  33. Wang, Y.; Yu, J.H.; Xi, C.Y.; Ling, L.S.; Zhang, S.L.; Wang, J.R.; Xiong, Y.M.; Han, T.; Han, H.; Yang, J.; et al. Topological semimetal state and field-induced Fermi surface reconstruction in the antiferromagnetic monopnictide NdSb. Phys. Rev. B 2018, 97, 115133. [Google Scholar] [CrossRef]
Figure 1. Energies of configurations S + + , S + , S + and S as a function of the strength of the field, H : (a) the system does not meet conditions (10); (b) the system meets conditions (10). The solid lines mark the range of H , where a configuration is a minimum, the broken lines mark the range where a configuration is not a minimum. Here, we use the following notations: H + + = max ( | q 12 | p 1 q 22 p 2 , | q 12 | p 2 q 11 p 1 ) , H + = q 22 p 2 + | q 12 | p 1 , H = min ( q 22 p 2 | q 12 | p 1 , q 11 p 1 | q 12 | p 2 ) and H + = q 11 p 1 + | q 12 | p 2 .
Figure 1. Energies of configurations S + + , S + , S + and S as a function of the strength of the field, H : (a) the system does not meet conditions (10); (b) the system meets conditions (10). The solid lines mark the range of H , where a configuration is a minimum, the broken lines mark the range where a configuration is not a minimum. Here, we use the following notations: H + + = max ( | q 12 | p 1 q 22 p 2 , | q 12 | p 2 q 11 p 1 ) , H + = q 22 p 2 + | q 12 | p 1 , H = min ( q 22 p 2 | q 12 | p 1 , q 11 p 1 | q 12 | p 2 ) and H + = q 11 p 1 + | q 12 | p 2 .
Entropy 25 01428 g001
Figure 2. Condition q 12 * < 0 is met for parameters lying in the painted area, i.e., between the curve p 1 = 1 / 2 and the curve p 1 = q 22 / ( q 11 + q 22 ) .
Figure 2. Condition q 12 * < 0 is met for parameters lying in the painted area, i.e., between the curve p 1 = 1 / 2 and the curve p 1 = q 22 / ( q 11 + q 22 ) .
Entropy 25 01428 g002
Figure 3. Typical curves for H = 0 : (a) partial magnetizations m 1 , 2 = m 1 , 2 ( K ) and full magnetization M = M ( K ) ; (b) heat capacity C = C ( K ) . The curves are drawn for q 11 = 2 ,   q 22 = 1 ,   q 12 = 1 ,   p 1 = 0.6 ,   p 2 = 0.4 .
Figure 3. Typical curves for H = 0 : (a) partial magnetizations m 1 , 2 = m 1 , 2 ( K ) and full magnetization M = M ( K ) ; (b) heat capacity C = C ( K ) . The curves are drawn for q 11 = 2 ,   q 22 = 1 ,   q 12 = 1 ,   p 1 = 0.6 ,   p 2 = 0.4 .
Entropy 25 01428 g003
Figure 4. Changes in functions of physical quantities at the “critical” antiferromagnetic interaction: (a) susceptibility χ = χ ( K , 0 ) at H = 0 ; (b) H -dependences of partial magnetizations m 1 = m 1 ( K c , H ) and m 2 = m 2 ( K c , H ) at the critical temperature. The curves are drawn at q 11 = 1 ,   q 22 = 2 ,   p 1 = 0.6 ,   p 2 = 0.4 . In panel (b), the curves m 1 , 2 = m 1 , 2 ( K c , H ) merge when q 12 = q 12 * .
Figure 4. Changes in functions of physical quantities at the “critical” antiferromagnetic interaction: (a) susceptibility χ = χ ( K , 0 ) at H = 0 ; (b) H -dependences of partial magnetizations m 1 = m 1 ( K c , H ) and m 2 = m 2 ( K c , H ) at the critical temperature. The curves are drawn at q 11 = 1 ,   q 22 = 2 ,   p 1 = 0.6 ,   p 2 = 0.4 . In panel (b), the curves m 1 , 2 = m 1 , 2 ( K c , H ) merge when q 12 = q 12 * .
Entropy 25 01428 g004
Figure 5. Temperature dependences at H < H c . From top to bottom: curves m 1 , 2 ( K ) , M ( K ) , U ( K ) and C ( K ) . (a) q 11 = 2 ,   q 22 = 1 ,   q 12 = 1 ,   p 1 = 0.6 ,   p 2 = 0.4 , H = 0.2 (b) q 11 = 3 ,   q 22 = 1 ,   q 12 = 1 ,   p 1 = 0.4 ,   p 2 = 0.6 , H = 0.2 . As we see, there is no phase transition when p 1 > p 2 , and there is a first-order phase transition when p 1 < p 2 (panel (b)).
Figure 5. Temperature dependences at H < H c . From top to bottom: curves m 1 , 2 ( K ) , M ( K ) , U ( K ) and C ( K ) . (a) q 11 = 2 ,   q 22 = 1 ,   q 12 = 1 ,   p 1 = 0.6 ,   p 2 = 0.4 , H = 0.2 (b) q 11 = 3 ,   q 22 = 1 ,   q 12 = 1 ,   p 1 = 0.4 ,   p 2 = 0.6 , H = 0.2 . As we see, there is no phase transition when p 1 > p 2 , and there is a first-order phase transition when p 1 < p 2 (panel (b)).
Entropy 25 01428 g005
Figure 6. Temperature dependences with H > H c . From top to bottom: curves m 1 , 2 ( K ) , M ( K ) , и C ( K ) . (a) q 11 = 2 ,   q 22 = 1 ,   q 12 = 1 ,   p 1 = 0.6 ,   p 2 = 0.4 , H = 0.61 (b) q 11 = 1 ,   q 22 = 2 ,   q 12 = q 12 * = 1 ,   p 1 = 0.4 ,   p 2 = 0.6 , H = 0.5 .
Figure 6. Temperature dependences with H > H c . From top to bottom: curves m 1 , 2 ( K ) , M ( K ) , и C ( K ) . (a) q 11 = 2 ,   q 22 = 1 ,   q 12 = 1 ,   p 1 = 0.6 ,   p 2 = 0.4 , H = 0.61 (b) q 11 = 1 ,   q 22 = 2 ,   q 12 = q 12 * = 1 ,   p 1 = 0.4 ,   p 2 = 0.6 , H = 0.5 .
Entropy 25 01428 g006
Figure 7. Temperature behavior of magnetizations in the symmetric model. (a) H = 0.75 H c ; (b) H = 0.995 H c ; (c) H = 1 . 0075 H c ; (d) H = 0 . 82 H c ; (e) H = 1 . 025 H c ; (f) H =   0 . 9998 H c . In all cases | q 12 | = 2 q , except plots (d) | q 12 | = 0.75 q and (f) | q 12 | = 5 q . q in all plots is equal to 4.
Figure 7. Temperature behavior of magnetizations in the symmetric model. (a) H = 0.75 H c ; (b) H = 0.995 H c ; (c) H = 1 . 0075 H c ; (d) H = 0 . 82 H c ; (e) H = 1 . 025 H c ; (f) H =   0 . 9998 H c . In all cases | q 12 | = 2 q , except plots (d) | q 12 | = 0.75 q and (f) | q 12 | = 5 q . q in all plots is equal to 4.
Entropy 25 01428 g007aEntropy 25 01428 g007b
Figure 8. Graphical solution of the equation R ( m s ) = H . Depicted is an example when at a given H , Equation (36) has one solution at | q 12 | > q (the upper curve), two solutions at | q 12 | = q (the middle curve), and no solutions at | q 12 | < q (the lower curve).
Figure 8. Graphical solution of the equation R ( m s ) = H . Depicted is an example when at a given H , Equation (36) has one solution at | q 12 | > q (the upper curve), two solutions at | q 12 | = q (the middle curve), and no solutions at | q 12 | < q (the lower curve).
Entropy 25 01428 g008
Figure 9. The layered model.
Figure 9. The layered model.
Entropy 25 01428 g009
Figure 10. Curves m 1 = m 1 ( K ) , m 2 = m 2 ( K ) and M = M ( K ) at different magnitudes of the external field. (a) H = 0 , (b) H = 0.25 H c , (c) H c = 0.975 H , (d) H = 0.9975 H c . Everywhere, L = 24 and q 11 = q 22 = 8 ,   q 12 = 8 .
Figure 10. Curves m 1 = m 1 ( K ) , m 2 = m 2 ( K ) and M = M ( K ) at different magnitudes of the external field. (a) H = 0 , (b) H = 0.25 H c , (c) H c = 0.975 H , (d) H = 0.9975 H c . Everywhere, L = 24 and q 11 = q 22 = 8 ,   q 12 = 8 .
Entropy 25 01428 g010
Figure 11. Temperature dependences of susceptibility χ . (a) q 11 = q 22 = 8 ,   q 12 = 8 ; (b) q 11 = q 22 = 8 ,   q 12 = 0.4 ; (c) q 11 = 9.6 ,   q 22 = 8 ,   q 12 = 8 ; (d) q 11 = 9.6 ,   q 22 = 8 ,   q 12 = 0.4 . The form of the curves χ ( K ) is independent of the lattice dimensions L for (a,b).
Figure 11. Temperature dependences of susceptibility χ . (a) q 11 = q 22 = 8 ,   q 12 = 8 ; (b) q 11 = q 22 = 8 ,   q 12 = 0.4 ; (c) q 11 = 9.6 ,   q 22 = 8 ,   q 12 = 8 ; (d) q 11 = 9.6 ,   q 22 = 8 ,   q 12 = 0.4 . The form of the curves χ ( K ) is independent of the lattice dimensions L for (a,b).
Entropy 25 01428 g011aEntropy 25 01428 g011b
Figure 12. Binder cumulants (a) q 11 = q 22 = 8 ,   q 12 = 8 , (b) q 11 = 16 ,   q 22 = 8 ,   q 12 = 8 .
Figure 12. Binder cumulants (a) q 11 = q 22 = 8 ,   q 12 = 8 , (b) q 11 = 16 ,   q 22 = 8 ,   q 12 = 8 .
Entropy 25 01428 g012
Figure 13. ln M as the function of ln H . The black broken line is a linear approximation used in the calculations of δ . (a) q 11 = q 22 = 8 ,   q 12 = 0.4 , which corresponds to q 12 = q 12 * ; (b) q 11 = 24 ,   q 22 = 8 ,   q 12 = 8 , which corresponds to q 12 q 12 * .
Figure 13. ln M as the function of ln H . The black broken line is a linear approximation used in the calculations of δ . (a) q 11 = q 22 = 8 ,   q 12 = 0.4 , which corresponds to q 12 = q 12 * ; (b) q 11 = 24 ,   q 22 = 8 ,   q 12 = 8 , which corresponds to q 12 q 12 * .
Entropy 25 01428 g013
Figure 14. The dependence of the cumulants on g s the reciprocal temperature K . The parameters of the model are: q 11 = q 22 = 8 ,   q 12 = 8 , H = 1 .
Figure 14. The dependence of the cumulants on g s the reciprocal temperature K . The parameters of the model are: q 11 = q 22 = 8 ,   q 12 = 8 , H = 1 .
Entropy 25 01428 g014
Figure 15. Determining the critical exponent β s . The parameters of the model are: q 11 = q 22 = 8 ,   q 12 = 8 , H = 2 . (a) Dependencies of ln ( m s M ) on ln ( t s ) when t s < 0 . The broken line is our linear approximation, from which we calculated β s . (b) Dependencies of ln ( m s M ) on ln t s when t s > 0 .
Figure 15. Determining the critical exponent β s . The parameters of the model are: q 11 = q 22 = 8 ,   q 12 = 8 , H = 2 . (a) Dependencies of ln ( m s M ) on ln ( t s ) when t s < 0 . The broken line is our linear approximation, from which we calculated β s . (b) Dependencies of ln ( m s M ) on ln t s when t s > 0 .
Entropy 25 01428 g015
Figure 16. Temperature dependencies of the susceptibility near the critical point. The parameters of the model are: q 11 = q 22 = 8 ,   q 12 = 8 . The magnitude of the external field is (a) H = 0.1 , (b) H = 2 .
Figure 16. Temperature dependencies of the susceptibility near the critical point. The parameters of the model are: q 11 = q 22 = 8 ,   q 12 = 8 . The magnitude of the external field is (a) H = 0.1 , (b) H = 2 .
Entropy 25 01428 g016
Figure 17. Determining the critical exponents δ s and δ s . The parameters of the model are: q 11 = q 22 = 8 ,   q 12 = 8 , H s = 2 .
Figure 17. Determining the critical exponents δ s and δ s . The parameters of the model are: q 11 = q 22 = 8 ,   q 12 = 8 , H s = 2 .
Entropy 25 01428 g017
Figure 18. Magnetization versus the field H : from bottom to top, p 1 p 2 = 0.1 , 0.2 , 0.3 . The curves are plotted for the case of K = 4 K c with q 11 = 4.4 , q 22 = 4 , q 12 = 2 .
Figure 18. Magnetization versus the field H : from bottom to top, p 1 p 2 = 0.1 , 0.2 , 0.3 . The curves are plotted for the case of K = 4 K c with q 11 = 4.4 , q 22 = 4 , q 12 = 2 .
Entropy 25 01428 g018
Table 1. Ground and metastable states in the zero-temperature limit ( K ).
Table 1. Ground and metastable states in the zero-temperature limit ( K ).
ConditionGround State
(Global Minimum)
Metastable State
(Local Minimum)
H < H c , p 1 > p 2 m 1 = 1 , m 2 = 1 , M = p 1 p 2 m 1 = 1 , m 2 = 1 , M = 1
H < H c , p 1 < p 2 m 1 = 1 , m 2 = 1 , M = p 2 p 1 m 1 = 1 , m 2 = 1 , M = 1
H > H c , p 1 > p 2 m 1 = 1 , m 2 = 1 , M = 1 m 1 = 1 , m 2 = 1 , M = p 1 p 2
H > H c , p 1 < p 2 m 1 = 1 , m 2 = 1 , M = 1 m 1 = 1 , m 2 = 1 , M = p 2 p 1
Table 2. Critical exponents.
Table 2. Critical exponents.
Critical Exponent q 11 p 1 + q 12 p 2 q 22 p 2 + q 12 p 1
( q 12 q 12 * )
q 11 p 1 + q 12 p 2 = q 22 p 2 + q 12 p 1
( q 12 = q 12 * )
α 00
β 1/23/2
γ = γ 10
δ 31
Scaling hypothesisconfirmed?
Table 3. Evaluation of the exponent δ for the layered model.
Table 3. Evaluation of the exponent δ for the layered model.
Model Parameters K c δ
Balanced system ( q 12 = q 12 * )
q 11 = q 22 = 8 ,   q 12 = 0.4 0.35800.9982
q 11 = q 22 = 8 ,   q 12 = 2 0.2710 0.9921
q 11 = q 22 = 8 ,   q 12 = 4 0.22170.9974
q 11 = q 22 = 8 ,   q 12 = 8 0.17071.0085
q 11 = q 22 = 8 ,   q 12 = 12 0.14301.0012
Unbalanced system ( q 12 q 12 * )
q 11 = 4 ,   q 22 = 8 ,   q 12 = 8 0.19852.2929
q 11 = 16 ,   q 22 = 8 ,   q 12 = 4 0.16154.4033
q 11 = 16 ,   q 22 = 8 ,   q 12 = 8 0.13102.5688
q 11 = 16 ,   q 22 = 8 ,   q 12 = 12 0.11253.1631
q 11 = 24 ,   q 22 = 8 ,   q 12 = 8 0.10103.1175
Table 4. Critical exponents for the layered model in the splitting point ( q 11 = q 22 = 8 ).
Table 4. Critical exponents for the layered model in the splitting point ( q 11 = q 22 = 8 ).
q 12 K c /MFM β m 1 , 2 H s K s /MFM/Equation (33) m s /MFM β s δ s / δ s
−40.2217/0.16670.28260.50.2257/0.1694/0.22460.1145/0.12710.61061.0688/0.8623
10.2387/0.1798/0.23500.2381/0.27000.58141.0474/0.6352
−80.1707/0.12500.31280.10.1707/0.1250/0.17070.0110/0.01250.60771.0425/1.0502
0.50.1717/0.1255/0.17120.0550/0.06270.59221.0703/1.0083
10.1744/0.1270/0.17280.1105/0.12630.58141.1014/0.7923
20.1866/0.1342/0.18010.2287/0.26220.56501.1015/0.6304
−160.1240/0.08330.303220.1271/0.0847/0.12540.1051/0.12600.50491.1609/0.7473
40.1383/0.0893/0.13010.2168/0.25880.49061.1924/0.5753
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kryzhanovsky, B.; Egorov, V.; Litinskii, L. Characteristics of an Ising-like Model with Ferromagnetic and Antiferromagnetic Interactions. Entropy 2023, 25, 1428. https://doi.org/10.3390/e25101428

AMA Style

Kryzhanovsky B, Egorov V, Litinskii L. Characteristics of an Ising-like Model with Ferromagnetic and Antiferromagnetic Interactions. Entropy. 2023; 25(10):1428. https://doi.org/10.3390/e25101428

Chicago/Turabian Style

Kryzhanovsky, Boris, Vladislav Egorov, and Leonid Litinskii. 2023. "Characteristics of an Ising-like Model with Ferromagnetic and Antiferromagnetic Interactions" Entropy 25, no. 10: 1428. https://doi.org/10.3390/e25101428

APA Style

Kryzhanovsky, B., Egorov, V., & Litinskii, L. (2023). Characteristics of an Ising-like Model with Ferromagnetic and Antiferromagnetic Interactions. Entropy, 25(10), 1428. https://doi.org/10.3390/e25101428

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop