Removal of fluoride from polluted waters using active carbon derived from barks of Vitex negundo plant

Abstract

Background

Deleterious effects of fluoride contamination in ground waters on the health of human beings are well known and intensive research on developing de-fluoridation methods is globally pursued. Of the various methodologies, increasing interest is being envisaged in using the adsorption methods based on active carbons derived from plant material. In the present investigation, Nitric acid activated carbon derived from barks of Vitex negundo plant (NVNC) is probed for its de-fluoridation abilities.

Methods

The activated carbon is characterized adopting various physicochemical methods and surface morphological studies are carried out using FT-IR and SEM-EDX techniques. The effect of various parameters such as pH, sorbent dosage, agitation time, initial concentration of fluoride, temperature, particle size and presence of foreign ions on the extraction of the fluoride is studied adopting Batch methods.The adsorption process is analyzed with Freundlich, Langmuir, Temkin and Dubinin-Radushkevich (D-R) isotherms and kinetics of adsorption is studied using pseudo first-order, pseudo second-order, Weber and Morris intraparticle diffusion, Bangham’s pore diffusion and Elovich equations. The methodology developed is applied to real ground water samples.

Results

De-fluoridation is maximum at the pH: 7.0, adsorbent dosage: 4.0g/lit; equilibrium time: 50 min, Particle size: 45μ and temperature: 30 ± 1°C. The correlation coefficient values for the adsorption isotherms: Freundlich, Langmuir, Temkin and Dubinin-Radushkevich are 0.929, 0.998, 0.980 and 0.946 respectively and for kinetic models: pseudo-first-order, pseudo-second-order, Weber and Morris intraparticle diffusion, Bangham’s pore diffusion and Elovich equations are 0.989, 0.994, 0.874, 0.902 and 0.912 respectively. The Temkin heat of sorption, B, and the Dubinin-Radushkevich mean free energy, E, for the activated carbon adsorbent are 0.196 J/mol and 7.07 kJ/mol respectively.

Conclusions

Nitric acid activated carbon derived from barks of Vitex negundo (NVNC) plant is found to be an effective adsorbent for the de-fluoridation of waters. The adsorption process is satisfactorily fitted with Langmuir adsorption isotherm with good correlation coefficient value and it indicates monolayer adsorption. The adsorption kinetics is found to follow pseudo-second-order kinetics. The Dubinin-Radushkevich mean free energy and Temkin heat of sorption confirm the physisorption nature as these are lower than 20kJ/mol. The procedure developed is remarkably successful in de-fluoridation of real ground water samples.

Introduction

Fluoride in drinking water has both beneficial and harmful effects on human health (Maheshwari and Meenakshi 2006). Though fluoride in minute amounts is an essential component for bones and for the formation of dental enamel in animals and humans (Grynpas et al. 2000; Jackson et al. 1973; Fawell et al. 2006; Kumar and Moss 2008; Underwood 1997), its high concentration cause irreversible demineralization of bones and tooth tissues which is known as dental and skeletal fluorosis, damage to the brain, liver, and kidney, headache, skin rashes, bone cancer, and even death in extreme cases (Susheela 2001; Barbier et al. 2010; Gazzano et al. 2010; Ayoob and Gupta 2006; Chaturvedi et al. 1990; Wang and Reardon 2001; Lounici et al. 1997; Srimurali et al. 1998; Savinelli and Black 1958; Ganvir and Das 2011; Chinoy 1991). The maximum permissible limit of fluoride in water is 1.5 mg/l as per WHO standards (WHO 1984, 2004; BIS 1991).

A variety of treatment procedures have been reported for the removal of excess fluoride from polluted waters based on precipitation (Aldaco et al. 2007; Akbar et al. 2008; Reardon and Wang 2000; Cengeloglu et al. 2002; Yadav et al. 2006; Nawlakhe et al. 1975; Saha 1993), ion exchange (Meenakshi and Viswanathan 2007; Castel et al. 2000; Feng Shen et al. 2003; Chubar et al. 2005; Apambire et al. 1997; Singh et al. 1999), reverse osmosis (Sehn 2008; Simons 1993), Donnan dialysis (Tor 2007; Garmes et al. 2002; Hichour et al. 1999), electrodialysis (Lahnid et al. 2008; Menkouchi et al. 2007; Hichour et al. 2000; Adikari et al. 1989; Amer et al. 2001), nanofiltration (Liu et al. 2007), membrane-based methods (Dieye et al. 1998; Mjengera and Mkongo 2003; Lhassani et al. 2001; Mameri et al. 2001), electrocoagulation (Hu et al. 2003), and adsorption on to various adsorbents (Shihabudheen et al. 2006; Onyango et al. 2006; Tripathy et al. 2006; Mohapatra et al. 2004; Raichur and Jyoti Basu 2001). The choice of the method depends on conditions like area, concentration, availability of resources, etc. Among all these techniques, adsorption methods have more advantages because of their greater accessibility, economical, ease of operation, and effectiveness in removing fluoride from water to the maximum extent (Chauhan et al. 2007; Ayoob and Gupta 2008; Venkata Mohan et al. 2007).

Activated carbons prepared from various raw materials exhibit good fluoride uptake capacity (McKee and Jhonston 1934). Activated coconut shell carbon (Arulanantham et al. 1989), zirconium-impregnated ground nut shell carbon (Alagumuthu and Rajan 2010a), zirconium-impregnated coconut shell carbon (Sai sathish et al. 2007), Dolichos lablab carbon (Rao et al. 2009), zirconium-impregnated cashew nut shell carbon (Alagumuthu and Rajan 2010b), Phyllanthus emblica-activated carbon (Alagumuthu and Veeraputhiran 2011), Acacia farnesiana carbon (Kishore and Hanumantharao 2011), Moringa indica-activated carbon (Karthikeyan and Siva Elango 2007), Cynodon dactylon-activated carbon (Alagumuthu et al. 2011), Typha angustata-activated carbon (Hanumantharao et al. 2012a), Pine wood and Pine bark chars (Mohan et al. 2012), and activated carbon derived from steam pyrolysis of rice straw (Daifullah et al. 2007) have been explored for their adsorption nature towards fluoride. These techniques suffer from one or the other drawbacks, and a universally accepted, simple, eco-friendly and economical methods are still alluding researchers.

