Chun-lu Chang¹, Wei Wang^,1,∗, He Ma^,1,∗, Han Huang¹, Jin-cheng Liu², Rui-ze Geng²¹School of Science, Shenyang University of Technology, Shenyang 110870, China
²School of Electrical Engineering, Shenyang University of Technology, Shenyang 110870, China

First author contact: ^∗Author to whom any correspondence should be addressed.
Received:2021-03-4Revised:2021-05-31Accepted:2021-06-1Online:2021-08-05


Abstract
The magnetic properties and magnetocaloric effect of an antiferromagnetic/ferromagnetic (AFM/FM) BiFeO₃/Co bilayer with mixed-spin (5/2, 3/2) have been studied based on Monte Carlo simulation. The magnetization, susceptibility, and critical temperature are investigated under various exchange couplings and an external magnetic field. In particular, the influence of exchange couplings and an external magnetic field on the magnetic entropy change, adiabatic temperature change, and the relative cooling power (RCP) are studied. The simulation results indicated that the decrease of the exchange coupling and the increase of external magnetic fields can cause an increase of magnetic entropy change, adiabatic temperature change, and RCP. In addition, the hysteresis loops of the system are presented for different exchange couplings and temperatures.
Keywords： BiFeO₃/Co bilayer;magnetization;critical temperature;entropy change;hysteresis loops;Monte Carlo simulation


PDF (1482KB)MetadataMetricsRelated articlesExportEndNote|Ris|BibtexFavorite
Cite this article
Chun-lu Chang, Wei Wang, He Ma, Han Huang, Jin-cheng Liu, Rui-ze Geng. Monte Carlo study of the magnetic properties and magnetocaloric effect of an AFM/FM BiFeO₃/Co bilayer. Communications in Theoretical Physics, 2021, 73(9): 095702- doi:10.1088/1572-9494/ac06bc

1. Introduction

BiFeO₃ belongs to the family of multiferroic materials and also encompasses ceramic materials. It can play a significant role in the application of magnetoelectric read heads [1], magnetoelectric transformers [2], memory devices [3], and tunable microwave devices [4]. In practical applications, BiFeO₃-based components, such as the BiFeO₃/Co bilayer, which can be synthesized by rapid hot pressing [5], solid-state reaction methods [6], the Pechini method [7], and so on, have attracted much attention. In addition, much effort has been devoted to the study of the magnetic, as well as electrical properties of the BiFeO₃/Co bilayer based on diverse methods from theoretical and experimental aspects [8–10].

In recent years, research on the magnetocaloric effect (MCE), which was first proposed in 1881 by Warburg [11], is attractive because of its potential application in sustainable energy-related emerging technologies, for instance, magnetic refrigeration [12–15]. As far as we know, MCE is an intrinsic and exciting nature of the magnetic system, which can be defined by the fact that the temperature of a ferromagnet varies with the strength of the magnetic field in an adiabatic process [16, 17]. From a physical point of view, MCE is directly related to the coupling between magnetic atoms in the magnetic sub-lattice and the applied magnetic field [18]. It is worth mentioning that isothermal entropy change (−ΔS), adiabatic temperature change (ΔT_ad), and relative cooling power (RCP) are representative characteristics that can better describe MCE, and these parameters are necessary to obtain a high magnetic refrigeration efficiency [19]. There are many investigations on the experimental and theoretical study of MCE. Experimentally, there are many reports of MCE on diverse nanoparticles or compounds. To name a few, the Ln_0.67Sr_0.33MnO₃ (Ln = La, Pr, and Nd) nanoparticles were fabricated by the sol-gel method, and much work was devoted to studying their structural and magnetic properties [20]. Another conclusion is that all the samples showed considerable MCE in a wide temperature range. Some research results indicated that the physical properties of thin films are different from their bulk counterparts [21, 22]. A report suggested that MCE was successively reached at a low magnetic field [23], and the study also paid attention to the MCE of Ni-Mn-Ga thin films. Further, some other studies focused on the LaFe_13−xSi_x compounds [24], Heulser alloys Ni-Mn-X with (X = Ga, Sn,…) [25], [Gd/Tb] multilayers [26], and so on. The experimental research of BiFeO₃-based ceramics shows that the external magnetic and electric field have a significant influence on the magnetic properties of the system [27, 28]. In addition, the atom doping concentration can modify the MCE of the BiFeO₃-based system [29].

