In the last decades, the synthesis of iron oxide magnetic nanoparticles (MNPs) has gained much attention,[
Chinese Physics B, Volume. 29, Issue 10, (2020)
Evaluating physical changes of iron oxide nanoparticles due to surface modification with oleic acid
The physical characterization of a colloidal system of superficially modified magnetic nanoparticles (MNPs) is presented. The system consists of oleic acid-coated iron oxide nanoparticles (OAMNP) suspended in water. A structural analysis is carried out by using standard physical techniques to determine the diameter and shape of the MNPs and also the width of the coating shell. The colloidal stability and the polydispersity index of this ferrofluid are determined by using Zeta potential measurements. Additionally, the magnetic characterization is conducted by obtaining the DC magnetization loops, and the blocking temperatures are determined according to the ZFC–FC protocol. Finally, the values of power absorption density P of the ferrofluid are estimated by using a magneto-calorimetric procedure in a wide range of magnetic field amplitude H and frequency f. The experimental results exhibit spherical-like shape of OAMNP with (20 ± 4) nm in diameter. Due to the use of coating process, the parameters of the magnetization loops and the blocking temperatures are significantly modified. Hence, while the uncoated MNPs show a blocking state of the magnetization, the OAMNP are superparamagnetic above room temperature (300 K). Furthermore, the reached dependence P versus f and P versus H of the ferrofluid with coated MNPs are clearly fitted to linear and quadratic correlations, respectively, showing their accordance with the linear response theory.
1. Introduction
In the last decades, the synthesis of iron oxide magnetic nanoparticles (MNPs) has gained much attention,[
Ferrofluids of MNPs with controlled diameters and shapes, are typically produced by using the well-known chemical co-precipitation procedure, where a base solution is added to liquid saline solutions of Fe2+/Fe3+, under an inert environment with adjustable temperature and pH.[
Superficial modification of the MNPs produces many physical changes, which are exploited to reach a good colloidal stability and diminish the internal dipolar interactions, thus avoiding forming magnetic clusters. Nevertheless, the magnetic properties and monodispersity of the MNPs should not be dramatically modified after this functionalization process. Thus, the main aim of this work is to evaluate the induced changes in the structure and magnetic properties of Fe3O4 MNPs which are coated with oleic acid. Later, the power absorption density of the ferrofluid with OAMNP is analyzed to explore its expected performance in magnetic hyperthermia experiment.
2. Theoretical aspect
The interaction between alternant magnetic field and ferrofluid system can be modeled by using the linear response theory. In this approximation, a linear relationship between the induced magnetization M and H is assumed. Then, the drained power density P on a volume V of ferrofluid is given by Eq. (1).[
Regarding the effective relaxation time τ, this parameter takes into account both the Néel relaxation time (τN) and the Brownian relaxation times (τB) as expressed in Eq. (3). Thus, τN includes the effective magnetic anisotropy constant κe of the suspended MNPs through the Arrhenius law (see Eq. (4)), where τ0 is the characteristic length of time (typically 1 × 10−9 s). In addition, τB includes the viscosity η of the liquid medium and the hydrodynamic diameter DH of the MNPs, following the Stokes–Einstein law (Eq. (5))
3. Materials and methods
3.1. Synthesis and preparation of ferrofluid
The ferrofluid of uncoated MNPs is prepared following a co-precipitation methodology. First, 30 ml of deionized water (previously deoxygenated with N2) is deposited in a glass test tube and 15-ml deionized water in another one, both are isolated from the atmosphere to avoid absorbing oxygen. A Schlenk flask is depressurized to 500 mm of Hg within 8 min and then, it is filled with N2. Afterwards, 0.4867 g of FeCl3 (Sigma Aldrich) is dissolved in the first test tube and 0.417 g of FeSO4 ⋅ 7H2O (Fermont) in the second tube. Subsequently, both solutions are deposited in the flask. The blend is heated at 30 °C under continuous stirring at 400 rpm. At this point a chronometer is started in order to control the addition sequence of the reactants and heating periods. Right after 5 min, 3 ml of 5-M NH4OH is added by using a peristaltic pump at a volumetric rate of 1.5 ml/min. After 15 min, a syringe needle is inserted into the septum of the flask. 30 min later, the synthesis is finished. A blackish suspension is obtained, indicating the formation of MNPs. The suspension of MNPs is washed four times with deionized water to eliminate the remaining by-products of the synthesis. Finally, a ferrofluid of MNP dispersed in water is obtained.
In the same sense, the procedure to obtain the ferrofluid of OAMNPs includes the aforementioned steps, followed by the addition of oleic acid. After 10 min, the dispersion is maintained at 80 °C for 10 min, and then, at 90 °C for 10 min. Then, 300 μl of oleic acid (Fluka Analytics, reagent grade) is added into the dispersion and kept stirring for 20 min. Afterwards, the stirring is stopped, and a colloidal dispersion of OAMNPs is obtained. The excess of oleic acid is eliminated by centrifugation at 1.2 × 104 rpm for 1 h and 10 °C. The supernatant is removed and the OAMNPs are resuspended in deionized water. The centrifugation and resuspension in deionized water steps are repeated four times, and finally, the obtained ferrofluid of OAMNP is dialyzed by using a cellulose membrane (Sigma-Aldrich, MWCO 12 kDa, 1 Da = 1.66054 × 10−27 kg).
3.2. Structural and magnetic characterization of ferrofluid
The core diameter of the MNPs is estimated through x-ray diffraction (XRD) measurements by using a Panalytical Empyrean system and its goniometer covers the angular displacement 5° ≤ 2θ ≤ 80°. Then, Δθ = 0.02° is the programmed angular increment and t = 30 s is the sampling time for data acquisition. For this purpose, dried samples with 50 mg of coated and uncoated MNPs are placed on specials glass sample holders for XRD. Additionally, dried samples are analyzed through Fourier transform infrared (FTIR) spectroscopy measurements by using a Thermo Scientific Nicolet iS5 spectrometer in an ATR mode.
The scanning electron microscopy (SEM) measurements are obtained by using an MIRA3 (LMU) device of TESCAN company, with 1 nm of resolution capability. Samples are prepared by diluting 10 μl of ferrofluid with coated and uncoated MNPs in 1 ml of deionized water; subsequently, 10-μl aliquot is placed on a carbon conductive tape (5 mm × 20 mm), which is adhered to the special sample holder for SEM. During the next 35 min, this sample is dried at 45 °C inside an oven and then they are covered via gold sputtering.
In order to analyze the colloidal stability of water suspended coated and uncoated MNPs, Zeta potentials are determined at different pH values by means of a Zetasizer SZ90 (Malvern) through using polystyrene folded capillary disposable cells. First, a 0.01-M NaCl solution is prepared with deionized water, then, 100 μl of MNP at 15 mg/ml is added into 10 ml of NaCl solution. Afterwards, the pH of the dilutions is adjusted to a desired value by adding NaOH or HCl. Dispersions with pH values ranging from 2 to 11 are measured three times to determine the mean values and standard deviation. Additionally, the corresponding Z-average hydrodynamic diameter (DH) and polydispersity index (PDI) are measured by using polystyrene disposable cuvettes (10-mm path length, four clear sides).
To determine the fraction of oleic acid covering the surface of the MNPs, thermo-gravimetric (TGA) measurements are realized with the help of a high-resolution Cahn Versa Therm Analyzer and two sample containers with 1.5-g capacity. A container is filled with 60 mg of dried coated MNPs while the second container is filled with the same mass of uncoated MNPs. Then, both samples are separately heated by sweeping the temperature interval 40 °C < T < 850 °C, at 10 °C/min of heating rate. During the experiments, the uncoated MNPs are always heated under an inert N2-atmosphere (at 25 ml/min of flux) and the coated samples are heated up to 750 °C under N2-atmosphere, but in the last interval (750 °C < T < 850 °C) the N2 is replaced by oxygen.
The magnetization loops (at room temperature 300 K) of the MNPs are obtained using a VersaLab vibrating sample magnetometer (VSM) of quantum design. Hence, 5 mg of dried coated and uncoated MNPs are placed into diamagnetic containers and uniform magnetic fields from –30 kOe–30 kOe (1 Oe = 79.5775 A/m) are applied. Magnetization traces following the ZFC–FC protocol are registered establishing H = 100 Oe, while the temperature covers the interval 50 K < T < 400 K.
To determine the power absorption density of the OAMNP ferrofluid, a previously reported induction heater is used.[
4. Results
The typical XRD spectra obtained from uncoated (dark line) and OAMNP (gray line) samples are displayed in Fig. 1(a), and both measurements show the characteristic peaks of the magnetite structure. Thus, the Miller indices (2 2 0), (3 1 1), (4 0 0), (4 2 2), (3 3 3), and (4 4 0) are clearly distinguished, which correspond to the respective phases 2θ = 30.16°, 35.52°, 43.17°, 53.56°, 57.10°, and 62.70° (JCPDS file, No. 19-0629). Using the Scherrer formula (in radian units) σ = kλ(β(2θ)cos(θ))−1, the estimated average crystallite sizes of uncoated σuc and coated σOAMNP samples are σuc=(13 ± 1) nm and σOAMNP = (15 ± 1) nm respectively. For this estimation a shape factor k = 0.94 is considered, the x-ray wavelength is λ = 1.54 × 10−10 m, and the determined parameters β(2θ) and θ are obtained from Fig. 1(b), which represent the close-up of the main peaks displayed in Fig. 1(a). In order to determine the standard deviation of σuc and σOAMNP, an experiment is performed three times, hence an approximated increase of 13% is observed in the magnetic core of OAMNP, these discrepancies can be attributed to random effects of the preparation processes.
Figure 1.(a) Typical XRD spectra of uncoated (dark lines) and OAMNP (gray lines) samples, (b) close-up of the main peaks, (c) SEM micrograph of the OAMNP, (d) size distribution plot fitted to a normal data regression, (e) SEM of uncoated MNPs, and (f) corresponding plot including the data regression.
Other spurious uncertainties of the average σ can be intrinsically involved in the Scherrer equation when the sample has a polydispersed diameter below 10 nm and a crystal structure with superficial strain.[
In relation to the superficial modification of the MNPs, the FTIR spectra of pure oleic acid, uncoated and coated MNP are displayed in Fig. 2. In the dark waveform, the presence of Fe3O4 is corroborated due to the vibrations observed at 577 cm−1 which corresponds to the Fe–O bond, as well as those at 1622 cm−1 and 3209 cm−1 of the O–H groups.[
Figure 2.FTIR spectra of uncoated Fe3O4 nanoparticles (dark), pure oleic acid (blue), and OAMNP (gray).
On the other hand, the relative mass dependent temperatures of uncoated (dark line) and OAMNP (gray line) are shown in Fig. 3(a). These measurements are obtained via TGA analysis and two clear different behaviors are observed. After a careful observation of the mass registers, the uncoated MNP sample has 2.87 % of mass loss at 39 °C < T < 296 °C and also has 0.45% of mass loss at temperatures above of 296 °C. This decrease is caused by evaporation of water adsorbed on the particle’s surfaces. Nevertheless, the sample of OAMNP decreases by 4.4% at 30 °C < T < 165 °C, followed by a 54.4-% mass loss at 165 °C < T < 337 °C. Also, 6.5 % of mass loss at 338 °C < T < 576 °C is observed, and finally, a 3-% decrease above 577 °C. At the end, 31.7% of iron oxide is reached above 750 °C. In the same sense, the derivative of the relative mass loss is shown in Fig. 3(b), where the characteristic peaks of the mass loss rate are highlighted. Hence, the peak at 77 °C of the OAMNP is associated with the water evaporation (4.4% of mass loss), the peaks at 340 °C and 541 °C are related to the boiling point of the oleic acid (60.9% of mass loss) and the last peak at 758 °C to the oxidation of iron oxide (34.7% of relative mass plus resultant material). With the change of atmosphere at T > 750 °C, clearly can be corroborated the expected phase transition from magnetite to hematite.[
Figure 3.(a) Relative mass
The measurements of Zeta potential over the pH interval 2.2 < pH < 11 are displayed in Fig. 4(a). For OAMNP, an initial negative decrease of the Zeta potential is observed to reach a minimum value (–80 mV) at pH = 8; this value increases with the systematic increase of pH and it is maintained below –45 mV for 9 < pH < 11. Thus, for pH > 4 a good colloidal stability of the OAMNP ferrofluid is achieved. In contrast, the ferrofluid with uncoated MNPs shows a low stability behavior for pH > 4, an incipient stability at 6 < pH < 7 and a good stability only for pH > 9. Additionally, the dependence of hydrodynamic diameters DH on pH is shown in Fig. 4(b). In general, the higher DH is observed in uncoated MNPs. Indeed, the largest diameter DHmax = 5720 nm is registered at pH = 8, but the Zetasizer cannot determine the corresponding DH for pH > 8, even on its higher resolution scale. When the ferrofluid of OAMNP is analyzed the behavior notably changes, because this parameter systematically decreases with the increase of pH, reaching up to 140 nm at pH = 11. The remarkable differences between the σOAMNP via SEM and the σOAMNP via DH are explained by the formation of islands of OAMNP in the ferrofluid, which are probably due to electrostatic and Van Der Waals interactions. A similar trend is observed on the PDI, where the higher values are observed in uncoated MNPs (see Fig. 4(c)). Indeed, for uncoated MNPs, the lowest value PDImin = 0.6 is reached at pH = 2, which indicates a completely polydisperse behavior over the studied pH range. For OAMNP, the measured PDI indicates some polydispersity for pH > 8 and pH = 6. Nevertheless, an acceptable monodispersity is observed at pH < 6 and pH = 7.4; indeed, the best stability values are reached for pH = 4 (PDI = 0.13 7#x00B1; 0.2) and pH = 7.4 (PDI = 0.18 ± 0.2).
Figure 4.The pH-dependent (a) Zeta potential, (b) their corresponding hydrodynamic diameters, and (c) polydispersity index of water-suspended coated (open circles) and uncoated (black circles) MNPs.
Magnetic properties of uncoated and coated MNPs under the action of static magnetic field are discussed now. In Fig. 5(a), the curves of magnetization M versus magnetic field of the uncoated MNP (dark line) and coated MNP (gray line) are plotted, which are obtained at room temperature. The magnetic saturation of the OAMNP is clearly diminished, reaching Ms = 24 emu/g, which is 2.92 times lower than the same parameter obtained for the uncoated sample. Regarding to the magnetic cohercivity Hc, a similar result is observed, reaching 11 Oe and 25 Oe for coated and uncoated samples respectively. Also, the low magnetic remanence Mr = 0.25 emu/g of OAMNP reaches up to 6.8 times less than the measured in uncoated samples. Hence, those significant differences in the magnetization loop are caused by the lower concentration of magnetite in the OAMNP sample (only 31.7%), in addition, the coating shell helps to diminish the internal dipolar interactions between neighboring particles. Likewise, the dependence of M on temperature following the ZFC–FC protocol (at H = 100 Oe) is shown in Fig. 5(b). Notable differences between both curves can be clearly seen, but the ZFC trace of uncoated samples does not exhibit any maximum value with null slope, indicating the magnetic blocking state of M. Whereas the magnetization of OAMNP exhibits a blocking temperature TB = 303 K, evidencing a superparamagnetic ordering. Despite a low Hc observed on coated samples, the inset of Fig. 5(a) shows some of ferrimagnetic ordering, but the obtained TB indicates a predominate superparamagnetic behavior above room temperature. These mixed magnetic properties are the consequence of the polydispersity of the samples, which coincides with the measured PDI of Fig. 4(c). Several studies have reported the modulation of TB as a consequence of a decrease in the dipolar interactions,[
Figure 5.(a) Magnetization loops at room temperature of uncoated and coated MNPs, and (b) their corresponding ZFC–FC graphs using
Using the core size σOAMNP obtained via XRD and the observed blocking temperature TB, we can obtain κe ≈ 9.1 × 10−4 J/cm3 from Eq. (6), which is the estimated effective anisotropy constant of OAMNP. In consequence, τN ≈ 88 s is the effective Néel relaxation time, which is computed by substituting κe into Eq. (4). Therefore, the reached Brownian relaxation time τB ≈ 2.7 ms is computed from Eq. (5) under the assumption of a water viscosity of η = 8.9 ° 104 N ⋅ s/m2 and DH ≈ 200 nm. This value is four orders of magnitude less than τN, and then τ ≈ τB. Thus, the effective relaxation time τ described by Eq. (3) is dominated by the rotation of the OAMNP.
The induced heating on the OAMNP ferrofluid is shown in Fig. 6(a). In this experiment, H = 25 mT is the constant amplitude applied during 2 min, while f has been increased in five steps from 182 kHz to 530 kHz. The five traces of the temperature increment over time exhibit an almost linear behavior, where the slope d T / d t systematically grows at each frequency. Indeed, those slopes are computed using data linear regressions in the time interval 20 s < t < 40 s, and also they are used to determine the dependence Pversus time at each value of f as displayed in Fig. 6(b). Moreover, the parameter P is determined first by computing the corresponding SAR values from Eq. (8), which later is introduced into Eq. (7). For these purposes, are used are the specific heat Cv = 4.186 kJ/(K ⋅ kg), density ρ = 1000 kg/cm3, and the suspended fraction of OAMNP mnp = 0.1. The delineated solid line on the plot of Fig. 6(b) is the best linear data fit, which reaches the slope mf = 0.047 ± 0.004 W/cm3 per kHz and it is characterized by the quality factor χ2 = 0.94. This statistical parameter evidences good accordance with P corresponding to the linear response model Eq. (1). Then, P ≈ 22 W is the maximum power dissipated by OAMNP to 1 cm3 of water when H = 25 mT and f = 530 kHz. This magnitude is a very important parameter for future applications in the induction heating of biologic agents.
Figure 6.Time dependent (a) temperature increment during 2 min applying
In other assays, the OAMNP ferrofluid is heated using the constant frequency f = 330 kHz, whereas the amplitude H is increased from 10 mT to 30 mT in steps of 5 mT. In Fig. 7(a), the five traces of the temperature increment during 2 min are shown and the gradual increase of d T / d t with H is evidenced. The same procedure described above (Fig. 6(a)) is used to determine the dependence of P on H and it is displayed in Fig. 7(b). There, the solid line is the best quadratic data fit estimated using the formula P = mh ⋅ H2, where
Figure 7.(a) Temperature increment during 2 min applying
The linearity condition for P is given by the expression μ0HVpMs < kBT, and this inequality can be computed using the employed parameters
5. Conclusions
In this work, the physical changes of a ferro-colloid system of superficially modified iron oxide MNPs are analyzed. Three principal changes must be highlighted as a consequence of the oleic acid coating, i.e. the high colloidal stability related to a low polydispersity index, the modulation of the magnetic properties which leads the magnetization to go up to a superparamagnetic ordering, and the good accordance of its power absorption density with the linear response theory.
[1] A K Gupta, M Gupta. Biomaterials, 26, 3995(2005).
[2] W Wu, Z Wu, T Yu, C Jiang, W S Kim. Science and Technology of Advanced Materials, 16(2015).
[3] L Babes, B Denizot, G Tanguy, J J Le Jeune, P Jallet. Journal of Colloid and Interface Science, 212, 482(1999).
[4] A Moore, E Marecos, A Bogdanov, R Weissleder. Radiology, 214, 568(2000).
[5] A Jordan, R Scholz, P Wust, H Fähling, R Felix. J. Magn. Magn. Mater., 201, 413(1999).
[6] A Jordan, R Scholz, P Wust, H Fähling, J Krause, W Wlodarczyk, B Sander, T Vogl, R Felix. International Journal of Hyperthermia, 13, 587(1997).
[7] R A Weissleder, D D Stark, B L Engelstad, B R Bacon, C C Compton, D L White, P Jacobs, J Lewis. American Journal of Roentgenology, 152, 167(1989).
[8] A P Philipse, M P Van Bruggen, C Pathmamanoharan. Langmuir, 10, 92(1994).
[9] L Shen, P E Laibinis, T A Hatton. Langmuir, 15, 447(1999).
[10] J De Vicente, A V Delgado, R C Plaza, J D Durán, F González-Caballero. Langmuir, 16, 7954(2000).
[11] P A Dresco, V S Zaitsev, R J Gambino, B Chu. Langmuir, 15, 1945(1999).
[12] L Shen, Y Qiao, Y Guo, S Meng, G Yang, M Wu, J Zhao. Ceram. Int., 40, 1519(2014).
[13] M Filippousi, M Angelakeris, M Katsikini, E Paloura, I Efthimiopoulos, Y Wang, D Zamboulis, G Van Tendeloo. J. Phys. Chem. C, 118(2014).
[14] K Petcharoen, A Sirivat. Mater. Sci. Eng. B, 177, 421(2012).
[15] T Ahn, J H Kim, H M Yang, J W Lee, J D Kim. J. Phys. Chem. C, 116, 6069(2012).
[16] F M Kievit, Z R Stephen, O Veiseh, H Arami, T Wang, V P Lai, J O Park, R G Ellenbogen, M L Disis, M Zhang. ACS Nano, 6, 2591(2012).
[17] A Shkilnyy, E Munnier, K Hervé, M Soucé, R Benoit, S Cohen-Jonathan, P Limelette, M L Saboungi, P Dubois, I Chourpa. J. Phys. Chem. C, 114, 5850(2010).
[18] Y V Kolen’ko, M Bañbre-Loóez, C Rodríguez-Abreu, E Carbó-Argibay, A Sailsman, Y Piñiro-Redondo, M F Cerqueira, D Y Petrovykh, K Kovnir, O I Lebedev, J Rivas. J. Phys. Chem. C, 118, 8691(2014).
[19] P I P Soares, C A T Laia, A Carvalho, L C J Pereira, J T Coutinho, I M M Ferreira, C M M Novo, J P Borges. Appl. Surf. Sci., 383, 240(2016).
[20] K Yang, H Peng, Y Wen, N Li. Appl. Surf. Sci., 256, 3093(2010).
[21] D Maity, D C Agrawal. J. Magn. Magn. Mater, 308, 46(2007).
[22] M Mahdavi, M B Ahmad, M J Haron, F Namvar, B Nadi, M Z A Rahman, J Amin. Molecules, 18, 7533(2013).
[23] Y Li, F Ma, X Su, L Shi, B Pan, Z Sun, Y Hou. Industrial & Engineering Chemistry Research, 53, 6718(2014).
[24] R E Rosensweig. J. Magn. Magn. Mater., 252, 370(2002).
[25] J Carrey, B Mehdaoui, M Respaud. J. Appl. Phys., 109(2011).
[26] J A Dearing, P M Bird, R J Dann, S F Benjamin. Geophys. J. Int., 13, 727(1997).
[27] E E Mazon, E Villa-Martínez, A Hernández-Sámano, T Córdova-Fraga, J J Ibarra-Sánchez, H A Calleja, J A Leyva Cruz, A Barrera, J C Estrada, J A Paz, L H Quintero. Rev. Sci. Instrum., 88(2017).
[28] E E Mazon, A H Sámano, H Calleja, L H Quintero, J A Paz, M E Cano. Measurement Science and Technology, 28(2017).
[29] D W Armitage, H H Le Veen, R Pethig. Phys. Med. Biol., 28, 31(1983).
[30] B Ingham. Crystallography Reviews, 21, 229(2015).
[31] G A Dorofeev, A N Streletskii, I V Povstugar, A V Protasov, E P Elsukov. Colloid Journal, 74, 675(2012).
[32] L Y Zhang, Y H Dou, L Zhang, H C Gu. Chin. Phys. Lett., 24, 483(2007).
[33] H Zheng, Y Yang, F S Wen, H B Yi, D Zhou, F S Li. Chin. Phys. Lett., 26(2009).
[34] Z L Wang, H Ma, F Wang, M Li, L G Zhang, X H Xu. Chin. Phys. Lett., 33(2016).
[35] N Wu, L Fu, M Su, M Aslam, K C Wong, V P Dravid. Nano Lett., 4, 383(2004).
[36] W F Nor, S K Soh, A A Azmi, M S Yusof, M Shamsuddin. Malaysian Journal of Analytical Sciences, 2, 768(2017).
[37] M El-Hilo, R W Chantrell, Y K O’Grady. J. Appl. Phys., 5114(1998).
[38] J L Dormann, L Bessais, D Fiorani. J. Phys. C: Solid State Phys., 21, 2015(1988).
[39] M Knobel, L M Socolovsky, J M Vargas. Rev. Mex. Fís., 50, 8(2004).
[40] W C Nunes, F Cebollada, M Knobel, D Zanchet. J. Appl. Phys., 99(2006).
Get Citation
Copy Citation Text
S Rosales, N Casillas, A Topete, O Cervantes, G González, J A Paz, M E Cano. Evaluating physical changes of iron oxide nanoparticles due to surface modification with oleic acid[J]. Chinese Physics B, 2020, 29(10):
Received: Feb. 1, 2020
Accepted: --
Published Online: Apr. 21, 2021
The Author Email: M E Cano (meduardo2001@hotmail.com)