BioMed Research International

BioMed Research International / 2015 / Article

Research Article | Open Access

Volume 2015 |Article ID 126298 |

Md. Sayedur Rahman, Kathiresan V. Sathasivam, "Heavy Metal Adsorption onto Kappaphycus sp. from Aqueous Solutions: The Use of Error Functions for Validation of Isotherm and Kinetics Models", BioMed Research International, vol. 2015, Article ID 126298, 13 pages, 2015.

Heavy Metal Adsorption onto Kappaphycus sp. from Aqueous Solutions: The Use of Error Functions for Validation of Isotherm and Kinetics Models

Academic Editor: José L. Campos
Received17 Mar 2015
Revised29 Jun 2015
Accepted30 Jun 2015
Published30 Jul 2015


Biosorption process is a promising technology for the removal of heavy metals from industrial wastes and effluents using low-cost and effective biosorbents. In the present study, adsorption of Pb2+, Cu2+, Fe2+, and Zn2+ onto dried biomass of red seaweed Kappaphycus sp. was investigated as a function of pH, contact time, initial metal ion concentration, and temperature. The experimental data were evaluated by four isotherm models (Langmuir, Freundlich, Temkin, and Dubinin-Radushkevich) and four kinetic models (pseudo-first-order, pseudo-second-order, Elovich, and intraparticle diffusion models). The adsorption process was feasible, spontaneous, and endothermic in nature. Functional groups in the biomass involved in metal adsorption process were revealed as carboxylic and sulfonic acids and sulfonate by Fourier transform infrared analysis. A total of nine error functions were applied to validate the models. We strongly suggest the analysis of error functions for validating adsorption isotherm and kinetic models using linear methods. The present work shows that the red seaweed Kappaphycus sp. can be used as a potentially low-cost biosorbent for the removal of heavy metal ions from aqueous solutions. Further study is warranted to evaluate its feasibility for the removal of heavy metals from the real environment.

1. Introduction

Heavy metal pollution due to rapid urbanization and industrialization is one of the most significant environmental problems. Heavy metals are released into the aquatic environment from several domestic (automobile exhaust, smelting processes, burning of fossil fuels, incineration of wastes, landfill leaches, use of sewage sludge, municipal wastewater, and urban runoff) and industrial processes (electroplating, refining ore, mining, electronic and metal-finishing industries, fertilizer industry, tanneries, painting, paper industries, and pesticides) [1]. Heavy metals have become a global issue of environment and public health concern due to their toxicities, bioaccumulation in human body and food chain, carcinogenicities, and mutagenesis in various living organisms [24].

Numerous methods such as chemical precipitation, ion exchange, coagulation-flocculation, flotation, membrane filtration, electrochemical treatment, magnetic separation and purification, biosorption, and nanotechnology are being used to treat or remove heavy metals from water and wastewater [1, 5, 6]. Among them, biosorption has been regarded as a promising cost-effective, sustainable, and ecofriendly technology for the removal of different types of organic and inorganic pollutants from water and wastewater [7]. Moreover, this process offers a number of advantages in comparison to the conventional methods [8].

A wide range of commercial and potentially low-cost adsorbents including living or dead microorganisms, seaweeds, plant materials, industrial and agricultural wastes, natural residues, and inorganic precursors including red mud, clays, blast furnace slags, zeolites, chitosan, and peat has been reported in literature [812]. Seaweeds are widely distributed in marine, freshwater, and terrestrial ecosystems, which can serve as good biosorbents due to their abundance, cost-effectiveness, reusability, and high metal sorption capacities [7, 10, 11]. Despite that fact that the red algae constitute carrageenan that provides different binding sites (e.g., hydroxyl, carboxyl, amino, and sulfhydryl) responsible for the adsorption for heavy metals [9], they are the least focused group [13]. Therefore, further research studies are warranted on the selectivity of algal species [9].

The red seaweed Kappaphycus sp. is one of the most important commercial sources of kappa-carrageenan, which also has different medicinal and industrial applications [14]. Malaysia produced 331,490 tonnes of Kappaphycus sp., being 17.039% of the total world production in 2012 [15]. Recent studies suggest that both the living biomass and the waste biomass of Kappaphycus alvarezii are a good biosorbent for the removal of nutrients [13] and heavy metals from the aqueous environment [16, 17]. However, there is no available literature report on the biosorption of heavy metals using the dry biomass of Kappaphycus sp. Hence, the present study was investigated to study the performance of the dried biomass of Kappaphycus sp. for the removal of Zn2+, Cu2+, Pb2+, and Fe2+ from aqueous solutions in batch system at laboratory scale under different parameters like pH of solution, contact time, temperature, and initial metal ion concentrations. The Fourier transform infrared (FTIR) spectral analysis was made to identify the main functional groups involved in the biosorption process of those metal ions.

In general, mechanistic or empirical equations are used to express heavy metal adsorption capacities of different types of biosorbents using batch or column method [18]. Available literature reports confirm that nearly two dozens of empirical models involving 2, 3, 4, or even 5 parameters have been used to fit batch equilibrium isotherm curves to biosorbents [18, 19]. Besides, kinetic models have been described by several authors elsewhere [2022]. The equilibrium and kinetic models are often validated on the basis of coefficient of regression () of the experimental data. In the present study, some error functions have been used for validating the experimental data along with an insight into the usual measures of model inferences.

2. Materials and Methods

2.1. Collection, Identification, and Preparation of the Biosorbent

The red seaweed Kappaphycus sp. was collected from the Semporna coast of Sabah, Malaysia, in April 2013. The alga was washed for several times with running water and subsequently with deionized water to remove epiphytes and salts. The washed biomass was then dried in an oven at 60°C for 48 h until a constant weight was attained. The dried biomass was then crushed with an analytical mill, sieved (250 μm size), and stored in polypropylene bottles until use.

The living biomass of the species was preliminarily identified following systemic morphological features [23]. It was then subjected to 28S DNA based molecular identification [24]. The species was identified as Kappaphycus sp. and the gene sequence of the nucleotide was submitted in the NCBI GenBank (accession number KM229320).

2.2. Chemicals and Reagents

All the chemicals and reagents used in this study were of analytic reagent grade. The working solutions of different concentrations (10–200 mg L−1) of the heavy metals (Zn, Cu, Pb, and Fe) were prepared by diluting the stand solutions ( mg L−1) of the metals (Merck, Germany) in double distilled deionized water. Different initial pH of the solutions was obtained by adding 0.1 N HCl (Sigma-Aldrich, USA) or 0.1 N NaOH (Merck, Germany).