On the other hand, more efforts have been devoted to the theoretical investigations of MCE for nanomaterials, such as nanoparticles, nano-compounds, and thin film by various methods. A study was implemented by Hu et al to study the MCE in anisotropic zinc ferrite nanoparticles based on Monte Carlo (MC) simulation [30]. The magnetic entropy change and adiabatic temperature change were predicted in detail, and the emphasis was on the effect of anisotropy on these properties. In order to foster favorable materials for magnetic refrigeration devices, a study was carried out on the magnetic properties and MCE of the MnFe₄Si₃ system based on the MC method [31]. The structure was composed of a hexagonal lattice with Fe (S = 3/2) and Mn (S = 3/2). The magnetic entropy change, adiabatic temperature change, RCP, as well as magnetization, were obtained for different parameters and discussed in detail. In addition, the second-order phase transition was predicted at the Curie temperature. The results in this study provided good proof that MnFe₄Si₃ can be a suitable material for magnetic refrigeration. At present, many researchers are dedicated to the investigation of MCE for thin films by various methods. Giant magnetic entropy change is mainly rooted in the magnetic state change, which was concluded from a study of magnetic properties and MCE of multifunctional α-FeRh thin films [32]. MC simulation, as one of the most conventional and successful methods, has also been extended to study the MCE of thin films. An investigation of MCE for La_1−xCa_xMnO₃ thin films based on this method revealed that the external magnetic field has an important impact on MCE for such thin films [33]. The surface effect on the magnetocaloric properties of perovskite ferromagnetic thin films was studied by Elyacoubi et al based on MC simulation [34]. Some results were discussed in detail, such as the maximal entropy change and that RCP increases when enlarging the external magnetic field, while RCP decreases when the value of the surface exchange coupling increases. In the field of theoretical calculation of magnetic nanomaterials, the mixed-spin Ising model is a well-known model which can help in understanding the magnetic nature of various magnetic nano-systems and could be taken into account to investigate the MCE. Albayrak studied the entropy and isothermal entropy change of the mixed-spin(1,1/2) Ising model [35]. The influence of the crystal field and external magnetic field on the variation of the peak value of isothermal entropy change was discussed emphatically. Apart from the research mentioned above, there are many other investigations on the MCE of thin films, for instance AlCMn₃ [36], NdMnO₃ [37], Ho₃Pd₂[38], Sr₂CrMoO₆ [39], Sr₂CrIrO₆ [40], and SmFe_1−xMn_xO₃ [41]. All these studies reveal that exchange coupling and an applied external field play a significant role in the variation of the MCE, and this is exactly what our present research focuses on.

In our previous works, the MC simulation and the Ising mixed-spin model were significantly applied to investigate the magnetic, electric, as well as thermal dynamic, properties of various nanosystems [42–48]. We have ever devoted ourselves to the study of the dielectric properties of the BiFeO₃/FM bilayer system [49]. It should be mentioned that there have been many experimental [50, 51] and theoretical [52, 53] works which have studied the MCE behavior of multiferroic materials, and a study has released that BiFeO₃ has great potential applications because of its higher Curie temperature (T_C = 1100 K) and Nell temperature (T_N = 650 K) [54, 55]. However, few theoretical works have been conducted on the MCE in an AFM/FM BiFeO₃/FM bilayer system. This further stimulated our interest to conduct such a study to discuss, in depth, the effect of exchange coupling, the external magnetic field on MCE, and the magnetic properties of the BiFeO₃/FM bilayer system. The rest of this paper is arranged as follows: in the next section, the model and the method applied to this study are presented. In the third section, we show the detailed results of the simulation and provide some discussion. The last section includes a brief summary.

2. Model and method

An AFM/FM BiFeO₃/Co bilayer with mixed-spin (5/2, 3/2) consisting of 2N² × L spins on a three-dimensional honeycomb lattice, as plotted in figure 1, is constructed. N and L are defined as the number of spins in each layer and the thickness of the bilayer system, respectively. For such a bilayer system, we adopt L = 2.

Figure 1.

New window|Download| PPT slide
Figure 1.The schematic of an AFM/FM mixed-spin (5/2, 3/2) Ising model with a BiFeO₃/FM bilayer structure.


The Hamiltonian of the system is considered as follows:(1)$\begin{eqnarray}\begin{array}{rcl}H & = & -{J}_{{ab}}\sum _{\langle i,j\rangle }{\sigma }_{{ia}}^{z}{\sigma }_{{jb}}^{z}-{J}_{{cd}}\sum _{\langle m,n\rangle }{S}_{{mc}}^{z}{S}_{{nd}}^{z}-{J}_{z}\\ & & \times \sum _{\langle i(j),m(n)\rangle }{\sigma }_{{ia}({jb})}^{z}{S}_{{mc}({nd})}^{z}\\ & & -h(\sum _{i(j)}{\sigma }_{{ia}({jb})}^{z}+\sum _{m(n)}{S}_{{mc}({nd})}^{z}),\end{array}\end{eqnarray}$where the first item J_ab (<0) represents the intralayer exchange coupling between sublattices a and b in the upper layer, and J_cd (>0) refers to the intralayer exchange coupling between sublattices c and d in the lower layer. J_z denotes the interlayer exchange coupling which is positive between sublattices a and c, but negative between sublattices b and d. In the whole simulation process, J_ab is considered to be the unit of energy and temperature, thus ∣J_ab∣ takes the value of 1. ⟨i, j⟩, ⟨m, n⟩, and ⟨i, m⟩ stand for the summation over all the first nearest neighbor sites at sublattices a, b, sublattices c, d, and between sublattices a(b) and c(d), respectively. ${\sigma }_{{ia}}^{z}$,${\sigma }_{{jb}}^{z}$ are the z-component spins of spin-5/2 with the values ±5/2, ±3/2, ±1/2 and ${S}_{{mc}}^{z}$, ${S}_{{nd}}^{z}$ are the z-component spins of spin-3/2 with the values ±3/2, ±1/2. The last item h is the applied external magnetic field.

