Next Article in Journal
Experimental Investigation of the Stress–Strain Behavior and Strength Characterization of Rubberized Reinforced Concrete
Next Article in Special Issue
Titanium(IV) Oxo-Complex with Acetylsalicylic Acid Ligand and Its Polymer Composites: Synthesis, Structure, Spectroscopic Characterization, and Photocatalytic Activity
Previous Article in Journal
Research Progress on Controlled Low-Strength Materials: Metallurgical Waste Slag as Cementitious Materials
Previous Article in Special Issue
Electrochemical Mechanism of Recovery of Nickel Metal from Waste Lithium Ion Batteries by Molten Salt Electrolysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Magnetite Impregnated Lignocellulosic Biomass for Zn(II) Removal

by
Christopher Asimbaya
1,
Nelly Maria Rosas-Laverde
1,
Salome Galeas
1,
Alexis Debut
2,
Victor H. Guerrero
1,* and
Alina Pruna
3,4,*
1
Department of Materials, Escuela Politécnica Nacional, Quito 170524, Ecuador
2
Centro de Nanociencia y Nanotecnología, Universidad de las Fuerzas Armadas (ESPE), Sangolqui 171103, Ecuador
3
Institute of Materials Technology, Universitat Politècnica de València, 46022 Valencia, Spain
4
Center for Surface Science and Nanotechnology, University Politehnica of Bucharest, 313 Splaiul Independentei, 060042 Bucharest, Romania
*
Authors to whom correspondence should be addressed.
Materials 2022, 15(3), 728; https://doi.org/10.3390/ma15030728
Submission received: 18 October 2021 / Revised: 24 November 2021 / Accepted: 26 November 2021 / Published: 19 January 2022

Abstract

:
Magnetic composites obtained by impregnation of lignocellulosic biomass with magnetite nanoparticles were used for zinc(II) removal from aqueous synthetic solutions. Laurel, canelo and eucalyptus sawdust, with a particle size between 74 and 150 µm were used as support. Structural and morphological examinations of the composites confirmed the presence of magnetite nanoparticles in the lignocellulosic support. Transmission Electron Microscopy showed nanoparticles with diameters of about 20 nm. The maximum removal efficiencies for 7 g L−1 of modified adsorbent were increased to 98.9, 98.8 and 97.6% for laurel, canelo and eucalyptus magnetic composites, respectively, in comparison to 60.9, 46.0 and 33.3%, for corresponding unmodified adsorbents. Adsorption data was analyzed using pseudo-first, pseudo-second order and intra-particle diffusion kinetic models and various isotherm models. The results determined that Freundlich isotherm fits the Zn ions adsorption on magnetite modified adsorbents while chemisorption and boundary diffusion were dominating the process.

1. Introduction

The increasing industrial activity during the last decades has resulted in a worldwide concern problem: heavy metals contamination of liquid effluents. These metals have become a serious problem because their ions are not biodegradable and accumulate in animal tissues, damaging human health and water quality [1,2,3]. One of the most common heavy metals in industrial discharges is zinc, and the main sources of pollution are the automotive, electroplating, iron and steel industries [2,4,5,6].
Nowadays, several techniques are used to remove heavy metals from water, such as: sedimentation, filtration, electrocoagulation, ion exchange and reverse osmosis, however, disadvantages of these methods are low efficiency, high energy requirements, incompatibility with pollutants and toxic residues generation, that require subsequent treatments [6,7]. For these reasons, alternative treatments have been developed and optimized in the few last decades. One of the most efficient, low cost, that is widely used is adsorption [8,9].
The most common adsorbent used for heavy metal removal is activated carbon, however it is expensive. Recent studies focus on alternative materials: inexpensive, abundant in nature, with minimal processing requirement and in some cases by-products of other industries [9]. These materials include waste wood, mainly sawdust, which due to its porosity has high free surface that results in removal of heavy metals through the formation of bonds between the metal ions and the oxygen of ionizable functional groups [10].
Therefore, nanotechnology has allowed the development of more efficient adsorbent materials for environmental remediation. For example, iron nano-oxides (magnetite, hematite and maghemite nanoparticles) have a magnetic core that improve the adsorption process. In addition, they allow an easy recovery, handling and less contamination [11]. This, combined with the adsorptive capacity of the sawdust allows the removal of metal ions from liquid effluents to be more effective and the recovery of the residue to be easier and less expensive [12].
Numerous studies performed during the last decades have reported lignocellulosic materials impregnated with magnetic materials that have been used successfully to remove metals from contaminated water. Moafi et al. [13] used a wood sawdust/hematite (α-Fe2O3) nanocomposite to adsorb arsenic (III) from water with a removal capacity of 58.8 mg g−1. Panneerselvam et al. [5] obtained a nickel removal of 38.3 mg g−1 using tea residues impregnated with magnetic nanoparticles. Additionally, Gupta and Nayak [14] used orange peel residues impregnated with magnetite nanoparticles to treat cadmium-contaminated effluents and reported an adsorption capacity of 76.9 mg g−1.
Wood residues are widely available around the globe, are cheap and have a particularly favorable composition, structural and morphological characteristics for adsorption. For instance, Ecuador has become one of the most biodiverse countries around the world due to the variety of climates which resulted in high availability of plant species. According to National Financial Corporation [15] statistics for 2016, forest area of Ecuador was around 54,000 ha, and the annual production of wood was about 421,000 metric tons. Additionally, wood industry is well known for generating big amounts of waste, up to 70% of the wood raw materials can become chippings, sawdust or off-cuts [16,17]. Therefore, wood residues can be advantageously used to remove a wide variety of contaminants from aqueous solutions [18]. Although the removal of heavy metals using magnetic composite bioadsorbents has been investigated by several previous authors, the use of supporting materials derived from species such as Ecuadorian laurel (Cordia alliodora), canelo (Ocotea javitensis) and eucalyptus (Eucalyptus glodulus) has not been sufficiently analyzed. In several studies, the improvement of the adsorptive performance obtained by using particulated materials to impart magnetic properties to the adsorbents is not well defined and understood. Moreover, in some cases the study of the kinetics and equilibrium processes associated to the adsorptive removal of pollutants using this type of materials is limited to fit the experimental data obtained to a particular set of models without verifying the validity and applicability of the underlying assumptions [19]. For this reason, in this work, we studied the absorptive removal of zinc(II) ions from synthetic wastewaters by using the previously mentioned wood residues to obtain magnetic adsorbents by magnetite nanoparticles impregnation. The obtained adsorbents were characterized by using several analytical techniques and the influence of the magnetite nanoparticles on the removal performance of the composites was determined. The experimental data obtained in batch tests was fitted to a series of isotherm and kinetics models, aiming to elucidate the nature and the factors that affect the ability of the composites for capturing the metal ions in the solution.