Hence, in this work, we searched for active carbons derived from plant materials belonging to different classes of plant kingdom. Our primary investigations indicated that there is strong affinity between fluoride and active carbons derived from barks of Vitex negundo plant. So, this work is devoted to study in depth the sorption characteristics of the said active carbon towards fluoride with respect to various physicochemical parameters such as adsorbent dosage, particle size, temperature, contact time, effect of pH, presence of foreign ions, and initial fluoride concentration. Further, the mechanism of sorption is probed in order to establish theoretical grounds for the observed sorption phenomenon by making morphological studies adopting such techniques like Fourier transform infrared spectroscopy (FTIR), scanning electron microscopy (SEM), and energy-dispersive X-ray spectroscopy (EDX). Kinetics of the adsorption process has been also investigated. Further, the adoptability of the methodologies developed in this work has been tested with diverse groundwater samples collected from different places in Guntur District of Andhra Pradesh.

Methods

Plant description

Of the various classes of plant materials tested for their sorption abilities towards fluoride, it has been noted that the active carbon derived from the barks of V. negundo plant shows affinity for fluoride.

The V. negundo plant (Figure 1), known as five-leaved chaste tree, or Vavili or Nalla-vavili belongs to the Lamiaceae family in the plant kingdom. It is a small tree growing from 2 to 8 m in height with reddish-brown barks and is commonly found near bodies of water, grasslands, and mixed open forests. All parts of the plant like leaves, roots, bark, fruits, flowers, and seeds can be used medicinally in the form of powder, decoction, juice, oil, tincture, sugar/water/honey paste, dry extract.

Figure 1

Vitex negundo plant.

Preparation of active carbon

The barks of V. negundo plant were collected in bulk, crushed into small pieces, washed with fluoride-free water, and completely dried under sunlight for 2 days. The dried material was carbonized in muffle furnace in the absence of air at 500°C for about 4 h. After carbonization, the carbon was washed with fluoride-free water; the process was repeated for several times and then filtered. After that, the carbon was dried in air oven at 110°C, and it was sieved into desired particle sizes and then the carbon was subjected to liquid phase oxidation by mixing with 0.1 N HNO₃ and boiled for 2 to 3 h on flame. Then, the carbons were washed with double distilled water to remove the excess acid and dried at 150°C for 12 h. This activated carbon prepared from V. negundo was named as NVNC.

Reagents and chemicals

All the chemicals used were of analytical reagent grade purchased from Merck. India Pvt. Ltd. (Bengaluru, India) and Sd. Fine Chemicals (Mumbai, Maharashtra, India), and all solutions were prepared by using double distilled water throughout this study. Stock solution of fluoride (100 ppm) was prepared and was suitably diluted as per the need. SPADNS solution, Zirconyl acid reagent, Acid Zirconyl-SPADNS reagent, and reference solution were prepared as per the literature.

Experimental procedure

Batch mode adsorption studies were adopted. Test solution of 5 mg/l of fluoride was prepared by diluting appropriate quantity of standard fluoride (100 mg/l) solution with double distilled water because the maximum concentration of fluoride reported in groundwater of most of the fluoride affected areas is around 5 mg/l. From this 5 mg/l of fluoride simulated solution, 100 ml of solution was pipetted out into a 250-ml conical flask at room temperature 30°C ± 1°C, and to it, weighed quantity of the prepared active carbon adsorbent was added. Then, the conical flask along with test solution and adsorbent was shaken in horizontal shaker at 120 rpm. At the end of the desired contact time, the conical flask was removed from shaker and allowed to stand for 2 min for settling the adsorbent, and the adsorbent was filtered using Whatman No.42 filter paper. The filtrate was analyzed for residual fluoride concentration by SPADNS method using UV-visible spectrophometer (Model No: Elico U.V-2600, ELICO, Hyderabad, India) as described in Standard Methods of Water and Waste Water Analysis (APHA (American Public Health Association) 1985) at λ _max 570 nm.

Fluoride ion analysis

The percentage removal of fluoride ion and amount adsorbed (in mg/g) were calculated using the following equations:

$$ \%\ \mathrm{Removal}\ \left(\%R\right) = \kern0.5em \frac{\left({C}_{\mathrm{i}}\hbox{-} {C}_{\mathrm{e}}\right)}{C_{\mathrm{i}}}\times 100 $$

$$ \mathrm{Amount}\ \mathrm{adsorbed}\ \left({q}_{\mathrm{e}}\right) = \kern0.5em \frac{\left({C}_{\mathrm{i}}\hbox{-} {C}_{\mathrm{e}}\right)}{m}V $$

where C _i = initial concentration of the fluoride solution in mg/l, C _e = equilibrium concentration of the fluoride solution in mg/l, m = mass of the adsorbent in grams, and V = volume of the fluoride test solution in liters.

The same procedure has been adopted for the experiments carried out by varying the physicochemical parameters such as adsorbent dosage, pH of the fluoride solution, initial concentration of the standard fluoride solution, particle size, temperature, presence of foreign ions, and agitation time. Further, the same method was adopted in monitoring the concentration of fluoride in the characterization and adsorption studies.

Methodologies adopted for the adsorbent characterization

Physicochemical parameters

Various properties of the bio-sorbent, NVNC, were explored using standard methods, and the most important features of the results obtained were presented in Table 1. The pH for the activated carbon adsorbent was determined using the Elico pH meter, model LI-120, and the pH_ZPC was determined using the pH equilibrium method (Kadirvelu et al. 2000; Marsh and Rodriguez-Reinoso 2006; Newcombe et al. 1993). Particle size was determined using American Standard Test Method (ASTM) sieves (El-Hendawy et al. 2001). Iodine number (ASTM D4607-94 2006; Hill and Marsh 1968), decolorizing power (Girgis and El-Hendawy 2002; Rozada et al. 2005), and other parameters such as apparent density, moisture, loss on ignition, ash, water soluble matter, and acid soluble matter were analyzed by using standard test methods (Namasivayam and Kadirvelu 1997; BIS (Bureau of Indian Standards) 1989).