The method used for this simulation to calculate the BiFeO₃/Co bilayer is MC simulation based upon the standard Metropolis algorithm [56]. For the present bilayer, the periodic boundary condition is applied in the xy plane, and the free boundary condition is considered in the z-direction. In order to determine the finite-size effect caused by the number of each layer 2 × N² for such a bilayer system, we present in figure 2 the temperature dependence of susceptibility χ when 2 × N² is changed from 2 × 5² to 2 × 50² with fixed ∣J_z∣/∣J_ab∣ = 0.8, J_cd/∣J_ab∣ = 1.2, T/∣J_ab∣ = 0.05. From the figure, we can see every χ curve has a peak that corresponds to a temperature defined as critical temperature (T_C). It is an important parameter for better understanding magnetic phase transition from order to disorder for magnetic nanostructures [28, 34, 35]. It is found that the χ curves intersect each other at the same T_C when 2 × N² is larger than 2 × 20², whereas the T_C appears earlier for 2 × N² smaller than 2 × 20², suggesting that the system can exhibit an obvious finite-size effect for small size. Consequently, in order to save computation time, we choose 2 × N² = 2 × 20² for the following simulation. In order to balance the system, the first 2 × 10⁵ MC steps (MCS) among the total was abandoned in our simulation and we used the rest 3 × 10⁵ MCS to calculate the average values of the thermodynamic quantities as follows:

Figure 2.

New window|Download| PPT slide
Figure 2.The temperature dependence of χ for different system sizes 2 × N² with ∣J_z∣/∣J_ab∣ = 0.8, J_cd/∣J_ab∣ = 1.2.


The magnetization of sublattices a, b, c, and d are calculated as:(2)$\begin{eqnarray}{M}_{a}=\displaystyle \frac{1}{{N}^{2}}\sum _{i}{\sigma }_{{ia}}^{z},\end{eqnarray}$(3)$\begin{eqnarray}{M}_{b}=\displaystyle \frac{1}{{N}^{2}}\sum _{j}{\sigma }_{{jb}}^{z},\end{eqnarray}$(4)$\begin{eqnarray}{M}_{c}=\displaystyle \frac{1}{{N}^{2}}\sum _{m}{S}_{{mc}}^{z},\end{eqnarray}$(5)$\begin{eqnarray}{M}_{d}=\displaystyle \frac{1}{{N}^{2}}\sum _{n}{S}_{{nd}}^{z}.\end{eqnarray}$Thus, the total magnetization of the system is shown as:(6)$\begin{eqnarray}M=\displaystyle \frac{1}{4}({M}_{a}+{M}_{b}+{M}_{c}+{M}_{d}).\end{eqnarray}$The total susceptibility of the bilayer system is:(7)$\begin{eqnarray}\chi =\beta (\langle {M}^{2}\rangle -\langle M{\rangle }^{2}),\end{eqnarray}$where $\beta =\tfrac{1}{{k}_{{\rm{B}}}T}$, k_B is the Boltzmann constant and T represents the absolute temperature. In this study, k_B is assigned a value of 1 just to simplify the calculation procedure.

The characterization of the MCE can be described by the following formulas:

The magnetic entropy per site is obtained by:(8)$\begin{eqnarray}{S}_{m}(T,h)={\int }_{0}^{T}\displaystyle \frac{{C}_{m}(T,h)}{T}{\rm{d}}T.\end{eqnarray}$

The magnetic entropy change is presented by:(9)$\begin{eqnarray}{\rm{\Delta }}{S}_{m}(T,h)={\int }_{0}^{h}{\left(\displaystyle \frac{\partial M}{\partial T}\right)}_{{h}_{i}}{\rm{d}}{h}_{i},\end{eqnarray}$where ${\left(\tfrac{\partial M}{\partial T}\right)}_{{h}_{i}}$ denotes themal magnetization at a fixed magnetic field h_i.

The adiabatic temperature change is deduced by:(10)$\begin{eqnarray}{\rm{\Delta }}{T}_{{ad}}(T,h)=-T\displaystyle \frac{{\rm{\Delta }}{S}_{m}(T,h)}{C(T,h)}.\end{eqnarray}$

The C appearing in the above equation is the representation of the total magnetic specific heat. The lattice term, electronic term, and magnetic term are taken into account. On the strength of the Debye model [57], the approximated value of it can be described as:(11)$\begin{eqnarray}C=7.378T+0.115{T}^{3}+{C}_{m}.\end{eqnarray}$

Finally, the RCP of the system is given by:(12)$\begin{eqnarray}\mathrm{RCP}={\int }_{{T}_{1}}^{{T}_{2}}{\rm{\Delta }}{S}_{m}(T){\rm{d}}T,\end{eqnarray}$T₁ and T₂ are the cold and hot temperatures corresponding to both ends of the 1/2-maximum value of ${\rm{\Delta }}{S}_{m}^{\max }$, respectively.

3. Result and discussion

In this section, some significant results of the effects of multifarious physical parameters, for instance, exchange couplings (J_cd, ∣J_z∣) and external magnetic field h on the magnetic properties, the MCE as well as the hysteresis loops of a mixed-spin (5/2, 3/2) Ising bilayer, are presented in detail and discussed.

3.1. Magnetization and susceptibility

