Original Research
Adsorption mechanism of malachite green onto activated phosphate rock: a kinetics and theoretical study

Nouar Sofiane Labidi* , Nour Elhouda Kacemi


Department of Materials Sciences, Institute of Science and Technologies. University Centre of Tamanrasset, Algeria


The author(s) declared that no grants were involved in supporting this work.



Adsorption kinetics of malachite green onto the Algerian activated phosphate rock was studied for better removal of the dye from wastewater. The prepared sorbent displayed à good surface area of 42.2 m²/g. The adsorption process appeared to be of physisorption nature and it took less than 60 min to get equilibrium whereas the kinetics indicated that the adsorption is likely a second order reaction, which is further proved with the high R2 value. The intraparticle diffusion model confirms an adsorption mechanism limited on two steps, i.e., (1) surface adsorption, and (2) pore diffusion with a diffusion parameter of Di=10-18 cm2 /s. Besides, semi-empirical theoretical calculations provide a new insight into adsorption mechanism as a principle of hydrogen bonding and ionic interaction.

Corresponding author: Nouar Sofiane Labidi *,

How to cite: Labidi, N.S, and Kacemi, N.E, 2016: Adsorption mechanism of malachite green onto activated phosphate rock: a kinetics and theoretical study. Bulletin of Environmental Studies 1(3): 69-74.


Copyright: © 2017 MN Publishers. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

Competing interests: The authors declare no conflict of interest.

Edited by: Muhammad Arslan (UFZ, Germany)

Reviewed by: Naeem Shahid (RWTH Aachen, Germany) & Mujtaba Baqar (GCU Lahore, Pakistan)

Received: 24/01/2016

Accepted: 22/07/2016

Published Online: 22/07/2016


Today, synthetic dyes are found in all spheres of our daily lives as more than 100,000 of commercial dyes have been produced as a factor of 7x105 to 106 tons per year (McMullan et al , 2001). Among them, 10-15% of are estimated to be discharged as effluents by textile industry (Ncibi et al., 2007). Even in smaller amounts, xenobiotic nature of dyes makes them resistant to deterioration and hence are persistent in flowing water resulting into inhibition of aquatic biota (Patil et al., 2011). The compounds are therefore required to be removed from the aquatic system. Several methods have been employed for remediation of dye-containing effluents including biological, physicochemical, and chemical techniques (Forgacs et al., 2004).

For the removal of dyes from wastewater, adsorption has been proven to be one of the effective procedure due to its high efficiency, adsorbent recyclability and less harmful effects on environment as compared to others conventional remediation techniques (Ertaş et al., 2010). Since commercially available synthetic adsorbent are expensive and are difficult to regenerate (Rubin et al.,2005), research is focusing on use of natural and low cost adsorbents especially phosphate containing compounds which possess great potential to remove dyes from solutions (Keleş et al.,2010; Elouear et al.,2008; Young et al., 2005; Aklil et al., 2004; Baraka et al.,2009; Alzaydien, 2009; Achkoun et al., 2012).

In order to develop an accurate understanding of the adsorption process, theoretical modeling of the physicochemical properties resulting from the interaction between adsorbent and adsorbate in terms of various computations quantum chemistry approaches is therefore, an important issue. Recently, quantum chemical methods have profoundly changed the science of designing, developing and interpreting structure–property relationships of different organic and inorganic materials (Pang and Xu, 2012; Wick et al., 2010; Vysotsky et al., 2015; Chernyshova et al., 2011). For experimentalists, exploitation of deep insights provided by the computational quantum chemistry may present an effective guide for the adsorption studies in terms of comprehension and interpretation. 

The aim of the current study is to study the mechanistic nature of malachite green adsorption onto low-cost Algerian activated phosphate rock. For this purpose, a combination of theoretical and kinetic model is used, i.e., general adsorption, adsorption diffusion, and semi-empirical modeling.

Materials and Methods


The phosphate used in this study was obtained from the region of Djebel onk (FERPHOS Mining Company, Algeria); whose basin is located about 100 km south of Tebessa city and 20 km from the Algerian-Tunisian border. The mountainous massive of Djebel Onk forms a set of 20 km length, which rises to 1198 m altitude in Djebel Tarfaya. The geographical coordinates are recorded as 34.4300° N, 8.0000° E (Mezghache et al., 2002).

As phosphate rock was received in the form of brown-coloured powder, prior purification was undergone to remove all impurities associated with the phosphate material. Particles fraction, having a particle size less than 80 μm, was abundantly washed with distilled water and then dried in an oven overnight at 105 °C. Subsequently, the fraction of the phosphate obtained was subjected to chemical activation (acid/base) in order to obtain an activated phosphate rock-support according to the protocol described by Elouearet al., (2008). Spectrometric and gravimetric methods were carried out for the chemical analysis and the specific surface area was measured by Sear’s method (Sears., 1956).