Table 1 Physicochemical properties of NVNC

Boehm titration

The surface functional groups of oxygen were determined according to Boehm titration (Meldrum and Rochester 1990a; Boehm 1994; Bandosz et al. 1992). One gram of carbon sample was placed in 50 ml of the solution containing 0.05 N of sodium hydroxide, 0.05 N of sodium carbonate, and 0.05 N of sodium bicarbonate. The bottles were sealed and shaken for 24 h, and the mixture was filtered consequently. The excess base was titrated with 0.05 N HCl solution. The number of acidic sites was determined under the assumptions that NaOH neutralizes carboxylic, lactonic, and phenolic groups; that Na₂CO₃ neutralizes carboxylic and lactonic groups; and that NaHCO₃ neutralizes only carboxylic groups. The number of basic sites was calculated from the amount of hydrochloric acid reacted with the carbon.

BET surface area

The Brunauer-Emmett-Teller (BET) surface area was determined using computer-controlled nitrogen gas adsorption analyzer at 77 K by Quantachrome NovaWin - Data Acquisition and Reduction for NOVA instruments version 10.01 (Quantachrome Instruments, Boynton Beach, FL, USA). The BET-N₂ surface area was obtained by applying the BET equation to the adsorption data (Kadirvelu et al. 2000; Hashim 1994; Brunauer et al. 1938).

FTIR analysis

The surface functional groups on activated carbon adsorbent were examined using FTIR. The spectra were measured from 4,000 to 500 cm⁻¹ on a BRUKER VERTEX 80/80v FTIR spectrometer (Bruker AXS, Inc., Madison, WI, USA), optical resolution of <0.06 cm⁻¹, with automatic and vacuum compatible beam splitter changer (BMS-c) option. Anhydrous KBr was used as a pellet material.

SEM

The SEM is one of the most versatile instruments available for the examination and analysis of the microstructure morphology and chemical composition characterizations. The microphotographs of these carbons were recorded using LEO 1420 VP compact variable pressure digital SEM, manufactured by Leo Electron Microscopy Ltd. (Cambridge, UK; beam voltage 500 to 2,000 V, magnifications ×250 to ×65,000, resolution 3 nm at 1,000 V).

EDS

EDS, EDX, or XEDS is an analytical technique used for the elemental analysis to identify what those particular elements are and their relative proportions (atomic %) of a sample. Each element has a unique atomic structure allowing unique set of peaks on its X-ray spectrum and was recorded using BRUKER EDX two-dimensional VANTEC-500 detector.

Results and discussions

The effect of various parameters such as adsorbent dosage, pH, contact time, initial fluoride ion concentration, particle size, and presence of foreign ions on the adsorption of the fluoride ion from aqueous solution by activated carbon adsorbent, NVNC, have been represented by plots as shown in Figure 2.

Figure 2

Effect of various parameters. Adsorbent dosage, pH, contact time, initial fluoride ion concentration, particle size, presence of foreign ions on the adsorption of fluoride ion (left to right).

Effect of adsorbent dosage

The percentage removal of the fluoride ion by adsorption onto activated carbon adsorbent was studied by varying the dosage of the activated carbon adsorbent, NVNC, in the range from1.0 to 9.0 g/l, with a solution of 5 mg/l of fluoride ion concentration, at agitation time of 50 min, at a fixed pH = 7 and at a temperature of 30°C ± 1°C. The results obtained were plotted as percentage removal of the fluoride ion versus adsorbent dosage as shown in Figure 2.

As can be inferred from the graphs that the percentage removal of the fluoride ion increases rapidly with increase in adsorbent dosage and after optimum dosage of 4.0 g/l, there has been no significant change in the percentage removal of the fluoride ion.

Effect of pH

The effect of pH of the fluoride ion solution on the removal of the fluoride ion from aqueous solution was studied by varying the pH from 2 to 11 under optimum parameters, i.e., with a solution of 5 mg/l of fluoride ion concentration, contact times of 50 min, dose of 4.0 g/l for NVNC, and at a temperature of 30°C ± 1°C. The results obtained were plotted as percentage removal of the fluoride ion versus pH as shown in Figure 2.

As pH is less than pH_ZPC (Table 1), the net charge on surface of the solid activated carbon adsorbents is positive due to adsorption of excess H⁺, which favors adsorption due to columbic attraction. So, it is expected that acidic conditions favor the fluoride adsorption on to the surface of the solid activated carbon adsorbent. But according to fluoride speciation, neutral hydrofluoride which has less affinity towards active carbon adsorbents is predominant in the pH at less than 3 (Lagergren 1898), and so, less adsorption is observed at low pH conditions. At low pH values, the fluoride ion converts into neutral HF, and thereby anion exchanging nature is lost.

When pH increases from 2 to 11, the percentage removal of the fluoride ion increases up to 7, and then onwards, the percentage removal of the fluoride ion decreases. The optimum pH range has been found to be 6 to 7 and below, and above this range, the percentage removal of the fluoride ion is less (Figure 2).

Effect of agitation time

In adsorption system, contact time plays an important role, irrespective of the other experimental parameters affecting the adsorption kinetics. In order to study the effect of contact time on kinetics of adsorption of the fluoride ion, the adsorption experiments were conducted and the extent of removal of the fluoride ion was known by varying the contact time from 10 to 100 min, with a solution of 5 mg/l of fluoride ion concentration, at a fixed pH = 7 and at a temperature of 30°C ± 1°C. The results obtained were plotted as percentage removal of the fluoride ion versus contact time (min) as shown in Figure 2.

As contact time increases, initially, the percentage removal of the fluoride ion is increased rapidly and after a certain time, approached an almost constant value indicating an attainment of equilibrium condition at which the rate of adsorption of fluoride onto the surface of the sorbent is equal to the rate of desorption. The rate of removal of the fluoride ion with time is higher at initial stages because of the availability of more active sites on the surface of the adsorbent, and with increase in contact time, the availability of active sites on the surface of the adsorbent decreases and this result in the decrease of the fluoride ion removal rate by the adsorbent. The decreased removal rate indicates the possible monolayer of the fluoride ions on the outer surface and pores of the adsorbent leading to pore diffusion onto inner surface of adsorbent particles (Namasivayam and Kadirvelu 1994; Yadav et al. 2006). From the observed results, the optimum contact time of activated carbon adsorbent NVNC is 50 min.