Figure 3 shows the magnetization (M), susceptibility (χ), and the critical temperature (T_C) of the bilayer system under the influence of the exchange interaction J_cd with other fixed parameters ∣J_z∣/∣J_ab∣ = 0.8, h/∣J_ab∣ = 0.5. Figure 3(a) shows the variation of magnetization with J_cd. Just one saturation is observed, so a conclusion can be drawn that M is not sensitive to the variation of interlayer exchange interaction J_cd. The M curve holds a constant when J_cd/∣J_ab∣ = 2.0 for T is lower (0 ∼ 1), while with the decreasing of the value of J_cd/∣J_ab∣, this phenomenon gradually becomes inconspicuous. In the perspective of theory, the magnetization is presented under the interaction of exchange coupling, temperature, and external magnetic field. A larger exchange coupling value is more advantageous in competition with temperature, which causes the spins in the sublattice to align in its direction, thus showing a stable magnetization externally. However, when the exchange coupling decreases or the temperature increases, the spins have more chance to flit, so the magnetization changes for the response. It gradually decreases with the increase in temperature, and at the same temperature, it decreases with the decrease of J_cd/∣J_ab∣. A similar phenomenon can also be seen in [38, 58]. The total susceptibility of the system varying with J_cd is displayed in figure 3(b). The six χ curves have the same shape. As the temperature rises, χ first increases and then gradually decreases after reaching a peak. The increase of J_cd/∣J_ab∣ can cause the position of the peak to shift to a higher temperature region. Physically, with the increase of J_cd/∣J_ab∣, more thermal energy is needed to overcome the forces from exchange coupling to break the magnetic order, thus T_C increases correspondingly. For a deeper understanding of the variation of T_C with J_cd, the T_C curve is plotted in the (T_C − J_cd/∣J_ab∣) plane. In figure 3(c), a smooth curve can be seen. With the increase of J_cd/∣J_ab∣, the curve gradually extends upward. In addition, the T_C curve splits the entire phase plane into two parts, with the ferromagnetic phase below T_C and the paramagnetic phase above T_C. A similar phenomenon can be seen in [47, 48].

Figure 3.

New window|Download| PPT slide
Figure 3.The influence of J_cd on a) M, b) χ, and c) T_C with fixed ∣J_z∣/∣J_ab∣ = 0.8, h/∣J_ab∣ = 0.5.


Similarly, as shown in figure 4, the influence of intralayer exchange coupling ∣J_z∣/∣J_ab∣ on M, χ, T_C when J_cd/∣J_ab∣ =1.2, h/∣J_ab∣ = 0.5 is studied. The M curve plotted in figure 4(a) is slightly different from that in figure 3(a). There is just one saturation as well, but as the temperature rises, the six M curves finally coincide with each other, which is different from the phenomenon in figure 3(a). Figure 4(b) shows the change of T_C with ∣J_z∣. Obviously, the peak of the curves shifts to the right as ∣J_z∣ increases. The curve is shown in figure 4(c), depicted in the (T_C − ∣J_z∣/∣J_ab∣) plane. It is clear that as ∣J_z∣/∣J_ab∣ increases, the T_C goes up in response. The phase plane is also divided into two parts by T_C—the lower and upper parts correspond to the magnetic phase and paramagnetic phase, respectively. Taking interlayer exchange coupling J_cd and intralayer exchange coupling ∣J_z∣ into consideration, we can conclude that both have a similar impact on χ and T_C. Similarly, these results are in good accordance with that in [46, 48, 58, 59].

Figure 4.

New window|Download| PPT slide
Figure 4.The influence of ∣J_z∣ on a) M, b) χ, and c) T_C with fixed J_cd/∣J_ab∣ = 1.2, h/∣J_ab∣ = 0.5.

3.2. Magnetocaloric effect

Figure 5 shows the magnetic entropy change (−ΔS) and adiabatic temperature change (ΔT_ad) of the bilayer system for different values of J_cd/∣J_ab∣ for fixed values of Δh/∣J_ab∣ = 0.5, ∣J_z∣/∣J_ab∣ = 0.8. From figure 5(a), a phenomenon can be seen whereby −ΔS initially increases with the increase of T/∣J_ab∣, and then decreases with the continuous increase of T/∣J_ab∣ after reaching the maximum value. We have to mention that the temperature corresponding to the maximum is near the critical temperature. This is consistent with that shown in [14, 37, 40], which present that the maximum of the −ΔS appears in the vicinity of the critical temperature. Thus, the prominent increase in the variation of −ΔS results from the heat that the system releases during the phase transition. Turning our attention to figure 5(b), similar curves to those shown in figure 4(a) are seen again. The biggest difference between them is the order of magnitude. When J_cd/∣J_ab∣ takes the value of 0.2, the −ΔS is observed to be 0.033 75, while the ΔT_ad is observed to be 0.004 25. Moreover, the ΔT_ad curve tends to decline with the increase of T/∣J_ab∣ more obviously than the −ΔS curve.

Figure 5.

New window|Download| PPT slide
Figure 5.The a) −ΔS and b) ΔT_ad as a function of T for various J_cd/∣J_ab∣ with fixed ∣J_z∣/∣J_ab∣ = 0.8 and h/∣J_ab∣ = 0.5.


The −ΔS and ΔT_ad of the bilayer system under the effect of ∣J_z∣/∣J_ab∣ for fixed J_cd/∣J_ab∣ = 1.2, Δh/∣J_ab∣ = 0.5 is shown in figure 6. In figure 6(a), each curve has a maximum, and the temperature corresponding to it is in the vicinity of T_C. The location of the maximum shifts to the high temperature region as ∣J_z∣ increases. The value of it can be seen as 0.026 when ∣J_z∣/∣J_ab∣ = 0.2. Figure 6(b) presents the maximum moves to the high temperature region with the increase of ∣J_z∣/∣J_ab∣, and it can reach 0.0028 for ∣J_z∣/∣J_ab∣ = 0.2.

