删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

Monte Carlo study of the magnetic properties and magnetocaloric effect of an AFM/FM BiFeO/Co bilayer

本站小编 Free考研考试/2022-01-02

Chun-lu Chang1, Wei Wang,1,, He Ma,1,, Han Huang1, Jin-cheng Liu2, Rui-ze Geng21School of Science, Shenyang University of Technology, Shenyang 110870, China
2School 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) BiFeO3/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: BiFeO3/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 BiFeO3/Co bilayer. Communications in Theoretical Physics, 2021, 73(9): 095702- doi:10.1088/1572-9494/ac06bc

1. Introduction

BiFeO3 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, BiFeO3-based components, such as the BiFeO3/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 BiFeO3/Co bilayer based on diverse methods from theoretical and experimental aspects [810].

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 [1215]. 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 (ΔTad), 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 Ln0.67Sr0.33MnO3 (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 LaFe13−xSix compounds [24], Heulser alloys Ni-Mn-X with (X = Ga, Sn,…) [25], [Gd/Tb] multilayers [26], and so on. The experimental research of BiFeO3-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 BiFeO3-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 MnFe4Si3 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 MnFe4Si3 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 La1−xCaxMnO3 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 AlCMn3 [36], NdMnO3 [37], Ho3Pd2[38], Sr2CrMoO6 [39], Sr2CrIrO6 [40], and SmFe1−xMnxO3 [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 [4248]. We have ever devoted ourselves to the study of the dielectric properties of the BiFeO3/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 BiFeO3 has great potential applications because of its higher Curie temperature (TC = 1100 K) and Nell temperature (TN = 650 K) [54, 55]. However, few theoretical works have been conducted on the MCE in an AFM/FM BiFeO3/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 BiFeO3/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 BiFeO3/Co bilayer with mixed-spin (5/2, 3/2) consisting of 2N2 × 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 BiFeO3/FM bilayer structure.


The Hamiltonian of the system is considered as follows:$\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 Jab (<0) represents the intralayer exchange coupling between sublattices a and b in the upper layer, and Jcd (>0) refers to the intralayer exchange coupling between sublattices c and d in the lower layer. Jz 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, Jab is considered to be the unit of energy and temperature, thus ∣Jab∣ 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 BiFeO3/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 × N2 for such a bilayer system, we present in figure 2 the temperature dependence of susceptibility χ when 2 × N2 is changed from 2 × 52 to 2 × 502 with fixed ∣Jz∣/∣Jab∣ = 0.8, Jcd/∣Jab∣ = 1.2, T/∣Jab∣ = 0.05. From the figure, we can see every χ curve has a peak that corresponds to a temperature defined as critical temperature (TC). 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 TC when 2 × N2 is larger than 2 × 202, whereas the TC appears earlier for 2 × N2 smaller than 2 × 202, suggesting that the system can exhibit an obvious finite-size effect for small size. Consequently, in order to save computation time, we choose 2 × N2 = 2 × 202 for the following simulation. In order to balance the system, the first 2 × 105 MC steps (MCS) among the total was abandoned in our simulation and we used the rest 3 × 105 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 × N2 with ∣Jz∣/∣Jab∣ = 0.8, Jcd/∣Jab∣ = 1.2.


The magnetization of sublattices a, b, c, and d are calculated as:$\begin{eqnarray}{M}_{a}=\displaystyle \frac{1}{{N}^{2}}\sum _{i}{\sigma }_{{ia}}^{z},\end{eqnarray}$$\begin{eqnarray}{M}_{b}=\displaystyle \frac{1}{{N}^{2}}\sum _{j}{\sigma }_{{jb}}^{z},\end{eqnarray}$$\begin{eqnarray}{M}_{c}=\displaystyle \frac{1}{{N}^{2}}\sum _{m}{S}_{{mc}}^{z},\end{eqnarray}$$\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:$\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:$\begin{eqnarray}\chi =\beta (\langle {M}^{2}\rangle -\langle M{\rangle }^{2}),\end{eqnarray}$where $\beta =\tfrac{1}{{k}_{{\rm{B}}}T}$, kB is the Boltzmann constant and T represents the absolute temperature. In this study, kB 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:$\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:$\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 hi.

The adiabatic temperature change is deduced by:$\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:$\begin{eqnarray}C=7.378T+0.115{T}^{3}+{C}_{m}.\end{eqnarray}$

Finally, the RCP of the system is given by:$\begin{eqnarray}\mathrm{RCP}={\int }_{{T}_{1}}^{{T}_{2}}{\rm{\Delta }}{S}_{m}(T){\rm{d}}T,\end{eqnarray}$T1 and T2 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 (Jcd, ∣Jz∣) 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 (TC) of the bilayer system under the influence of the exchange interaction Jcd with other fixed parameters ∣Jz∣/∣Jab∣ = 0.8, h/∣Jab∣ = 0.5. Figure 3(a) shows the variation of magnetization with Jcd. Just one saturation is observed, so a conclusion can be drawn that M is not sensitive to the variation of interlayer exchange interaction Jcd. The M curve holds a constant when Jcd/∣Jab∣ = 2.0 for T is lower (0 ∼ 1), while with the decreasing of the value of Jcd/∣Jab∣, 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 Jcd/∣Jab∣. A similar phenomenon can also be seen in [38, 58]. The total susceptibility of the system varying with Jcd 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 Jcd/∣Jab∣ can cause the position of the peak to shift to a higher temperature region. Physically, with the increase of Jcd/∣Jab∣, more thermal energy is needed to overcome the forces from exchange coupling to break the magnetic order, thus TC increases correspondingly. For a deeper understanding of the variation of TC with Jcd, the TC curve is plotted in the (TCJcd/∣Jab∣) plane. In figure 3(c), a smooth curve can be seen. With the increase of Jcd/∣Jab∣, the curve gradually extends upward. In addition, the TC curve splits the entire phase plane into two parts, with the ferromagnetic phase below TC and the paramagnetic phase above TC. A similar phenomenon can be seen in [47, 48].

Figure 3.

New window|Download| PPT slide
Figure 3.The influence of Jcd on a) M, b) χ, and c) TC with fixed ∣Jz∣/∣Jab∣ = 0.8, h/∣Jab∣ = 0.5.