Effect of initial concentration

The effect of initial concentration of the fluoride ion solution on the extent of removal of the fluoride ion from aqueous solution was studied by varying the initial concentration of the fluoride ion solution from 1 to 12 mg/l under constant parameters, i.e., at a fixed pH = 7, contact times of 50 min with a dose of 4.0 g/l, and at a temperature of 30°C ± 1°C. The results obtained were plotted as percentage removal of the fluoride ion versus initial concentration of the fluoride ion solution as shown in Figure 2.

With an increase in initial concentration of the fluoride ion solution, the percentage removal of the fluoride ion has been decreased due to insufficient number of active sites that are available on the activated carbon adsorbents to adsorb the fluoride ions from highly concentrated solution of the fluoride ions. At low concentrated solution of the fluoride ions, sufficient numbers of active sites are available on activated carbon adsorbents, and hence, most of the fluoride ions interact with the active sites on the activated carbon adsorbent, and thus, percentage removal of fluoride is more.

With an increase in the initial concentration of the fluoride ion solution from 1 to 12 mg/l, the percentage removal (%R) of the fluoride ion is decreased from 99.2% to 49.9%, and adsorption capacity (q _e) increased from 0.248 to 1.497 mg/g for NVNC.

Effect of particle size

The effect of particle size on fluoride removal by NVNC was explored to have a better understanding of the adsorption process with 45 to 150 μ mesh sized particles under constant parameters, i.e., with a solution of 5 mg/l of fluoride ion concentration, at a fixed pH = 7, contact time of 50 min with a dose of 4.0 g/l, and at a temperature of 30°C ± 1°C. The results obtained were plotted as percentage removal of the fluoride ion vs. particle size as shown in Figure 2.

With an increase in particle size of activated carbon adsorbent, the percentage removal of the fluoride ion has been decreased due to insufficient number of surface active sites that are available on activated carbon adsorbent to adsorb fluoride ions. The rate of fluoride adsorption on smaller particles was significantly greater than that on larger sized particles. This is attributed to the fact that the lesser the particle size, the more will be the surface area and the more will be the number of active sites available for adsorption processes for a given amount of adsorbent. The percentage removal (%R) of the fluoride ion increases with decrease in particle size of activated carbon adsorbent, and hence, 45 μ is fixed as the optimum size of the particles of activated carbon adsorbent.

Interfering ions

The effect of interfering ions on fluoride removal by NVNC was explored to have a better understanding of the adsorption process with 50 mg/l concentration of interfering ions such as chloride, nitrate, sulfate, bicarbonate, and phosphate under constant parameters, i.e., with a solution of 5 mg/l of fluoride ion concentration, at a fixed pH = 7, contact time of 50 min with a dose of 4.0 g/l, and at a temperature of 30°C ± 1°C. The results obtained were presented in Table 2 and plotted as shown in Figure 2.

Table 2 Effect of interfering ions on fluoride ion removal from aqueous solution by the activated carbon: NVNC

The impact of interfering ions present in water on fluoride adsorption by the activated carbon adsorbent, NVNC, follows the order ${\mathrm{PO}}_4^{3\hbox{-} }>{\mathrm{HCO}}_3^{\hbox{-} }>{\mathrm{SO}}_4^{2\hbox{-} }>{\mathrm{NO}}_3^{\hbox{-} }>{\mathrm{Cl}}^{\hbox{-} }$

Previous research indicates that chloride and nitrate ions form outer sphere surface complexes while sulfate ions form both outer sphere and inner sphere surface complexes (Onyango et al. 2004). Hence, chloride ions as well as nitrate ions have less interference with fluoride removal while sulfate ions have some significant effect on fluoride removal efficiency. The most important factor affecting fluoride removal efficiency was the bicarbonate alkalinity of the water, and it reduces the positive charge on the active sites of the active carbon and thereby reduces the affinity of the active sites of the adsorbent for fluoride adsorption. This results in the decrease in the percentage removal of the fluoride ion from water.

Phosphate ion which is having high negative charge compared to other anions needs three close surface groups and adsorbed on adsorbents as inner sphere surface complex. This inner spherically adsorbing phosphate ion (Goldberg and Sposito 1984a, 1984b; Zhang and Spark 1990) can significantly interfere with the fluoride ion and hence decrease in the percentage removal of the fluoride ion from water. In the presence of interfering ions chloride, nitrate, sulfate, bicarbonate, and phosphate, the percentage removal of the fluoride ion decreases from 92% to 90.9%, 89.4%, 84.2%, 81.5%, and 79.3%, respectively.

Effect of temperature

The effect of solution temperature was studied by conducting the experiment at different temperatures 303, 313, and 323 K and at optimum conditions, i.e., with a solution of 5 mg/l of fluoride ion concentration, at a fixed pH = 7, contact time of 50 min, dosage of 4.0 g/l for NVNC, and results obtained were plotted as ln K _d vs. 1/T as shown in Figure 3. Thermodynamic parameters of the adsorption process such as change in free energy (ΔG) (kJ/mole), change in enthalpy (ΔH) (kJ/mole), and change in entropy (ΔS) (kJ/mole) were determined at different temperatures by using the equations (Alagumuthu and Rajan 2010a; Karthikeyan and Siva Elango 2007)

Figure 3

Effect of temperature on the adsorption of fluoride ion.

$$ \Delta G = \hbox{-} RT\ \ln\ {K}_{\mathrm{d}} $$

$$ \ln\ {K}_{\mathrm{d}} = \Delta S/R\ \hbox{-}\ \Delta H/RT $$

$$ {K}_{\mathrm{d}} = {q}_{\mathrm{e}}/{C}_{\mathrm{e}}\ \mathrm{and} $$

$$ \Delta G = \Delta H\ \hbox{-}\ T\Delta S $$