Figure 6.

New window|Download| PPT slide
Figure 6.The a) −ΔS and b) ΔT_ad as a function of T for various ∣J_z∣/∣J_ab∣ with fixed J_cd/∣J_ab∣ = 1.2 and h/∣J_ab∣ = 0.5.


All the results regarding the properties of the bilayer system are simulated in the presence of an external magnetic field, so it is of great necessity to investigate the dependence of −ΔS and ΔT_ad on the Δh, which is displayed in figure 7. Figure 7(a) gives the −ΔS as a function of Δh and T. The other parameters are selected as J_cd/∣J_ab∣ = 1.2, ∣J_z∣/∣J_ab∣ = 0.8, respectively, and this figure shows us an intuitive phenomenon that when T changes from 6.0 to 9.8, there exists the maximum and it increases as Δh changes from 3.25 to 4.0. Figure 7(b) presents the ΔT_ad as a function of Δh and T. The maximum of the curves can be seen when T changes from 4.5 to 8.8, and it increases with the enlargement of Δh. Some comparisons may be implemented between these results and other studies, which investigated the influence of an external field on the −ΔS and ΔT_ad of the system [14, 31, 34, 35, 37, 38].

Figure 7.

New window|Download| PPT slide
Figure 7.The phase diagram of the a) −ΔS and b) ΔT_ad under the effects of Δh and T for J_cd/∣J_ab∣ = 1.2, ∣J_z∣/∣J_ab∣ = 0.8.


The relative cooling power (RCP) is a significant physical parameter to elucidate the character of the MCE. As given in figure 8, the Δh dependence of RCP corresponds to different values of exchange couplings—figure 8(a) for J_cd/∣J_ab∣ = 0.6, 1.2, and figure 8(b) for ∣J_z∣/∣J_ab∣ = 0.2, 0.8. It can be seen from figure 8(a) that RCP increases monotonously and smoothly as the field Δh increases, and can reach values of 2.00 and 1.75 for J_cd/∣J_ab∣ = 0.6 and 1.2, respectively. Similar results can be seen in [31, 36, 38], and RCP decreases gradually with the increase of J_cd/∣J_ab∣ for a certain Δh. A similar phenomenon in figure 8(a) can also be seen in figure 8(b)—that is, RCP goes smoothly with the increase of Δh, and the large ∣J_z∣/∣J_ab∣ corresponds to the small RCP. Therefore, a conclusion can be drawn that weakening the exchange coupling and heightening the external magnetic field can contribute to magnifying the RCP of the system and make the material a better candidate for a magnetic refrigerant. This conclusion is in accordance with [34, 37, 40].

Figure 8.

New window|Download| PPT slide
Figure 8.The influence of h on the RCP (a) for various J_cd/∣J_ab∣ = 0.6, 1.2 with fixed ∣J_z∣/∣J_ab∣; (b) for various ∣J_z∣/∣J_ab∣ = 0.2, 0.8 with fixed J_cd/∣J_ab∣ = 1.2.

3.3. Hysteresis loops

In this subsection, we discuss the hysteresis loops of the system under the influence of exchange couplings, and the temperature dependence of the hysteresis loop is also presented. Figure 9 presents the magnetization as a function of the external magnetic field h for different values of J_cd/∣J_ab∣ with ∣J_z∣/∣J_ab∣ = 0.8, T = 0.05. Only one loop can be seen and the area of the loops gets wider with J_cd/∣J_ab∣ changing from 0.1 to 2.1. It is reasonable that the coercive field and the residual magnetization augment with the increase of exchange coupling. This action originates from the enlargement of the exchange interaction energy. An analysis of these results shows that the present AFM/FM BiFeO₃/Co can exhibit some properties similar to those found in ferromagnetic materials [60, 61].

Figure 9.

New window|Download| PPT slide
Figure 9.The hysteresis loops of the system for different values of J_cd/∣J_ab∣ with ∣J_z∣/∣J_ab∣ = 0.8, T/∣J_ab∣ = 0.05, a) J_cd/∣J_ab∣ = 0.1, b) J_cd/∣J_ab∣ = 0.8, c) J_cd/∣J_ab∣ = 1.4, d) J_cd/∣J_ab∣ = 2.1.


Analogically, figure 10 displays the magnetization as a function of the external magnetic field h for various values of ∣J_z∣/∣J_ab∣ with J_cd/∣J_ab∣ = 1.2, T = 0.05. Similarly, the increase of ∣J_z∣/∣J_ab∣ is not advantageous for the number of hysteresis loops, but it is beneficial for the increase of the area for the loops.

Figure 10.

New window|Download| PPT slide
Figure 10.The hysteresis loops of the system for different values of ∣J_z∣/∣J_ab∣ with J_cd/∣J_ab∣ = 1.2, T/∣J_ab∣ = 0.05, a)∣J_z∣/∣J_ab∣ = 0.1, b) ∣J_z∣/∣J^ab∣ = 0.6, c) ∣J_z∣/∣J^ab∣ = 1.3, d) ∣J_z∣/∣J_ab∣ = 1.8.


Finally, the impact of temperature on the hysteresis loop with the parameter J_cd/∣J_ab∣ = 1.2, ∣J_z∣/∣J_ab∣ = 0.8 is shown in figure 11. From the figure, we can see that the remaining magnetization and the coercive field decrease gradually as T changes from a lower value to a larger; eventually, the loop vanishes at T = 6.5, corresponding to the critical temperature. Similar phenomena have also been found in the investigations of nanoislands [62, 63], graphene-like quantum dots [64], ladder-type graphene nanoribbons [65], a graphene monolayer [66] and triple-layer [67] as well as other magnetic nanosystems [68–72].

