1. Introduction
Alginate is a polysaccharide extracted from the cell walls of marine brown algae; it is widely used in various biomedical applications, healthcare, and food products. Detailed information on alginate production, properties, and applications can be found in numerous books and reviews, for instance, see references [
1,
2,
3,
4,
5]. Polyelectrolyte behavior is the main feature distinguishing alginate from other natural polysaccharides like agarose, starch, guar, etc. Each monomeric unit in an alginate polymeric chain holds a carboxyl residue normally neutralized by a cation of a metal. If the polymeric alginate salt is formed by a monovalent cation such as Na+ or K+, alginates can dissociate in water providing highly viscous solutions. The alginates of divalent cations such as Ca++, Mg++, etc., however, are insoluble and form elastic hydrogels. In this case, the divalent cation is coordinated by several monomeric units of adjacent macromolecules and forms a crosslink of the gel networking [
1,
2,
3,
4,
5].
The mono/di valent ionic exchange provides the versatility of alginate applications. Thus, the fluidic water solution of sodium alginate can be readily transformed to a certain shape of an elastic gel of calcium alginate by exposure to the solution of simple calcium salts [
2,
3,
4,
5,
6]. This feature is widely used in the encapsulation of different bioactive substances in alginate beads for their further in vivo transportation in various medical applications [
7,
8].
Alginate gels are extensively tested as scaffolds for cell proliferation and growth [
9,
10]. In these applications, the fluid-to-gel transition of alginate driven by the mono/di valent ionic exchange is especially useful for the in vivo tissue regeneration and engineering [
11]. Fluidic sodium alginate solution can easily be injected into a living organism and the gelation might then be provided in the desired place of the interior. In this regard, the modification of alginate gels and solutions with dispersed magnetic particles gives extra options to control the rheo-elastic properties both for the fluidic and gel forms of alginate [
12]. Several studies on the structure and magneto-rheological properties of alginate ferrogels had been reported [
13,
14,
15,
16,
17,
18]. Thus [
14,
15], ferrogels filled with iron micron-sized magnetic particles were prepared by the Na+/Ca++ ionic exchange in magnetic suspensions in alginate solution according to the “internal gelation” procedure, and their mechanical properties were tested using dynamic mechanical analysis (DMA). It was shown that Ca-alginate ferrogels had substantially higher levels of the storage (elastic) modulus compared to the loss modulus within the linear viscoelastic region, which is the typical behavior of the cross-linked polymeric networks. The values of viscoelastic moduli intensively enhanced if the concentration of magnetic particles increased by up to 30% (vol). The application of the magnetic field with up to 300 kA/
v intensity considerably strengthened the ferrogels: for instance, the storage modulus increased by seven-fold approximately. DMA study was also performed on Ca-alginate ferrogels with rod-like composite particles with magnetite shells [
17] and the same magnetorheological effect for the influence of the applied magnetic field on the viscoelastic moduli was revealed. The enforcement of the gel network by the embedded particles was reported for Ca-alginate ferrogels with super-paramagnetic iron oxide particles [
16].
The properties and behavior of alginate ferrogels are the implicit functions of their internal structure which is governed by multiple types of molecular non-valent interactions among the constituents, which include the polymer network, the particles, and the water medium. Whereas the interaction of alginate with water might be considered the same for all alginate-based ferrogels, the interaction between alginate and embedded magnetic particles might be different depending on the chemical nature of the particles. It puts forward a question of whether interaction among particles and alginate would notably affect the viscoelastic properties of alginate ferrogels. In the present study, we have examined this opportunity for the viscoelastic suspensions of magnetic particles in sodium alginate solutions.
We have taken four different types of submicron-sized magnetic particles: magnetite, metallic iron, metallic nickel, and metallic nickel with deposited carbon shell and used them for the preparation of magneto-rheological suspensions (MRS) based on sodium alginate solution in water. Viscoelastic moduli were determined for these suspensions as a function of the content of the embedded magnetic particles. Using the thermodynamic approach based on experimental calorimetry studies we have evaluated the enthalpy of interaction between alginate polymer and these types of magnetic particles. It will be shown below that interaction between alginate and particles strongly affects the rheological behavior of the suspension and can even cause its gelation.
3. Results
As the first step in the dynamic mechanical analysis of alginate-based magneto-rheological suspensions, we tested water solutions of sodium alginate in concentrations 5, 10, and 15% (wt.). Frequency dependence experiments in the frequency range 0.1–100 Hz at 5 Pa stress showed that
G′ and
G″ were linearly increasing in double log scale, with G″ > G′ and their cross-over at approximately 80 Hz. Based on these results, the frequency of 1 Hz was chosen as the level for the study on stress-dependent rheological properties for MRSs. Typically, the process of gelation in polymer solutions is examined at a low frequency of the oscillatory stress and we took 1 Hz as a close value to a lower limit in frequency measurements (0.1 Hz).
Figure 2 presents the shear stress dependence of storage modulus (
G′), loss modulus (
G″), and shift angle (
δ) between the applied stress and the strain in DMA measurement at the frequency of oscillations 1 Hz for two concentrations of Na-alginate (5 and 10%). The dependences for the 15% solution were very much like that for the 10% solution but absolute values of moduli were higher by approximately an order of magnitude.
In the range of shear stress, 1–50 Pa the values of dynamic moduli were constant for both concentrations, except the storage modulus for 5% solution. In that case,
G′ kept constant in the range of shear stress 1–8 Pa and gradually diminished at higher values. According to common settlement [
24] the plateau on the shear stress dependence of modulus is an intrinsic feature of the gel-like behavior of polymeric systems. In the case of Na-alginate it, in general, indicated the presence of physical gelation due to the formation of H-bonds.
Meanwhile, for all concentrations of Na-alginate, the values of loss modulus (
G″) were substantially higher than the values of storage modulus (
G′). It is clear from
Figure 3, which presents the dependences of
G′,
G″,
δ on Na-alginate concentration at 1 Hz frequency and selected shear stress
τ = 1 Pa. The values at 1 Pa were taken as characteristics for the modulus plateau.
The prevalence of G″ over G′ means that the balance in the viscoelastic behavior of the Na-alginate solution is biased to fluidic properties. It is clearly indicated by the value of the shift angle. For the idealized elastic material, the strain is in-phase with the stress and the shift angle is 0, while in the idealized fluid, it is 90 degrees. The shift angle was 85 degrees for the solution with 5% Na-alginate and it indicated that the viscous flow dominated the dynamic deformation. As the shift angle went down to 72 degrees for 10% solution and to 57 degrees for 15% solution, it meant that the elastic contribution to deformation progressively increased with the concentration of Na-alginate. It looks reasonable, as the upsurge in concentration had to result in the increase of the probability of inter-chain contacts, which were the origin of physical networking due to H-bonding.
The dependences for
G′ and
G″ were found linear in the double-logarithmic scale. It means that
G′ and
G″ had a power dependence on the concentration of Na-alginate in water (
C). It led to the following scaling equations:
The power exponents in Equation (1) are substantially larger than the classical exponent 1/3 in Flory-Rehner theory for the elasticity of polymeric networks [
25]. Most likely it is the consequence of the dominance of the fluidity of the solutions. Although the elastic contribution to the deformation of Na-alginate solutions is clearly present, it is still less pronounced than the contribution from the viscous flow. In other words, the solutions of Na-alginate are more “liquids” than “elastomers”. At the molecular level, it probably means that the H-bonds, which presumably act as the cross-links in the physical networking of the solution, are unstable and they constantly appear and dissolve in a stochastic fashion. It is also likely that the polyelectrolyte nature of Na-alginate contributes considerably to the fluidity of the solution. Basically, interaction among negatively charged ionized base-units of the Na-alginate chain is repulsive and it hampers the formation of stable cross-links in the physical network of inter-chain H-bonds.
Let us now turn to the dynamics of magneto-rheological suspensions (MRS) based on Na-alginate solutions with dispersed magnetic particles (MPs). For technical reasons, MRS were prepared using a 5% solution of Na-alginate only. The viscosity of 10% and 15% solutions was too high to provide a uniform mixing with MNPs.
Figure 4 presents the shear stress dependences of storage modulus (
G′), loss modulus (
G″), and shift angle (
δ) at 1 Hz for MRS based on 5% Na-alginate solution with 2.8% (wt.) of Fe
3O
4 MPs. The dependence of all these dynamic characteristics was essentially the same as that presented in
Figure 2 for the 5% Na-alginate solution without MPs. The storage modulus had a plateau at a low level of shear stress and moderately diminished with its increase. The loss modulus preserved an almost constant value and was substantially higher than the storage modulus. The shift angle was almost constant, and it was close to 90 degrees. It all meant that MRS with 2.8% (wt.) of Fe
3O
4 MPs demonstrated the same liquid-type behavior as the carrier solution of Na-alginate.
Similar features were observed for suspensions with all other MPs, which contained 2.8% (wt.) of MPs. The dependences of G′, G″, δ were qualitatively very much like that for the initial 5% Na-alginate solution. However, the situation changed drastically if the content of MPs in suspensions was enlarged.
Figure 5 presents the dependence of the dynamic moduli (
G′ and
G″) on the volume fraction of MPs in suspension. The values of
G′ and
G′′ were averaged over the range of their constancy, which was 0.1–1 Pa.
In the case of Fe
3O
4 MPs, both dynamic moduli moderately increased with particle volume fraction. At any concentration, the loss modulus was higher than the storage modulus. In the case of Fe, Ni, and Ni@C MPs there was a steep increase in both moduli if the content of particles exceeded 3–4% (vol.) Note that the ordinate in
Figure 5 is logarithmic and the growth in the moduli was several orders of magnitude. For instance, the storage modulus of 12% (1.7% vol.) suspension of Fe MPs was
G′ = 3.1 Pa and for 30% (5.2% vol.) it had lifted to 42.3 kPa, which was exactly 4 orders of magnitude boost. According to the degree of elevation of
G′ with concentration the MPs fell in a sequence: Fe
3O
4–Ni@C–Ni–Fe.
The loss modulus of suspensions also enlarged with MPs concentration. The influence of the nature of MPs on the degree of elevation was the same as for G′. At the same time the increase in G″ was not as drastic as the increase in G′. As a result, there was a cross-over in the G′ vs. G″ relation at high content of Ni@C, Ni, Fe MPs in the suspension. While at low concentrations of these MPs G″ was larger than G′, at high concentration this relation had reversed to G″ < G′. It meant that at high MPs content (except Fe3O4 MPs) elastic deformation prevailed over the viscous dissipation of mechanical energy.
It can be illustrated by
Figure 6 which gives the dependence of the shift angle on MPs concentration in the suspensions.
It is clear from
Figure 6 that at low concentrations of MPs the dynamic deformation of suspensions was predominantly the viscous flow, as the shift angle was above 80 degrees close to the value of 90 degrees for the ideal fluid. In the case of suspensions with Fe
3O
4 the fluidic nature of deformation preserved up to the highest content of MPs (30% wt., 8.3% vol.). Meanwhile, for the suspensions with other MPs, the turn-over occurred for the high concentration of MPs. The shift angle for the suspensions with Ni@C, Ni, and Fe MPs dropped down to values below 10 degrees which were close to the characteristic shift angle
δ = 0 for the ideal solid. It meant that magneto-rheological suspensions with Ni@C, Ni, and Fe MPs in a concentration above 3–4% (vol.) became elastic gels rather than viscous fluids. Note that even in the fluidic state the effect of the particle concentration on the moduli is much higher than that predicted by the classical Einstein theory of effective viscosity of suspensions.
According to the degree of the descent of the shift angle, the suspensions form a sequence Fe3O4–Ni@C–Ni–Fe, which is identical to that marked above concerning G′ data.
What is the underlying reason for the transition of fluid to elastic gel to occur in the case of suspensions with Ni@C, Ni, and Fe MPs, and not to happen in the case of suspension with Fe
3O
4? Our hypothesis was that its background is in different interfacial interactions in these suspensions. To check this supposition, we studied the interfacial adhesion of Na-alginate polymer to the surface of MPs using the thermodynamic approach, which was successfully elaborated for this purpose in earlier studies [
26,
27,
28].
From the viewpoint of basic thermodynamics, the intensity of the adhesion of polymer at the interface with solid might be characterized by the value of the enthalpy of adhesion (Δ
Hadh), which is the gain in enthalpy if the polymer is put in contact with a solid surface. The thermodynamic equation for this process is:
The heat effect of this process cannot, certainly, be measured directly as both interacting substances are solids. In this case the appropriate thermodynamic cycle is elaborated, which includes the sequence of measurable steps from the initial state to the final state. In the case of composite formation, these steps are based on the dissolution of polymer and composite in a good solvent. Then, water taken as a good solvent for Na-alginate and its composites, the steps of the thermodynamic cycle which corresponds to Equation (2) are:
All these steps can be performed in a calorimeter and the resulting combination of the enthalpy values reads:
It is worthwhile noting that values ΔH1 and ΔH2 relate to individual components but the value ΔH4 for the dissolution of a composite is a function of the polymer content in it. Typically, the value ΔH3 related to the mixing of a solution with a suspension of MPs (Equation (5)) is small, it falls within the limits of the experimental error for the other values in the cycle, and it is usually neglected.
Composites based on Na-alginate polymer with embedded MPs in different content were prepared separately as described above in the Materials and Methods part.
Figure 7 presents the concentration dependence of the enthalpy of dissolution (Δ
H4) of these composites in water at 25 °C measured in a direct experiment in calorimeter. Ni and Ni@C MPs were taken as examples.
Figure 7 contains all the necessary information to calculate the enthalpy of adhesion using the thermodynamic cycle in Equation (7). The value at the left-side ordinate axis corresponds to the composite with 0% content of MPs, which is pure Na-alginate. It is the Δ
H1 value in Equation (3). The value at the right-side axis is Δ
H2 in Equation (4) as it corresponds to the composite with 100% of MPs. Other points in the plot correspond to the values of the enthalpy of dissolution of composites with progressively increasing MPs content (Δ
H4 in Equation (6)).
In fact, any composite is a mixture of components and even if there is no interaction between components, the enthalpy of dissolution would depend on the proportion of the components in a mixture as the values of Δ
H1 and Δ
H2 are different. Theoretically, the enthalpy of dissolution of an additive mixture of components would be a linear function between the values for the constituents. It is noticeable in
Figure 7 that it was the case of Na-alginate composites with embedded Ni MPs. The experimental points for the enthalpy of dissolution in this system are close to a straight line between Δ
H1 and Δ
H2. Meanwhile, the values of the enthalpy of dissolution for Na-alginate composites with Ni@C deviated upward from the straight line of an additive mixture. It meant that there was the interaction between components in composites with Ni@C particles, unlike the system with Ni MPs. Most likely it was the result of the surface modification of Ni with a carbon layer.
Experimental data presented in
Figure 7 and similar data for Na-alginate composites with Fe
3O
4 and Fe MPs were used for the calculation of the enthalpy of adhesion of Na-alginate to the surface of MPs using Equation (7). The results are plotted in
Figure 8, which gives concentration dependence on the enthalpy of adhesion. The values of the enthalpy of adhesion are related to the unit area (m
2) of the surface of the embedded MPs.
The dependences in
Figure 8 are cupola curves, which pass through zero at both ends of the abscissa axis, that is at 0% and 100% of MPs content. The values of Δ
Hadh were negative in the case of Na-alginate composites with Fe
3O
4 and Ni@C composites, they were very close to zero for the composites with Ni and were slightly positive for the composites with Fe. It meant that the interaction of Na-alginate with Fe
3O
4 and Ni@C MPs is energetically favorable and there is a gain in enthalpy in this process. There is almost no interaction of Na-alginate with Ni MPs, and interaction with Fe MPs is slightly unfavorable.
It is noticeable that the values of the enthalpy of adhesion for Na-alginate increase in a sequence Fe
3O
4–Ni@C–Ni–Fe, which is the same sequence as that pointed out above concerning
G′,
G″, and
δ for magneto-rheological suspensions. Combining these two results we may conclude that good interaction of Na-alginate with Fe
3O
4 MPs resulted in fluid-type behavior of suspensions (see
Figure 6). On the contrary, Fe MPs poorly interacted with Na-alginate, the enthalpy of adhesion was positive, and the suspensions with Fe MPs exhibited dominance of elasticity at their volume fraction above 3% (vol.).
The pair Ni–Ni@C showed the same feature. Ni MPs were almost non-interacting with Na-alginate and it resulted in the distinct elasticity of suspensions with a volume fraction of Ni above 3% (see
Figure 6). Deposition of the carbon layer had improved the adhesion of Na-alginate to the surface of MPs but it resulted in a drawback in the elasticity of the suspensions. As one may see in
Figure 6, the suspensions with 3.5% (vol.) of Ni@C MPs remained fluids, while the suspensions with Ni MPs at the same volume fraction were dominantly elastic.
4. Discussion
From the basic physical point of view, any measurable property of the physicochemical system takes its origin in the balance of all molecular interactions of the interior. If there are no chemical reactions in the system as in the case of alginate-based MRSs then all the features of the macroscopic properties are governed by the following types of interaction: hydration of alginate, hydration of solid magnetic particle, interfacial adhesion of alginate to the surface of particles, and magnetic forces among particles. The study was focused on the influence of interfacial adhesion on the rheological properties of alginate MRSs. Presented experimental results revealed the correlation between the viscoelastic properties of alginate-based suspensions and the enthalpy of adhesion of alginate polymer to the surface of magnetic particles in composites. If the values of the enthalpy of adhesion of Na-alginate to the solid surface of MPs were negative, which meant the energetically favorable interaction, then it led to a fluidic type of deformation of magneto-rheological suspensions at the same time. In other words, good interfacial adhesion at the interface depressed the elevation of the dynamic moduli and the descent of the shift angle with the increase in MPs content. Thus, the enhancement of polymer-solid interaction might be favorable to preserve the fluidity of the suspension. Oppositely, if the interaction of the polymer and the solid at the interface was poor, it provided the boost of the dynamical moduli and the dropdown of the shift angle at a relatively low volume fraction of MPs. Thus, based on the obtained results we may conclude that to shift the properties of MRS to elasticity rather than fluidity and to transform the fluid to a gel phase one needs to block polymer-solid interaction somehow.
At the molecular level, we may suppose the following explanation of the observed dependence of the DMA behavior of MRS on the interaction at the interfaces. If the interaction of Na-alginate with the surface of MPs is energetically favorable, then polymer chains adsorb at the surface and form a steric protective layer which prevents a close contact between magnetic particles. Based on the extended Derjaguin-Landau-Verwey-Overbeek theory it was shown earlier [
26], that magnetic interactions between Fe MPs in colloids dominated over van der Waals and electrostatic forces. Therefore, thick steric protective polymeric layers on the surface of MPs are needed to prevent aggregation of MPs due to the magnetic forces and to provide colloidal stability of the suspension.
In the present study, no special alignment of particles in the preparation of MRSs was provided. Therefore, the distribution of magnetic particles in MRS was macroscopically random. Meanwhile, the possibility of aggregation at the microscopic scale could not be completely excluded. The analysis of particle aggregation in condensed phases is very complicated and to a large extent an unresolved problem. In the present study, we could not evaluate the aggregation as the content of magnetic particles in suspensions was too large for optical microscopy. In this respect, the following model consideration might refer not only to individual MPs but too small aggregates of MPs as well.
The caliper diameter of MPs in the present study (see
Table 1) was above the threshold to a single domain state which is about 128 nm for Fe
3O
4, 14–30 nm for Fe, and 55 nm for Ni [
29]. Thus, the MPs under study consisted of several domains each. In this situation magnetic interaction between the nearest domains in the closely situated particle can be not fully compensated by the interaction between the other domains in these particles and the total magnetic interaction of these “bare” particles can be quite significant, provoking their agglomeration and gelation of the suspension. The ensemble of individual MPs covered with the polymeric protective layer and with partially compensated domain magnetic moments would likely favor the independent mobility of particles under the applied mechanical stress and consequently the fluid-type deformation. Schematically this scenario is given in
Figure 9a. Apparently, it occurred in the case of Fe
3O
4 and Ni@C MPs.
Figure 9b schematically illustrates the other situation if there is no interfacial adhesion of Na-alginate at the surface of MPs. In this case, polymer does not adsorb at the surface and there is no polymeric layer at the surface which prevents close contact between particles. At low MPs concentration if the probability of inter-particle contacts is small this case does not differ from the previous one. The situation changes if the content of particles increases, and the probability of their contact becomes substantial. As MPs are not covered with steric polymeric layers, they may establish close contact and aggregate with each other. In this aggregation, magnetic forces may appear between domains that are located at the neighboring particles. Although the cumulative magnetic moment of the multi-domain particle is close to zero, domains at the surface have non-zero moments. Then, the magnetic interaction of local magnetic moments at the surfaces of neighboring particles may take place at their close contact. As a result of this interaction, a network structuring of magnetic particles may appear, which causes the gelation of the suspension. The evidence in favor of this scenario was reported in reference [
30] for ferrocolloids with Co nanoparticles, which formed chain-like structures, and in reference [
31] concerning ferrogels with embedded strontium hexaferrite magnetic particles which showed network-type aggregation of multi-domain magnetic particles. The networking of magnetic particles due to their orientation by attractive magnetic interaction apparently would lead to the prevalence of the elastic deformation of MRS.
It is worthwhile noting that the presented specific behavior of MRS suspensions was observed at a low frequency (1 Hz) of deformation. Low-frequency mode is modeling slow motions in biomimetic materials, for instance, in their application as scaffolds for tissue regeneration or as working bodies in actuators. Meanwhile, there can also be a high-frequency mode for such materials, for example in drug delivery, if the gel drug carrier is activated by ultrasound. In this respect, it might be interesting to examine the mechanical moduli of the magneto-rheological suspensions and the influence of the interfacial adhesion using ultrasound. An example of such an approach to the properties of gel materials with volume phase transition was demonstrated in ref. [
32].