where K _d is the distribution coefficient for the adsorption, q _e is the amount of fluoride ion adsorbed on the activated carbon adsorbent per liter of solution at equilibrium, C _e is the equilibrium concentration of fluoride ion solution, T is the absolute temperature in Kelvin, R is the gas constant. ΔG is the change in free energy, ΔH is the change in enthalpy, and ΔS is the change in entropy.

The values of ΔH and ΔS were obtained from the slope and intercept of a plot between ln K _d and 1/T and ΔG values were obtained from the equation ΔG = ΔH – TΔS and tabulated (Horsfall and Spiff 2005; Viswanathan and Meenakshi 2010).

It is observed that with an increase in the temperature from 303 to 323 K (30°C to 50°C), the percentage removal of fluoride ion increases from 92% to 95.8% for NVNC.

As the temperature increases, the thickness of the outer surface of the activated carbon adsorbent decreases and the kinetic energy of the fluoride ion increases, and hence, the rate of diffusion of the fluoride ion increases across the external boundary layer and internal pores of the activated carbon adsorbent.

As can be inferred from the Table 3, the values of ΔH are positive, which indicates the physisorption and endothermic nature of adsorption (Bouberka et al. 2005). The R ² values close to one also indicates that adsorption process is endothermic nature. The positive values of ΔS indicate the increased disorder and randomness at the solid solution interface of the fluoride ion with the adsorbent (Sairam Sundaram et al. 2009). The negative values of ΔG indicate the spontaneous nature of adsorption process.

Table 3 Thermodynamic parameters of fluoride ion adsorption on NVNC

Adsorption isotherms

The adsorption isotherms are one of the most significant methods for representing the adsorption capacity of the adsorbent and the mechanism of the adsorption system. The purpose of an adsorption isotherm is to evaluate the relation between the fluoride concentrations remaining in the bulk solution to the amount of fluoride adsorbed at the solid/solution interface. Four well-known models, Freundlich, Langmuir, Dubinin-Radushkevich, and Temkin adsorption isotherms have been selected for describing adsorption isotherms at a constant temperature for water and wastewater treatment application.

Linear form of Freundlich isotherm equation is as follows:

$$ \mathbf{log}\ \left({q}_{\mathrm{e}}\right) = \mathbf{log}\ {k}_{\mathrm{f} +}\left(\kern0.3em \frac{1}{n}\right)\ \mathbf{log}\ {C}_{\mathrm{e}} $$

Linear form of Langmuir isotherm equation is as follows:

$$ \left({C}_{\mathrm{e}}/{q}_{\mathrm{e}}\right) = \left({a}_{\mathrm{L}}/{k}_{\mathrm{L}}\right){C}_{\mathbf{e}} + \mathbf{1}/{k}_{\mathrm{L}} $$

The significant features of the Langmuir isotherm model can be defined by the dimensionless constant separation factor R _L which is expressed by the following equation:

$$ {R}_{\mathrm{L}} = \mathbf{1}/\ \left(\mathbf{1} + {a}_{\mathrm{L}}{C}_{\mathrm{i}}\right) $$

where, k _F and 1/n are the Freundlich constants, C _i is the initial fluoride ion concentration, q _e (mg/g) is the amount of fluoride ions adsorbed per unit weight of the adsorbent (mg/g), and k _L and a _L are the Langmuir constants related to capacity and energy of adsorption, respectively.

In defluoridation of water using activated carbon adsorbents, linear plot of log (q _e) vs. log (C _e) at different fluoride ion concentrations were found to be linear as shown in Figure 4 and confirmed the applicability of Freundlich isotherm model. When C _e/q _e is plotted against C _e, a straight line was observed as shown in Figure 4. The slope, a _L/k _L, and intercept, 1/k _L, of the straight line showed that the adsorption followed the Langmuir isotherm.

Figure 4

The Freundlich, Langmuir, Dubinin-Radushkevich, and Temkin adsorption isotherm (clock wise direction).

The Freundlich constants k _F and 1/n were calculated from the intercept and slope of the plots, respectively. The Langmuir constants a _L and k _L were calculated from the slope and intercept of the straight line, respectively. The observed linear relationships were statistically important as evidenced from the correlation coefficients (R ² values) close to unity, which indicated that the applicability of these two adsorption isotherms confirmed the heterogeneous surface of the adsorbent and the monolayer coverage of the fluoride ion on the active carbon surface.

According to Hall et al. (1966), the separation factor R _L indicates the isotherm’s shape and the nature of the adsorption process as unfavorable (R _L > 1), linear (R _L = 1), favorable (0 < R _L < 1), and irreversible (R _L = 0). The correlation coefficient (R ²) value was higher for the activated carbon adsorbent NVNC in the Langmuir isotherm than in the Freundlich isotherm, and the value of dimensionless separation factor (R _L), 0.0299, which was found to be fraction in the range of 0 to 1 indicates the favorability of the Langmuir isotherm than the Freundlich isotherm.

Linear form of Temkin isotherm equation is as follows:

q _e = B ln C _e + B ln A where RT/b = B

Linear form of Dubinin-Radushkevich isotherm equation is as follows:

ln q _e = − Βε ² + ln Q _m where ε = RT ln(1 + 1/C _e) where B is the Temkin constant related to heat of sorption (J/mol), A is the Temkin isotherm constant (L/g), b is the Temkin isotherm constant, q _m (mol/g) is the Dubinin-Radushkevich monolayer adsorption capacity, β is a constant related to energy, E = 1√2β, ε is the polanyi potential, R is a gas constant (8.314 J/mol K), T is the absolute temperature, C _e is the equilibrium concentration of the fluoride solution in mg/l, and q _e is the amount of fluoride ions adsorbed per unit weight of the adsorbent (mg/g).

In defluoridation of water using activated carbon adsorbents, linear plot of (q _e) vs. ln (C _e) at different fluoride ion concentrations were found to be linear as shown in Figure 4 and confirmed the applicability of Temkin (R ² = 0.980) isotherm model. When ln (q _e) (mol/g) is plotted against ε ², a straight line with slope β and intercepts ln q _m showed that the adsorption followed the Dubinin-Radushkevich (R ² = 0.946) isotherm as shown in Figure 4.