2.3. Batch Biosorption Experiments

All the experiments were conducted in a batch system using 150 mL Erlenmeyer flasks in a thermostatic shaker (25°C, 180 rpm), unless otherwise stated. Each flask was filled with 50 mL of solution and biosorbent as appropriate. The influence of several operational parameters on the biosorption characteristics of the metals such as pH of the aqueous solution (2–7), contact time (0−120 min), initial metal ion concentration (25−200 mg L−1), and temperature (25−50°C) were assessed using a constant biomass dosage (4 g L−1). Competitive adsorption of the four metal ions under mixed condition was also evaluated.

The adsorption studies were conducted with 50 mL of the metal solutions at an initial concentration of 10 mg L−1. For the kinetic studies sample solutions were withdrawn at regular intervals and the residual concentration of the heavy metals in the aqueous phase was analyzed after filtration as stated above.

The amount of the metal ions remaining in the solutions was measured by using Atomic Absorption Spectrometer (AAnalyst700, Perkin-Elmer, USA) after separation of the biosorbent by filtration through Whatman Filter number 1.

The amount of metal adsorbed per gram of the biosorbent at equilibrium, (mg g−1), was calculated from the difference of the metal concentration in the aqueous phase before and after biosorption as follows:where and are the initial and equilibrium concentration of metal ions in the solution (mg/L), respectively, is the volume of metal solution (L), and is the mass of the dry biosorbent (g).

The percentage of metal removal (, %) from the solution was calculated as follows:

Each experiment was done in triplicate and the data were expressed as the mean of the triplicate results. Statistical analyses were performed using Microsoft Office Excel 2007 (Microsoft Corp., USA).

2.4. Application of Adsorption Models

In the present experiment, four two-parameter isotherm models: Langmuir, Freundlich, Temkin, and Dubinin-Radushkevich (D–R); four two-parameter kinetic models: pseudo-first-order (PFO), pseudo-second-order (PSO), Elovich, and intraparticle diffusion (IpD); and kinetic model were applied to describe the sorption behaviour of the adsorbent. The equations of these models are given in Table 1.











(mg L−1): adsorbate initial concentration, (mg L−1): adsorbate equilibrium concentration, (mg g−1): observed biosorption capacity at equilibrium, (mg g−1): maximum biosorption capacity, (L mg−1): Langmuir constant related to the energy of adsorption, (): a dimensionless constant, known as separation factor, (mg g−1) (L mg−1): Freundlich isotherm constant related to the sorption capacity, : a constant which gives an idea of the grade of heterogeneity, (8.314 J mo−1): universal gas constant, (°K): absolute temperature, (L mg−1): equilibrium binding constant corresponding to the maximum binding energy, (J mol−1): Temkin constant related to heat of sorption, (mol2 kJ−2): Dubinin-Radushkevich isotherm constant, : Polanyi potential related to the equilibrium concentration, (kJ mol−1): mean free energy of biosorption, (mg g−1): equilibrium adsorption uptake at time, , (min−1): pseudo-first-order rate constant of adsorption, (g mg−1 min−1): pseudo-second-order rate constant of adsorption, (mg g−1 min−1): initial adsorption rate, (min): half-adsorption time, (mg g−1 min−1): initial adsorption rate constant, β (g mg−1): desorption constant, (mg g−1): boundary layer diffusion effect, (mg g−1 min−0.5): rate constant for intraparticle diffusion, (kJ mol−1): change in Gibbs free energy, (kJ mol−1): change in enthalpy, (kJ mol−1 K−1): change in entropy, and : thermodynamic equilibrium constant.
2.5. Error Function Analysis

In order to evaluate the suitability of the equation to the experimental data error function is the best optimization procedure. Apart from the regression coefficient (), nine error functions such as sum of square error (SSE), average relative error (ARE), hybrid functional error (HYBRID), sum of absolute error (EABS), Marquardt’s percent standard deviation (MPSD), normalized standard deviation ((%)), coefficient of determination (), nonlinear chi-square test (), and residual root mean square error (RMSE) were calculated to evaluate the best fit of the modeled equation to the experimental data. The equations of the error functions are given in Table 2.











(mg g−1): value obtained from the batch experiment, (mg g−1): calculated value from the isotherm for corresponding , (mg g−1): mean of , : number of observations in the experimental isotherm, and : number of parameters in the respective model.
2.6. FTIR Analysis

FTIR spectral analysis was carried out to determine the possible functional groups present in the dried biomass of Kappaphycus sp. Infrared spectra of the raw and metal-loaded biomass were obtained using a Fourier transform infrared (FTIR) spectrometer (Spectrum GX, Perkin-Elmer, USA).

3. Results and Discussion

3.1. Effect of Solution pH

In the adsorption process of metal ions from aqueous solutions, pH of the solution plays an important role. It is apparent from the results represented graphically in Figure 1(a) that, with the increase in pH, the biosorption increased gradually. The maximum biosorption (54.13%, 81.84%, 84.17%, and 20.94% for Zn2+, Cu2+, Pb2+, and Fe2+, resp.) was observed at pH 5. At lower pH (2–4), biosorption of metal ions was inhibited greatly. This can be explained on the basis that cell wall of the Kappaphycus sp. contains various functional groups (as described in Section 3.9). The positively charged functional groups increase competition between protons and metal cations for binding active sites of biomass, resulting in decreasing the metal cations adsorption on the biomass surfaces [48]. At higher pH values (6–8), the biosorption efficiency of metal ions was significantly decreased (Figure 1(a)), which may be attributed to the formation of anionic hydroxide complexes that decrease the dissolved metal concentration in solution and their competition with the active sites [49]. Therefore, all the rest biosorption experiments were carried out at pH 5.

3.2. Effect of Contact Time

As shown in Figure 1(b), biosorption of metal ions on the adsorbent increased with an increase in contact time and the equilibrium biosorption was attained within 90–120 min reflecting rapid initial biosorption. Maximum uptake of Cu2+, Pb2+, Fe2+, and Zn2+ was reached up to 83.88, 85.89, 21.27, and 54.13%, respectively, within 90 min. A decrease in the biosorption was noticed during the subsequent time of incubation indicating the maximum adsorption level as a saturation point of biosorption. The rapid kinetic mechanism can be attributed to the formation of exterior surface complexes neglecting intraparticle diffusion, which is very advantageous in biotechnological processes for wastewater treatment [48]. In general, heavy metal biosorption efficiency of seaweeds attained a maximum level within 30 and 90 min [50]. Hence, a contact time of 120 min was selected for further experiments ensuring attainment of equilibrium conditions.