Figure 11.

New window|Download| PPT slide
Figure 11.The hysteresis loops of the system for different values of T/∣J_ab∣ with J_cd/∣J_ab∣ = 1.2, ∣J_z∣/∣J_ab∣ = 0.8, a) T/∣J_ab∣ = 0.05, b) T/∣J_ab∣= 0.5, c) T/∣J_ab∣ = 1.6, d) T/∣J_ab∣ = 3.5, e) T/∣J_ab∣ = 4.8, f) T/∣J_ab∣ = 6.5.

4. Conclusion

The magnetic properties, magnetocaloric effect as well as hysteresis behaviors have been investigated and predicted for a mixed-spin (5/2, 3/2) AFM/FM Ising BiFeO₃/Co bilayer applying the Monte Carlo method. The results show that either exchange coupling or an external magnetic field has an important influence on the magnetic properties and the magnetocaloric effects. The presence of J_cd, ∣J_z∣ and h causes the system to undergo a rough process to accomplish phase transition, while the decrease of J_cd, ∣J_z∣, and the increase of h has a significant advantage in increasing the values of Δh, ΔT_ad, and RCP. Furthermore, the larger J_cd and ∣J_z∣ are beneficial for the enlargement of the area for the hysteresis loops. These results indicate that the BiFeO₃/Co bilayer has appropriate properties to be a suitable candidate for use as a magnetic refrigerant material. We hope the study that we have implemented can provide direction for experimental research and theoretical reference for other works. It can be remarked that, although much effort has gone into the magnetic and electric properties of the BiFeO₃/Co bilayer, there is still much more to do, such as dynamic magnetic properties as well as the effect of atom doping, which will be investigated in detail in the future.

Reference View Option By original order
By published year
By cited within times
By Impact factor

[1]Zhang Y Li Z Deng C Ma J Lin Y Nan C W 2008 Appl. Phys. Lett. 92 152510
DOI:10.1063/1.2912032 [Cited within: 1]

[2]Roy S Majumder S B 2012 J. Alloys Compd. 538 153
DOI:10.1016/j.jallcom.2012.05.125 [Cited within: 1]

[3]Thomas R Scott J F Bose D N Katiyar R S 2010 J. Phys.: Condens. Matter 22 423201
DOI:10.1088/0953-8984/22/42/423201 [Cited within: 1]

[4]Scott J F 2007 Nat. Mater. 6 256
DOI:10.1038/nmat1868 [Cited within: 1]

[5]Wang C G Yang L Li Z W Zeng M Zhang A H Qin M H Lu X B Gao X S Gao J W Lam K H 2017 Ceram. Int. 43 12764
DOI:10.1016/j.ceramint.2017.06.163 [Cited within: 1]

[6]Liu N T Liang R H Zhao X B Zhang Y Y Zhou Z Y Tang X D Dong X L 2018 J. Alloys Compd. 740 470
DOI:10.1016/j.jallcom.2018.01.036 [Cited within: 1]

[7]Huo S X Yuan S L Qiu Y Ma Z Z Wang C H 2012 Mater. Lett. 68 8
DOI:10.1016/j.matlet.2011.09.081 [Cited within: 1]

[8]Zhang Y Q et al. 2016 Appl. Surf. Sci. 367 418
DOI:10.1016/j.apsusc.2016.01.210 [Cited within: 1]

[9]Wang K Si N Zhang Y L Zhang F Guo A B Jiang W 2019 Vacuum 165 105
DOI:10.1016/j.vacuum.2019.04.009

[10]Zhang H X Cheng S Li B Cheng X W Cheng Q F 2018 Sep. Purif. Technol. 202 242
DOI:10.1016/j.seppur.2018.03.072 [Cited within: 1]

[11]Warburg E 1881 Ann. Phys. 249 141
DOI:10.1002/andp.18812490510 [Cited within: 1]

[12]Mañosa L Planesa A Acet M 2013 J. Mater. Chem. 1 4925
DOI:10.1039/c3ta01289a [Cited within: 1]

[13]Franco V Blázquez J S Ingale B Conde A 2012 Annu. Rev. Mater. Res. 42 305
DOI:10.1146/annurev-matsci-062910-100356

[14]Benhouria Y Khossossi N Houmad M Essaoudi I Ainane A Ahuja R 2019 Phys. E 105 139
DOI:10.1016/j.physe.2018.09.008 [Cited within: 2]

[15]Balli M Jandl S Fournier P Kedous-Lebouc A 2017 Appl. Phys. Rev. 4 021305
DOI:10.1063/1.4983612 [Cited within: 1]

[16]Pecharsky V K Gschneidner K AJr 1999 J. Magn. Magn. Mater. 200 44
DOI:10.1016/S0304-8853(99)00397-2 [Cited within: 1]

[17]Laghrissi A Zahraouy E M S H E Benyoussef A 2016 J. Magn. Magn. Mater. 412 259
DOI:10.1016/j.jmmm.2016.04.009 [Cited within: 1]