The Temkin constants A and B were calculated from the intercept and slope of the plots. The Dubinin-Radushkevich constant, β, can be calculated from the slope of the straight line. The observed linear relationships were statistically important as evidenced from the correlation coefficients (R ² values) close to unity, which indicated the applicability of these two adsorption isotherms and confirmed the heterogeneous surface of the adsorbent. And the Dubinin-Radushkevich mean free energy, E, can be calculated from the relation E =1√2β. The mean free energy, E, was found to be 7.07 kJ/mol for the activated carbon adsorbent NVNC. This indicates that the mechanism of adsorption is ‘physisorption’. When the E < 8 kJ/mol, it is an indication of physisorption (Monika et al. 2009) (nonspecific adsorption) dominating the chemisorptions and ion exchange. Physisorption is also called nonspecific adsorption which occurs as a result of long-range weak van der Waals forces between adsorbates and adsorbents. According to Atkins (Atkins 1999) also, if the free energy, E and Temkin heat of sorption (B) are less than 20 kJ/mol, the physisorption is predominant.

The values of the Freundlich and Langmuir Temkin and Dubinin-Radushkevich adsorption isotherm constants together with the correlation coefficients values were presented in Table 4.

Table 4 Adsorption and kinetic parameters

Kinetic study of adsorption

In the present work, some kinetic models, namely, pseudo-first-order model, pseudo-second-order model, Weber and Morris intraparticle diffusion model, Bangham’s pore diffusion model, and Elovich equations are discussed to study the rate and kinetics of adsorption of the fluoride ion onto the activated carbon adsorbent NVNC. The kinetics of adsorption describes the solute uptake rate, which in turn governs the residence time of adsorption reaction.

The pseudo-first-order model

On adsorption of the fluoride ion onto the activated carbon adsorbent, linear plot of log (q _e − q _t) vs. t at different contact times for the extent of removal of the fluoride ion from water is applied to confirm the applicability of pseudo-first-order model. The pseudo-first-order equation is log (q _e–q _t) = log q _e – k ₁ t/2.303

It was found that the plot of log (q _e − q _t) vs. t should give a linear relationship from which values of the k ₁ and q _e could be determined from the slope and intercept of the plot, respectively.

The pseudo-second-order model

On adsorption of fluoride ion on to the activated carbon adsorbent, linear plot of (t / q _t) vs. t at different contact times for the extent of removal of the fluoride ion from water is applied to confirm the applicability of pseudo-second-order model. The pseudo-first-order equation is $t/{q}_{\mathrm{t}}=1/{k}_2{q}_{{}_{\mathrm{e}}}^2-\left(1/{q}_{\mathrm{e}}\right)t$

The plot of (t / q _t) vs. t should give a linear relationship from which the values of the k ₂ and q _e could be determined from the intercept and slope of the plot, respectively.

Weber and Morris intraparticle diffusion model

On adsorption of the fluoride ion onto the activated carbon adsorbent, linear plot of (q _t) vs. t ^1/2 at different contact times for the extent of removal of fluoride ion from water is applied to confirm the applicability of Weber and Morris intraparticle diffusion model. The Weber and Morris intraparticle diffusion equation is as follows:

$$ {q}_{\mathrm{t}} = {k}_{\mathrm{ip}}{t}^{1/2} + c $$

The plot of (q _t) vs. t ^1/2 should give a linear relationship from which the value of the Weber and Morris intraparticle diffusion rate constant, k _ip, could be determined from the slope of the plot.

Bangham’s pore diffusion model

On adsorption of the fluoride ion onto the activated carbon adsorbent, linear plot of log [log (C _i/(C _i − q _t m))] vs. log(t) at different contact times for the extent of removal of the fluoride ion from water is applied to confirm the applicability of Bangham’s pore diffusion model. Bangham’s pore diffusion equation is as follows:

$$ \log\ \left[ \log \left({C}_{\mathrm{i}}/{C}_{\mathrm{i}}-{q}_{\mathrm{t}}m\right)\right] = \kern0.5em \log\ \left({k}_{\mathrm{o}}/2.303\ V\right) + \alpha\ \log (t) $$

The linear plots of log [log (C _i/(C _i − q_tm))] vs. log (t) should give a linear relationship from which the value of the Bangham’s pore diffusion rate constant, k _o, and a constant, α, could be determined from the intercept and slope of the plot, respectively.

Elovich equation

On adsorption of the fluoride ion onto the activated carbon adsorbent, linear plot of (q _t) vs. ln(t) at different contact times for the extent of removal of the fluoride ion from water is applied to confirm the applicability of Elovich equation. The Elovich equation is as follows:

$$ {q}_{\mathrm{t}} = 1/\beta\ \ln \left(\alpha \beta \right) + 1/\beta\ \ln (t) $$

The linear plot of (q _t) vs. ln(t) should give a linear relationship from which the value of the constants α and β could be determined from the intercept and slope of the plot, respectively.

The plots of all these five kinetic models are as shown in Figure 5.

Figure 5

Pseudo-first-order, pseudo-second-order, Bangham’s pore diffusion, Elovich and Weber Morris intraparticle diffusion models (left to right).

The values of pseudo-first-order, pseudo-second-order, Weber and Morris intraparticle diffusion, Bangham’s pore diffusion, and Elovich model constants together with the correlation coefficients values were presented in Table 4.

The experimental data revealed that of the five kinetic models, namely, pseudo-first-order, pseudo-second-order, Weber and Morris intraparticle diffusion, Bangham’s pore diffusion, and Elovich models when correlated with the linear forms, the correlation coefficient value of R ² = 0.994 for the pseudo-second-order model is greater than the other kinetic models, and this indicates that the pseudo-second-order model is the best fit to the experimental data of the present studied adsorption system. The next to follow the order is the pseudo-first-order model with R ² = 0.989, Elovich model with R ² = 0.912, Bangham’s pore diffusion model with R ² = 0.902, and the least is Weber and Morris intraparticle diffusion model (R ² = 0.874).