Similarly, as shown in figure 4, the influence of intralayer exchange coupling ∣Jz∣/∣Jab∣ on M, χ, TC when Jcd/∣Jab∣ =1.2, h/∣Jab∣ = 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 TC with ∣Jz∣. Obviously, the peak of the curves shifts to the right as ∣Jz∣ increases. The curve is shown in figure 4(c), depicted in the (TC − ∣Jz∣/∣Jab∣) plane. It is clear that as ∣Jz∣/∣Jab∣ increases, the TC goes up in response. The phase plane is also divided into two parts by TC—the lower and upper parts correspond to the magnetic phase and paramagnetic phase, respectively. Taking interlayer exchange coupling Jcd and intralayer exchange coupling ∣Jz∣ into consideration, we can conclude that both have a similar impact on χ and TC. 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 ∣Jz∣ on a) M, b) χ, and c) TC with fixed Jcd/∣Jab∣ = 1.2, h/∣Jab∣ = 0.5.


3.2. Magnetocaloric effect

Figure 5 shows the magnetic entropy change (−ΔS) and adiabatic temperature change (ΔTad) of the bilayer system for different values of Jcd/∣Jab∣ for fixed values of Δh/∣Jab∣ = 0.5, ∣Jz∣/∣Jab∣ = 0.8. From figure 5(a), a phenomenon can be seen whereby −ΔS initially increases with the increase of T/∣Jab∣, and then decreases with the continuous increase of T/∣Jab∣ 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 Jcd/∣Jab∣ takes the value of 0.2, the −ΔS is observed to be 0.033 75, while the ΔTad is observed to be 0.004 25. Moreover, the ΔTad curve tends to decline with the increase of T/∣Jab∣ more obviously than the −ΔS curve.

Figure 5.

New window|Download| PPT slide
Figure 5.The a) −ΔS and b) ΔTad as a function of T for various Jcd/∣Jab∣ with fixed ∣Jz∣/∣Jab∣ = 0.8 and h/∣Jab∣ = 0.5.


The −ΔS and ΔTad of the bilayer system under the effect of ∣Jz∣/∣Jab∣ for fixed Jcd/∣Jab∣ = 1.2, Δh/∣Jab∣ = 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 TC. The location of the maximum shifts to the high temperature region as ∣Jz∣ increases. The value of it can be seen as 0.026 when ∣Jz∣/∣Jab∣ = 0.2. Figure 6(b) presents the maximum moves to the high temperature region with the increase of ∣Jz∣/∣Jab∣, and it can reach 0.0028 for ∣Jz∣/∣Jab∣ = 0.2.