[18]Omari L H Lekdadri A Hlil E K Lassri H Moubah R Tozri A Dhahri E Kerouad M Zaim A 2020 Mater. Today Commun. 23 100857
DOI:10.1016/j.mtcomm.2019.100857 [Cited within: 1]

[19]Bhattacharyya A Giri S Majumdar S 2011 J. Magn. Magn. Mater. 323 1484
DOI:10.1016/j.jmmm.2011.01.004 [Cited within: 1]

[20]Wang G F Li L R Zhao Z R Yu X Q Zhang X F 2014 Ceram. Int. 40 16449
DOI:10.1016/j.ceramint.2014.07.154 [Cited within: 1]

[21]Dunand D C Müllner P 2011 Adv. Mater. 23 216
DOI:10.1002/adma.201002753 [Cited within: 1]

[22]Shinde K P Sinha B B Oh S S Kim H S Ha H S Baik S K Chung K C Kim D S Jeong S 2015 J. Magn. Magn. Mater. 374 144
DOI:10.1016/j.jmmm.2014.08.007 [Cited within: 1]

[23]Oumezzine M Amaral J S Mompean F J Hernandez M G Oumezzine M 2016 J. RSC Adv. 6 32193
DOI:10.1039/C6RA03328H [Cited within: 1]

[24]Kuzmin M D Richter M 2007 Phys. Rev. B 76 092401
DOI:10.1103/PhysRevB.76.092401 [Cited within: 1]

[25]Dubenko I Khan M Pathak A K Gautam B R Stadler S Ali N 2009 J. Magn. Magn. Mater. 321 754
DOI:10.1016/j.jmmm.2008.11.043 [Cited within: 1]

[26]Lambert C H El Hadri M S Hamedoun M Benyoussef A Mounkachi O Mangin S 2017 J. Magn. Magn. Mater. 443 1
DOI:10.1016/j.jmmm.2017.06.115 [Cited within: 1]

[27]Rahman A Wang Z P Dai R C Zhang Z M 2018 J. Magn. Magn. Mater. 462 230 235
DOI:10.1016/j.jmmm.2018.05.005 [Cited within: 1]

[28]Pan Q Fang C Luo H S Cun B J 2019 J. Eur. Ceram. Soc. 39 1057 1064
DOI:10.1016/j.jeurceramsoc.2018.12.045 [Cited within: 2]

[29]Makoed I I Amirov A A Liedienov N A Pashchenko A V Yanushkevich K I 2019 Solid State Sci. 95 105920
DOI:10.1016/j.solidstatesciences.2019.06.009 [Cited within: 1]

[30]Hu T Y Chi X D Lu Q Yu L Li R J Liu Y Du A Li Z B Shi F Hu Y 2019 J. Alloys Compd. 801 465
DOI:10.1016/j.jallcom.2019.06.155 [Cited within: 1]

[31]Bouachraoui R Ziat Y Sbai Y Rhazouani O E Goumrhar F Bahmad L 2019 J. Alloys Compd. 809 151785
DOI:10.1016/j.jallcom.2019.151785 [Cited within: 3]

[32]Zhou T J Cher M K Shen L Hu J F Yuan Z M 2013 Phys. Lett. A 377 3052
DOI:10.1016/j.physleta.2013.09.027 [Cited within: 1]

[33]Restrepo-Parra E Ramos-Rivera L Londoño-Navarro J 2014 J. Magn. Magn. Mater. 351 65
DOI:10.1016/j.jmmm.2013.09.011 [Cited within: 1]

[34]Elyacoubi A S E Masrour R JAbar A 2018 Appl. Surf. Sci. 459 537
DOI:10.1016/j.apsusc.2018.08.020 [Cited within: 4]

[35]Albayrak E 2020 Phys. A 559 125079
DOI:10.1016/j.physa.2020.125079 [Cited within: 3]

[36]Ennassiri N Tahiri N El Bounagui O Ez-Zahraouy H Benyoussef A 2018 J. Alloys Compd. 741 1196
DOI:10.1016/j.jallcom.2018.01.223 [Cited within: 2]

[37]Masrour R JAbar A Benyoussef A Hamedoun M Hlil E K 2016 J. Magn. Magn. Mater. 401 91
DOI:10.1016/j.jmmm.2015.10.019 [Cited within: 4]

[38]Kadim G Masrour R JAbar A 2020 J. Magn. Magn. Mater. 499 166263
DOI:10.1016/j.jmmm.2019.166263 [Cited within: 4]

[39]Rhazouani O E Slassi A 2017 Comput. Condens. Matter 11 55
DOI:10.1016/j.cocom.2017.05.002 [Cited within: 1]

[40]Rhazouani O E Slassi A Ziat Y Benyoussef A 2017 J. Phys. Chem. Solids 104 32
DOI:10.1016/j.jpcs.2016.12.033 [Cited within: 3]

[41]Elyacoubi A S E Masrour R Jabar A 2018 Solid State Commun. 271 39
DOI:10.1016/j.ssc.2017.12.015 [Cited within: 1]

[42]Wang W Chang C L Li Q Xue F L Sun L Huang T 2019 Superlattices Microstruct. 136 106325
DOI:10.1016/j.spmi.2019.106325 [Cited within: 1]

[43]Wang W Wu H J Xie P N Lv J H Peng Z 2020 Eur. Phys. J. Plus 135 605
DOI:10.1140/epjp/s13360-020-00627-x

[44]Lv D Yang Y Jiang W Wang F Gao Z Y Tian M 2019 Phys. A 514 319
DOI:10.1016/j.physa.2018.09.089