2. Materials and Methods

2.1. Materials

Iron(III) chloride hexahydrate (FeCl3∙6H2O) (LOBA Chemie, Mumbai, India), iron(II) sulfate heptahydrate (FeSO4∙7H2O) (Panreac, Barcelona, Spain) and ammonium hydroxide solution (NH4OH, 29%) (Fisher Chemicals, New York, NY, USA). These reagents were used to obtain magnetic nanoparticles. Zinc chloride (ZnCl2) (Merck, Kenilworth, NJ, USA) was used to prepare Zn2+ synthetic solutions for removal tests. All reagents were of analytical grade. The solutions were prepared with pure water (Type 2). Lignocellulosic residues (laurel, yellow canelo and eucalyptus wood chips) were obtained from a local sawmill in Quito, Ecuador.

2.2. Lignocellulosic Biomass Conditioning and Characterization

Wood chips of laurel, canelo and eucalyptus (separately) were washed several times with water and dried at 100 °C for 12 h, then reduced up to 0.5 mm with a knife-mill (Thomas Scientific, model 3379-K05, Philadelphia, PA, USA), and classified by sieving (vibratory sieve shaker, Retsch, AS 200 basic, Llanera, Spain) for 30 min in order to retain the fraction between 74 and 150 µm to be used as support of the magnetic nanocomposites.
Chemical composition of lignocellulosic residues was determined using ASTM standards methods: D1106-96 (lignin), D1109-84 (hemicellulose), ASTM D1110-84 and D1107-96 (extractive components). Cellulose content was determined by difference. Moisture, ash and volatile matter content were measured by ASTM standards (D4442-07, D1102-84 and ASTM E872–82, respectively, ASTM International, Lima, Peru). Fixed carbon content was determined by difference. Functional groups of lignocellulosic residues were identified by Fourier Transform Infrared Spectroscopy (FTIR, Perkin Elmer, São Paulo, Brazil, Spectrum 100/Spotlight 200, range: 4000 to 500 cm−1). Surface morphology and particle size were determined by Scanning Electron Microscopy (SEM, Tescan, Mira 3, Quito, Ecuador) and Transmission Electron Microscopy (TEM, FEI, G2 Spirit Twin, Madrid, Spain), respectively.

2.3. Magnetite Impregnation of Lignocellulosic Biomass and Characterization

Magnetite nanoparticles were synthesized by controlled precipitation method. FeCl3∙6H2O (6.1 g) and FeSO4∙7H2O (4.2 g) were added to 100 mL of water, heated to 35 °C in an argon atmosphere with stirring (1250 rpm). Next, 10 mL of 29% ammonia solution was added and the reaction was carried out for 30 min. Finally, the lignocellulosic residue was added to the resultant solution keeping it at 35 °C and constant stirring for another 30 min (magnetite/residue ratio: 1:1). The obtained black precipitate was naturally cooled down to room temperature and thoroughly rinsed with pure water until the unreacted reagents were removed and a neutral pH of the solution was reached. The precipitate was magnetically separated from the solution and dried at 70 °C for 12 h. This procedure was performed for each type of wood.
Magnetite impregnated sorbents were characterized by X-ray Diffraction (Panalytical, Empyrean XRD with a Cu Kα source, 40 kV, 40 mA, 2θ range: 10 to 90°, Quito, Ecuador), FTIR microscopy in ATR mode (Perkin Elmer, Spectrum 100 Spectrometer with Spotlight 200 microscope, São Paulo, Brazil), SEM (Tescan, Mira 3, Quito, Ecuador) and TEM (FEI, G2 Spirit Twin, Madrid, Spain) microscopy.

2.4. Adsorption Studies

