Pedro A.
Sánchez
*ab,
Thomas
Gundermann
c,
Alla
Dobroserdova
b,
Sofia S.
Kantorovich
ab and
Stefan
Odenbach
c
aUniversity of Vienna, Sensengasse 8, Vienna, 1090, Austria. E-mail: pedro.sanchez@univie.ac.at
bUral Federal University, Lenin av. 51, Ekaterinburg, 620000, Russia
cTechnische Universität Dresden, Institute of Fluid Mechanics, Dresden, D-01062, Germany
First published on 12th February 2018
Being able to predict and understand the behaviour of soft magnetic materials paves the way to their technological applications. In this study we analyse the magnetic response of soft magnetic elastomers (SMEs) with magnetically hard particles. We present experimental evidence of a difference between the first and next magnetisation loops exhibited by these SMEs, which depends non-monotonically on the interplay between the rigidity of the polymer matrix, its mechanical coupling with the particles, and the magnetic interactions in the system. In order to explain the microstructural mechanism behind this behaviour, we used a minimal computer simulation model whose results evidence the importance of irreversible matrix deformations due to both translations and rotations of the particles. To confirm the simulation findings, computed tomography (CT) was used. We conclude that the initial exposure to the field triggers the inelastic matrix relaxation in the SMEs, as particles attempt to reorient. However, once the necessary degree of freedom is achieved, both the rotations and the magnetisation behaviour become stationary. We expect this scenario not only to be limited to the materials studied here, but also to apply to a broader class of hybrid SMEs.
Current experimental techniques for the synthesis of magnetic gels and elastomers offer different mechanisms for the embedding of the magnetic particles within the polymer matrix. In one of the main cases, the magnetic particles simply become physically trapped into hollow ‘pockets’ of the matrix during the synthesis process. Under these conditions, the particles experience a translational constraint imposed by the ‘walls’ formed by the surrounding mesh of polymers, while keeping a higher orientational freedom in case the volume of the cavity is larger than those of the magnetic particles.1 In the second case, the surface of the magnetic particles is functionalised, allowing the formation of chemical bonds with the polymers. In this way, the particles play the role of network crosslinkers and not only their translations, but also their rotations are constrained by a mechanical coupling with the polymer matrix.7–9 Finally, in the third case there is no chemical binding between the particles and the matrix but the latter is so compact that both translations and rotations might be hindered, especially when the particles are large and/or anisotropic. In all cases, the mechanical coupling between the magnetic particles and the polymer matrix makes the field-induced microscopic displacements and reorientations of the former responsible for the macroscopic deformations of the material.
The macroscopic response produced by the field-induced self-assembly of magnetic micro- and nanoparticles is a phenomenon extensively studied in another class of soft magnetic materials: the dispersions of magnetic particles in inert carrier liquids known as magnetic fluids. It is well known that magnetic particles in dispersion tend to form chain-like arrangements under applied external fields. These chains tend to follow the field lines by keeping a head-to-tail orientation of the magnetic moments of their particles.26,27 The same type of field-induced rearrangement takes place for the magnetic particles embedded in magnetic gels and elastomers,3,28,29 but in this case the mechanical constraints imposed by the polymer matrix may hinder to some extent the particle rearrangements. Basically, the local displacements and reorientations of the particles induced by external fields tend to produce deformations of the polymer matrix with an associated mechanical energy penalty. This phenomenon, known as magnetodipolar striction (MDS),30 is fundamental in the behaviour of these materials. In particular, as long as the elastic limit of the polymer matrix is not exceeded, such deformations tend to be spontaneously reverted after field removal. In this way, these materials can exhibit important shape memory properties.2,28,31
The reversibility of the changes in the shape and mechanical properties of magnetic gels and elastomers is crucial for many prospective applications. For this reason, the study of their behaviour in elastic regimes has been the subject of most research studies to date. However, the characterisation of their irreversible deformations under strong enough external fields is also essential from both perspectives, their fundamental understanding and the design of technological applications. Related to this important aspect, very recent experimental measurements have shown that a novel type of hybrid magnetic elastomer, obtained by embedding a mixture of magnetically soft and hard particles within the polymer matrix, may exhibit certain unconventional features: among other properties that provide important technological advantages, these materials can exhibit a significant nonelastic behaviour when still freshly synthesised. In particular, they display a substantial change in their magnetic response during their first exposure to external fields.32 This is experimentally evidenced by consecutive measurements of conventional magnetisation loops: in some cases a notable difference in the shape of the first hysteresis loop can be observed, particularly in the second half of the backward curve, when compared to the next ones (see, for instance, Fig. 4(c) in ref. 32). In addition, in these magnetically hybrid materials the magnetisation loops tend to show some bias, with different values of the forward and backward coercivities.32,33 A reasonable explanation for this behaviour is that it corresponds to the signature of irreversible internal changes in the material induced by the first-time application of an external field, but the nature of such changes is still unclear. Another important effect might be the complex interplay between the magnetically hard and soft particles. Some efforts to understand this phenomenon that do not take into account inelastic processes have already been made: a very recent theoretical study presented a microscale two-dimensional model of SMEs, based on a magnetomechanical continuum formulation, which predicts a large difference in the first magnetisation loop as a consequence of the strain energy stored in the elastic matrix in the initial magnetisation curve.34 Even this model naturally includes the detailed magnetic interactions between close particles, the predicted difference is much larger than the experimental one, which suggests that the relevant underlying physical mechanisms might not be fully captured by this approach. Therefore, the contribution of eventual irreversible deformations of the matrix to the experimental behaviour remains to date as an open question. Unfortunately, even though it is already possible to observe the deformations of the matrix by means of confocal microscopy techniques,35 the available experimental methods are not accurate enough to provide a direct observation of its irreversible microstructural changes. This makes the study of such irreversible deformations and their effects particularly challenging. Another complication for the study of hybrid SMEs is the limited resolution in the observation of the magnetic particles. Current X-ray tomography techniques are able to image relatively large solid particles, such as the magnetically hard ones used in these materials. The imaging of the smaller, magnetically soft particles used in the aforementioned hybrid SMEs, with typical sizes not larger than 5 μm, requires an increase of their resolution power that is still under development. Alternatively, one could also create hybrid SMEs using magnetically soft particles large enough to be observable with available experimental methods, but in this case it would be hard to distinguish them from the magnetically hard ones.
The difficulties discussed above suggest the use of a reductionist strategy in order to unveil the role of inelastic matrix deformations in the distinctive first magnetisation loop originally observed for hybrid SMEs. With this goal, here we study experimentally and theoretically a simpler SME material, composed of only relatively large magnetically hard particles. In this way, we avoid the complex magnetic interactions between heterogeneous particles to focus on the analysis of elastic and inelastic microstructural deformations induced by the external field. This analysis includes several results. First, we provide the first experimental evidence of simple SMEs also showing a difference in their first magnetisation loop, obtained from measurements using a vibrating sample magnetometer (VSM), proving that this behaviour is not exclusive of hybrid materials. Second, we study by means of computer simulations with a minimal model of SMEs if irreversible deformations of the matrix alone can reproduce the experimental loops, showing that a good qualitative agreement is obtained when such deformations increase both the translational and rotational freedom of the magnetic particles. Finally, we characterise experimentally the movements of the particles as a response to the applied field using micro-computed tomography imaging (Xμ-CT), showing that they are consistent with the theoretical assumptions.
This paper is organized as follows. First, in Section 2.1, the experimental approach is presented. Section 2.2 describes the theoretical modeling approach. The results are presented in Section 3. Finally, conclusions and outlook are discussed in Section 4.
Considering the SME Material Xμ-CT is a tool to analyze particle structures36,37 up to a single particle detection29,35,38 inside such SME materials. This leads to statements about particle distribution39 and particle motion under the influence of an external magnetic field.29,35,38
Of great importance when considering such SME materials is the effect of the microstructural changes on macroscopic behaviours such as the mechanical response and the magnetic properties of the material. To establish this relationship, magnetic measurements of the samples using a vibrating sample magnetometer (VSM) have been carried out with respect to the CT-scans to enable a physical description of the material. Particular hysteresis effects which were already found in a previous study32 are to be explained by using this method. For this purpose, a VSM Lake Shore 7404 was used to perform the magnetic measurements.
Three different samples containing 5 vol% of magnetically hard NdFeB microparticles dispersed isotropically in an elastomeric matrix ELASTOSIL RT 623 A/B by Wacker Corp. Germany were produced. The samples were of a cylindrical shape with a diameter of 4 mm and a height of 4 mm. A silicon oil with a viscosity of 100000 Pa s was used to adjust the elastic modules and to avoid sedimentation. Thereby a sample with a tensile module of E = 30 kPa (SME1), a sample with a tensile module of E = 80 kPa (SME2) and a sample with a tensile module of E = 900 kPa (SME3) were produced. In contrast to the work in ref. 32, no magnetically soft particles were added to the samples.
To carry out the CT investigations we used a home build Xμ-CT System based on a nano-focus tube with a maximum acceleration voltage of 160 kV and a detector with a photodiode array of 2940 × 2304 (horizontal × vertical) pixels. The sample was placed on a movable sample holder between the X-ray tube and detector. The sample holder consists of two cylindrical permanent magnets which were mounted displaceable and could generate a homogeneous field parallel to the direction of gravity in the range of B = 0–400 mT. This enables investigations of cylindrical samples with a height of 4 mm and a diameter of 4 mm within a magnetic field with a homogeneity of more than 95% in a radial and axial direction.40 The temperature during the CT investigations was kept constant at 20 °C and projection images were generated by rotating the sample with a 0.25° angular increment for a tube current of 170 μA and an acceleration voltage of 90 kV. The exposure time was varied between 0.8 s and 1.0 s to achieve the highest image quality. The magnification was 20, which results in a resolution of 1 pixel = 2.5 μm. The reconstruction of the CT scans was carried out using home-made software. The evaluation of the reconstructed images was carried out by using MATLAB 7.10.0 software within the toolbox DIPimage. In order to separate the particles from the surrounding matrix, a threshold was set with respect to the known particle volume concentration in the sample. After separation, the particles were identified individually by using a watershed algorithm.41 This separation of clustered particles leads to the possibility to calculate the size, surface area, position and main axes, including their directions, for each single particle. An example of the evaluation process results is shown as an image of the separated particles in a reconstructed CT measurement in Fig. 1(a). One can see that the particles have a flake like shape and a high polydispersity. The mean size of the particles was determined by sieving and has resulted in a range between 10 μm and 100 μm. Due to the possibility of calculating the size of the particles, a size distribution could be determined by an approximation of the particles by a sphere. The size distribution is shown in Fig. 1(b) and it confirms the size-range that was estimated from the sieving. The image processing allows us to calculate the direction of the major axes relative to a given vector. In our CT setup, the magnetic field is always vertical. The rotation of the particles is given by the angle between the longest major axis and the direction of the field. In order to analyse the reorientations of the particles during magnetisation–demagnetisation, CT scans at different points of the hysteresis loops were performed.
Here, we model a simple magnetic elastomer system filled with magnetically hard particles, incorporating a novel minimal representation of the effects of inelastic deformations of the polymer matrix on the magnetic response.
(1) |
In case a uniform external magnetic field, , is applied, every particle i also experiences the effect of the Zeeman interaction:
UH(i) = −i· | (2) |
The steric repulsions between the magnetic particles are modelled by a truncated and shifted Lennard-Jones (LJ) interaction, also known as Weeks–Chandler–Andersen (WCA) potential.59 This is defined as:
(3) |
Finally, we chose a minimal implicit representation for the constraints imposed by the polymer matrix on these magnetic particles: for any system formed by N magnetic beads, we assume every one of them to be trapped by the matrix and mechanically forced to remain within a limited region of the space around its initial position, corresponding to a local elastic equilibrium. For simplicity, such initial positions are placed randomly in space with a given number density, ρ. At this point, in order to complete the description of our modelling approach, we have to distinguish between the two aforementioned mechanisms of coupling between the particles and the matrix.
(4) |
So far, we have defined the models for the elastic behaviour of the mechanical constraints affecting the magnetic particles in a elastomer. The last ingredient in our minimal modelling approach is to take into account the inelastic deformations induced by the application of external fields.
For M1, every time the distance between the centre of a particle and its anchoring point, r, becomes larger than a given threshold, r > rth, a shift of the anchoring point towards the particle centre is performed with a certain arbitrary probability,
PshiftM1 = ps, | (5) |
(6) |
Fig. 3 Inelastic deformations for models M1 (a) and M2 (b), represented as shifts of the anchoring points of the springs that reduce their stretching. |
Even based on the same strategy, inelastic relaxations of the rotational constraints in model M2 require a slightly more sophisticated approach. In this case, the criterium to attempt an irreversible change in the constraints of a given particle is a continuous distribution that grows with the misalignment of the dipole with respect to the orientation corresponding to the mechanical equilibirum of the springs. This misalignment is measured by θ, that is defined as the angle between the dipole moment of the particle, , and the displacement vector between the anchoring points of the springs corresponding to the tail and the head of the dipole, = head − tail, as depicted in the scheme of Fig. 3(b). In general, angle θ will tend to grow statistically with the strength of the external field, H = ‖‖, due to the torque produced by the Zeeman interaction on the magnetic dipole of the particles, = × . We take the modulus of this torque, nomalised with the saturation field strength, Hmax, as the main part of the probability distribution for performing an irreversible relaxation of the constraints, with an added arbitrary prefactor, ps, as free tuning parameter. Therefore, the final expression is
(7) |
(8) |
In molecular dynamics simulations with coarse-grained models like these, it is usual to use a Langevin thermostat in order to approximate statistically the effects of the degrees of freedom neglected by the coarse-grained approach,62 particularly the thermal fluctuations of the background fluid. Such fluctuations are very relevant in magnetic fluids and gels, in which the background fluid is a liquid. In elastomers, that are dry materials, thermal fluctuations at room temperature can be considered negligible in front of the magnetoelastic interactions. Since the latter are the only focus of our study, we basically need to perform athermal simulations. However, it might be still useful to keep some little thermal fluctuations in the simulation of these systems to favour their mechanical relaxation, because the strongly simplified representation of the effects of the polymer matrix, that ignores any long range mechanical coupling, makes them prone to kinetic trapping. Therefore, we chose to use a Langevin thermostat with a relatively low temperature.
In our simulations, we used an arbitrary system of reduced units. Note that the qualitative nature of our analysis makes the values we assign to the different parameters irrelevant when considered individually. Only their ratios, that determine the relative strengths of the interactions, are relevant for our purpose. In particular, we chose to measure lengths in units of the diameter of the particles, σ = 1. The energy scale of the magnetic interactions is defined by an arbitrary value of the minimum sampled Zeeman energy of each particle, defined by (2), that corresponds to the maximum absolute value of the saturation fields we apply during the magnetisation loops and the magnetic moment assigned to the dipoles. The chosen values for such magnitudes are, respectivley, UminH = −60, Hmax = 30 and μ = 2. For all sampled systems a random elastic constant, Ki, is assigned to each individual (M1) or pair (M2) of springs that anchor every particle i according to (4). This random value is taken from a Normal distribution with mean α and standard deviation s2, N(α,s2), that is defined to span over a fixed arbitrary interval [Kmin, Kmax]. This is carried out by taking α = (Kmin + Kmax)/2, s = (Kmax − Kmin)/6 and rejecting all values that do not lie within the chosen interval during the random assignation procedure. In any case the minimum and maximum values of Ki we sampled are, respectively, Kmin ≥ 1 and Kmax ≤ 30. For the case of the model M2 systems, we chose the initial separation between anchoring points, defined for expressions (7) and (8), to be || = 5. We chose the thermal energy in our simulated systems to be much lower than those associated with the magnetoelastic interactions by taking a reduced temperature of T = 10−4. Finally, for the sake of simplicity, we chose the prefactor of the soft core steric interactions defined by (3) to be ε = 1. This latter choice is loosely connected to the rest of the energy scales we defined above, since the purpose of interaction (3) here is only to keep a minimum centre-to-centre distance between particles not much lower than their diameters, σ = 1. Therefore, any higher value for ε would also serve this purpose.
The Langevin dynamics method adds to the Newtonian equations of motion stochastic and friction terms that satisfy the usual fluctuation–dissipation relationships. In particular, the Langevin translational and rotational equations of motion governing the dynamics of a particle i are
mi(di/dt) = i − ΓTi + i,T | (9) |
i·(di/dt) = i − ΓRi + i,R, | (10) |
Fig. 4 shows the experimental magnetisation loops obtained for three different samples of magnetic elastomer filled with a 5% in volume of NdFeB microparticles. The measurements were performed right after their synthesis. Fig. 4(a) corresponds to a sample with a tensile modulus of E = 30 kPa (SME1), whereas one of the samples in Fig. 4(b) (SME2) is E = 80 kPa and one of the samples in Fig. 4(c) (SME3) is E = 900 kPa. In order to ease the comparison with the simulation results, fields and magnetisations in these figures are normalised with their corresponding values at saturation, Hmax and Mmax respectively. For all samples the first magnitude was approximately Hmax ≈ 1660 kA m−1, whereas the saturation magnetisation of samples SME1 and SME3, Mmax ≈ 38 kA m−1, was slightly lower than that corresponding to sample SME2, Mmax ≈ 50 kA m−1.
In the case of SME1 and SME2 the area enclosed by the hysteresis loops is rather small due to the moderate stiffness of the matrix. For SME3 the area is larger due to the higher stiffness of the matrix. It can also be observed that the coercive fields corresponding to the upward and backward curves are not symmetric. This is caused by the internal interaction field due to the magnetisation of the particles.68 Importantly, all samples show a significant difference in their first magnetisation loop with respect to the next ones, being particularly clear in the second half of the backward curves.
Also interesting is the fact that such a difference is more clear for sample SME2. The reason for this is the interplay between the particle–matrix coupling and the stiffness of the matrix itself. If the stiffness of the matrix is very low (SME1) the restoring forces of the matrix acting on the particles are low. The particles are able to rotate almost freely when an external magnetic field is applied. The latter results in a weaker asymmetric behaviour. If the stiffness of the matrix is too high (SME3), the particles are strongly constrained by the matrix and the rotations are hindered, causing a disappearance of the asymmetric behaviour. As a consequence, the largest asymmetry is observed for sample SME2, for which strong field induced rotations are possible. In this case the rotations might require a partial destruction of the coupling between the particles and the matrix.
For comparing the experimental results of the magnetisation and the simulation results, only SME1 and SME2 have been used due to the fact that almost no asymmetry was found for sample SME3.
First, we present the simulation results corresponding to the case when no inelastic constraint relaxations are allowed, i.e., when a fully elastic regime is simulated. Under these conditions, no difference between the first and next magnetisation loops has been observed in any simulated system. Inelastic deformations will be discussed afterwards.
Finally, some examples of the magnetisation measurements obtained for model M2 under different concentrations and spring rigidities are shown in Fig. 6. In this case, despite the maximum strength of the reduced field we sampled, Hmax = 15 was much larger than the one applied to M1 systems, the effect of these elastic constraints is so strong that a significant hysteresis is only observed for M2 systems with the lowest sampled rigidities, K ∈ [1,10].
After this overview of the elastic response of elastomer models M1 and M2, we can now discuss the effects of inelastic deformations in the next section.
Fig. 7 shows the simulation results of the first three normalised magnetisation loops obtained for M1 systems with K ∈ [1,10], ps = 0.2 and four different values of ρ. One can see that a notable difference in the first magnetisation loop with respect to the next ones only occurs for the lowest concentration: as ρ is increased, such differences tend to vanish and the loop profiles tend to have a large enclosed area, as corresponds to the low rigidity of the constraints. Significant differences for the first loop were not observed in any of the sampled M1 systems with higher rigidities. Moreover, a simple visual comparison of the difference observed at low ρ and K (Fig. 7(a)) shows that it is rather distinct from the ones obtained in the experiments (Fig. 4). Therefore, we can conclude that the model M1 does not provide a reasonably good qualitative representation of the experimental behaviour, at least within the broad ranges of the parameters we sampled.
The magnetic response of model M2 turns out to be much more affected by the inelastic deformations than that of M1. Fig. 8 and 9 show four selected results of the first three normalised magnetisation loops obtained from the M2 samples with ps = 0.002 and different combinations of ρ and K. In contrast to Fig. 6, here in all cases clear hysteresis loops are obtained. Importantly, a clear difference in the first magnetisation loop with respect to the rest can be seen. This difference follows qualitatively quite well the experimental behaviour: the backward curve of the first loop shows moderately higher values of magnetisation than the subsequent ones along its lower half. Results for higher irreversible shift probabilities (not shown) indicate that the difference in the first loop tends to decrease as ps increases. One can also observe that the area enclosed by the second and third loops, as expected, tends to decrease as the rigidity of the springs increases, but in this case the latter does not hinder the difference in the profile of the first loop. Three configuration snapshots included for the system with ρ = 0.1, K ∈ [20,30] (Fig. 9(a)) illustrate the evolution of the constraining springs during the first loop: the initial configuration (top left) shows the random orientation of particles and springs and the initial stretching length of the latter; at the saturation field of the initial magnetisation curve (top right), the distance between the anchoring points of most particles has decreased significantly due to their irreversible shifts, allowing the particles to highly align with the field and form chain-like aggregates; finally, after the first reverse saturation field is reached and the irreversible shifts are not allowed anymore, almost all pairs of anchoring points have largely reduced their separation to practically zero (bottom left), showing a configuration very similar to that observed for M1 under saturation. This evolution is better seen in Fig. 9(b), that shows an example of the change in the average distance between the anchoring points of the springs of each particle, 〈‖‖〉, as a function of the measured steps along the first magnetisation loop. Here we can see that during the first magnetisation of the system 〈‖‖〉 decreases notably, in a moderate way at the beginning and rather strongly afterwards, so that at the first saturation point its value is approximately 10% of the initial one. Finally, 〈‖‖〉 keeps decreasing during the first field inversion, reaching an almost zero value at the first inverse saturation point. After this, the system has become essentially identical to an M1 sample. According to what we found for M1 systems, the existence of eventual further shifts would not have any impact on the profile of subsequent magnetisation loops. A similar behaviour has been observed for any system that reproduces well the experimental difference in the first magnetisation loop.
Our results show that, despite its simplicity, model M2 with irreversible deformations represents rather well the experimental behaviour. This suggests that the translational and rotational constraints imposed by the polymer matrix on the embedded magnetic particles, as well as their irreversible relaxation due to field-induced deformations, play an important role in the difference observed experimentally for the first magnetisation loop of magnetic elastomers.
Below we verify this conclusion experimentally by means of direct observations of microstructural changes.
Fig. 10 Schematic representation of the different stages in which tomographic investigations of the sample have been undertaken. The different stages are highlighted in different colors with respect to the detailed views of the CT-scans in Fig. 11. |
Fig. 11 Detailed views of CT Scans at different points of magnetisation which are explained in Fig. 10. |
To ease the visualisation of the field-induced microstuctural changes, in Fig. 11 only a small part of the sample is shown. Applying a field of B = 400 mT and comparing this to the CT scan of the unmagnetised sample, the position of the particles changes only slightly, as seen in Fig. 11(a). This evidences a very strong initial mechanical coupling between the matrix and the magnetic particles. However, after the first positive saturation of the sample, as one can see in Fig. 11(b), most of the particles manage to align their longest major axis perpendicular to the magnetic field. In other words, the saturation field is strong enough to reduce the aforementioned coupling. It is worth mentioning that earlier X-ray diffraction measurements35 have shown that the easy axis of these magnetically hard particles is the shortest major axis, so that the particles tend to rotate to orient their longest major axis perpendicular to the magnetic field. Particle rotations during the first inversion of the field, are again hindered by the matrix, that has to further deform in order to allow reorientation of the particles. From this point on, the particles are unconstrained enough to follow the field without the need of further irreversible deformations of the surrounding matrix, i.e., field-induced rotations will only produce elastic deformations.
In order to provide a quantitative analysis of the reorientation occurring between points 2–6 and 3–7, visualised in Fig. 11(h) and (i), we calculate the distributions of angles between particle long major axis and the field for all CT scans, averaged over approximately 45000 particles, and plot them in Fig. 12. In our case, the angle of 90° corresponds to the alignment of particle magnetisation and the external field. One can see, that the initial distribution is basically uniform, with a slight bias towards perpendicular orientation of particles with respect to the vertical axis. Had the magnetisation process been totally reversible, the distributions of the angles for point 2 (green pointing up triangles in Fig. 12) and 6 (pink pointing right triangles in Fig. 12) and also for point 3 (yellow pointing down triangles in Fig. 12) and point 7 (teal circles in Fig. 12) would have been identical. However, we notice a clear difference between 2–6 and 3–7. These discrepancies can be attributed to the fact that the interaction forces between the matrix and particles during the first hysteresis loop are very high, so that some of the connections are destroyed. After that, the magnetisation will follow the same final hysteresis loop. In order to confirm this, CT-scans during the third and fourth magnetisation loops were carried out and no further evolution of the angular distribution was found. This additionally supports the protocol used in our computer simulations, in which the irreversible deformations were allowed only along the first field cycle.
Fig. 12 Distribution of the angles between the major axis of the flake like NdFeB particles and the magnetic field direction for the different magnetisation states shown in Fig. 10. |
In order to understand the microstructural mechanism behind this behaviour, we put forward a minimal simulation model of these materials that takes into account the rotational and/or translational constraints imposed by the matrix on the movement of the particles, as well as their relaxation due to field-induced irreversible deformations. Our simulations show that only when both types of constraints are acting and allowed to relax during the first magnetisation loop, the experimental difference is reproduced. This suggests that irreversible deformations of the matrix play a major part in this behaviour.
In order to confirm the conclusions drawn from the theoretical approach, we performed CT scans at different field strengths and various magnetisation loops. Direct observations of the changes in the microstructure allowed us to quantify the angular distributions of the particles with respect to the field direction at different stages of the magnetisation measurements. These distributions revealed a clear difference in the rotation of the particles for the first magnetisation loop: initially, particles turned to be strongly constrained by the matrix; after the first field inversion, only partial reorientation of the particles was observed; after the second field inversion and subsequent saturation of the samples, the rotations induced by further field cycles become fully reversible and the profile of the magnetisation loops does not change anymore. In other words, even though CT scans cannot directly provide the information about irreversible deformations of the matrix, the indications found from angular distributions and magnetisation measurements, combined with computer simulations, provide a clear picture of their role.
This study is the first step towards the profound understanding of the complex magnetic response of soft magnetic elastomers, that is totally essential for their technological applications. Besides that, being able to predict particle–particle and particle–matrix interactions in these systems, depending on composition and external stimuli, will allow the creation of materials with tunable macroscopic behaviours. That is why our next step is to clarify the part of the unavoidable particle polydispersity and shape anisotropy in the properties of SMEs by extending our experimental and theoretical approaches.
This journal is © The Royal Society of Chemistry 2018 |