[45]Sun L Wang W Li Q Wang F Wu H J 2020 Superlattices Microstruct. 147 106701
DOI:10.1016/j.spmi.2020.106701

[46]Yang M Wang W Li B C Wu H J Yang S Q Yang J 2020 Phys. A 539 122932
DOI:10.1016/j.physa.2019.122932 [Cited within: 1]

[47]Wu H J Wang W Lv D Chang C L Li B C Tian M 2020 J. Magn. Magn. Mater. 515 167306
DOI:10.1016/j.jmmm.2020.167306 [Cited within: 1]

[48]Wu H J Wang W Wang F Li B C Li Q Xu J H 2020 J. Phys. Chem. Solids 136 109174
DOI:10.1016/j.jpcs.2019.109174 [Cited within: 3]

[49]Chang C L Wang W Lv D Geng R Z Liu J C Huang H 2020 Ceram. Int. 46 22907
DOI:10.1016/j.ceramint.2020.06.064 [Cited within: 1]

[50]Dey K Indra A Majumdar S Giri S 2017 J. Magn. Magn. Mater. 435 15
DOI:10.1016/j.jmmm.2017.03.068 [Cited within: 1]

[51]Ursic H Bobnar V Malic B Filipic C Vrabelj M Drnovsek S Jo Y Wencka M Kutnjak Zd 2016 Sci. Rep. UK 6 26629
DOI:10.1038/srep26629 [Cited within: 1]

[52]Liu Y Zhang G Z Qi L Laurent B Scott J F Dkhil B Wang Q 2016 Phys. Rev. B94 214113
DOI:10.1103/PhysRevB.94.214113 [Cited within: 1]

[53]Bai G Xie Q Y Qin X S Xu J Yan X B Gao C F 2018 Ferroelectrics 531 186
DOI:10.1080/00150193.2018.1497414 [Cited within: 1]

[54]Cao X S 2019 J. Alloys Compd. 806 1008
DOI:10.1016/j.jallcom.2019.07.330 [Cited within: 1]

[55]Ramachandran B Ramachandra Rao M S 2009 Appl. Phys. Lett. 95 142505
DOI:10.1063/1.3242411 [Cited within: 1]

[56]Metropolis N Rosenbluth A W Rosenbluth M N Teller A H 1953 J. Chem. Phys. 21 1087
DOI:10.1063/1.1699114 [Cited within: 1]

[57]Debye P 1912 Ann. Phys. 39 789
DOI:10.1002/andp.19123441404 [Cited within: 1]

[58]Benhouria Y Essaoudi I Ainane A Ahuja R Dujardin F 2017 Ferroelectrics 507 58
DOI:10.1080/00150193.2017.1283170 [Cited within: 2]

[59]Sun L Wang W 2020 Commun. Theor. Phys. 72 115703
DOI:10.1088/1572-9494/abb7d0 [Cited within: 1]

[60]Wang W Li Q Wang M Z Ma Y Guo A B Huang T 2019 Physica E 111 63 74
DOI:10.1016/j.physe.2019.02.028 [Cited within: 1]

[61]Lv D Ma Y Luo X H Jiang W Wang F Li Q 2020 Physica E 116 113721
DOI:10.1016/j.physe.2019.113721 [Cited within: 1]

[62]Yang Y Wang W Ma H Li Q Gao Z Y Huang T 2019 Physica E 108 358 371
DOI:10.1016/j.physe.2018.11.038 [Cited within: 1]

[63]Yang S Q Wang W Wang F Li B C Wu H J Yang M Xu J H 2019 J. Phys. Chem. Solids 135 109110
DOI:10.1016/j.jpcs.2019.109110 [Cited within: 1]

[64]Sun L Wang W Lv D Gao Z Y Li Q Li B C 2021 J. Magn. Magn. Mater. 528 167820
DOI:10.1016/j.jmmm.2021.167820 [Cited within: 1]

[65]Gao Z Y Lv D Wang W Yu J 2021 Polymer 223 123678
DOI:10.1016/j.polymer.2021.123678 [Cited within: 1]

[66]Wei W Sun L Li Q Lv D Gao Z Y Haung T 2021 J. Magn. Magn. Mater. 527 167692
DOI:10.1016/j.jmmm.2020.167692 [Cited within: 1]

[67]Sun L Lv D Wang W Gao Z Y Li B C 2021 Phys. Scr. 151 075809
DOI:10.1088/1402-4896/abfc80 [Cited within: 1]

[68]Chen J Y Zhang X D Yang L M Wang F 2021 Commun. Theor. Phys. 73 045702
DOI:10.1088/1572-9494/abe367 [Cited within: 1]

[69]Mai Z Y Zhang X D Liu Y T Yu H Wang F 2020 Vacuum 177 109427
DOI:10.1016/j.vacuum.2020.109427

[70]Chang C L Wang W Lv D Liu Z Y Tian M 2021 Eur. Phys. J. Plus 136 290
DOI:10.1140/epjp/s13360-021-01280-8

[71]Li Q Li R D Wang W Geng R Z Huang H Zheng S J 2020 Phys. A 555 124741
DOI:10.1016/j.physa.2020.124741

[72]Lv D Zhang D Z Yang M Wang F Yu J 2021 Superlattices Microstruct. 96 106833
DOI:10.1016/j.spmi.2021.106833 [Cited within: 1]