Sorbate Dye

Malachite green dye was purchased from Sigma-Aldrich. A stock solution of dye was prepared by dissolving 1 gram dye in 1000 ml distilled water. Afterwards, the solution diluted using distilled water as per the required concentration range of 1-100 mg/L.

Kinetics study

The adsorption kinetics of malachite green was performed using 100 mg/L dye concentration and an activated phosphate rock mass of 1 gram. The agitation speed was kept constant at 250 rpm for 85 min at 30 °C. At programmed intervals of time, malachite green concentrations were measured using the UV-visible spectrometric method at a wavelength of 674 nm. Blanks containing distilled water were used as negative controls.

Theoretical calculation (modeling approach)

The computations were made with a P4 Duo Proc X64 (8G0 RAM) microcomputer in Windows 8.0 environment. All calculations were made using HyperChem v7 software (http:// www.hyperchem.com). The geometry of malachite green was optimized at the MM+ level of theory (r.m.s = 0.0001 Kcal/Å). Mulliken atomic charges and electrostatic potential were assessed by the semi-empirical parametric method 3 (PM3) being derived from the Hartree-Fock theory (Stewart., 2004).

We studied possible interaction sites through computation of electronic charge distribution and electrostatic potential. The interactions between phosphate groups (PO4-3 ) of activated phosphate rock and malachite green dye is considered as background knowledge due to (a) acid–base interaction due to hydrogen bonding; (b) adsorption by dispersion forces as a principle of electrostatic interactions; and (c) ion pairing by adsorption onto oppositely unoccupied charged sites.

Results and Discussion

Phosphate rock characterization

The chemical analysis revealed that the main components in phosphate rock are CaO, P2O5, and CO2 followed by flour and quartz (Table 1). The phosphate rock powders showed different chemical composition before and after the activation due to acid/base treatment. A reduction in the mass fraction of CO2, CaO and F2 after chemical treatment is observed, which was further decreased from 9.83%, 42.8% and 2.8 % (w/w) to 3.83%, 42.72% and 2.04 % (w/ w), respectively. However, the percentage of P2O5 content increased from 23.05% to 23.89%.


The surface area of activated phosphate rock is observed to be 42.20 m²/g, which is greater than the natural phosphate rock whose initial surface area was of 27.67 m²/g. This value indicates that the phosphate has low porosity.

The pH measurements are an important factor in adsorption study, as it indicates surface charge properties of phosphate particles and in turn can suggest the ability of the compound to get adsorb, penetrate and adhere onto the surface of the adsorbent. In distilled water, phosphate particles had an isoelectric point for pH ranging between 4-6, while the surface charge is appeared to be negative above this pH. The activated phosphate rock used in this study had a pH of 8.58 that confirms the negative charge of activated phosphate rock over the surface area (Figure 1).


Adsorption equilibrium – effect of contact time

The adsorption experiments were performed using two systems with natural and activated phosphate rock. The adsorption capacity of malachite green in both systems was measured using 1 g of the prepared phosphate combined with an initial solution of 20 mg/L at adjusted pH of 8. The solutions were continuously agitated at 250 rpm and a temperature of 30 °C.

Figure 2 illustrates that the adsorption rate was fast in the beginning of the process which became slower during the stirring time, and reached to the equilibrium state in less than an hour. It is observed that the adsorption capacity of the activated phosphate rock was larger than that of natural phosphate rock. In addition, as adsorption rate was faster with activated phosphate rock compared to the phosphate rock; the difference could be due to the nature of the interactions of each phosphate surface with the malachite green dye and to the changes in phosphate surface area (Elouear et al., 2008).

Adsorption kinetics


The adsorption reaction and adsorption diffusion models are well-established methods to describe the adsorption kinetics. Hereby, we considered the kinetic laws of the pseudo-first-order (Ho and McKay, 1999), and the pseudo-second-order (Ho et al., 2000) for intraparticle diffusions as described earlier (Manlin et al., 2014).

Pseudo first-order kinetic model: The pseudo first-order equation, proposed by Lagergren (Ho and McKay, 1999), can be expressed as follows

log (qe - qt) = log(qe) - (K1 / 2.303) x t ........... (1)

where qe and qt are the quantities of dye adsorbed at “equilibrium” and at time “t”, respectively (mg/g) while K1 is the adsorption rate constant (min-1 ).

The straight line obtained by plotting log (qe-qt) versus t gives log qe as the intercept and K1/2.303 as the slope. Thus, the amount of solute adsorbed at equilibrium (qe) and the first order kinetic constant (K1) can be obtained from the intercept and the slope.

The values of the kinetic model parameters and the correlation coefficient are presented in Table 2 & Figure 3 showed the pseudofirst order kinetics modelling of malachite green dye adsorption onto activated phosphate rock. The equilibrium rate constant of pseudo-first order sorption (K1) was recorded to be of 2.543x10-2 min-1 . The correlation coefficients R2 was below than 0.930 and the calculated amount of solute adsorbed qe did not appear to be good in agreement with the experimental data. Therefore, it is confirmed that the adsorption of malachite green onto activated phosphate rock is not a first order kinetics.


Pseudo-second-order kinetic model: The pseudo-second-order kinetic equation is expressed as follows (Ho et al., 2000);

t / qt = 1 / K2qe2 + (1/qe) t ................ (2)

where K2 (g/mg·min) is the second order kinetic constant for the adsorption reaction.

The results in Figure 4 showed a linear relationship with a correlation coefficient closer to the unity (R2 =0.997), hence indicating an adequate model to correctly describe the adsorption kinetics. Hereby, calculated qe values found to agree the experimental data set and therefore it can be concluded that the adsorption of malachite green dye onto activated phosphate rock was likely a second order reaction. Similar results were found in the adsorption of malachite green onto different adsorbents such as a chitosan bead, degreased coffee bean, modified rice husk and maize cob (Bekçi et al., 2008; Baek et al., 2010; Zou et al., 2011; Sonawane et al., 2009).


Intraparticle diffusion model: The intraparticle diffusion model is presented by the following equation (Manlin et al., 2014);

qt = Ktt1/2................... (3)

where Kt is the intraparticle diffusion rate constant (g/mg.min1/2), which is derived from the slope of the linear part of the equation for this model. The qt is the quantity of the dye adsorbed at “equilibrium” and at time “t”.

The values of the diffusion parameter (Di) are deduced from the following formula;

Di = 0.03.ro / t1/2 .............. (4)

where t1/2 is a half-time reaction (s) (adsorption equilibrium reaction), ro is the particle diameter of the adsorbent in (cm), and D is intraparticle diffusion parameter (cm2 /s).

The introduction of the equation 3 to the experimental data indicated that the plot of qt versus t1/2 exhibit two separate regions with two straight lines, implying the existence of several steps that may have controlled the adsorption reaction process. The intraparticle diffusion constants observed to have decreasing values, i.e., Kt1=1.42 (g/mg.min1/2) and Kt2=0.023 (g/mg.min1/2), successively; which are directly related to the nature of the internal porosity of activated phosphates (Zou et al., 2011).

In Figure 5, two plateau are observed, indicating that the diffusion of malachite green onto the activated phosphate rock surface is carried out in two stages. In addition, deviations of the straight lines from the origin and the calculated values of the diffusion parameter are estimated to be of Di= 10-18 cm2 /s, hence, suggesting that the intraparticle diffusion was not the only rate limiting mechanism in the adsorption mechanism.


The behavior may be due to the dimensionality of malachite green dye since the malachite green is a voluminous molecule with a calculated polarizability value of 42Å3 . It can be deduced from here that the dye diffused slowly in the macro-pores of phosphate particles during the adsorption reaction. The phenomenon, therefore, confirms that the mechanism of malachite green adsorption onto the activated phosphate rock could have been due to the contribution of both surface adsorption and intraparticle diffusion. Similar results were reported for slow rate adsorption in case of big molecules such as for methylene blue and malachite green dyes, onto fly ash and activated carbon respectively (Kumar et al., 2005; Hameed and El-Khaiary, 2008).

Theoretical calculations

Mulliken atomic charges calculated by the semi-empirical PM3 theoretical method are presented in Figure 6. It shows that, for malachite green dye, negative charges were situated onto carbon atoms while positive charges on nitrogen and hydrogen atoms. In addition, malachite green has three extremity sites, i.e., two alkyl amine groups -N(CH3)2 and =N(CH3)2 + with positive charges +0.105e, +0.103e, respectively; and ring with low positive charge (+6x10-3 e).


From the calculated charges and the electrostatic potential model, it is clear that the negative charge of the phosphate (PO4-3 ) onto the surface of activated phosphate rock interacts with malachite greenalkyl amine groups existing in solution by electrostatic attraction (Figure 7). In addition, the negative charges of phosphate groups (PO4-3 ) facilitated the formation of hydrogen bonding between the positively charged hydrogen atoms in alkyl amine groups of malachite green dye. Therefore, the cooperative contribution of both hydrogen bonding and electrostatic attraction between alkyl amine groups in malachite green and phosphate groups (PO4-3 ) leads to adsorption and enhancement of its capacity. This suggests that the main mechanism for the adsorption behavior of malachite green onto activated phosphate rock is physisorption in nature.


To confirm the physical adsorption nature, Dubinin–Radushkevich isotherm is applied and the mean free energy of adsorption (E, kJ/ mol) is calculated as described earlier (Manlin et al., 2014). Following expression used for this purpose;

E = [1/√-2B] .............. (5)

where B is the constant of sorption energy (mol²/kJ²).

It has been established that, adsorption free energy ranging from 1 to 7 kJ/mol indicate a physisorption process (Li et al., 2004). In this work, we have calculated an adsorption free energy of 0.5 kJ/ mol, attesting a physical adsorption reaction (Figure 8). Similar results have been reportedfor physical adsorption reaction with a free energy E=1kJ/mol using layered double hydroxide (HDL) carbon dots composite (Manlin et al., 2014). A previous study mentioned that the hydrogen bonding took an important role in physical adsorption process (Renault et al., 2008).



Mechanism study based on the theoretical calculations revealed that the malachite green dye can be adsorbed onto activated phosphate rock by the cooperative contribution of hydrogen bonding and ionic interaction formed between the alkyl amine groups -N (CH3)2, =N(CH3)2 + and the negative surface charges of phosphate groups, these double interactions are beneficial to the improvement of adsorption capacities. The calculated Dubinin–Radushkevich means sorption energy (E) gives a value of 0.50 kJ/mol, indicating a physical adsorption in nature and further, confirms the proposed mechanism. The kinetics study reveals that the adsorption is described by the pseudo-second-order kinetic model and taking place in two successive steps, firstly surface adsorption and then slow pore diffusion characterized by a small value of effective diffusion parameter.

Compliance with ethical standards 

Conflict of Interest

The authors declare that they have no conflict of interests


Achkoun, A., Naja, J., M’Hamdi, R., 2012. Elimination of cationic and anionic dyes by natural phosphate. Journal of Chemistry and Chemical Engineering. 6(8), 721-725. Aklil,

Aklil, A.,Mouflih, M., Sebti, S., 2004. Removal of heavy metal ions from water by using calcined phosphate as a new adsorbent. Journal of Hazardous Materials. 112,183-19.

Alzaydien, A.S., 2009. Adsorption of methylene blue from Aqueous Solution onto a Low-Cost Natural Jordanian Tripoli. American Journal of Applied Sciences. (1), 197-208.

Baek, M.H., Ijagbemi, C.O., Se-Jin, O., Kim, D.S., 2010.Removal of malachite Green from aqueous solution using degreased coffee bean. Journal of Hazardous Materials.176, 820-828.

Baraka, N., Assabbane, A.,Nounah, A., Laanab, L., Aît-Ichou, Y.,2009.Removal of textile dyes from aqueous solutions by natural phosphate as a new adsorbent. Desalination 235, 264- 275.

Bekçi, Z., Özveri, C., Seki, Y., Yurdakoç, K., 2008.Sorption of malachite green on chitosan bead. Journal of Hazardous Materials.154 ,254-261.

Chernyshova, I.V., Ponnurangam, S., Somasundaran, P., 2011.Adsorption of Fatty Acids on Iron (Hydr)oxides from Aqueous Solutions. Langmuir.27(16), 10007-10018.

Elouear, Z., Bouzid, J., Boujelben, N., Feki, M., Jamoussi, F., Montiel, A., 2008. Heavy metal removal from aqueous solutions by activated phosphate rock. Journal of Hazardous Materials. 156, 412 - 420.

Ertaş, M., Acemioğlu, B., Alma, M.H., Usta, M., 2010. Removal of methylene blue from aqueous solution using cotton stalk, cotton waste and cotton dust. Journal of Hazardous Materials.183, 421-427.

Forgacs, E., Cserhati, T., Oros, G., 2004. Removal of synthetic dyes from wastewaters: a review. Environment international.30, 953-971.

Hameed, B.H., El-Khaiary, M.I., 2008. Equilibrium, kinetics and mechanism of malachite green adsorption on activated carbon prepared from bamboo by K2CO3 activation and subsequent gasification with CO2. Journal of Hazardous Materials.157, 344-351.

Ho, Y.S., Ng, J.C.Y., McKay, G., 2000. Kinetics of pollutant sorption by biosorbents. review. Separation & Purification Reviews.29 (2), 189-232.

Ho, Y.S., McKay, G., 1999. Pseudo-second-order model for sorption process. Process Biochemistry. 34 (5), 451-465.

HyperChem (Molecular Modeling System)., 2009. Hypercube, Inc, 1115 NW, 4th Street, Gaines ville, FL 32601; USA, (http:// www.hyperchem.com).

Keleş, E., Özer, A. K., Yörük, S., 2010. Removal of Pb2+ from aqueous solutions by phosphate rock (low-grade). Desalination.253, 124–128.

Kumar, K.V., Ramamurthi, V., Sivanesan, S.,2005.Modeling the mechanism involved during the sorption of methylene blue onto fly ash. Journal of Colloid and Interface Science. 284(1), 14-21.

Li, C., Chen, H., Li, Z., 2004. Adsorptive removal of Cr(VI) by Fe -modified steam exploded wheat straw. Process Biochemistry. 39(5), 541–545.

Manlin, Z., Qingfeng, Y., Chao, L., Zenghe, L., Wenxing, W., 2014.Layered double hydroxide carbon dots composite: high performance adsorbent for removal of anionic organic dye.ACS Applied Materials & Interfaces. 6(22), 20225- 20233.

Mezghache, H., Hani, A., 2002. Typologie chimique des phosphates dugisement de Djemi-Djema, bassin de Djebel Onk (Algérie orientale). Géologie Méditerranéenne, pp.95-106.

Marchant, R., Smyth, W., 2001. Microbial decolourisation and degradation of textile dyes. Applied Microbiology and Biotechnology. 56, 81-87.

Ncibi, M.C., Mahjoub, B.,Seffen, M., 2007.Kinetic and equilibrium studies of methylene blue biosorption by posidonia oceanica (L) fibres. Journal of Hazardous Materials.139, 280-285.

Pang, J., Xu, G., 2012.Comparison of the influence of fluorocarbon and hydrocarbon surfactants on the adsorptions of SDS, DTAB and C12E8 at the air/water interface by MD simulation. Chemical Physics Letters. 537, 118-125.

Patil, S., Renukdas, S., and Patel, N.,2011. Removal of methylene blue, a basic dye from aqueous solutions by adsorption using teak tree bark powder. International Journal of Environmental Sciences.1(5), 711-726.

Renault, F., Morin-Crini, N., Gimbert, F., Badot, P.M., Crini, G., 2008.Cationized starch-based material as a new ionexchanger adsorbent for the removal of C.I. acid blue 25 from aqueous solutions. Bioresource Technology.99(16), 7573- 7586.

Rubin, E., Rodrigues, P., Herrero, R., Cremades, J., Barbara, I., Sastre de Vicente, M., 2005. Removal of methylene blue from aqueous solutions using as biosorbent sargassum muticum. Journal of Chemical Technology and Biotechnology. 80(3), 291–298.

Sears, G., 1956. Determination of specific surface area of colloidal silica by titration with sodium hydroxide. Analytical Chemistry.28 (12), 1981-1983.

Sonawane, G.H., Shrivastava, V.S., 2009.Kinetics of decolourization of malachite green from aqueous medium by maize cob (Zea maize): An agricultural solid waste. Desalination. 247, 430-441.

Stewart, J.J.P., 2004. Optimization of parameters for semiempirical methods IV: extension of MNDO, AM1, and PM3 to more main group elements. Journal of Molecular Modeling. 10 (2), 155–164.

Vysotsky, Y.B., Kartashynska, E. S., Vollhardt, D., 2015. Theoretical description of 2D-cluster formation of nonionic surfactants at the air/water interface. Colloid and Polymer Science.293(11),3065-3089. Wick, C.D., Chen, B., Valsaraj, K.T., 2010.Computational investigation of the influence of surfactants on the air-water interfacial behavior of

Wick, C.D., Chen, B., Valsaraj, K.T., 2010.Computational investigation of the influence of surfactants on the air-water interfacial behavior of polycylic aromatic hydrocarbons. Journal of Physical Chemistry C.114 (34), 14520-14527. Young, J.L., Evert, J.E., Richard, J.R., 2005. Sorption mechanisms of zinc on hydroxyapatite, systematic uptake

Young, J.L., Evert, J.E., Richard, J.R., 2005. Sorption mechanisms of zinc on hydroxyapatite, systematic uptake studies and EXAFS spectroscopy analysis. Environmental Science & Technology. 9 (11), 4042- 4048.

Zou, W.,Li, K., Bai, H.,Shi, X., Han, R., 2011.Enhanced Cationic Dyes Removal from Aqueous Solution by Oxalic Acid Modified Rice Husk. Journal of Chemical & Engineering Data.56 (5), 1882-1891.