Figure 6.

New window|Download| PPT slide
Figure 6.The a) −ΔS and b) ΔTad as a function of T for various ∣Jz∣/∣Jab∣ with fixed Jcd/∣Jab∣ = 1.2 and h/∣Jab∣ = 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 ΔTad 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 Jcd/∣Jab∣ = 1.2, ∣Jz∣/∣Jab∣ = 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 ΔTad 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 ΔTad 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) ΔTad under the effects of Δh and T for Jcd/∣Jab∣ = 1.2, ∣Jz∣/∣Jab∣ = 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 Jcd/∣Jab∣ = 0.6, 1.2, and figure 8(b) for ∣Jz∣/∣Jab∣ = 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 Jcd/∣Jab∣ = 0.6 and 1.2, respectively. Similar results can be seen in [31, 36, 38], and RCP decreases gradually with the increase of Jcd/∣Jab∣ 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 ∣Jz∣/∣Jab∣ 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 Jcd/∣Jab∣ = 0.6, 1.2 with fixed ∣Jz∣/∣Jab∣; (b) for various ∣Jz∣/∣Jab∣ = 0.2, 0.8 with fixed Jcd/∣Jab∣ = 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 Jcd/∣Jab∣ with ∣Jz∣/∣Jab∣ = 0.8, T = 0.05. Only one loop can be seen and the area of the loops gets wider with Jcd/∣Jab∣ 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 BiFeO3/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 Jcd/∣Jab∣ with ∣Jz∣/∣Jab∣ = 0.8, T/∣Jab∣ = 0.05, a) Jcd/∣Jab∣ = 0.1, b) Jcd/∣Jab∣ = 0.8, c) Jcd/∣Jab∣ = 1.4, d) Jcd/∣Jab∣ = 2.1.


Analogically, figure 10 displays the magnetization as a function of the external magnetic field h for various values of ∣Jz∣/∣Jab∣ with Jcd/∣Jab∣ = 1.2, T = 0.05. Similarly, the increase of ∣Jz∣/∣Jab∣ 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 ∣Jz∣/∣Jab∣ with Jcd/∣Jab∣ = 1.2, T/∣Jab∣ = 0.05, a)∣Jz∣/∣Jab∣ = 0.1, b) ∣Jz∣/∣Jab∣ = 0.6, c) ∣Jz∣/∣Jab∣ = 1.3, d) ∣Jz∣/∣Jab∣ = 1.8.


Finally, the impact of temperature on the hysteresis loop with the parameter Jcd/∣Jab∣ = 1.2, ∣Jz∣/∣Jab∣ = 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 [6872].

Figure 11.

New window|Download| PPT slide
Figure 11.The hysteresis loops of the system for different values of T/∣Jab∣ with Jcd/∣Jab∣ = 1.2, ∣Jz∣/∣Jab∣ = 0.8, a) T/∣Jab∣ = 0.05, b) T/∣Jab∣= 0.5, c) T/∣Jab∣ = 1.6, d) T/∣Jab∣ = 3.5, e) T/∣Jab∣ = 4.8, f) T/∣Jab∣ = 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 BiFeO3/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 Jcd, ∣Jz∣ and h causes the system to undergo a rough process to accomplish phase transition, while the decrease of Jcd, ∣Jz∣, and the increase of h has a significant advantage in increasing the values of Δh, ΔTad, and RCP. Furthermore, the larger Jcd and ∣Jz∣ are beneficial for the enlargement of the area for the hysteresis loops. These results indicate that the BiFeO3/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 BiFeO3/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 By original order
By published year
By cited within times
By Impact factor

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]

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

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]

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

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]

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]

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]

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

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

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]

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

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

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

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]

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

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

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]

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]

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

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]

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

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]

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]

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

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]

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]

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]

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]

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]

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]

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]

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]

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]

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]

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

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]

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]

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

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

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]

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]

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]

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

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

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

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]

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]

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]

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]

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]

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]

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]

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]

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

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

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]

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

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

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

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]

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]

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]

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]

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]

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

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]

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]

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]

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

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

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

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]

相关话题/Monte Carlo study