Batch removal tests were performed at 298 K in order to optimize experimental conditions. Each of the obtained nanocomposites were added to zinc chloride synthetic solutions (ZnCl2) at constant stirring (150 rpm, orbital shaker KS 130, IKA, Staufen, Germany) and pH was adjusted at 6 [20], by addition of 0.1 M NaOH or HCl solutions. Zinc concentration after adsorption was analyzed by visible spectroscopy (HACH, DR1900, Loveland, CO, USA) as established in the APHA 3500—Zn method. All the tests were performed in duplicate [21].
The optimal contact time was determined using 50 mg L−1 zinc precursor solutions with adsorbent dosage of 5 g L−1. Samples were collected at 30, 45, 60, 120, 180, 240 and 300 min in order to determine evolution of removal process. Then, using the optimal contact time (3 h), the effect of adsorbent dosage on zinc removal was studied. The magnetite modified adsorbents were added to 50 mg L−1 zinc solutions at a dosage of 3, 5, 7 and 10 g L−1, at pH value of 6. Non-modified adsorbents were employed as reference samples, in the same conditions.
Adsorption isotherms were obtained using the optimal adsorbent dosage and contact time. Zinc chloride solutions of 25, 50, 100, 150, 200 and 300 mg L−1 [22] were used. The experiments were carried out at 298 K, constant stirring at 150 rpm and pH 6.
For the adsorption experiments, the following equation was employed to calculate the equilibrium adsorption capacity, q e (mg g−1) of adsorbent:
q e = V ( C o C e ) / m ,
where V (L) is the solution volume, C o and C e (mg L−1) are the initial and equilibrium concentrations of Zn solutions, respectively, and m (g) is the mass of dry sorbent. The removal efficiency is calculated as a ratio between adsorbed value at time t and the initial zinc concentration (mg L−1).
The Equations (2)–(4) for pseudo first-order (PFO) model (Equation (2)), pseudo second-order (PSO) model (Equation (3)) and intra-particle diffusion model [23] (Equation (3)):
q ( t ) = q e ( 1 e k 1 t ) ,
t / q ( t ) = 1 / ( k 2 q e 2 ) + t / q e
q ( t ) = k p t 0.5 + C
were used to analyze the kinetic data, where q(t) (mg g−1) is the amount of Zn adsorbed on the sorbent at any time t (h), k1 (min−1), k2 (g (mg min)−1) and kp (g mg−1 min−0.5) are the PFO, PSO and intraparticle diffusion kinetic model constants, respectively.
The adsorption isotherm data was fitted with Langmuir, Freundlich and Temkin models in their linearized forms as described by Equations (5)–(7) and the non-linear forms as given in Equations (8)–(10) [24], by using Origin 8.5 software (OriginLab Corp., Northampton, MA, USA):
C e q e = 1 ( q m × K L ) + 1 q m   C e
log q e = 1 n log C e + log K F
q e = B   lnK T + B   lnC e
q e = C e   q m   K L 1 + K L C e
q e = K F   C e 1 / n
q e = q m ln ( K T C e )
where qe is the amount of solute adsorbed at equilibrium per adsorbent unit weight (mg·g−1), Ce represents the concentration of solute at equilibrium (mg·L−1), qm represents the maximum adsorption capacity (mg g−1), KL is the Langmuir isotherm constant (L mg−1), KF is the Freundlich isotherm constant (L g−1), n is the coefficient of the Freundlich equation related to heterogeneity (measure of adsorption linearity), B is related to sorption heat and KT (L mg−1) is Temkin isotherm equilibrium binding constant. The sum of squares error (SSE) is employed to estimate the goodness of non-linear fit. The best model is characterized by the smallest SSE value.
The Dubinin–Radushkevich isotherm model was also evaluated by using the Equation (11):
ln   q e = ln q m β   ε 2 ,
where
ε = RT   ln ( 1 + 1 / C e ) ,
where qm (mg Zn g−1) is the maximum adsorption capacity, β (mol2 kJ−2) is a coefficient related to the mean sorption energy, ε is the Polanyi potential, R (J mol−1 K−1) is the gas constant and T (K) is the absolute temperature. The sorption energy E (kJ mol−1) is given by:
E = 1 / 2 β
The fittings of the experimental data collected to three different kinetics models aimed to study the adsorption rate and the rate-controlling step during the removal of the heavy metal ions. The first two models considered, PFO and PSO, are reaction kinetics models that are usually adequate for cases in which there is a high and low concentration of the pollutant in the solution, respectively. In the PSO model, chemisorption is supposed to be the rate-controlling step in the process. In the third model considered, the rate-controlling step is the intraparticle diffusion of the pollutant into the pores of the adsorbent. This model is usually adequate for porous adsorbents [25,26].

3. Results

3.1. Adsorbents Characterization

3.1.1. Chemical Composition of Lignocellulosic Biomass

In wood, the major component is cellulose and lignin [19]. Studies performed during the last decades [19,20], showed that lignin plays a significant role in the removal of heavy metals due to the presence of phenolic groups. However, this is not the predominant factor, since there are other components such as hemicellulose that have ionizable carboxylic groups in their composition that contribute to the adsorption of metal ions. As shown in Table 1, the values in this study, vary between 63 and 74% for cellulose, while lignin content varied from 18 to 32%.

3.1.2. Infrared Spectroscopy

Figure 1 shows FTIR spectra of the lignocellulosic residues and the magnetic composites. The laurel, canelo and eucalyptus residues showed similar spectra to literature [27] with bands located at 3332, 3341 and 3322 cm−1, respectively, which correspond to O–H bond stretches of phenol and alcohol groups. The bands at 2895, 2895 and 2903 cm−1 attributed to the C–H bonds stretches of the aliphatic chains, while 1733, 1734 and 1736 cm−1 bands correspond to C=O groups [28]. Characteristic bands of lignin for laurel, canelo and eucalyptus were seen at 1502, 1505 and 1502 cm−1, respectively, which correspond to the benzene ring vibration. Bands at 1733, 1734 and 1736 cm−1 are attributed to the carbonyl group vibration of cellulose and hemicellulose. Bands between 1154 and 1024 cm−1 represent the vibrational stretching of the C–OH groups. Although very weak, the bands at 577, 591 and 571 cm−1, respectively in the magnetic composites of laurel, canelo and eucalyptus could be attributed to the Fe–O bond in magnetite [5,12]. Such weak bands could be attributed to the region of the sample studied in ATR mode. Nevertheless, their presence could be also retrieved by a slight shift in the biomass bands upon modification. Moreover, as the microscopy and diffraction measurements further show, the presence of magnetite is confirmed.

3.1.3. Scanning and Transmission Electron Microscopy

Figure 2 shows SEM images of lignocellulosic residues and composite materials. Non-homogeneous surfaces with smooth sections are seen in the lignocellulosic residues. Magnetic composites show rough surfaces due to the presence of magnetite nanoparticles deposited on the lignocellulosic residue surfaces. Agglomerates of different size are formed due to the magnetic properties of the nanoparticles and the high surface area [29]. TEM images of magnetite particles (a) and magnetic composite of laurel (b) are shown in Figure 3. These pictures reveal nanoscaled quasi-spherical particles with particle diameter of 20 nm, approximately. In previous studies [6,30,31], for similar synthesis methods of magnetite nanoparticles, particles with size between 10 and 50 nm and a similar geometry to this study were obtained [32]. Furthermore, according to Shahid and Choi, the dissolution-recrystallization, occurring during coprecipitation process is responsible for the spherical shape of the magnetite nanoparticles [32]. An aggregation phenomenon is evident in the magnetite nanoparticles, which could be attributed to the small size and their magnetic properties. This phenomenon seems to be reduced in the composites as the particles are distributed on the lignocellulosic matrix surface.

3.1.4. X-ray Diffraction

Figure 4 shows the X-ray diffraction patterns of magnetite and magnetic composites. Peaks at 30.0°, 35.5°, 43.0°, 53.5°, 57.0° and 62.5° corresponding to indices (220), (311), (400), (422), (511) and (440), were observed. These peaks correspond to a cubic spinel structure of magnetite [29]. Based on these results, the presence of magnetite in the composite material can be verified. The same behavior was observed for canelo and eucalyptus magnetic composites. According to the Scherrer equation mentioned elsewhere [33], the crystallite size calculated for magnetite nanoparticles, and magnetite modified laurel, canelo and eucalyptus composites were 10 nm, 15.3 nm, 19.5 nm and 13.5 nm, respectively. These results are comparable to the sizes observed with TEM and reported scale values [32].