3.3. Effect of Temperature

The solution temperature plays a vital role on the metal ions biosorption, which was found to increase with the increase of solution temperature. The rate of Zn2+, Cu2+, Pb2+, and Fe2+ biosorption by the dried biomass of Kappaphycus sp. was rapid reaching a maximum of 61.74, 89.12, 86.68, and 40.06%, respectively, at 50°C. This phenomenon indicates that the biosorption process of the metal ions onto the biomass is endothermic. It can be attributed that, at the higher temperatures, the activation of the biosorbent surfaces is enlarged facilitating more active sites for biosorption of the metal ions. Moreover, an easy mobility and enhanced accessibility of metal ions from the bulk solution to the biomass active sites could also be the possible reason for the maximum biosorption of metal ions at higher temperatures [51].

3.4. Effect of Initial Metal Ion Concentration

Biosorption capacity of the biomass was found to increase with increasing initial concentration of the metal ions. This phenomenon can be attributed to an increase in electrostatic interactions involving sites of progressively lower affinity for the metal ions up to the point of saturation [51, 52]. It was further observed that the percentage removal of the metal ions decreased markedly from 77.52% to 30.58% for Zn2+, 87.52% to 36.09% for Cu2+, 87.12% to 40.04% for Pb2+, and 80.72% to 31.52% for Fe2+ with an increase in the initial concentration of the metal ions from 25 to 200 mg L−1. This might be due to the rapid saturation of all metal binding active sites of the biosorbent at a certain concentration of the metal ions [52, 53] and an equilibrium state between adsorbate and biosorbent was attained.

3.5. Biosorption Isotherm Studies

The equilibrium adsorption isotherms are essential data source to design, understand, and optimize the biosorption process. The data express the intrinsic properties of the biosorbent and interaction between adsorbate and adsorbent. The data can be used to compare the biosorptive capacities of the biosorbent for different pollutants.

3.5.1. Langmuir Isotherm Model

The model isotherm parameters together with regression coefficient are represented in Table 3. As shown in Figure 2(a), approximately linear relationship () exists in the adsorption isotherms for Kappaphycus sp. The maximum Langmuir monolayer adsorption capacity of the sorbent, (mg g−1), for the experimental metal ions followed an increasing order: Pb2+ (22.27) > Cu2+ (19.46) > Fe2+ (17.09) > Zn2+ (16.78), suggesting that Pb2+ has a preferential uptake compared to the other metals, which can be attributed to its low tendency in forming strong complex [46]. Another reason might be attributed to carboxylate polysaccharides in seaweeds that show preferential binding of cations with large ionic radii [48]. However, the preferential sorption order of the metal ions in the present study can be explained by Pauling’s electronegativity [54]: Pb2+ (2.33) > Cu2+ (1.190) > Fe2+ (1.83) > Zn2+ (1.65). This implies the fact that the higher the ion’s electronegativity the higher the attraction for its electrons, and the attraction becomes stronger to the negative charge of the biomass ligands [55]. Furthermore, the separation factor, , for the metal ions (Table 3) falls within the range of suggesting that adsorption of the experimental ions is favourable at all the concentrations investigated [26]. Hence, the Kappaphycus sp. is a suitable biosorbent for the sorption of the experimental metal ions from aqueous solutions.

ModelParameterMetal ion

Langmuir (mg g−1)22.1719.4916.9216.23
(L mg−1)0.06760.07280.0730.0561