The correlation coefficient of R ² = 0.912 value for the Elovich equation suggest that the diffusion of the fluoride ion follows the Elovich kinetics pattern and the rate determining step is diffusion in nature, and this diffusional rate limiting is more prominent in fluoride ion adsorption onto the activated carbon adsorbent. The correlation coefficient of R ² = 0.902 value indicates that the adsorption kinetics follows the Bangham’s equation and the adsorption of fluoride ion onto the activated carbon adsorbent is Bangham’s pore diffusion controlled.

The low correlation coefficient value of R ² = 0.874 for the Weber and Morris intraparticle diffusion model indicates that the intraparticle diffusion model does not apply to the present studied adsorption system.

Characterization of the adsorbent NVNC

The different physicochemical properties of the active carbon, NVNC, were presented in Table 1. The surface characteristics make it fit for adsorption of the fluoride ion. Zero point charge (pH_ZPC) of the adsorbent is more important for adsorption processes. In the present study, adsorption of anion is favored because for the NVNC, the obtained result is pH < pH _ZPC. The BET surface area before defluoridation was 262.6 m²/g, and the value reduces to 194.8 m²/g after defluoridation. The presence of surface functional groups of oxygen like phenol, carbonyl, hydroxyl, and lactones were determined according to Boehm titration, and from Table 1, it was clear that for activated carbon, the total basic groups were greater than the total acidic groups. The presence of these groups was confirmed by FTIR results.

FTIR spectroscopy of NVNC before and after defluoridation

It gives confirmation for the presence of specific functional groups on the surface of carbon materials. Several characteristic bands were observed in the FTIR spectrum of NVNC (vide Table 5 and Figure 6), and each of the bands can be assigned to specific functional group based on the previous assignments made in literature.

Table 5 Bands assigned to the surface functional groups of NVNC before and after

Figure 6

FTIR spectra of NVNC before and after fluoride adsorption.

Even though a cluster of functional groups were present on the carbon surface, the prominent among them was a sharp and intense band centered around 1,691.51 cm⁻¹ which was attributed to the carbonyl (-C=O) stretching vibration of quinine or quinone or conjugated ketone (Ji et al. 2007; Biniak et al. 1997; Yu et al. 2008; Ishizaki and Marti 1981; Shin et al. 1997; Moreno-Castilla et al. 1998; Starsinic et al. 1983; Zawadzki 1989). The carbon surface is oxidized by treatment with concentrated HNO₃ leading to the generation of such quinone type carbon functional groups which bear significance in the redox chemistry of carbon materials. Carbonyl functional groups are known to be pronounced in oxidized carbon materials rather than the original parent carbon material (Budinova et al. 2006). In addition, stretching vibration bands of surface and hydrogen bonded -O-H group of alcohols, phenols, and chemisorbed water (Daifullah et al. 2003; Ibrahim et al. 1980; Yang and Lua 2003; Puziy et al. 2003), the peaks pertaining to asymmetric -C-H stretching vibration of aliphatic -CH₃ or -CH₂ groups (Biniak et al. 1997; Yu et al. 2008; Puziy et al. 2003; Rajeshwari et al. 2001); peaks due to the in plane bending vibration of -C-H of methylene group (Budinova et al. 2006; Rajeshwari et al. 2001; Ozgul et al. 2007); peaks due to -C-O stretching in alcohols, phenols, ethers, esters, acids, epoxides, lactones, and carboxylic anhydrides (Shin et al. 1997; Budinova et al. 2006; Rajeshwari et al. 2001; Gomez-Serrano et al. 1994; Figueiredo et al. 1999; El-Hendawy 2003; Park et al. 1997; Attia et al. 2006; Lapuente et al. 1998); the peaks pertaining to the -C=O stretching in carbonyl and carboxyl groups and in lactones (Zawadzki 1989; Nageswara Rao et al. 2011; Fanning and Vannice 1993; Painter et al. 1985; Zhuang et al. 1994); peaks due to the out-of-plane deformation vibrations of -C-H group in aromatic structures (Nageswara Rao et al. 2011; Meldrum and Rochester 1990b, 1990c) have been noted. The absence of specific peak pertains to -C-F, suggesting that the adsorption process was physisorption but not chemisorptions. The shifts or decrease in the percentage of transmittance in FTIR spectra before and after deflouidation in the range 4,000 to 500 cm⁻¹ (Figure 6 and Table 5) indicates that the sorption of fluoride on the surface of the adsorbent is through physisorption.

SEM analysis of NVNC before and after defluoridation

SEM micrographs were studied for surface and morphological characteristics of activated carbon material (vide Figure 7). The examination of the SEM micrographs of the activated carbon material showed that in micrographs, dark areas indicated pores and grey areas indicated the carbon matrix and showed rough surface of the adsorbent that provided large surface area for adsorption. Overall, a well-developed porous surface was observed at higher magnification, and further, randomly distributed pore size was observed in all micrographs. The grey surface area of every micrograph contained smaller microparticles (nm to μm) which might be due to the activated sites or surface functional groups of the carbon. The adsorption of fluoride might be due to the presence of pores and active groups on the surface. SEM micrographs proved the fluoride adsorption.

Figure 7

SEM analysis of NVNC before (left) and after (right) defluoridation at ×5,000 and ×6,000 magnifications.

EDX elemental analysis

The EDX elemental analysis highlighted the presence of carbon and oxygen ions in untreated and fluoride treated sample (Figure 8 and Table 6). Nitrogen is also present but in less mass amount compared with carbon and oxygen. The intensity of the fluoride signals was higher in samples treated with fluoride and was below the limit of detection on the untreated samples. It can clearly conclude that the oxygen element concentration has been reduced after adsorption process, but simultaneously, fluoride concentration is increased in the adsorbent and hence the fluoride ion may replace the ion-containing oxygen atom (OH⁻) on the surface of adsorbent.

Figure 8

The EDX spectra of NVNC before and after defluoridation.

Table 6 Elemental analysis of NVNC by EDX spectra before and after defluoridation

Comparison with other carbon adsorbents

The efficiency of the adsorbent developed in this work, viz., NVNC has been compared with the already existing adsorbents reported in the literature based on the adsorbent mass needed to reduce the fluoride concentration from 4 mg/l (average of fluoride concentration in potable water) to 1 mg/l which is a reasonable concentration for health standards. In this concern, the amount of adsorbent has been estimated using a mass balance equation (Ramos et al. 1999):

