A.
Omelyanchik
ab,
S.
Villa
a,
M.
Vasilakaki
c,
G.
Singh
d,
A. M.
Ferretti
e,
A.
Ponti
f,
F.
Canepa
a,
G.
Margaris
c,
K. N.
Trohidou
c and
D.
Peddis
*ag
aDepartment of Chemistry and Industrial Chemistry (DCIC), University of Genova, Genova, Italy. E-mail: davide.peddis@unige.it
bImmanuel Kant Baltic Federal University, Kaliningrad, Russia
cInstitute of Nanoscience and Nanotechnology, National Center for Scientific Research Demokritos, Athens 15310, Greece
dEngineering School of Biomedical Engineering, Sydney Nano Institute, The University of Sydney, Sydney, Australia
eIstituto di Scienze e Tecnologie Chimiche “Giulio Natta”, Via G. Fantoli 16/15, 20138 Milano, Italy
fIstituto di Scienze e Tecnologie Chimiche “Giulio Natta”, Via C. Golgi 19, 20133 Milano, Italy
gIstituto di Struttura Della Materia, CNR, 00015 Monterotondo Scalo, RM, Italy
First published on 4th October 2021
The synthesis strategy and magnetic characterisation of two systems consisting of nanoparticles with core/shell morphology are presented: an assembly of hard/soft nanoparticles with cores consisting of magnetically hard cobalt ferrite covered by a magnetically soft nickel ferrite shell, and the inverse system of almost the same size and shape. We have successfully designed these nanoparticle systems by gradually varying the magnetic anisotropy resulting in this way in the modulation of the magnetic dipolar interactions between particles. Both nanoparticle systems exhibit high saturation magnetisation and display superparamagnetic behaviour at room temperature. We have shown strong exchange coupling at the core/shell interface of these nanoparticles systems which was also confirmed by mesoscopic modelling. Our results demonstrate the possibility of modulating magnetic anisotropy not only by chemical composition but also by adopting the proper nano-architecture.
Precise control of morpho-structural features (i.e. size of the core, thickness of the shell) and chemical composition of the core/shell MNPs adds the necessary degrees of freedom to control magnetic properties for specific applications. The fine tuning of coercivity and blocking temperature has been achieved by the modulation of the morphology (i.e., core or shell) of soft and hard materials.7,16–20 Although the general trend of the simple linear dependence (additivity) of the coercivity on the volume ratio of materials holds for the hard/soft system, it was observed that it breaks for the inverted soft/hard system.17 As was reported by Song and Zhang,17 the linear trend is oversimplified and does not take into account magnetisation reversal processes which are expected to be different when the order of shell and core materials changes. Moreover, the significant contribution of the surface spins in the magnetic anisotropy of single magnetic phase CoFe2O4 nanoparticles has been observed.21 For core/shell systems, a phenomenological model considering both surface and interface effects was developed by Trohidou et al. which estimated Keff for a two-phase exchange-coupled nanoparticle system:22
(1) |
In core/shell MNPs assemblies, a very important but not efficiently explored issue is the interplay between intra- and interparticle interactions and their influence on the macroscopic magnetic properties of the system. Ming Da Yang et al.,23 demonstrated that modulating interparticle interactions in core/shell systems results in the stabilisation of suspension properties, thus increasing the heating performance in hyperthermia as well as R2 relaxivity for contrast enhancement in MRI. Samuel D. Oberdick et al. employed theoretical calculations and showed that the Dzyaloshinskii–Moriya interactions were responsible for the experimentally observed spin canting in bi-magnetic Fe3O4/MnxFe1−xO4 (x > 0.8) particles.24 Moreover, interparticle interactions lead to the correlation of spin canting angle between neighbouring particles. Thus, the macroscopic magnetic properties of this system are determined by the strong interplay of intra- and interparticle interactions. Theoretical description of interacting bi-magnetic core/shell nanoparticle systems is a complex task due to the need to consider both long-range and short-range interactions. Margaris et al. have developed a mesoscopic scale model to accomplish this task in the case of dipolarly interacting Co nanoparticles surrounded by a thin CoO shell.25 In this study, they showed that the interplay of intraparticle and dipolar and exchange interparticle interactions defined the observed magnetic behavior.
The synthesis and characterisation of bi-magnetic core/shell ferrite nanoparticles with conventional (hard/soft) and inverse (soft/hard) structures are reported and compared with single-phase systems. The ferrites CoFe2O4 (CFO) and NiFe2O4 (NFO) were chosen due to the increasing interest in their applications (e.g., biomedicine, catalysis, etc.26,27) as well as the significant difference between the magnetocrystalline anisotropy, KNFOV = −6.2 × 103 J m−3 and KCFOV = 2 × 105 J m−3, and the relatively high value of their magnetisation Ms (MCFOs = 88 A m2 kg−1 and MNFOs = 55 A m2 kg−1).28,29 While most works deal with iron oxides or Mn/Zn substituted ferrites as a soft phase, only a few works are devoted to nickel ferrite in core/shell systems,4,14,30 although nickel ferrite being a promising material for catalysis31 and spintronics32 applications.
In this context, the paper focuses on the interplay between inter- and intraparticle interactions in core/shell systems. We approach this issue by studying the magnetization reversal mechanisms in terms of the switching field distribution (SFD),33 and the anisotropy energy barrier distribution.34–36 Finally, Monte Carlo simulations were performed on interacting assemblies of these core/shell nanoparticles taking into account intraparticle (i.e. interface and surface effects) and interparticle interactions. A mesoscopic model has been developed considering each nanoparticle with the minimum number of spins necessary to describe its morphology. This model includes for each particle four spins and it considers the actual size of the magnetic moments for each region, the actual contribution of the volume anisotropies in each region and the interparticle interactions. To the best of our knowledge, it is the first attempt to apply a mesoscopic theoretical model for core/shell magnetic nanoparticles system including interface and surface effects to investigate the proximity effects and how these effects modulate the interparticle interactions. We compare the findings of the model with the experimental results.
In the present work, all the syntheses were conducted under oxygen-free conditions in a Schlenk line. The synthesis of core/shell nanoparticles was carried out using a two-step procedure in which preformed Co or Ni ferrite nanoparticles were used as seeds for the subsequent growth of Co or Ni ferrites. The core was synthesised using a modified procedure reported elsewhere.37 For example, to obtain the CoFe2O4, Fe(acac)3 (2 mmol, 0.706 g), Co(acac)2 (1 mmol, 0.267 g), 1,2-hexadecanediol (10 mmol, 2.58 g), oleic acid (6 mmol, 1.69 g) and oleylamine (6 mmol, 1.60 g), were dissolved in benzyl ether (20 mL) and magnetically stirred under a flow of nitrogen. The mixture was heated (5 °C min−1) up to 200 °C for 100 min and then heated (5 °C min−1) to reflux (≈300 °C) for 60 min. The black-coloured mixture was left to cool to room temperature (RT ≈ 25 °C) and washed with toluene and isopropanol. Next, the precipitate was washed with ethanol and collected by centrifugation (4000 rpm, 10 min). The washing with ethanol was repeated three times in total.
A seed-mediated growth at high temperature was used to achieve core/shell nanoparticles. For example, for CoFe2O4/NiFe2O4 MNPs, cobalt ferrite seeds were sonicated in 5 mL benzyl ether in the presence of oleic acid (600 μL): then, the suspension was added to the flask containing Ni(acac)2 (0.33 mmol, 0.0848 g), Fe(acac)3 (0.67 mmol, 0.237 g) and oleic acid (600 μL) were dissolved benzyl ether (15 mL). After 1 h under vacuum at RT, the solution was heated to reflux at 290 °C (5 °C min−1) under argon atmosphere and kept at this temperature for 30 min. The black-coloured mixture was left to cool to RT and magnetically washed overnight with toluene and isopropanol. Next, the precipitate was washed with ethanol and collected by centrifugation (4000 rpm, 10 min). The washing with ethanol was repeated three times in total.
(2) |
Scanning Transmission Electron Microscopy (STEM) was used to investigate the size distribution and morphology of the particles by means of a Hitachi S5500 STEM operating at 30 kV. Finely ground samples were dispersed in ethanol and submitted to an ultrasonic bath. The suspensions were then dropped on carbon-coated copper grids for the STEM observations. A log-normal function38 was fitted to the particle size distribution obtained by STEM image analysis:
(3) |
The TEM images and EFTEM (Energy Filtered TEM) maps used to identify elemental distribution with nanometric resolution inside the core/shell MNPs were collected by a ZEISS LIBRA 200FE-HR TEM, operating at 200 kV and equipped with a second generation in-column Ω filter. Elemental maps of the elements were recorded at 708 eV, corresponding to the iron L3 white line, at 779 eV corresponding to the cobalt L3 white line, and at 855 eV, corresponding to the nickel L3 white line, using an energy window of 12 eV. The images were processed by means of the TEM Imaging Platform software (Olympus). The samples were prepared by dropping 7 μL of core/shell MNP suspension on a copper grid, coated with a carbon ultrathin film (3–5 nm thick), and let it dry overnight.
The DC magnetisation measurements were performed by a quantum design superconducting quantum interference device (SQUID) magnetometer (μ0Hmax = 5 T). Since after synthesis some surfactant (oleic acid, oleylamine) remains bound to the surface of the particles, the mass of the inorganic component of the samples was measured by thermogravimetric analysis and the magnetisation was calculated with respect to the inorganic component. The sample in the form of a powder was immobilised in epoxy resin to avoid any movement of the nanoparticles during the measurements, which can distort the shape of the hysteresis loop and eliminate remanence effects following elaboration of data.39,40 Hysteresis loops were obtained in a −5/+5 T applied magnetic field at temperatures from 5 to 300 K. The value of saturation magnetization (Ms) was estimated by law of approach to saturation (i.e., H: M = Ms(1 − A/H − B/H2), where A and B are free parameters) to high values of fields.41
The magnetic study of the interparticle interactions was performed at 5 K by field dependent remanent magnetization.42 The isothermal remanent magnetisation (IRM) was obtained starting from a totally demagnetised state. A positive magnetic field (μ0HR) was applied for 10 s, then removed and the remanence MIRM(H) was measured. The process was repeated, increasing the field up to 5 T. In the direct current demagnetisation (DCD) measurements, the sample has been saturated at −5 T; then positive magnetic field was applied and, after 10 s, it was switched off and the remanent magnetization was measured. This was repeated, increasing the field up to 5 T. Then the ΔM-plot was plotted and analysed.36,42–44 Qualitative analysis of interparticle interactions is based on the calculation of the interaction field μ0HINT as half the difference in the position of the maxima (μ0HCr) of the derivatives dMDCD(H)/dH and dMIRM(H)/dH.
The thermal dependence of magnetization was investigated by the zero field cooled (ZFC) and field cooled (FC) protocols. In ZFC protocol the sample was first cooled down from 300 to 5 K in zero magnetic field, then a static magnetic field of 2.5 mT was applied and MZFC was measured during warming up. In FC protocol the sample was cooled down to 5 K under the same magnetic field and MFC was measured during the cooling.
Each nanoparticle in the assembly consists of four regions as in an atomic scale model:45 the core, the interface (including both the core and shell interface), the shell and the surface. The four regions are represented by a set of four classical macro spins vectors: core, IF, shell, surf. The size of these spins in each assembly is given in Table 1 and it is calculated as sregion = Ms,region·Nregion/Ms,core·Ncore. Here Nregion is the number of spins in each of the above particle regions and Ms,core, Ms,shell the bulk magnetisation of the core and the shell in the CFO/NFO and NFO/CFO systems. In our model for pure CoFe2O4 or NiFe2O4 nanoparticles, we consider two regions, the core and the surface, consequently, we use two macro spins, one for the core and one for the surface to describe them.
Sample | Size of macro spins | Anisotropy constant | ||||||
---|---|---|---|---|---|---|---|---|
s core | s IF | s shell | s surf | k core | k IF | k shell | k surf | |
CFO | 0.796 | 0.0 | 0.0 | 0.204 | 0.8 | 0.0 | 0.0 | 2.6 |
CFO/NFO | 0.796 | 0.396 | 2.642 | 0.473 | 0.8 | 2.73 | 0.087 | 3.12 |
NFO | 0.796 | 0.0 | 0.0 | 0.204 | 0.019 | 0.0 | 0.0 | 0.095 |
NFO/CFO | 0.796 | 0.509 | 4.912 | 0.88 | 0.019 | 2.974 | 3.6024 | 0.774 |
In all cases, the number of spins in each region is calculated from an atomic scale model for each nanoparticle, reproducing the geometry of the inverse spinel structure and taking into account the dimensions of the nanoparticle in the assembly. We consider a random anisotropy assembly, so that each particle is assigned a randomly oriented uniaxial easy axis, common to all four spins in the particle.
The effective anisotropy constants of the four spins have been calculated as kregion = Kregion·Vregion/Kcore·Vcore where Kregion and Vregion are the anisotropy energy density and the volume of the corresponding region. Kcore, and Kshell are taken to be the values of the bulk for CoFe2O4 and NiFe2O4 in each case,46 the values of Kcore,IF and Kshell,IF are 13 times larger than the core bulk values due to interface effects (proximity effects, local frustration) expressing the strong coupling between the two phases. In the case of pure CoFe2O4 and NiFe2O4 nanoparticles, Ksurf is taken to be 13 and 20 times larger than the corresponding bulk values of the core material respectively due to surface effects, in agreement with our experimental estimations of the enhanced effective anisotropy, Keff, of the samples. Accordingly, for the NiFe2O4 core/CoFe2O4 shell, the Ksurf is taken to be 1.2 times larger than the hard Kshell. For the CoFe2O4 core/NiFe2O4 shell, Ksurf is taken to be 200 times larger than the soft KV,shell because of the proximity effect with the hard core which induces a higher anisotropy in the nickel ferrite when it grows on more magnetically hard cobalt ferrite seeds, justifying the high experimental coercive field value. The values of the anisotropy constants in the above regions are also listed in Table 1 for the corresponding systems.
In our simulations, we divide the energy terms of the four systems with the energy term kcr = Kcore·Vcore, where Kcore is the bulk CoFe2O4 core anisotropy, so our energy terms are dimensionless. In our model, in each nanoparticle, Heisenberg exchange interactions have been considered between nearest neighbour regions, namely: (a) between the bulk-core spin and the interface spin (jcore-IF), (b) between the interface spin and the bulk-shell spin (jIF-sh) and (c) between the bulk-shell spin and the surface spin (jsh-surf). We assume that these intraparticle exchange strengths are proportional to the surface area of contact between the two regions and not to their volume. Since in the literature there is no calculation for the interface exchange coupling, we consider the exchange coupling jcore-IF as a free parameter and we set it equal to two for CFO and equal to 1.25 for the NFO system. We calculate the other two exchange parameters by multiplying each of them by the ratio of the shell and the surface respectively to the core-interface region. Thus jIF-sh = 2.22 and jsh-surf = 6.34.
Moreover, since the particles are coated with an oleic acid surfactant, they interact in the assembly only via long-range dipolar forces. Each spin in a particle interacts with all the spins in all the other particles in the assembly. The reduced dipolar interaction strength is calculated using the expression g = μ0Ms,core2/24Kcore (Dcore/α)3 where Ms,core is the core saturation magnetisation MCFOs,core = 0.45 MA m−1 and MNFOs,core = 0.33 MA m−1.46 For a random monodisperse dense assembly α ≈ dparticle so for the pure CoFe2O4, CoFe2O4/NiFe2O4 core/shell, pure NiFe2O4 and NiFe2O4/CoFe2O4 core/shell the calculated dipolar energy values are 0.0236, 0.0044, 0.0197, 0.0036 respectively.
The external magnetic field is defined as h = (μ0Ms,core·Ncore/Kcore·Vcore)·H and the thermal energy is denoted by t = kBT/Kcore·Vcore. Including the energy parameters of the anisotropy, the nearest neighbour Heisenberg exchange interactions between the spins in each particle, the dipolar energy and the Zeeman energy terms the total energy of the system is given by the equation:
(4) |
It has been demonstrated in the literature that the saturation magnetisation of the ferrite nanoparticles is reduced due to the spin canting effect. This spin canting decreases with the increase of the temperature47 or the increase of the external applied field24 and depends on the size and the method of production of the nanoparticles.47–51 Thus, in our calculations of the low temperature hysteresis loops, spin canting effects are introduced by estimating a 15% and 30% reduction of the surface magnetization for the pure CoFe2O4 and NiFe2O4 nanoparticles respectively. Accordingly, in the case of CoFe2O4 core/NiFe2O4 shell particles, we consider a reduction in magnetisation of 15% at the interface and 30% at the surface. However, in the case of NiFe2O4 core/CoFe2O4 shell particles, we assume a very small magnetisation reduction at the interface (0.5%) and the surface (1%) because the large hard shell (82% volume ratio) dominates over the spin canting effects at the interface with the soft core material and at the surface.17
Monte Carlo simulations using the standard Metropolis algorithm52 were performed to calculate the hysteresis loops, the demagnetization DCD (mDCD) and isothermal remanence IRM (mIRM) curves at a low temperature (t = 0.01) following the same procedure as in the experimental situation as described in Section 2.2. For the dipolar energy calculation, the Ewald summation technique has been implemented taking into account the long-range character of the dipolar interactions, using periodic boundaries in all directions and Dij is the dipolar interaction tensor.25,53 The Monte Carlo simulations result for a given temperature and applied field were averaged over 10 samples with various spin configurations, easy-axis distribution and different spatial configurations for the nanoparticles.
Sample | Composition | d C, nm | a, nm | V hard/Vd | ||
---|---|---|---|---|---|---|
STEMb | SDb | XRDc | ||||
a The error in the last digit is reported in parenthesis. b Median diameter of the particle (dC) and standard deviation (SD) of best-fit log-normal size distribution. c Crystallite size calculated using Scherrer formula for [400] reflections. d The volume of phases was calculated by assuming spherical particles. | ||||||
CFO | CoFe2O4 | 6.9(1) | 0.9 | 6.2(2) | 0.840(1) | 1 |
CFO/NFO | CoFe2O4@NiFe2O4 | 11 | 3 | 9.2(4) | 0.838(1) | 0.3 |
NFO | NiFe2O4 | 8.4(2) | 1.8 | 6.2(4) | 0.838(1) | 0 |
NFO/CFO | NiFe2O4@CoFe2O4 | 12 | 2 | 9.0(2) | 0.837(1) | 0.6 |
The analysis of STEM images confirms the presence of the rounded particles with high crystallinity (Fig. 2). Median particle diameter (dc) and standard deviations (SD), calculated by the log-normal fitting of diameter distribution with eqn (3) (Fig. 2e, f and Table 2), indicate, within experimental error, that the sizes in both the single phases and in both core/shell structures are comparable. It is worth underlining that both core/shell samples are about 2 nm larger than the corresponding seeds. In conjunction with the XRD results above, the STEM data confirm the growth of a ∼2 nm shell around the seeds. Additional confirmation of the core/shell structure is provided by TEM energy-filtered elemental maps (Fig. S3†). Also, an increase in SD was observed after the treatment of CFO to form CFO/NFO MNPs and, on the contrary, a reduction of SD was observed in the formation of the inverted system, i.e. from NFO to NFO/CFO MNPs.
The field dependence of magnetisation recorded at 5 K (Fig. 3) shows hysteretic behaviour. All the extracted magnetic parameters (i.e. coercivity, μ0Hc; saturation magnetization, Ms; and reduced remanent magnetisation, Mr/Ms) are reported in Table 3. The observed values of Ms (i.e., ∼98 and ∼75 A m2 kg−1 for CFO and NFO) are higher than the bulk values (88 and 55 A m2 kg−1 at 5 K for bulk cobalt ferrite and nickel ferrite, respectively) probably due to some difference in magnetic structure (i.e., cationic distribution and surface effect) at the nanoscale. It is interesting to underline that, for both core/shell samples, the Ms values are very close to the mean of the magnetisation of two materials weighted by the volume composition. The monotonic shape of the hysteresis loops of the core/shell CFO/NFO MNPs, without any anomalies around zero field,55 suggests a strong coupling between the hard core and the soft shell phases.17 In all samples, the reduced remanent magnetisation increases going from the single phase to the core/shell system. CFO show Mr/Ms ≅ 0.47 and an evident increase is observed in CFO/NFO system (Mr/Ms ≅ 0.61) approaching the theoretical value of ∼0.83 expected for single-domain non-interacting particles with cubic magnetic anisotropy.41 The lower value of Mr/Ms of ∼0.21 for NFO (theoretical value of ∼0.5 for uniaxial particles41) can be ascribed to the presence of a fraction of still relaxing particles at 5 K, due to the presence of inter-particle interactions or to some effect due to a demagnetising field.56 A significant increase in Mr/Ms is observed in NFO/CFO core/shell systems. Curiously, the largest Mr/Ms (i.e., the tendency to cubic anisotropy) has been observed for the system with a higher volume of nickel ferrite (sample CFO/NFO), which is expected to have the lower magnetocrystalline anisotropy. Indeed, usually cobalt ions induce dominating cubic-type effective magnetic anisotropy because of strong magnetocrystalline anisotropy57). This can be ascribed to the proximity effects that increase anisotropy in nickel ferrite when it grows on magnetically harder cobalt ferrite seeds due to the induced strong pinning of interfacial spins.58,59 In contrast, a smaller Mr/Ms was observed in the system with the predominant cobalt ferrite volume.
Sample | M s,b A m2 kg−1 | μ 0 H c, T | μ 0 H irr, T | M r/Ms | T B, K | K eff,c ×105 J m−3 | E dip/kB,dK | H INT/Hc (%) |
---|---|---|---|---|---|---|---|---|
a The error in the last digit is reported in parenthesis. b The value of Ms obtained at 5 K with the fitting of high field region with equation M = Ms·(1 − a/H − bH2), the value in brackets is the fitting error. c The Keff is calculated following eqn (5). d The estimated error for the Edip is about 20%. | ||||||||
CFO | 98(5) | 1.3(2) | 3.2(2) | 0.47 | 175(5) | 8.3(2) | 44 | −6 |
CFO/NFO | 83(3) | 0.83(1) | 1.8(1) | 0.61 | 174(4) | 3.9(2) | 198 | −2 |
NFO | 75(3) | 0.025(1) | 0.3(1) | 0.21 | 28(3) | 0.6(1) | 57 | −33 |
NFO/CFO | 86(3) | 0.30(1) | 1.3(1) | 0.45 | 173(4) | 2.9(1) | 295 | −2 |
The coercive field (μ0Hc) and the irreversibility field, (μ0Hirr), (i.e., the point at which the difference between the magnetising and demagnetising branches is below 3% of their maximum value60), of the CFO sample were reduced when it was covered with a softer nickel ferrite shell (CFO/NFO sample). On the contrary, μ0Hc and μ0Hirr of NFO/CFO were increased with respect to NFO seeds (Table 3). Notably, the design of the magnetic nanoarchitecture (i.e. relative position (core or shell) of soft and hard materials) plays a crucial role in their effective magnetic anisotropy. Indeed anisotropy estimated from eqn (1) without consideration of surface and interface effects predicts a different trend.
Fig. 3 shows the Monte Carlo simulation results for the low temperature (t = 0.01) hysteresis loops for the assemblies of CFO and CFO/NFO (Fig. 3c) and the NFO and NFO/CFO (Fig. 3d) systems. The simulated hysteresis loops of core/shell particles are also smooth, indicating a strong exchange coupling between the two phases in agreement with the experimental findings. The proximity effect at the interface between the CFO core and NFO shell enhances the anisotropy of the shell component resulting in a coercive field value higher than that of the NFO/CFO system. On the other hand, the proximity effect of the CFO shell and NFO core enhances the coercive field of the pure NFO sample by a factor of almost 13. Our Monte Carlo simulations are in good qualitative agreement with the experimental results, indicating that the magnetic properties core/shell systems are tuned via proximity effects.
M–H curves measured at 300 K (inset Fig. 3a and b) show completely reversible behaviour (i.e., Mr = 0, Hc = 0) suggesting that all the samples are in the superparamagnetic state at room temperature. To better understand the magnetic properties of these systems, their magnetothermal behavior has been also investigated.
The ZFC and FC magnetisation curves for both single phase and core/shell samples (Fig. S4†) showed typical behaviour of interacting mono-domain particles with different anisotropy (see details in the ESI†). An estimation of the distribution of anisotropy energy (DAE) barriers can be obtained by the −dMFC-ZFC(T)/dT curve reported in Fig. 4 (see ESI† for a detailed explanation of the method).21,61 In Fig. 4 the filled areas under the −dMFC-ZFC(T)/dT curves represent the cumulative DAE (eqn (S1)†). Considering the mean blocking temperature (TB) as the temperature at which half of the particles are in the superparamagnetic state, the values of TB for all the samples are reported in Table 3. For the NFO and CFO samples, the different values of the blocking temperature (∼28 and ∼175 K, respectively) can be ascribed to the different magnetocrystalline anisotropy of the NFO and CFO. On the other hand, a significant increase of TB is observed in NFO after the coverage with the CoFe2O4, the same value of TB within an experimental error was also observed for the CFO and CFO/NFO samples.62
Fig. 4 shows the temperature dependence of the reduced remanent magnetisation estimated from hysteresis loops recorded at different temperatures. The Mr/Ms trend for all the samples substantially confirms the landscape obtained by DAE. The Mr/Ms vanishes at the upper limit of the DAE where the entire system is in the superparamagnetic state for T ≥ Tirr. The only exception is the NFO sample for which Mr/Ms vanishes at lower temperatures, closer to TB.
Following the Stoner and Wohlfarth model63 and assuming negligible magnetic interparticle interactions,36,60,64 the anisotropy constant can be defined as:
Keff = μ0HKMs/2. | (5) |
It can be considered that μ0HK ≈ μ0Hirr, an approximation first suggested by Kodama et al. for NiFe2O4 nanoparticles,65 which has also been used for other similar nanoparticle systems.36,66,67 The trend in anisotropy constants' values follows qualitatively the trend of the coercive and irreversibility fields' values, confirming the efficiency of interface coupling in tuning the anisotropy.
As it is expected, the HCr (Fig. 5a) has the same trend of μ0Hc and μ0Hirr and Keff. It is interesting to note that SFD, normalised in X-axis by corresponding HCr and in Y-axis by Mr (Fig. S6†), shows that the seed-mediated growth of the shell does not generate new factors which would broaden the distribution of magnetic properties. Indeed, the standard deviation of the SFD for both core/shell systems does not exceed those of the corresponding seeds, Table S1†). In addition, it appears that the magnetically softer samples NFO seeds and NFO/CFO core/shell particles have a higher dispersion than that of the CFO sample. This is consistent with the observed higher SD from TEM measurements and the expected stronger influence of the surface effects in magnetically softer nanoparticle samples. On the other hand, by comparing both core/shell systems, one can note that the SFD is broader in the NFO/CFO sample than in the CFO/NFO sample, while the SD from TEM analysis has the opposite trend. This behavior can be explained by the dominant role of the magnetically harder counterpart (CFO core), governing the magnetisation processes in the CFO/NFO system.
According to the Stoner–Wohlfarth model,42,70 for an assembly of non-interacting uniaxial single-domain particles, the remanent magnetisation recorded by DCD and IRM protocols are linearly related (mDCD(H) = 1 − 2MIRM(H)). Two factors may cause deviations from this trend: (i) interparticle interactions and/or (ii) multiaxial anisotropy.36,44 For the analysis, we have used the Henkel plots modified by Kelly (ΔM-plot):71
ΔM(H) = mDCD(H) − 1 + 2 MIRM(H) | (6) |
A negative contribution to the ΔM-plot indicates the dominant demagnetizing character of dipolar interactions between particles; a positive contribution, on the contrary, is caused by the dominant “positive” magnetizing character of the exchange interactions.44 The ΔM-plots (Fig. 5b) in all samples show negative values. For both core/shell systems and CFO samples, the absolute value of the |ΔM| dip is in the range 0.2–0.3 while for NFO samples this value is almost twice as large reaching ∼0.45 indicating that for NFO samples interparticle interactions play a significant role. The interaction field μ0HINT (i.e., (μ0HDCDCr − μ0HIRMCr)/2) shows negative values for all seed and core/shell systems confirming the dipolar interparticle interactions role. The ratio HINT/Hc shows the effect of the interparticle interactions compared with the coercive field (Table 3). These values are relatively low except in the case of the magnetically softer NFO seeds, where the intrinsic anisotropy is very small. The absolute value of dipolar interaction energy (Table 3) has been roughly estimated using Edip ≈ μ0μ2/4πd3 (where d is the average distance, which at first approximation is the diameter of the particles plus the thickness of two monolayers of oleic acid (one on each particle) each of them has thickness ∼2 nm).57 This calculated Edip cannot explain the picture of the interparticle interactions described by the ΔM-plot. We attribute the observed magnetic behavior to the complex interplay between intraparticle and interparticle interactions. This is further confirmed by the calculated SFD and ΔM-plots (Fig. 5c and d), showing very good qualitative agreement with the experimental results (Fig. 5a and b). To determine the relative contribution of the dipolar interparticle interactions and the interface exchange coupling, ΔM-plot in the absence of the dipolar interparticle interactions (g = 0) has been calculated. The results are shown in Fig. 6 for g = 0 (dotted lines) and g ≠ 0 (solid lines) for the CFO/NFO (black) and the NFO/CFO (red) core/shell nanoparticle systems. The results for the g = 0 case, where a small but significant negative ΔM still exists, show the contribution of intraparticle interactions (core/shell structure) in the ΔM-plots in agreement with the findings of ref. 72 and refs therein.
A new 4-spin mesoscopic theoretical model has been implemented to describe interacting core/shell magnetic nanoparticles including surface and interface effects in order to investigate the interplay between intraparticle and interparticle interactions in these systems (Fig. 1). By comparing experimental and Monte Carlo simulations results, we have shown that the magnetic properties of bi-magnetic core/shell nanoparticles assemblies depend both on the intraparticle interfacial exchange coupling and interparticle dipolar interactions and their interplay. The interfacial exchange coupling is the signature of true core/shell nanoparticles and it is manifested in the dramatic changes of the magnetic properties, such as the Mr/Ms and the coercivity Hc, and the smooth hysteresis loops observed after shell growth. The interparticle dipolar interactions affect the switching fields of core/shell nanoparticles in the intermediate range (0.5 to 3 T) much stronger than the interfacial exchange coupling. Interestingly, our study suggests that by tuning the effective anisotropy in bi-magnetic nanoparticle systems via proximity effects, a control over the magnetic interactions between particles can also be achieved.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d1na00312g |
This journal is © The Royal Society of Chemistry 2021 |