3.2. Adsorption Studies

3.2.1. Optimal Time and Adsorbent Dosage

The effect of contact time on adsorption of zinc was studied. As shown in Figure 5, the three magnetic composites have the same behavior. The adsorption rate was fast for the first minutes, after 60 min this rate decreases significantly. According to Zhang et al. [34], during the first few minutes, the adsorbent has many free active sites, and as the adsorption process progresses they are occupied by zinc ions. Furthermore, they generate repulsive forces towards the ions still present in the treated solution, which makes it difficult to occupy all the active sites on the adsorbent surface.
It was observed that after 60 min the Zn(II) removal exceeds 95% for the three composite materials and the equilibrium was achieved after 180 min. This time was considered as the optimal time for the next tests.
Figure 6 shows the effect of adsorbent dosage on zinc removal from synthetic solutions. It is observed for the three lignocellulosic residues (left) that the dosage has a direct influence in the adsorption efficiency. Laurel sawdust showed the maximum removal percentage (60.9%) of the three residues, which could be attributed to the lignin content of the lignocellulosic residue. For the magnetic composites (right), it is observed a maximum removal of 98.9, 98.8 and 97.6% for laurel, canelo and eucalyptus magnetic composites, respectively, when using 7 g L−1 of adsorbent. By using higher dosages of adsorbent (10 g L−1), the removal percentage does not increase which is attributed to a higher dispute over available active sites. Thus, the number of agglomerates increases, which causes the diffusion path to be longer and the adsorption process gets hindered [13].
When comparing the results obtained in the zinc removal, it is observed that the addition of magnetite nanoparticles to lignocellulosic residues improve notably the removal percentage. For laurel, the zinc removal increases from 60.9 to 98.9%. This behavior occurs due to the presence of nanoparticles, which considerably increases the surface area. Furthermore, the Fe and O atoms of the nanoparticles adsorb OH radicals from the solution, causing their surface to become rich in OH groups. This increases the interactions between these groups and the metal ions in the solution, which makes the removal process more efficient [31].

3.2.2. Adsorption Kinetics

As shown in Figure 7a–c, the Zn adsorption rates on the magnetite modified canelo, laurel and eucalyptus waste-based adsorbents increased rapidly in the initial stage and then gradually decreased reaching an equilibrium state. This behavior could be explained by the existence of numerous accessible active sites on the adsorbent surface and thus, when they get occupied, the adsorption rate decreases. The kinetics of the removal of Zn(II) ions using magnetite modified adsorbents was fitted by both pseudo-first order model (PFO) and pseudo-second order model (PSO) and the corresponding plots are given in Figure 7a–c. As it is clearly observed, the fitted values of q were closer to the experimental ones for PSO kinetic model for magnetite modified laurel and canelo adsorbents, which is indicative of the suitability of PSO model for the adsorption kinetics of Zn on the obtained bio-sorbents.
The parameters corresponding to the fitting of kinetic models are further given in the Table 2. The correlation coefficient, R2, obtained from the linear regression analysis was used to analyze the linear forms for PFO and PSO models. Thus, the higher values of R2 for PSO with respect to the PFO ones (above 0.963), confirm a better fit of the adsorption process by PSO model, independently of the employed adsorbent. This indicates that Zn adsorption process on each magnetic adsorbent could be controlled by chemisorption [35,36], which is also supported by the results of the adsorbent characterization (presence of active sites/functional groups). The fitted adsorption capacities are similar for the modified adsorbents, reaching a value of about 10 mg g−1. On the other hand, the parameter h increases from eucalyptus to canelo and laurel magnetite-modified adsorbents, indicating the initial adsorption has the highest rate for laurel-based adsorbent. Such results could be attributed to the biomass structure and functional groups [37]. The adsorption rate constants obtained from the plots were decreasing from 0.0141, 0.0116 and 0.0041 g mg−1 min−1 for laurel, canelo and eucalyptus, respectively.
Despite the good fit of PSO for laurel and canelo based sorbents, the low R2 values obtained for eucalyptus sorbent for both PFO and PSO models suggest that a different model could fit the adsorption behavior of the obtained materials, as the sorbent supports have similar characteristics. Thus, intraparticle diffusion model as proposed by Weber and Morris was employed. This model assumes that the rate is not limited by mass transfer from the bulk solution to the external surface of adsorbent as the solution is stirred during the sorption period [38]. Thus, the corresponding plot would be linear [39]. Figure 7d shows multi-linear plots for all sorbents, where the first step is attributed to fast (instantaneous) adsorption of Zn on the sorbent surface, the second step is characterized by a slower sorption due to diffusion of Zn ions into the inner adsorption sites of the adsorbent material, that would be the rate-limiting stage and a third step which corresponds to the equilibrium when the sorbent gets saturated with the ions. This behavior is similar to that of other adsorbents in literature [40,41]. As it can be observed in Figure 7c (kp, mg g−1 min−0.5), the intra-particle diffusion rate constant for the first step is the highest for laurel and followed by canelo and eucalyptus with similar values, that is initial adsorption process is faster on laurel-based sorbent. The second and third steps are characterized by similar trend of the constants. Since the constant value was markedly higher for the first step than the other two, the boundary diffusion is indicated as dominating the process followed by the intra-particle diffusion.

3.2.3. Adsorption Isotherms