Freundlich (mg g−1) (L mg−13.8364.0003.4913.164

Temkin (L mg−1)1.0371.3971.2281.056
(kJ mol−1)596.703720.350800.999853.158

D-R (mg g−1)15.5814.3313.2212.16
(mol2 kJ−2)2.2842.2244.7735.768
(kJ mol−1)0.4680.4740.3240.294

A comparative study on maximum heavy metal adsorption capacity of different low-cost adsorbents has been given in Table 4. The study shows that the dried biomass of Kappaphycus sp. is more promising than some other low-cost adsorbents for the removal of the metal ions.

Low-cost sorbentAdsorption capacity (mg g−1)Reference

Activated carbon from coconut4.56[41]
Activated carbon from seed hull of the palm tree3.58[41]
Epichlorohydrin-crosslinked chitosan34.1335.4610.21[42]
Hazelnut husk13.056.645[43]
Natural muscovite0.630.618[44]
Modified orange peel73.5315.27[46]
Coconut tree sawdust25.003.8923.81[47]
Sugarcane bagasse21.283.6540.00[47]
Kappaphycus sp.22.2719.4617.0916.78Present study

3.5.2. Freundlich Isotherm Model

The linearized Freundlich isotherm model is shown in Figure 2(b); Freundlich constants and are represented in Table 3. The results suggest that the biosorption of Pb2+ can be moderately described by the Freundlich model (). The magnitude of Freundlich isotherm constant, , suggests that the sorption capacity of the experimental metal ions was in the order of Cu2+ > Pb2+ > Fe2+ > Zn2+. The values of suggest heterogeneity of the biomass surface, and the metal ions are favourably and intensively biosorbed by the dried biomass of Kappaphycus sp. under the experimental conditions.

3.5.3. Temkin Isotherm Model

The Temkin isotherm model indicates the adsorption potentials of the adsorbent for adsorbates. The Temkin isotherm plots (Figure 2(c)) and parameters (Table 3) indicate that the model fits the experimental data well () for describing the metal ions (Pb2+, Cu2+, and Fe2+) adsorption. The lower values of the Temkin adsorption potential, (L mg−1), in the range of 0.769 to 1.455 indicate a lower sorbent-metal ion potential. Furthermore, the lower values (0.585–0.792) of the Temkin constant (kJ mol−1) indicate a weak sorbate-sorbent interaction [45].

3.5.4. D-R Isotherm Model

The D-R isotherm parameters (Table 3) indicate that the D-R model does not fit the experimental data well () for describing the metal ions biosorption, suggesting the involvement of metal sorption mechanisms other than van der Waals force [56]. The mean free energy of biosorption, (0.297−0.489 kJ·mol−1), for the metal ions suggests that the sorption process is physisorption [57] and corroborative to the earlier reports in literature [46, 48]. The positive values of indicate the endothermic nature of the sorption process [46]. Furthermore, the values of (<16 kJ mol−1) suggest that the mechanism of the ion exchange process is film-diffusion controlled [58].

3.5.5. Competitive Adsorption

Competitive adsorption of the metal ions under quaternary system shows the adsorption preference of Pb2+ > Cu2+ > Fe2+ > Zn2+ with the rate of metal removal as 90.39, 90.00, 75.00, and 58.38%, respectively. The results suggest that the potentiality of the adsorbent in the quaternary system remains the same as that in the single metal system, which proves its unique adsorption quality.

3.6. Biosorption Kinetic Studies

The kinetic data are essential to understand the rate and nature of adsorption onto the adsorbents. The data can be used to compare the kinetics of the biosorbent for different pollutants.

3.6.1. Pseudo-First-Order Kinetic Studies

As shown in Table 5 and Figure 3(a), the regression coefficient () of the pseudo-first-order model suggests that the experimental data accurately support the PFO model to describe adsorption kinetics of the metal ions. But the differences between the experimental values, , were higher than the modelled values, . It refers to the fact that both the metal ions and adsorbent were involved in the adsorption process [52]. Therefore, it is suggested that the pseudo-first-order model is not suitable to explain the kinetic sorption of the experimental metal ions onto the dried biomass of Kappaphycus sp. over the range of experimental time and metal ion concentrations. Similar results have been reported for the sorption kinetic of different metal ions onto different adsorbents including seaweeds in the literature [51, 52, 59, 60].

ModelParameterMetal ion

(mg g−1)2.1062.04671.35540.5255

Pseudo-first-order (min−1)0.0350.04350.04540.0329
(mg g−1)0.14560.11620.530.1262

Pseudo-second-order (g mg−1 min−1)0.61330.90470.17710.5968
(mg g−1)2.11592.05591.39990.535
(mg g−1 min−1)2.74573.82410.34710.1708

Elovich (mg g−1 min−1)5.67E + 211.62E + 2262.47353276.307
(g mg−1)27.624328.81847.961827.0124

Intraparticle diffusion (mg g−1 min−0.5)0.00540.00990.03110.0112
(mg g−1)2.04571.96361.05970.4142

3.6.2. Pseudo-Second-Order Kinetic Model

The values of the regression coefficient of the linearized PSO kinetic model as shown in Figure 3(b) were the highest () among the studied kinetic models, and the experimental values matched well with the calculated data (Table 5). Therefore, it can be suggested that the experimental data accurately support the best fit of the PSO model for the adsorption of the metal ions. Hence, chemisorption is the rate-limiting step which involves valence forces through the sharing or exchange of electrons between the metal ions and different functional groups in the sorbent [13, 60].

The pseudo-second-order rate constant, (g mg−1 min−1), was found in the range of 0.1874 to 0.9548, which supports that the metal ions uptake onto the sorbent from aqueous solution was more rapid and favourable. As shown in Figure 3(b), adsorption kinetic of the metal ions on Kappaphycus sp. occurred in two steps: a fast initial uptake rate, (0.17–3.82 mg g−1 min−1), in the first 30 min, where more than 85% of the total metal adsorption occurred, followed by a slower uptake rate leading to the equilibrium state (~120 min). Similar observation was reported in literature [51, 52]. This phenomenon supports that the diffusion is the rate-controlling step in the sorption process [60]. The half-adsorption time (min) defined as the time required for the adsorption to take up half amount of the equilibrium metal ions was found within the range of 0.54 to 4.03 indicating high affinity between the adsorbate and adsorbent molecules [61].

3.6.3. Elovich Model

The values of the regression coefficient ( = 0.79–0.98) of the Elovich kinetic model (Table 5, Figure 3(c)) suggest that kinetic data did not follow the Elovich model. However, the higher values of the Elovich constants, (mg g−1 min−1) and (g mg−1), as shown in Table 5 are suggestive of an increased rate of chemisorption [13].

3.6.4. Intraparticle Diffusion Model

The nonlinear regression data of versus plots as shown in Figure 3(d) for different heavy metal ions suggests multilinearity (two phases in Pb2+ adsorption and three phases in Cu2+, Fe2+, and Zn2+ adsorption). The intraparticle diffusion rate constant () as shown in Table 5 was obtained from the slope of the second linear portions of the plot of versus for the metal ions. Apparently intraparticle diffusion plays a significant role in the adsorption of Pb2+, Cu2+, and Zn2+ () onto the dried biomass of Kappaphycus sp. suggesting the fact that there is a significant relationship between and for the metal ions at the experimental conditions. However, versus plots did not pass through the origin () in any of the cases, suggesting that even though the adsorption process involved intraparticle diffusion, it was not the only rate-controlling step [13, 62], and external mass transfer had also played an important role in the metal ions sorption by the dried biomass of Kappaphycus sp. [13].

3.7. Thermodynamic Studies

The values of the thermodynamic parameters are shown in Table 6. The linearized Van’t Hoff plots of ln () versus are represented in Figure 4. The negative values of indicate that the thermodynamic process was spontaneous and feasible for all the tested metal ions [35]. Moreover, the increase in negative values with an increase in temperature shows an increased feasibility of adsorption at higher temperature, which is corroborative to the earlier reports [48, 52, 59].

Metal ion (kJ mol−1)
(kJ mol−1)(kJ mol−1 K−1)298° K303° K313° K323° K


The positive values of enthalpy change () suggest endothermic nature of the metal adsorption process [13, 51, 59, 63, 64]. In addition, the extent of enthalpy value gives indicative information on the type of biosorption, which can be either physical or chemical. The enthalpy change () in the range of 2.1–20.9, 20.9–80.0, and 80.0–418.4 kJ mol−1 is indicative of physisorption, physisorption together with chemisorptions, and chemisorptions, respectively [48]. Based on the values of , it can be presumed that the biosorption process took place physically for all the tested metal ions. This was also supported by D-R isotherm results with the (< 8 kJ mol−1) values of the metal ions (Table 3). Further, positive values of entropy change () are suggestive of increased randomness at the solid-solution interface during the biosorption process of the metal ions on the active sites of the biosorbent [59].

3.8. Error Function

In the real-world, data samples from each experiment in a series of experiments differ due to measurement error affecting data precision. In order to ensure accurate measurement results, statistical error function is the measure to compensate data errors [65]. Hence, the isotherm and kinetic data were further analyzed using nine error functions in order to test the fitness of the models. Lower value of SSE, ARE, HYBRID, EABS, MPSD, , , and RMSE and higher value of indicate the best fit of the model.

The correlation of regression () for the adsorption isotherm models (Table 3) suggests that Pb2+, Fe2+, and Zn2+ follow the Langmuir model while Cu2+ follows the Temkin model accurately. The error functions of the isotherm data (Table 7) suggest that the Temkin model provides the best fit to the experimental data. Again, the correlation of regression () for the kinetic models (Table 5) shows that PSO is the best fit model. But the error functions of the kinetic data (Table 8) suggest that the best fit of the kinetic models is intraparticle diffusion. It is, therefore, strongly suggested that the regression coefficient () is not an appropriate method for comparing the best fitting of the isotherm and kinetic models; rather some forms of error analysis could be a better criterion for avoiding data errors.

Metal ionIsotherm modelError function





Metal ionKinetic modelError function

Elovich9.30E − 050.160200.224280.023120.300640.984564.52E − 050.191280.00431
Intraparticle diffusion2.39E − 060.038990.116980.002450.106720.984371.14E − 060.052100.00155

Pseudo-second-order0.000150.160290.224400.022320.387120.978697.51E − 050.251510.00544
Intraparticle diffusion3.29E − 060.050870.152620.003100.127200.987181.62E − 060.063110.00181

Intraparticle diffusion0.000310.667251.334490.034361.089370.964720.000240.784270.01236

Intraparticle diffusion3.87E − 060.181190.362390.003620.195970.996457.68E − 060.225670.00139

3.9. FTIR Spectral Analysis

The FTIR spectra of Kappaphycus sp. (Figure 5) consist of a number of absorption peaks which indicate complex nature of the biomass. The strong broad peak observed at 3358.1 cm−1 in the raw biomass corresponds to O–H group from cellulose and N–H groups from proteins in the seaweed [66]. In the spectra other dominant peaks were observed at wavenumbers (cm−1) 2917.1, 1636.5, 1375.0, 1220.1, 1155.1, 1035.2, 924.0, and 842.9 which are characterized to the asymmetric C–H stretching vibrations of the aliphatic groups [67], C=O stretching vibration of carboxylate groups [48, 60], asymmetric stretching of – bonds in sulfonic acid [48], C=O stretching vibration of carboxylate groups [48, 60], symmetric stretching of – bonds in sulfonic acid [48], C–O stretching vibration of carboxyl groups [60], S–O stretching [60], and S=O stretching bands of sulfonate groups [48], respectively.

After biosorption of Pb2+ the peaks were shifted to 3355.9, 2917.8, 1638.8, 1370.4, 1222.0, 1154.1, 1033.5, 924.8, and 844.5 cm−1, respectively. After biosorption of Cu2+ the peaks were changed to 3324.7, 2919.9, 1638.3, 1370.5, 1216.0, 1153.3, 1032.7, 925.3, and 845.9 cm−1, respectively. The peaks after Fe2+ biosorption were changed to 3351.9, 2916.8, 1637.2, 1369.2, 1223.9, 1154.5, 1032.1, 925.2, and 844.2 cm−1, respectively. After biosorption of Zn2+ the peaks were shifted to 3328.4, 2918.1, 1636.9, 1370.3, 1221.1, 1154.2, 1030.9, 924.7, and 845.3 cm−1, respectively. In the quaternary system, the peaks after biosorption were shifted to 3348.3, 2918.0, 1636.1, 1355.1, 1224.0, 1154.9, 1033.1, 924.9, and 845.1 cm−1, respectively. The significant change in the intensity of the peaks shows interaction between the metal ions and the functional groups. Because intensity depends on change in dipole moment and total number of functional groups present on biosorbent surface. Therefore, it can be concluded that the carboxylic, sulfonic acid, and sulfonate groups of Kappaphycus sp. dried biomass are involved in the biosorption of the metal ions.

4. Conclusion

In the present study, we examined adsorption of four heavy metal ions such as Pb2+, Cu2+, Fe2+, and Zn2+ onto the dried biomass of the red seaweed Kappaphycus sp. from Malaysia. The adsorption isotherm data for the metal ions fitted well with the Temkin model followed. Kinetic data for all the metal ions can be best described by the intraparticle diffusion model. Adsorption process was feasible, spontaneous, and endothermic in nature. We strongly suggest that analysis of error functions is a better criterion for validating isotherm and kinetic models in order to evaluate adsorptive behaviour of a typical adsorbent using linear method.

Heavy metal adsorption process onto the dried biomass of Kappaphycus sp. was the complex one involving more than one mechanism. Both homogeneous and heterogeneous active sites were found to exist in the dried biomass. The FTIR study revealed the presence of carboxylic, sulfonic acid, and sulfonate groups in the cell wall matrix of the biomass that was involved in the adsorption of the metal ions. The dried biomass of Kappaphycus sp. may be used as a low-cost biosorbent for removal of heavy metal ions from aqueous solutions. Further study is warranted to evaluate the potentiality of the biosorbent for heavy metal removal from the real environment.

Conflict of Interests

The authors hereby declare no conflict of interests.


The authors greatly acknowledge the financial support of the AIMST University (AURGC/1/FAS/2013) for the research.


  1. R. K. Gautam, S. K. Sharma, S. Mahiya, and M. C. Chattopadhyaya, “Contamination of heavy metals in aquatic media: transport, toxicity and technologies for remediation,” in Heavy Metals In Water: Presence, Removal and Safety, S. K. Sharma, Ed., pp. 1–24, The Royal Society of Chemistry, Cambridge, UK, 2015. View at: Google Scholar
  2. S.-L. Wang, X.-R. Xu, Y.-X. Sun, J.-L. Liu, and H.-B. Li, “Heavy metal pollution in coastal areas of South China: a review,” Marine Pollution Bulletin, vol. 76, no. 1-2, pp. 7–15, 2013. View at: Publisher Site | Google Scholar
  3. A. Sarkar, J. Bhagat, and S. Sarker, “Evaluation of impairment of DNA in marine gastropod, Morula granulata as a biomarker of marine pollution,” Ecotoxicology and Environmental Safety, vol. 106, pp. 253–261, 2014. View at: Publisher Site | Google Scholar
  4. P. Chowdhury, A. Elkamel, and A. K. Ray, “Photocatalytic processes for the removal of toxic metal ions,” in Heavy Metals in Water: Presence, Removal and Safety, S. K. Sharma, Ed., pp. 25–43, The Royal Society of Chemistry, Cambridge, UK, 2015. View at: Google Scholar
  5. M. Zhang, B. Gao, J. Jin et al., “Use of nanotechnology against heavy metals present in water,” in Heavy Metals in Water: Presence, Removal and Safety, S. K. Sharma, Ed., pp. 177–192, The Royal Society of Chemistry, Cambridge, UK, 2015. View at: Google Scholar
  6. S. Majumder, S. Gupta, and S. Raghuvanshi, “Removal of dissolved metals by bioremediation,” in Heavy Metals in Water: Presence, Removal and Safety, S. K. Sharma, Ed., pp. 44–56, The Royal Society of Chemistry, Cambridge, UK, 2015. View at: Google Scholar
  7. J. He and J. P. Chen, “A comprehensive review on biosorption of heavy metals by algal biomass: materials, performances, chemistry, and modeling simulation tools,” Bioresource Technology, vol. 160, pp. 67–78, 2014. View at: Publisher Site | Google Scholar
  8. T. Macek and M. Mackova, “Potential of biosorption technology,” in Microbial Biosorption of Metals, P. Kotrba, M. Mackova, and T. Macek, Eds., pp. 7–17, Springer Netherlands, 2011. View at: Google Scholar
  9. M. Grassi, G. Kaykioglu, V. Belgiorno, and G. Lofrano, “Removal of emerging contaminants from water and wastewater by adsorption process,” in Emerging Compounds Removal from Wastewater, G. Lofrano, Ed., pp. 15–37, Springer Netherlands, 2012. View at: Google Scholar
  10. M. Bilal, J. A. Shah, T. Ashfaq et al., “Waste biomass adsorbents for copper removal from industrial wastewater-a review,” Journal of Hazardous Materials, vol. 263, part 2, pp. 322–333, 2013. View at: Publisher Site | Google Scholar
  11. I. Michalak, K. Chojnacka, and A. Witek-Krowiak, “State of the art for the biosorption process—a review,” Applied Biochemistry and Biotechnology, vol. 170, no. 6, pp. 1389–1416, 2013. View at: Publisher Site | Google Scholar
  12. M. Fomina and G. M. Gadd, “Biosorption: current perspectives on concept, definition and application,” Bioresource Technology, vol. 160, pp. 3–14, 2014. View at: Publisher Site | Google Scholar
  13. M. Rathod, K. Mody, and S. Basha, “Efficient removal of phosphate from aqueous solutions by red seaweed, Kappaphycus alverezii,” Journal of Cleaner Production, vol. 84, no. 1, pp. 484–493, 2014. View at: Publisher Site | Google Scholar
  14. V. Webber, S. M. de Carvalho, and P. L. M. Barreto, “Molecular and rheological characterization of carrageenan solutions extracted from Kappaphycus alvarezii,” Carbohydrate Polymers, vol. 90, no. 4, pp. 1744–1749, 2012. View at: Publisher Site | Google Scholar
  15. FAO, Global Aquaculture Production 1950–2012, 2015,
  16. O. L. Kang, N. Ramli, M. Said, M. Ahmad, S. M. Yasir, and A. Ariff, “Kappaphycus alvarezii waste biomass: a potential biosorbent for chromium ions removal,” Journal of Environmental Sciences, vol. 23, no. 6, pp. 918–922, 2011. View at: Publisher Site | Google Scholar
  17. K. O. Lee, N. Ramli, M. Said, M. Ahmad, S. M. Yasir, and A. Ariff, “Copper (II) and nickel (II) sorption onto seaweed (Kappaphycus alvarezii) waste biomass: equilibrium and mechanism studies,” Middle-East Journal of Scientific Research, vol. 9, no. 1, pp. 84–89, 2011. View at: Google Scholar
  18. D. Park, Y.-S. Yun, and J. M. Park, “The past, present, and future trends of biosorption,” Biotechnology and Bioprocess Engineering, vol. 15, no. 1, pp. 86–102, 2010. View at: Publisher Site | Google Scholar
  19. S. Basha, S. Jaiswar, and B. Jha, “On the biosorption, by brown seaweed, Lobophora variegata, of Ni(II) from aqueous solutions: equilibrium and thermodynamic studies,” Biodegradation, vol. 21, no. 5, pp. 661–680, 2010. View at: Publisher Site | Google Scholar
  20. H. Qiu, L. Lv, B.-C. Pan, Q.-J. Zhang, W.-M. Zhang, and Q.-X. Zhang, “Critical review in adsorption kinetic models,” Journal of Zhejiang University SCIENCE A, vol. 10, no. 5, pp. 716–724, 2009. View at: Publisher Site | Google Scholar
  21. S. Sen Gupta and K. G. Bhattacharyya, “Kinetics of adsorption of metal ions on inorganic materials: a review,” Advances in Colloid and Interface Science, vol. 162, no. 1-2, pp. 39–58, 2011. View at: Publisher Site | Google Scholar
  22. W. Plazinski and A. Plazinska, “Equilibrium and kinetic modeling of adsorption at solid/solution interface,” in Application of Adsorbents for Water Pollution Control, A. Bhatnagar, Ed., pp. 32–80, Bentham Science Publishers, AG Bussum, The Netherlands, 2012. View at: Google Scholar
  23. K. E. Carpenter and V. H. E. Niem, FAO Species Identification Guide for Fishery Purposes. The Living Marine Resources of the Western Central Pacific. Volume 1. Seaweeds, Corals, Bivalves and Gastropods, FAO, Rome, Italy, 1998.
  24. A. R. Sherwood, A. Kurihara, K. Y. Conklin, T. Sauvage, and G. G. Presting, “The Hawaiian Rhodophyta Biodiversity Survey (2006–2010): a summary of principal findings,” BMC Plant Biology, vol. 10, article 258, 2010. View at: Publisher Site | Google Scholar
  25. I. Langmuir, “The adsorption of gases on plane surfaces of glass, mica and platinum,” The Journal of the American Chemical Society, vol. 40, no. 9, pp. 1361–1403, 1918. View at: Publisher Site | Google Scholar
  26. T. W. Weber and R. K. Chakravorti, “Pore and solid diffusion models for fixed-bed adsorbers,” AIChE Journal, vol. 20, no. 2, pp. 228–238, 1974. View at: Publisher Site | Google Scholar
  27. H. M. F. Freundlich, “Über die adsorption in lösungen,” Zeitschrift für Physikalische Chemie (Leipzig), vol. 57A, pp. 385–470, 1906. View at: Google Scholar
  28. S. Z. Roginsky and Y. B. Zeldovich, “Die katalische oxidation von kohlenmonoxyd auf mangandioxyd,” Acta Physiochimica URSS, vol. 1, no. 3-4, pp. 554–594, 1934. View at: Google Scholar
  29. M. M. Dubinin and L. V. Radushkevich, “Equation of the characteristic curve of the activated charcoal,” Chemisches Zentralblatt, vol. 1, no. 1, pp. 875–890, 1947. View at: Google Scholar
  30. J. P. Hobson, “Physical adsorption isotherms extending from ultrahigh vacuum to vapor pressure,” The Journal of Physical Chemistry, vol. 73, no. 8, pp. 2720–2727, 1969. View at: Publisher Site | Google Scholar
  31. S. Y. Lagergren, “Zur theorie der sogenannten adsorption gelöster Stoffe,” Kongliga Svenska Vetenskaps-Akademiens Handlingar, vol. 24, no. 4, pp. 1–39, 1898. View at: Google Scholar
  32. Y.-S. Ho, Adsorption of heavy metals from waste streams by peat [Ph.D. thesis], University of Birmingham, Birmingham, UK, 1995.
  33. S. H. Chien and W. R. Clayton, “Application of elovich equation to the kinetics of phosphate release and sorption in soils1,” Soil Science Society of America Journal, vol. 44, no. 2, pp. 265–268, 1980. View at: Publisher Site | Google Scholar
  34. W. J. Weber and J. C. Morris, “Kinetics of adsorption on carbon from solution,” Journal of the Sanitary Engineering Division, vol. 89, no. 2, pp. 31–60, 1963. View at: Google Scholar
  35. J. W. Gibbs, “A method of geometrical representation of the thermodynamic properties of substances by means of surfaces,” Transactions of the Connecticut Academy of Arts and Sciences, vol. 2, pp. 382–404, 1873. View at: Google Scholar
  36. B. Subramanyam and A. Das, “Linearised and non-linearised isotherm models optimization analysis by error functions and statistical means,” Journal of Environmental Health Science & Engineering, vol. 12, no. 1, article 92, 2014. View at: Publisher Site | Google Scholar
  37. L. Wang, J. Zhang, R. Zhao, Y. Li, C. Li, and C. Zhang, “Adsorption of Pb(II) on activated carbon prepared from Polygonum orientale Linn.: kinetics, isotherms, pH, and ionic strength studies,” Bioresource Technology, vol. 101, no. 15, pp. 5808–5814, 2010. View at: Publisher Site | Google Scholar
  38. J. He, S. Hong, L. Zhang, F. Gan, and Y.-S. Ho, “Equilibrium and thermodynamic parameters of adsorption of methylene blue onto rectorite,” Fresenius Environmental Bulletin, vol. 19, no. 11, pp. 2651–2656, 2010. View at: Google Scholar
  39. Y.-S. Ho and A. E. Ofomaja, “Pseudo-second-order model for lead ion sorption from aqueous solutions onto palm kernel fiber,” Journal of Hazardous Materials, vol. 129, no. 1–3, pp. 137–142, 2006. View at: Publisher Site | Google Scholar
  40. K. Vijayaraghavan, T. V. N. Padmesh, K. Palanivelu, and M. Velan, “Biosorption of nickel(II) ions onto Sargassum wightii: application of two-parameter and three-parameter isotherm models,” Journal of Hazardous Materials, vol. 133, no. 1–3, pp. 304–308, 2006. View at: Publisher Site | Google Scholar
  41. S. Gueu, B. Yao, K. Adouby, and G. Ado, “Kinetics and thermodynamics study of lead adsorption on to activated carbons from coconut and seed hull of the palm tree,” International Journal of Environmental Science & Technology, vol. 4, no. 1, pp. 11–17, 2007. View at: Publisher Site | Google Scholar
  42. A.-H. Chen, S.-C. Liu, C.-Y. Chen, and C.-Y. Chen, “Comparative adsorption of Cu(II), Zn(II), and Pb(II) ions in aqueous solution on the crosslinked chitosan with epichlorohydrin,” Journal of Hazardous Materials, vol. 154, no. 1–3, pp. 184–191, 2008. View at: Publisher Site | Google Scholar
  43. M. Imamoglu and O. Tekir, “Removal of copper (II) and lead (II) ions from aqueous solutions by adsorption on activated carbon from a new precursor hazelnut husks,” Desalination, vol. 228, no. 1–3, pp. 108–113, 2008. View at: Publisher Site | Google Scholar
  44. J.-S. Yang, J. Y. Lee, Y.-T. Park, K. Baek, and J. Choi, “Adsorption of As(III), As(V), Cd(II), Cu(II), and Pb(II) from aqueous solutions by natural muscovite,” Separation Science and Technology, vol. 45, no. 6, pp. 814–823, 2010. View at: Publisher Site | Google Scholar
  45. S. Shahmohammadi-Kalalagh, H. Babazadeh, A. H. Nazemi, and M. Manshouri, “Isotherm and kinetic studies on adsorption of Pb, Zn and Cu by kaolinite,” Caspian Journal of Environmental Sciences, vol. 9, no. 2, pp. 243–255, 2011. View at: Google Scholar
  46. M. R. Lasheen, N. S. Ammar, and H. S. Ibrahim, “Adsorption/desorption of Cd(II), Cu(II) and Pb(II) using chemically modified orange peel: equilibrium and kinetic studies,” Solid State Sciences, vol. 14, no. 2, pp. 202–210, 2012. View at: Publisher Site | Google Scholar
  47. W. Pranata Putra, A. Kamari, S. Najiah Mohd Yusoff et al., “Biosorption of Cu(II), Pb(II) and Zn(II) ions from aqueous solutions using selected waste materials: Adsorption and characterisation studies,” Journal of Encapsulation and Adsorption Sciences, vol. 4, no. 1, pp. 25–35, 2014. View at: Publisher Site | Google Scholar
  48. S. Yalçin, “The mechanism of heavy metal biosorption on green marine macroalga Enteromorpha linza,” CLEAN—Soil, Air, Water, vol. 42, no. 3, pp. 251–259, 2014. View at: Publisher Site | Google Scholar
  49. F. B. de Souza, S. M. A. G. U. de Souza, A. A. U. de Souza et al., “Modeling of trivalent chromium speciation in binding sites of marine macroalgae Sargassum cymosum,” Clean Technologies and Environmental Policy, vol. 15, no. 6, pp. 987–997, 2013. View at: Publisher Site | Google Scholar
  50. W. M. Ibrahim, “Biosorption of heavy metal ions from aqueous solution by red macroalgae,” Journal of Hazardous Materials, vol. 192, no. 3, pp. 1827–1835, 2011. View at: Publisher Site | Google Scholar
  51. M. Arshadi, M. J. Amiri, and S. Mousavi, “Kinetic, equilibrium and thermodynamic investigations of Ni(II), Cd(II), Cu(II) and Co(II) adsorption on barley straw ash,” Water Resources and Industry, vol. 6, pp. 1–17, 2014. View at: Publisher Site | Google Scholar
  52. M. D. Meitei and M. N. V. Prasad, “Adsorption of Cu (II), Mn (II) and Zn (II) by Spirodela polyrhiza (L.) Schleiden: equilibrium, kinetic and thermodynamic studies,” Ecological Engineering, vol. 71, pp. 308–317, 2014. View at: Publisher Site | Google Scholar
  53. A. A. Al-Homaidan, H. J. Al-Houri, A. A. Al-Hazzani, G. Elgaaly, and N. M. S. Moubayed, “Biosorption of copper ions from aqueous solutions by Spirulina platensis biomass,” Arabian Journal of Chemistry, vol. 7, no. 1, pp. 57–62, 2014. View at: Publisher Site | Google Scholar
  54. L. Pauling, The Nature of the Chemical Bond: An Introduction to Modern Structural Chemistry, Cornell University Press, Ithaca, NY, USA, 3rd edition, 1960.
  55. J. M. Lezcano, F. González, A. Ballester, M. L. Blázquez, J. A. Muñoz, and C. García-Balboa, “Biosorption of Cd(II), Cu(II), Ni(II), Pb(II) and Zn(II) using different residual biomass,” Chemistry and Ecology, vol. 26, no. 1, pp. 1–17, 2010. View at: Publisher Site | Google Scholar
  56. S. Basha, Z. V. P. Murthy, and B. Jha, “Isotherm modeling for biosorption of Cu(ll) and Ni(ll) from wastewater onto brown seaweed, cystoseira Indica,” AIChE Journal, vol. 54, no. 12, pp. 3291–3302, 2008. View at: Publisher Site | Google Scholar
  57. M. M. Dubinin, “The potential theory of adsorption of gases and vapors for adsorbents with energetically nonuniform surfaces,” Chemical Reviews, vol. 60, no. 2, pp. 235–241, 1960. View at: Publisher Site | Google Scholar
  58. G. E. Boyd and B. A. Soldano, “Self-diffusion of cations in and through sulfonated polystyrene cation-exchange polymers,” Journal of the American Chemical Society, vol. 75, no. 24, pp. 6091–6099, 1953. View at: Publisher Site | Google Scholar
  59. C. L. Massocatto, E. C. Paschoal, N. Buzinaro et al., “Preparation and evaluation of kinetics and thermodynamics studies of lead adsorption onto chemically modified banana peels,” Desalination and Water Treatment, vol. 51, no. 28–30, pp. 5682–5691, 2013. View at: Publisher Site | Google Scholar
  60. J. Plaza Cazón, M. Viera, E. Donati, and E. Guibal, “Zinc and cadmium removal by biosorption on Undaria pinnatifida in batch and continuous processes,” Journal of Environmental Management, vol. 129, pp. 423–434, 2013. View at: Publisher Site | Google Scholar
  61. Q. Shi, A. Li, Z. Zhu, and B. Liu, “Adsorption of naphthalene onto a high-surface-area carbon from waste ion exchange resin,” Journal of Environmental Sciences, vol. 25, no. 1, pp. 188–194, 2013. View at: Publisher Site | Google Scholar
  62. D. D. Maksin, S. O. Kljajević, M. B. Dolić et al., “Kinetic modeling of heavy metal sorption by vinyl pyridine based copolymer,” Hemijska Industrija, vol. 66, no. 6, pp. 795–804, 2012. View at: Publisher Site | Google Scholar
  63. M. I. Din, M. L. Mirza, S. Ata, M. Athar, and I. U. Mohsin, “Thermodynamics of biosorption for removal of Co(II) ions by an efficient and ecofriendly biosorbent (Saccharum bengalense): kinetics and isotherm modeling,” Journal of Chemistry, vol. 2013, Article ID 528542, 11 pages, 2013. View at: Publisher Site | Google Scholar
  64. N. Rajamohan, M. Rajasimman, R. Rajeshkannan, and V. Saravanan, “Equilibrium, kinetic and thermodynamic studies on the removal of Aluminum by modified Eucalyptus camaldulensis barks,” Alexandria Engineering Journal, vol. 53, no. 2, pp. 409–415, 2014. View at: Publisher Site | Google Scholar
  65. W. R. Leo, Techniques for Nuclear and Particle Physics Experiments: A How-to Approach, Springer, Berlin, Germany, 2nd edition, 1994.
  66. A. Bhatnagar, V. J. P. Vilar, C. Ferreira, C. M. S. Botelho, and R. A. R. Boaventura, “Optimization of nickel biosorption by chemically modified brown macroalgae (Pelvetia canaliculata),” Chemical Engineering Journal, vol. 193-194, pp. 256–266, 2012. View at: Publisher Site | Google Scholar
  67. R. B. Nessim, A. R. Bassiouny, H. R. Zaki, M. N. Moawad, and K. M. Kandeel, “Biosorption of lead and cadmium using marine algae,” Chemistry and Ecology, vol. 27, no. 6, pp. 579–594, 2011. View at: Publisher Site | Google Scholar

Copyright © 2015 Md. Sayedur Rahman and Kathiresan V. Sathasivam. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.