$$ V{C}_{\mathrm{i}} = mq + V{C}_{\mathrm{e}} $$

because q and C _e must be in equilibrium, q for adsorbents can be calculated using Langmuir equation, and C _i, C _e, and V are 4, 1 mg/l, and l m³ (1,000 liters), respectively. The results and also the values pertaining to the other active carbons reported in literature have been presented in the Table 7.

Table 7 Comparison study with other carbon adsorbents

From the table, it may be noted that the active carbon developed in this work, namely, NVNC is more effective than commercially available F-4000 carbon and also than many other active carbons reported in literature such as A. farnesiana fruits carbon, T. angustata plant carbon, Lagenaria siceraria shell carbon, and aluminum-impregnated carbon. Though KMnO₄-modified carbon prepared from rice straw, aluminum-supported carbon nanotubes, and carbons loaded with specific chemical moieties pecan nut shells (CMPNS) seem to be scoring better than NVNC, they have their inherent disadvantages. In the case of KMnO₄-modified carbon from rice straw, the optimum pH condition is 2 which is a limiting factor for its application to defluoridation of drinking water. In view of the availability of plenty of raw material, the ease of deriving the active carbons and, further, the simplicity of the technique make the NVNC better than the rest two.

Further, it may be noted that the values of various physicochemical properties of NVNC have been found to be comparable with the other activated carbons, i.e., the moisture content value with the activated carbons of Euphorbia antiquorum L. (Palanisamy and Sivakumar 2008), loss on ignition value with the activated Kaza’s carbons (Sreenivasa Rao et al. 2011), ash content value with the activated carbon of bituminous coal (Cuhadaroglu and Aydemir Uygun 2008), pH value with the activated carbons of E. Antiquorum L. (Palanisamy and Sivakumar 2008) and coconut shell (Gimba and Turoti 2006), and the basicity values with the activated carbon prepared from rice straw (Daifullah et al. 2007).

Moreover, the calculated BET surface area value for the adsorbent, NVNC, has been found to be 262.6 m²/g which is more than some other active carbons such as L. siceraria shell carbon (198.5 m²/g) (Hanumantharao et al. 2012c), peanut hull (208 m²/g) (Periasamy and Namasivayam 1994), and cassava peel (200 m²/g) (Rajeshwari et al. 2001), and hence, NVNC shows more adsorption capacity than the said active carbons.

It may be inferred from this discussion that the prepared carbon, NVNC, exhibits considerably greater fluoride adsorption potential when compared with some of the less attractive and low-cost materials.

Applications

The adoptability of the methodology developed with the new bio-sorbent in this work for removing fluoride has been tried with some real water samples collected from groundwaters in fluoride-affected areas in Vinukonda Mandal of Guntur District of Andhra Pradesh.

The samples were subjected to extraction for fluoride using the said bio-sorbents developed in this work at optimum conditions of pH, equilibration time, and sorbent concentration. The results obtained were presented in Table 8.

Table 8 Fluoride ion concentration (before and after defluoridation) of ground water samples

From Table 8, it is evident that the concentration of fluoride in all groundwater samples collected from various villages of Vinukonda Mandal, Guntur District have been varied from 3.37 to 4.27 mg/l. The data indicates that most of the samples contain excess of fluoride beyond the permissible World Health Organization limit 1.5 mg/l (WHO 1984, 2004; BIS 1991). Hence, in the present work, the defluoridation studies were carried out on these particular samples using the adsorbent, NVNC, in order to reduce the fluoride content below the permissible limit. The concentrations of fluoride ion in these samples after defluoridation were analyzed and reported in Table 8.

It may be inferred that NVNC effectively decrease the fluoride content in groundwater samples to below permissible limits under optimum experimental conditions.

It can be inferred from Table 8 that the methodology developed in this research work using active carbon NVNC is remarkably successful.

Conclusions

A new activated carbon adsorbent has been developed for fluoride removal from aqueous solution in this study. The results indicated that the maximum fluoride adsorption takes place at the optimum pH of 7.0, the adsorbent dosage of 4.0 g/l, equilibration time of 50 min, and at a temperature of 30°C ± 1°C.

It was observed that the adsorption process satisfactorily fitted with Langmuir adsorption isotherm which had good correlation coefficient value indicating monolayer adsorption and also confirmed the heterogeneous surface of the adsorbent due to Freundlich isotherm. The Temkin and Dubinin-Radushkevich isotherms also established the linear relationship which indicated the applicability of these two adsorption isotherms and confirmed the heterogeneous surface of the adsorbent. The Dubinin-Radushkevich mean free energy, E, was found to be 7.07 kJ/mol for the activated carbon adsorbent NVNC. An energy value of <8 kJmol⁻¹ was an indication of physisorption. Hence, adsorption is not restricted to monolayer coverage as purposed for chemisorption.

In kinetic studies of defluoridation, pseudo-first-order, pseudo-second-order, Weber and Morris intraparticle diffusion model, Bangham’s pore diffusion model, and Elovich equations were applied to identify the rate and kinetics of adsorption process. The adsorption process had good correlation coefficient values with pseudo-second-order, Bangham’s pore diffusion, and Elovich equations which indicated that the process fitted with pseudo-second-order model and pore diffusion played a very important role in controlling the rate of the reaction. Even though the plots of Weber and Morris intraparticle diffusion made straight lines with correlation coefficients, they fail to pass through origin. This indicated that the process of the mechanism of adsorption was complex in nature with the more than one mechanism limiting the rate of adsorption, i.e., particle diffusion of fluoride adsorption was more towards the rate controlling step than intraparticle diffusion model. The FTIR studies indicated the participation of the surface sites of the adsorbent in the adsorbent interaction. Characterization of the activated carbon adsorbent through FTIR and SEM-EDX techniques confirmed the adsorption of the fluoride ion on the adsorbent surface. The active carbon developed has been found to be effectively decreasing the fluoride content in real groundwater samples below the permissible limits, and hence, the active carbon NVNC can be successfully applied in wastewater treatment technologies in controlling the fluorides.