In order to better understand the adsorption of Zn ions on the obtained magnetite modified adsorbents, linearized forms of Langmuir, Freundlich and Temkin models were used for fitting. The corresponding plots are showed in Figure 8, respectively. Langmuir model is employed for homogeneous adsorption which assumes monolayer adsorption in the absence of interaction or steric hindrance between adsorbed molecules. The D-R model considers on the other hand a porous structure of adsorbent with heterogeneous surface. The Freundlich model is suitable for the adsorption on heterogeneous surfaces. The Temkin model assumes uniform binding energies and that the sorption heat decreases linearly as the sorbent surface gets covered with the ions.
As it can be observed from Figure 8, the Langmuir, Freundlich and Temkin fittings resulted in correlation coefficient R2 above 0.9. The adsorption parameters from the linear fitting are further presented in Table 3.
By comparing the correlation coefficient values, one can observe larger values above 0.96 for Freundlich model which indicates better fitting with this model, independently of adsorbent. This suggests the Zn adsorption on the magnetite modified laurel, canelo and eucalyptus adsorbents takes place via multilayer adsorption on heterogeneous surfaces where active sites have different adsorption capacity. The fitting parameters show Freundlich constant (adsorption intensity) decreasing in the order eucalyptus > canelo > laurel for the magnetite modified adsorbents along with increasing n value, and thus indicating the Zn adsorption is more favorable on eucalyptus modified adsorbent.
From Langmuir fitting, it can be seen that qm slightly decreased from eucalyptus to laurel along with an increase in KL value. Temkin model can be employed in order to explain the interactions between Zn ions and the sorbent in regard to the heat of adsorption indicated by the B parameter. It was observed that KT was the highest for eucalyptus and laurel, while B decreased in the order eucalyptus, canelo, laurel magnetite-modified adsorbents, indicating physical adsorption process.
The order of the isotherm best fits in this work is Freundlich > Langmuir > Temkin. Although such isotherms are widely employed, they do not give information in regards to the mechanism of sorption. For this purpose, the Dubinin–Radushkevich (D-R) isotherm model employed for adsorption process taking place on both homogeneous and heterogeneous surfaces was further applied for modelling the adsorption process (see Table 3 for fitting parameters). The correlation fitting R2 value is below 0.8, indicating poor fitting with this model. The maximum adsorbed amount decreased from eucalyptus to laurel sorbent while and the sorption energy decreased in the same order and indicated a physical nature of the adsorption process [42]. The Xm value is slightly smaller than qm which could be attributed to the different definition of qm in the models.
On the other hand, given the errors introduced by the linearized forms of the isotherm models, non-linear regression is preferable and more relevant to estimate the isotherm model parameters. The plots are represented in Figure 9 showing the fitting curves for each model with respect to the experimental data. The fitting parameters are depicted in Table 4. It can be observed that the non-linear regression resulted in R2 values similar to those from linear regression. The SSE values decreased in the order Temkin < Langmuir < Freundlich for all sorbents. As the best fit is indicated by lower SSE values, it can be concluded that Freundlich non-linear model best fits the adsorption on the obtained materials.
To verify the suitability of the magnetite modified composite as adsorbents, a comparison with reports on other similar adsorbents in the literature was performed. Although direct comparison of the maximum adsorption capacity qm of studied lignocellulosic biomass-based adsorbents with literature is difficult, owing to the different experimental conditions, the maximum adsorption capacity obtained is in line with other reports on either same biomass sorbents (different analyte) [43] or adsorbents for Zn(II) obtained from various biomass such as Lagenaria vulgaris shell [44], modified orange peel [45] or [46] or magnetic water hyacinth [47]. Considering these findings and the cost-effective source of the obtained adsorbents, the laurel, canelo and eucalyptus biomass is considered a promising alternative to obtain highly effective adsorbents for water treatment.

4. Conclusions

Laurel, canelo and eucalyptus residues were used as support of magnetic composites to be used as adsorbents in zinc removal from synthetic solutions. XRD, FTIR, SEM and TEM analysis showed that magnetite nanoparticles with diameters of approximately 20 nm were obtained and that were effectively impregnated in lignocellulosic matrices. Impregnation of magnetite in the lignocellulosic matrices of laurel, canelo and eucalyptus increased the zinc removal in 79.7, 172.3 and 282.2%, respectively, in comparison to the lignocellulosic residues alone. The optimal contact time was 3 h and an adsorbent dosage of 7 g/L allowed the maximum removal of zinc, which was about 98.9, 98.8 and 97.6% for laurel, canelo and eucalyptus magnetic composites, respectively. The adsorption of zinc on magnetite modified adsorbents is best fit by the isotherm models: Freundlich > Langmuir > Temkin > Dubinin–Radushkevich, indicating multilayer adsorption on heterogeneous surfaces while kinetics models shown the adsorption process is dominated by both the chemisorption and diffusion.

Author Contributions

Conceptualization, N.M.R.-L. and C.A.; formal analysis, C.A. and S.G.; investigation, characterization A.D., C.A., N.M.R.-L., S.G. and V.H.G.; writing—original draft preparation, C.A. and S.G.; writing—review and editing N.M.R.-L., S.G., V.H.G. and A.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Escuela Politécnica Nacional, project grants PIMI-14-18 and PIGR-19-05.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors thank David Castro for useful discussion on adsorption analysis.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Singh, J.; Kalamdhad, A.S. Effects of Heavy Metals on Soil, Plants, Human Health and Aquatic Life Making bricks using variety of solid waste View project Anaerobic digestion View project. Int. J. Res. Chem. Environ. 2011, 1, 15–21. [Google Scholar]
  2. Ali, A.; Saeed, K.; Mabood, F. Removal of chromium (VI) from aqueous medium using chemically modified banana peels as efficient low-cost adsorbent. Alexandria Eng. J. 2016, 55, 2933–2942. [Google Scholar] [CrossRef] [Green Version]
  3. Setyono, D.; Valiyaveettil, S. Chemically modified sawdust as renewable adsorbent for arsenic removal from water. ACS Sustain. Chem. Eng. 2014, 2, 2722–2729. [Google Scholar] [CrossRef]
  4. Miretzky, P.; Cirelli, A.F. Cr(VI) and Cr(III) removal from aqueous solution by raw and modified lignocellulosic materials: A review. J. Hazard. Mater. 2010, 180, 1–19. [Google Scholar] [CrossRef]
  5. Panneerselvam, P.; Morad, N.; Tan, K.A. Magnetic nanoparticle (Fe3O4) impregnated onto tea waste for the removal of nickel(II) from aqueous solution. J. Hazard. Mater. 2011, 186, 160–168. [Google Scholar] [CrossRef] [PubMed]
  6. Sureshkumar, V.; Kiruba Daniel, S.C.G.; Ruckmani, K.; Sivakumar, M. Fabrication of chitosan–magnetite nanocomposite strip for chromium removal. Appl. Nanosci. 2016, 6, 277–285. [Google Scholar] [CrossRef] [Green Version]
  7. Abdolali, A.; Guo, W.S.; Ngo, H.H.; Chen, S.S.; Nguyen, N.C.; Tung, K.L. Typical lignocellulosic wastes and by-products for biosorption process in water and wastewater treatment: A critical review. Bioresour. Technol. 2014, 160, 57–66. [Google Scholar] [CrossRef] [PubMed]
  8. Owalude, S.O.; Tella, A.C. Removal of hexavalent chromium from aqueous solutions by adsorption on modified groundnut hull. Beni-Suef Univ. J. Basic Appl. Sci. 2016, 5, 377–388. [Google Scholar] [CrossRef] [Green Version]
  9. Gakwisiri, C.; Raut, N.; Al-Saadi, A.; Al-Aisri, S.; Al-Ajmi, A. A critical review of removal of zinc from wastewater. Lect. Notes Eng. Comput. Sci. 2012, 2197, 627–630. [Google Scholar]
  10. Hegazi, H.A. Removal of heavy metals from wastewater using agricultural and industrial wastes as adsorbents. HBRC J. 2013, 9, 276–282. [Google Scholar] [CrossRef] [Green Version]
  11. Hua, M.; Zhang, S.; Pan, B.; Zhang, W.; Lv, L.; Zhang, Q. Heavy metal removal from water/wastewater by nanosized metal oxides: A review. J. Hazard. Mater. 2012, 211–212, 317–331. [Google Scholar] [CrossRef]
  12. Cheng, Z.; Gao, Z.; Ma, W.; Sun, Q.; Wang, B.; Wang, X. Preparation of magnetic Fe3O4 particles modified sawdust as the adsorbent to remove strontium ions. Chem. Eng. J. 2012, 209, 451–457. [Google Scholar] [CrossRef]
  13. Moafi, H.F.; Ansari, R.; Sadeghnia, S. Preparation of wood sawdust/Fe2O3 nanocomposite and its application for arsenic (III) ion removal from aqueous solutions. Cellul. Chem. Technol. 2018, 52, 271–282. [Google Scholar]
  14. Gupta, V.K.; Nayak, A. Cadmium removal and recovery from aqueous solutions by novel adsorbents prepared from orange peel and Fe2O3 nanoparticles. Chem. Eng. J. 2012, 180, 81–90. [Google Scholar] [CrossRef]
  15. CFN. Ficha Sectorial: Explotación De Viveros Forestales y Madera En Pie; CFN: Quito, Ecuador, 2017. [Google Scholar]
  16. Daian, G.; Ozarska, B. Wood waste management practices and strategies to increase sustainability standards in the Australian wooden furniture manufacturing sector. J. Clean. Prod. 2009, 17, 1594–1602. [Google Scholar] [CrossRef]
  17. Chompu-Inwai, R.; Jaimjit, B.; Premsuriyanunt, P. A combination of Material Flow Cost Accounting and design of experiments techniques in an SME: The case of a wood products manufacturing company in northern Thailand. J. Clean. Prod. 2015, 108, 1352–1364. [Google Scholar] [CrossRef]
  18. Zamora-Ledezma, C.; Negrete-Bolagay, D.; Figueroa, F.; Zamora-Ledezma, E.; Ni, M.; Alexis, F.; Guerrero, V.H. Heavy metal water pollution: A fresh look about hazards, novel and conventional remediation methods. Press 2021, 22, 101504. [Google Scholar] [CrossRef]
  19. Hubbe, M.A.; Azizian, S.; Douven, S. Implications of Apparent Pseudo-Second-Order Adsorption Kinetics onto Cellulosic Materials: A Review. BioResources 2019, 14, 7582–7626. [Google Scholar] [CrossRef]
  20. Sciban, M.; Klasnja, M.; Skrbic, B. Modified hardwood sawdust as adsorbent of heavy metal ions from water. Wood Sci. Technol. 2006, 40, 217–227. [Google Scholar] [CrossRef]
  21. Šćiban, M.; Radetić, B.; Kevrešan, Ž.; Klašnja, M. Adsorption of heavy metals from electroplating wastewater by wood sawdust. Bioresour. Technol. 2007, 98, 402–409. [Google Scholar] [CrossRef]
  22. Bulut, Y.; Tez, Z. Removal of heavy metals from aqueous solution by sawdust adsorption. J. Environ. Sci. 2007, 19, 160–166. [Google Scholar] [CrossRef]
  23. Weber, W.J., Jr.; Morris, J.C. Kinetics of Adsorption on Carbon from Solution. J. Sanit. Eng. Div. 1963, 89, 31–59. [Google Scholar] [CrossRef]
  24. Okumu, F. A Novel Polyaniline Titanium Oxide Sawdust Composite Adsorbent for Polychlorinated Biphenyls. Sci. J. Chem. 2013, 1, 29. [Google Scholar] [CrossRef] [Green Version]
  25. McCoy, M.A.; Liapis, A.I. Evaluation of kinetic models for biospecific adsorption and its implications for finite bath and column performance. J. Chromatogr. A 1991, 548, 25–60. [Google Scholar] [CrossRef]
  26. Almeida-naranjo, C.E.; Belén, M.; Cabrera, G.; Guerrero, V.H. Caffeine removal from synthetic wastewater using magnetic fruit peel composites: Material characterization, isotherm and kinetic studies. Environ. Challenges 2021, 5, 100343. [Google Scholar] [CrossRef]
  27. Azzaz, A.A.; Jellali, S.; Assadi, A.A.; Bousselmi, L. Chemical treatment of orange tree sawdust for a cationic dye enhancement removal from aqueous solutions: Kinetic, equilibrium and thermodynamic studies. Desalin. Water Treat. 2016, 57, 22107–22119. [Google Scholar] [CrossRef]
  28. Md Salim, R.; Asik, J.; Sarjadi, M.S. Chemical functional groups of extractives, cellulose and lignin extracted from native Leucaena leucocephala bark. Wood Sci. Technol. 2021, 55, 295–313. [Google Scholar] [CrossRef]
  29. Petcharoen, K.; Sirivat, A. Synthesis and characterization of magnetite nanoparticles via the chemical co-precipitation method. Mater. Sci. Eng. B Solid-State Mater. Adv. Technol. 2012, 177, 421–427. [Google Scholar] [CrossRef]
  30. Tajabadi, M.; Khosroshahi, M.E. Effect of Alkaline Media Concentration and Modification of Temperature on Magnetite Synthesis Method Using FeSO4/NH4OH. Int. J. Chem. Eng. Appl. 2012, 3, 206–210. [Google Scholar] [CrossRef]
  31. Nethaji, S.; Sivasamy, A.; Mandal, A.B. Preparation and characterization of corn cob activated carbon coated with nano-sized magnetite particles for the removal of Cr(VI). Bioresour. Technol. 2013, 134, 94–100. [Google Scholar] [CrossRef]
  32. Shahid, M.K.; Choi, Y. Characterization and application of magnetite Particles, synthesized by reverse coprecipitation method in open air from mill scale. J. Magn. Magn. Mater. 2020, 495, 165823. [Google Scholar] [CrossRef]
  33. Fujimoto, K.; Oku, T.; Akiyama, T. Fabrication and Characterization of ZnO/Cu2O Solar Cells Prepared by Electrodeposition. Appl. Phys. Express 2013, 6, 086503. [Google Scholar] [CrossRef]
  34. Zhang, Y.; Zheng, R.; Zhao, J.; Zhang, Y.; Wong, P.K.; Ma, F. Biosorption of zinc from aqueous solution using chemically treated rice husk. Biomed Res. Int. 2013, 2013, 365163. [Google Scholar] [CrossRef] [PubMed]
  35. Feng, N.C.; Guo, X.Y.; Liang, S. Kinetic and thermodynamic studies on biosorption of Cu(II) by chemically modified orange peel. Trans. Nonferrous Met. Soc. China Engl. Ed. 2009, 19, 1365–1370. [Google Scholar] [CrossRef]
  36. Zhao, H.; Liu, X.; Cao, Z.; Zhan, Y.; Shi, X.; Yang, Y.; Zhou, J.; Xu, J. Adsorption behavior and mechanism of chloramphenicols, sulfonamides, and non-antibiotic pharmaceuticals on multi-walled carbon nanotubes. J. Hazard. Mater. 2016, 310, 235–245. [Google Scholar] [CrossRef]
  37. Ertugay, N.; Bayhan, Y.K. The removal of copper (II) ion by using mushroom biomass (Agaricus bisporus) and kinetic modelling. Desalination 2010, 255, 137–142. [Google Scholar] [CrossRef]
  38. Zhu, S.; Hou, H.; Xue, Y. RETRACTED: Kinetic and isothermal studies of lead ion adsorption onto bentonite. Appl. Clay Sci. 2008, 40, 171–178. [Google Scholar] [CrossRef]
  39. Nandi, B.K.; Goswami, A.; Purkait, M.K. Adsorption characteristics of brilliant green dye on kaolin. J. Hazard. Mater. 2009, 161, 387–395. [Google Scholar] [CrossRef]
  40. Behnamfard, A.; Salarirad, M.M. Equilibrium and kinetic studies on free cyanide adsorption from aqueous solution by activated carbon. J. Hazard. Mater. 2009, 170, 127–133. [Google Scholar] [CrossRef] [PubMed]
  41. Tan, I.A.W.; Ahmad, A.L.; Hameed, B.H. Adsorption of basic dye on high-surface-area activated carbon prepared from coconut husk: Equilibrium, kinetic and thermodynamic studies. J. Hazard. Mater. 2008, 154, 337–346. [Google Scholar] [CrossRef]
  42. Wang, C.C.; Juang, L.C.; Hsu, T.C.; Lee, C.K.; Lee, J.F.; Huang, F.C. Adsorption of basic dyes onto montmorillonite. J. Colloid Interface Sci. 2004, 273, 80–86. [Google Scholar] [CrossRef]
  43. León, G.; Aldás, M.; Guerrero, V.; Landázuri, A.; Almeida-Naranjo, C. Caffeine and irgasan removal from water using bamboo, laurel and moringa residues impregnated with commercial TiO2 nanoparticles. MRS Adv. 2019, 4, 3553–3567. [Google Scholar] [CrossRef]
  44. Mitic-Stojanovic, D.-L.; Bojic, D.; Mitrovic, J.; Andjelkovic, T.; Radovic, M.; Bojic, A. Equilibrium and kinetic studies of Pb(II), Cd(II) and Zn(II) sorption by Lagenaria vulgaris shell. Chem. Ind. Chem. Eng. Q. 2012, 18, 563–576. [Google Scholar] [CrossRef]
  45. Feng, N.C.; Guo, X.Y. Characterization of adsorptive capacity and mechanisms on adsorption of copper, lead and zinc by modified orange peel. Trans. Nonferrous Met. Soc. China Engl. Ed. 2012, 22, 1224–1231. [Google Scholar] [CrossRef]
  46. Castro, D.; Rosas-Laverde, N.M.; Aldás, M.B.; Almeida-naranjo, C.; Guerrero, V.H.; Pruna, A.I. Chemical modification of agro-industrial waste-based bioad- sorbents for enhanced removal of Zn(II) ions from aqueous solutions. Materials 2021, 14, 2134. [Google Scholar] [CrossRef]
  47. Nyamunda, B.C.; Chivhanga, T.; Guyo, U.; Chigondo, F. Removal of Zn (II) and Cu (II) ions from industrial wastewaters using magnetic biochar derived from water hyacinth. J. Eng. 2019, 2019, 5656983. [Google Scholar] [CrossRef] [Green Version]
Figure 1. FTIR spectra of lignocellulosic biomass before (x) and after impregnation with magnetite (x*): (a) (a* vs. a) laurel, (b) (b* vs. b) canelo and (c) (c* vs. c) eucalyptus.
Figure 1. FTIR spectra of lignocellulosic biomass before (x) and after impregnation with magnetite (x*): (a) (a* vs. a) laurel, (b) (b* vs. b) canelo and (c) (c* vs. c) eucalyptus.
Materials 15 00728 g001
Figure 2. SEM image of non-impregnated lignocellulosic biomass from: (a) eucalyptus, (b) canelo and (c) laurel, and their magnetite impregnated counterparts: (df).
Figure 2. SEM image of non-impregnated lignocellulosic biomass from: (a) eucalyptus, (b) canelo and (c) laurel, and their magnetite impregnated counterparts: (df).
Materials 15 00728 g002
Figure 3. TEM images of the magnetite nanoparticles (a) and magnetite modified sorbent: canelo (b), eucalyptus (c) and laurel (d).
Figure 3. TEM images of the magnetite nanoparticles (a) and magnetite modified sorbent: canelo (b), eucalyptus (c) and laurel (d).
Materials 15 00728 g003
Figure 4. X-ray diffraction pattern of magnetite nanoparticles and magnetic composites.
Figure 4. X-ray diffraction pattern of magnetite nanoparticles and magnetic composites.
Materials 15 00728 g004
Figure 5. Contact time effect on zinc removal using magnetic composites.
Figure 5. Contact time effect on zinc removal using magnetic composites.
Materials 15 00728 g005
Figure 6. Effect of adsorbent dosage on zinc removal using lignocellulosic biomass before and after magnetite impregnation.
Figure 6. Effect of adsorbent dosage on zinc removal using lignocellulosic biomass before and after magnetite impregnation.
Materials 15 00728 g006
Figure 7. Adsorption kinetics of Zn on magnetite modified lignocellulosic waste-based adsorbents: laurel (a), canelo (b) and eucalyptus (c). Weber–Morris plot for intraparticle diffusion kinetics of Zn adsorption on magnetite eucalyptus modified adsorbent (d).
Figure 7. Adsorption kinetics of Zn on magnetite modified lignocellulosic waste-based adsorbents: laurel (a), canelo (b) and eucalyptus (c). Weber–Morris plot for intraparticle diffusion kinetics of Zn adsorption on magnetite eucalyptus modified adsorbent (d).
Materials 15 00728 g007
Figure 8. Linear regression analysis of Zn ions adsorption of magnetite modified adsorbents by Langmuir (a); Freundlich (b); Temkin (c) and Dubinin–Radushkevich (d), (ε2 value multiplied with 103) models.
Figure 8. Linear regression analysis of Zn ions adsorption of magnetite modified adsorbents by Langmuir (a); Freundlich (b); Temkin (c) and Dubinin–Radushkevich (d), (ε2 value multiplied with 103) models.
Materials 15 00728 g008
Figure 9. Non-linear regression models of Zn ions adsorption on magnetite modified adsorbents: laurel (a), canelo (b) and eucalyptus (c).
Figure 9. Non-linear regression models of Zn ions adsorption on magnetite modified adsorbents: laurel (a), canelo (b) and eucalyptus (c).
Materials 15 00728 g009
Table 1. Chemical composition and proximate analysis of the lignocellulosic biomass.
Table 1. Chemical composition and proximate analysis of the lignocellulosic biomass.
ComponentLaurelCaneloEucalyptus
Chemical compositionCellulose [wt %]50.2053.4958.44
Lignin [wt %]31.8428.2818.28
Hemicellulose [wt %]12.7612.4516.01
Extractives [wt %]5.205.787.27
Proximate analysisVolatile matter [wt %]70.9478.7081.52
Fixed carbon [wt %]17.9912.9312.54
Ash content [wt %]1.380.160.25
Moisture [wt %]9.698.215.69
Table 2. Kinetics parameters for Zn(II) ions removal using the magnetite impregnated adsorbents.
Table 2. Kinetics parameters for Zn(II) ions removal using the magnetite impregnated adsorbents.
Lignocellulosic BiomassPseudo First-OrderPseudo Second-Order
q e   [ m g   Z n g ] k 1   [ m i n 1 ] R2 q e   [ m g   Z n g ] k 2   [ g m g   Z n × m i n ] h   [ m g   Z n g × m i n ] R2
Laurel8.6020.03230.87310.0130.01411.4120.998
Canelo6.9130.02960.94610.0520.01161.1710.999
Eucalyptus7.810.02750.94710.6130.00410.4560.963
Table 3. Linear adsorption isotherm fitting parameters on magnetic adsorbents.
Table 3. Linear adsorption isotherm fitting parameters on magnetic adsorbents.
Isotherm Model/AdsorbentParameterLaurelCaneloEucalyptus
Langmuirqm, mg g−10.110.150.31
KL (L mg−1)23.820.77.9
R20.9670.9130.971
Freundlichqm, mg g−12.312.612.18
KF, (L g−1)0.120.180.5
n1.251.392.31
R20.9890.9670.961
TemkinKT, (L mg−1)0.890.380.95
B0.690.840.58
R20.9140.9260.931
D-RXm, (mg g−1)1.361.681.78
E, (kJ mol−1)0.3530.370.452
R20.6160.6330.805
Table 4. Non-linear adsorption isotherm fitting parameters on magnetic adsorbents.
Table 4. Non-linear adsorption isotherm fitting parameters on magnetic adsorbents.
Isotherm Model/AdsorbentParameterLaurelCaneloEucalyptus
Langmuirqm, mg g−164.22.7
KL (L mg−1)0.0120.0350.1
R20.9710.9560.978
SSE0.0650.10.078
Freundlichqm, mg g−12.312.622.24
KF, (L g−1)0.160.210.454
n1.521.482.31
R20.9890.9720.966
SSE0.0460.080.039
TemkinKT, (L mg−1)0.340.540.951
B0.6930.720.585
R20.9140.9210.931
SSE0.1820.2480.081
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Asimbaya, C.; Rosas-Laverde, N.M.; Galeas, S.; Debut, A.; Guerrero, V.H.; Pruna, A. Magnetite Impregnated Lignocellulosic Biomass for Zn(II) Removal. Materials 2022, 15, 728. https://doi.org/10.3390/ma15030728

AMA Style

Asimbaya C, Rosas-Laverde NM, Galeas S, Debut A, Guerrero VH, Pruna A. Magnetite Impregnated Lignocellulosic Biomass for Zn(II) Removal. Materials. 2022; 15(3):728. https://doi.org/10.3390/ma15030728

Chicago/Turabian Style

Asimbaya, Christopher, Nelly Maria Rosas-Laverde, Salome Galeas, Alexis Debut, Victor H. Guerrero, and Alina Pruna. 2022. "Magnetite Impregnated Lignocellulosic Biomass for Zn(II) Removal" Materials 15, no. 3: 728. https://doi.org/10.3390/ma15030728

APA Style

Asimbaya, C., Rosas-Laverde, N. M., Galeas, S., Debut, A., Guerrero, V. H., & Pruna, A. (2022). Magnetite Impregnated Lignocellulosic Biomass for Zn(II) Removal. Materials, 15(3), 728. https://doi.org/10.3390/ma15030728

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop