Magnetic solid-phase extraction and Uv/vis Spectrophotometric determination of trace amount of Copper in vegetable and fruit samples after preconcentration of its Pentetate complex

Document Type : Original Research Paper


Department of Chemistry, Payame Noor University (PNU), Tehran, Iran


Core-shell structure Magnetite@polydopamine (Fe3O4@PDA) nanoparticles have been successfully synthesized and were characterized by FT-IR, TEM, FE-SEM, XRD, EDX, Raman Spectroscopy, VSM and Zeta potential. The effects of the different parameter on copperpentetate complex (CuDTPA) sorption on Fe3O4@PDA were investigated by using Box-Behnken design for Spectrophotometric determination of trace amount of Copper in vegetable and fruit samples after preconcentration of its Pentetic complex. Under optimal experimental conditions, L.R., LOD, and LOQ of Cu (II) were 2-700 µg.L-1, 0.0146 µg.L-1 and 0.0487 µg.L-1 respectively. The Preconcentration Factor 100 was obtained.  Low relative standard deviations (RSD=0.7 %) for 10 replicates are achieved. The method was applied to the determination of trace amount of Cu (II) in food samples and the results showed that no significant difference between the proposed method and Flame Atomic Absorption Spectroscopy.



Metals are normal ingredients in Earth’s crust. The ultra-trace amounts of poisonous metals are usually available in trace concentrations in environmental samples (1). Heavy metals perhaps available in foods naturally, or from contamination during protection, cooking and industrial processes or as an outcome of human processes (such as exhaust gases, industrial emissions, etc.). Among the attendance of heavy metals in foods and long-term effect on health, a relevance exists therefore, it is necessary to keep them at an acceptable low level of toxicity. (2).

Trace elements such as Copper are necessary for trace amount in human health and is straight involved in the Fe metabolism. Consumption of fundamental metal ions from food and beverages are the main origin for human (3). However, the higher level of fundamental trace metals in our body makes well-being danger, such as hepatic difficulty. But owing to these reasons, the precise and accurate determinations of trace metals for analytical chemists are significant (4).

FAAS and/or GFAAS, ICP-OES and/or ICP-MS and electro analytical techniques are the main instrumental techniques for the measuring trace amounts of heavy metal ions (5-8). One of the main difficulty of these instrumental techniques is the attendance of a trace amount of heavy metals lower than the limit of detection (9). Preconcentration techniques, for example, liquid-liquid extraction (LLE), membrane filtration, coprecipitation, cloud point extraction (CPE), solid-phase extraction (SPE) and electro deposition, are common to eliminate this difficulty are used by researcher around the world (10-13).

Different chelators like ethylenediamine tetra acetic acid (EDTA), triethylenetetramine hexaaceticacid (TTHA), diethylenetriamine pentaaceticacid (DTPA), 5, 10, 15, 20- tetrakis(penta fluorophenyl)porphyrin (H2TFPP)(14) and methyl-2-pyridylketoneoxime (MPKO) (15) have been widely studied, however EDTA and DTPA, are more than other common chelators have been used. EDTA is a hexa dentate chelator capable of complexing stoichiometrically with almost any transition metal in the periodic Table (16). The effectiveness of chelator for a specific metal ion by the stability constant (log K) is governed. The log K of EDTA-Cu (II) is 18.8 and the log K for DTPA-Cu (II) is 21.4 at 20 ºC in 0.1M ionic strength.

DTPA (LD50=0.59 g kg-1 for rat; λ max=237 nm) is one of the synthetic poly amino poly carboxylic acids that can form stable complexes with a large number of metal ions, like Cu 2+, Ni 2+, Co 2+, Cd 2+
Zn 2+, Mn 2+ and Ca 2+. Free DTPA contains 5 carboxylic acid and 3 amine groups, and 5 pKa values of 1.79, 2.56, 4.42, 8.76 and 10.42 have been reported (17). Like EDTA, DTPA chelates Cu ions at equimolar ratio over pH 2.5 to 11.0.

Magnetic-SPE (MSPE) which has received worthy attention in this decade, is based on magnetic sorbents with great extraction efficiency and comfortable operation, (18, 19). Magnetic nano particles (MNPs) have recently been the focus of continued research as the adsorbent for MSPE thanks to their high extraction capacity, high surface area, diverse functionalized possibilities and easy separation sufficiency (20). Between the Magnetic nano particles sorbents used, magnetite (Fe3O4) is most popular due to its low toxicity, unique magnetic property and low- price (18, 19). Extraction efficiency and selectivity of bare Fe3O4 MNPs could be improved via the surface modification (21, 22). Fe3O4 MNPs with organic modification like dopamine could be used as a sorbents for the extraction of low concentration of heavy metals from the food samples (23, 24).

Dopamine (DA) that includes amine and catechol functional groups is the first catecholamine. Polydopamine (PDA) films as multifunctional coatings can be easily grown and adhered to a large variety of organic and inorganic materials in pH>7 through self-polymerization of DA (22). Owing to the versatility and easy polymerization process, PDA used for surface modification of solid materials. The use of PDA as a coating agent for the construction of multifunctional nanocomposites can provide various advantages like enhancement of water dispersibility and stacking interaction between the analyte and sorbent (24, 25).

The multivariate design approach is a famous collection of statistical and calculation techniques for designing a set of experiments, expanding a mathematical model, measuring the efficacy of factors and obtaining the best condition of factors for appropriate responses. (26).These methods let the concurrent optimization of a number of variables, are quicker and more cost-effective than normal optimization strategy “one variable at a time (OVAT)”. (19, 27). The first step in multivariate optimization screening the factors that studied in order to get the significant efficacy of the analytical system. After determining the significant factors, in the second step, the optimum operating conditions are gained by using more complex experimental designs like Box-Behnken design (BBD). BBD as a type of multivariate design approach, is a class of rotatable second-order designs based on imperfect factorial designs and by the equation 2 k (k + 1) + cp; in this equation cp is the number of duplications at the central point and k is the number of factors; the number of experiments is obtained (19, 26, 28, 29).

The aim of this study is that stablishing a new, speed, effective and simple MSPE procedure by uv/vis for determination Cu(II) ions after preconcentration of its Pentetate complex in vegetable and fruit samples on Fe3O4@PDA.The magnetic sorbent is characterized by TEM, FT-IR, FE-SEM, XRD, EDX, Raman Spectroscopy, VSM and Zeta potential. The effects of different parameter on the sorption of complex CuDTPA on Fe3O4@PDA were investigated by using BBD.



Tris (hydroxy methyl) amino methane hydro chloride (H2NC (CH2OH)3. HCl) (PubChem CID: 93573), dopamine hydrochloride(PubChem CID: 65340), Penteticacid (diethylenetriamine pentaaceticacid or DTPA) (PubChem CID: 3053), Calcium chloride(PubChem CID: 24854), triethylamine (TEA) (PubChem CID: 8471), ferricchloridehexahydrate (FeCl3.6H2O) (PubChem CID: 24810), ferrouschloridetetrahydrate (FeCl2.4H2O) (PubChem CID: 16211588), sodium hydroxide, acetic acid, hydrochloric acid, potassium hydrogen phthalate, nitric acid, hydrogen peroxide, ethanol and copper(II) chloride (PubChem CID: 24014) were purchased from Merck.

A 100.0 mg L−1 stock solution of Cu2+ was prepared by dissolving 0.0212 gr of CuCl2 in a 100.0 mL volumetric flask and was diluted to the mark with doubly distilled water (DDW). DTPA solution was prepared by dissolving 1.967 gr DTPA, 13.1 mL TEA and 0.5 gr of CaCl2 in approximately 1000.0 mL of DDW and pH was adjusted to 7.3 with 1:1 HCl solution (30). Acetate (pH 5.0, 0.1 M) and/or Tris-HCl (pH 8.5, 0.1 M) buffer solution was prepared by dissolving an appropriate amount of CH3COOH or H2NC (CH2OH)3. HCl in DDW and pH was adjusted by adding appropriate amounts of 0.1 M NaOH or HCl solution.


The pH measurements were performed with a JENWAY model 3510 pHmeter. An ultrasonic water bath (SONICA 2200ETH S3) was used to disperse the nanoparticles in solution. The UV-Vis absorbance spectra were recorded using double beam UV-Visible spectrophotometer (PG Instrument, T80+), equipped with 10 mm quartz cuvettes. Flame atomic absorption spectrometer (Perkin–Elmer Analyst 800) equipped with Zeeman-effect background correction and air/acetylene burner (17 L. min -1) was used for the determination of Cu(II) as reference method. FT-IR spectra on a Perkin-Elmer spectrometer (Spectrum 65) were determined over a potassium bromide pellet in the range 4000-400 cm-1. TEM images at the accelerating voltage of 80 KV were taken with a Zeiss -EM10C. FESEM images were obtained by Hitachi S-4160-30KV scanning electron microscope on Au grid. The crystalline structures of MNPs were studied by XRD analysis on a Bruker D8 Advance diffractometer with Cu Kα radiation at 20 mA and 40 kV. Raman spectra were evaluated by WITec Confocal Raman instrument with a 514 nm laser wavelength. Magnetization measurements were carried out on a BHV-55 VSM. Malvern Zeta sizer ZS were applied for recording zeta potential of particles.

Synthesis of Fe3O4 and Fe3O4@PDA magnetic nano particles

Via the chemical co-precipitation method the Fe3O4 MNPs were prepared (19, 24, 31, 32). 10.4 g of FeCl3.6H2O, mixed with 4.0 g of FeCl2.4H2O, and 1.7 ml of HCl (concentrated) in 50 mL water and in next step the mixture was added drop by drop to 0.5 L NaOH (1.5 M) solution under N2 protection and vigorous stirring at 80°C for half hour. The precipitate were separated by external magnet and washed with 4×50.0 mL DDW, then the products at 50°C for 6 h were dried.

In next step, 1.0 g Fe3O4 nanoparticles were dispersed in 0.5 L Tris-HCl (pH 8.5, 0.1 M) buffer solution for 5 min by using ultrasonic, and then 1.0 g dopamine hydrochloride was added, then the mixture solution was mechanically stirred at 25 °C for 24 h. After reaction, the precipitate Fe3O4@PDA nanoparticles were separated by external magnet and washed with DDW and anhydrous ethanol (1:1) and finally dried at 60 °C for 4 h (32).

Procedure of Magnetic Solid Phase extraction

To prepare working solutions, transfer 2-700 μg L-1 of Cu2+ stock solution (100 mg. L-1) to a series of 150.0 mL beakers, then add 15.0 mL acetate (pH=5, 0.1M) buffer solution and finally dilute each to the mark with DTPA solution. Then, 28 mg of magnetic sorbent (Fe3O4@PDA) was added to the 150.0 mL of solution and the mixture stirred for 25 min minutes and finally, external disk magnet (1.5 T) was used to separate magnetic sorbent (Fe3O4@PDA) from the solution. The adsorbed complex onto the magnetic sorbent was desorbed by 1.5 mL of ethanol and collected for the determination Copper pentetate complex by UV/Vis in λ max=276 nm.

Preparation of vegetable and fruit samples

The vegetable (12 samples) and fruit (6 samples) samples were purchased from local markets in Kermanshah city (west of Iran). Each kind of product, if possible, is chosen from different manufacturers. Each sample was washed thoroughly with tap water and DDW. The shells of samples were robbed when necessary. Finally, the samples were dried for 24 h at 105 ºC and grounded.

One gram of each sample was put in a beaker, and 10.0 mL of HNO3 (concent.) was added to the beaker. The mixture was evaporated near to dryness for 4h at about 130 ºC on a hot plate. After cooling to 25 ºC, 3.0 mL of H2O2 (30%, w/w) was added. The mixture was again evaporated near to dryness. The resulting solution was diluted to 25.0 mL with distilled water (33). The filtration procedure by whatman filtration paper no.1 was made for some samples which were not completely dissolved. Then 2.5 mL of each sample transfer to 200 mL beaker, add to 20.0 mL acetate buffer solution (pH=5, 0.1M) and finally dilute each to the mark with DTPA solution (150.0 mL of this solution used for UV/Vis analysis and 50.0 mL used for FAAS analysis as a reference method). The MSPE technique that described above was used for the intended sample solutions.

Box-Behnken Design

A three-factor, three levels Box-Behnken design was employed in order to achieve the highest amount of the Cu (DTPA) adsorbed by the Fe3O4@PDA. Sorbent mass, pH and the uptake time were investigated as independent variables and Cu (DTPA) desorption was considered as the dependent variable. The experimental domains of these factors were established as preliminary experiments. The variables and levels of the BBD model are given in Table 1. The multinomial equation generated by this experimental design is as follows:

Y= b0 + b1X1 + b2X2 + b3X3 + b11X12 + b22X22

+ b33X3+ b12X1X2 + b13X1X3 + b23X2X3

Where Y is the estimated response (dependent variable), b1 to b33 are the regression coefficients, X1 to X3 are independent variables and b0 is the intercept (19).


Characterization studies

FT-IR spectroscopy was used to verify the functional groups of the synthesized Fe3O4 and Fe3O4@PDA nanoparticles. In the spectrum of bare Fe3O4 (Fig. 1a) the peaks at 443 and 588 cm-1 are respectively related to the octahedral bending and tetrahedral stretching vibration of the Fe-O function group and the peaks at 1625 and 3410 cm-1 are corresponding to the surface-adsorbed water and hydroxyl groups (31). For the Fe3O4@PDA, new broad and strong absorption bands can be seen in the range of 1000-1800 cm-1, which can be attributed to the C-O stretching of phenol compounds and the stretching vibration of the aromatic rings. Finally, result confirms the formation of the polydopamine shell on the surface of Fe3O4 NPs via both chemical and physical adsorption (24).

TEM measurement was used to confirm the nanostructure of Fe3O4@PDA. As shown in Fig. 1b TEM image of Fe3O4@PDA suggest an obvious PDA coat (about 6-10 nm-in-thickness) is immobilized on the surface of Fe3O4 nanoparticles, revealing the core/shell structure of Fe3O4@PDA.

The size and morphology of the Fe3O4@PDA were studied by field-emission scanning electron microscopy (FE-SEM). The FE-SEM image of Fe3O4@PDA was shown that the sorbent was formed in nanometer-sized particles as quasi-spherical with an average diameter of 20±5 nm (Fig. 1c).

In Fig. 1d the existence of PDA shell on the surface of magnetite NPs by the EDX detector coupled to the SEM which showed the presence of Fe, C, N and O were also confirmed, and as well as that figure shown, Au peak refer to grid that by EDX operator was used.

The crystal nano structure of the magnetic sorbent was investigated by powder XRD, as shown in Fig. 1e. The characteristic diffraction peaks in the Fe3O4@PDA at 2θ of 18.1 o, 30.2o, 35.6o,36.2 o, 43.4o, 53.8o, 57.2o ,and 62.8o have corresponded to the diffraction of (111), (220), (311), (222), (400), (422), (511), and (440) of the Fe3O4. The entire diffraction peaks match with the magnetic face centered cubic (FCC) structure of Fe3O4 (JCPDS card number 19-0262). No other characteristic peaks owing to the contamination of hydroxides or hematite were detected. Notably, the XRD pattern of the Fe3O4@PDA has similar diffraction peaks with Fe3O4, and the result have indicated that the PDA coating does not change the crystalline phase of Fe3O4 and polymer shell is formless in nature.

To evidently confirm the formation of PDA shell, Raman spectroscopy was used in investigating the changes of the interface of magnetic sorbent. Four strong and broad bands around 398, 489, 526 and 687 cm-1 can be observed in the Raman spectrum of bare magnetic cores, which correspond to the T2g and A1g modes of symmetry, As shown in Fig. 1f (34). However, after coating with the PDA polymer shell, the intensity of these characteristic bands decreased, whereas two new broad bands around 1577 and 1355 cm-1 are present in the Raman spectrum. They are characteristic bands of polydopamine that can be imputed to the deformation of the catechol group.

The magnetic evaluations were accomplished by VSM at 25 ºC. In the VSM magnetization curve of Fe3O4@PDA nanoparticles, there is no hysteresis, and the residue and coercivity is negligible, result indicating the super-para magnetism of these nanoparticles (Fig. 1g). The saturation magnetization (Ms) value of Fe3O4@PDA (30.4 emu g-1) is lower than Fe3O4 (61.2 emu g-1) (data was not shown) that owing to the existence of nonmagnetic PDA coating. However, the magnetic sensitivity of the prepared sorbent is strong enough to provide a quick, cheap, easy, effective and reusable way to separate the sorbent from samples.

Also, the successful grafting of dopamine molecules was confirmed by the zeta-potential analysis. As shown in Fig. 1h, the naked magnetic cores show positive values(8.89 ± 0.94 mV), while the Fe3O4@PDA nano sorbent because of the deprotonation of the phenolic group on the polydopamine shells, shows a negative zeta-potential (-28.42 ± 0.47 mV).

Optimization conditions for the solid phase extraction 
of Cu (DTPA).

Comparison of the preliminary results on the desorption of 100.0 mL of Cu(DTPA) solution (500.0 μg.L−1) by 0.2 g unmodified and modified Fe3O4 by 2.0 mL ethanol showed that the desorption of Cu (DTPA) was 0.382 ± 0.007 and 1.092 ± 0.006, respectively. Thus, bonded PDA on the surface of Fe3O4 had an effective role in increasing desorption of Cu (DTPA). In order to achieve the best performance, the conditions of desorption of Complex onto Fe3O4@PDA were optimized separately.

Optimization of the sorption step

The optimization step of the Cu (DTPA) sorption onto Fe3O4@PDA was performed using a BBD. The design matrix and the relevant results of BBD experiments to determine the influences of the three variables including pH, uptake time, and sorbent mass are shown in Table 2.

According to ANOVA analysis (Table 3), the model F-value was 98488.47. This shows that the terms in the model have a notable efficacy on the response. The determination coefficient (R2) and adjusted determination coefficient (R2adj) values were obtained 0.9999 and 0.9999, respectively. The value of R2 showed that the value of R2adj was very close to 1, which shows the high significance of the model. The P value for the lack of fit (0.00) also shows that multinomial model is statistically notable for the response. According to the model, X1 to X2X3 are significant model terms for the sorption of Cu (DTPA). Based on BBD, empirical relationships among the Cu (DTPA) sorption (Y) and significant independent variables were expressed by the following:

Y= 0.753500 - 0.096230 X1 + 0.003297 X2 + 0.009315 X3 
+ 0.002800 X1 2 - 0.000241X2 2 + 0.000013 X3 2 + 0.000467 X1X2 - 0.000920 X1X3 + 0.000213 X2X3

Whereas it is noted in Table 1, X1 to X3 denotes coded values of the pH of experiment, uptake time and utilized sorbent mass, respectively.

Fig. 2.demonstrates response surface plots for the relationship between pH, Uptake time and Sorbent mass on the sorption of Cu (DTPA) from the aqueous solution. The curvatures of these plots show the practicable interaction between the variables.

Fig. 2a represents the efficacy of Uptake time and Sorbent mass on the sorption of Cu (DTPA) at constant pH (5.0, 0.1 M). This plot shows that sorbent mass and uptake time have positive mutual effect on sorption of Cu (DTPA) complex. As seen in Fig. 2b by increasing the pH, the absorption of complex; due to hydrolysis of copper ion; is reduced. Thus it can be concluded that with adding time and keep solution in weak acidic pH can be achieved to the greatest absorption of copper complex. And finally, as shown in Fig. 2c, in weak acidic pH with increasing the mass of sorbent, the absorption of complex owing to more absorption capacity of nano-sorbent increased clearly. According to the results of the optimizing absorption step, pH 5.0 (acetate buffer solution, 0.1 M), 28 mg sorbent mass and 25min uptake-time were selected as the optimum conditions.

Optimization of the volume sample.

The efficacy of changes in the volume of sample solution was investigated by using different sample volumes (50.0, 100.0, 150.0, 200.0 and 250.0 mL) with the identical concentration of Cu (DTPA) in optimum pH, sorbent mass and uptake time to get the highest preconcentration factor. Result represented that 150.0 mL of volume sample has the highest sorption (Table 4). Therefore, the suitable volume (150.0 mL) was chosen as volume sample for next studies.

Optimization of type and volume of the eluent

In optimum condition, desorption of Cu (DTPA) complex from the Fe3O4@PDA sorbent was studied by using different eluents including chloroform, Ethyl acetate, Acetone, Aceto nitrile, Methanol, Ethanol, Acetic acid, Methanol-acetic acid (1:1) and ethanol-acetic acid (1:1) (eluent volume of 2.0 mL). Consequently, base on result obtained (Table 5A), Ethanol has been used as an eluent in subsequent studies. To choose the proper volume of eluent, the retained copperpentetate complex on the sorbent were stripped with varying volumes (0.5, 1.0, 1.5, 2.0, 2.5 mL) of Ethanol, therefore, the suitable volume of Ethanol (1.5 mL) was selected as eluent (Table 5B).

Analytical figures of merit.

By employing optimal experimental conditions, the developed sorbent presented good linearity (r = 0.998) in the concentration range of 2-700 μg L−1. The linear regression equation for this range was A = 8 × 10−3C − 0.027, where A is the absorbance and C is the concentration of Cu (DTPA) complex, μg L−1. The limit of detection (LOD) and quantification (LOQ) were 0.0146 and 0.0487 μg L−1, respectively. The relative standard deviation (RSD) resulting from the analysis of 10 replicates of 150 mL solutions containing 100 μg L−1 of Copper pentetate complex was 0.7 %. The preconcentration factor (P.F.) (the ratio of the volume aqueous phase/organic phase) of this method in optimum condition was 100.

Efficacy of diverse ions.

The efficacy of a potential interference caused by common species present in samples was examined. For this purpose diverse ions (up to a maximum tolerance ratio of 200) were added to solutions containing 100 µg L-1 of Copper pentetate complex, and the solutions submitted to the recommended procedure. An ion was considered as interfering when it caused a variation in the absorbance of the analyte greater than 3%. The results (Table 6) show that excess amounts of common cations and/or anions not interfere with the determination of trace quantities of Cu (DTPA) complex. The obtained results show that EDTA have seriously interfering effects, whereas it doesn’t exist in real sample matrix.

Reusability of sorbent.

To evaluate the reusability of sorbent on Copper complex, sorbent was used for 10 times. Sorbent with 3 ml of DDW after each absorption was washed and then dried again as a sorbent was used. This procedure was repeated up to 10 times. The results showed that sorbent can be used up to 5 times with high recovery up to 90 percent (Fig. 3).

Real sample analysis

To evaluate the applicability of this method Copper in vegetable and fruit samples was measured. Moreover, the accuracy of results FAAS measurements were done by official AOAC methods (35). Then, to verify the results and accuracy of the proposed approach and ensure the absence of systematic error T-test was used. Results measured 6 times in these samples (3 times for each method) (Table 7) , and a comparison of the mean concentration of Cu in some vegetable and fruit samples from this study with the results obtained from the some other country (Brazil, Turkey and Egypt) studies are summarized in Table 8.


The magnetic nano-adsorbent coated with PDA was successfully synthesized and used as a novel sorbent for SPE of trace amount of Copper by UV/Vis spectroscopy. The results showed that the nano-sorbent was successfully absorbed Cu (DTPA) and desorbed with the least amount of eluent. In optimal experimental condition (obtained by BBD), the linear range 2-700 µg L-1, limit of detection (LOD) 0.0146 µg L-1, limit of quantification (LOQ) 0.0487 µg L-1, low relative standard deviation (RSD) 0.7% (n=10) and preconcentration factor (P.F) (the ratio of the volume aqueous phase/organic phase) 100 was obtained. The method was used for the determination of trace amount of Cu in different real samples and the results showed that no significant difference between the proposed method and FAAS. Compare the results with other studies also showed that the proposed method due to high preconcentration factor, good accuracy, high sensitivity, ease of method, low cost, high selectivity and biocompatibility; and can be used successfully for determination of trace amounts of copper in different real samples.


The authors are thankful to the Payame Noor University (PNU) Research Council for financial support of this research.


The authors declare that there is no conflict of interests regarding the publication of this manuscript.


  1. L. Juszczak, R. Socha, J. Roznowski, T. Fortuna, K. Nalepka, Food Chem., 113, 538 (2009).
  2. F.L. Coco, P. Monotti, F. Cozzi, G. Adami, Food Control, 17, 966 (2006).
  3. H.A. Panahi, M. Karimi, E. Moniri, H. Soudi, Afr. J. Pure Appl. Chem., 2, 96 (2008).
  4. M.K. Jamali, T.G. Kazi, M.B. Arain, H.I. Afridi, N. Jalbani, R.A. Sarfraz, J.A. Baig, J. Hazard. Mater., 158, 644 (2008).
  5. A.Q. Shah, T.G. Kazi, M.B. Arain, M.K. Jamali, H.I. Afridi, N. Jalbani, G.A. Kandhro, J.A. Baig, R.A. Sarfraz, R. Ansari. Food Chem., 113, 1351 (2009).
  6. H. Gençcelep, Y. Uzun, Y. Tunçtürk, K. Demirel, Food Chem., 113, 1033 (2009).
  7. M. Diaz-Perez, M. Aboal-Somoza, P. Bermejo-Barrera, A. Bermejo-Barrera, Cent. Eur. J. Chem., 6, 520 (2008).
  8. T.G. Kazi, M.B. Arain, M.K. Jamali, H.I. Afridi, N. Jalbani, R.A. Sarfraz, R. Ansari, J. Hazard. Mater., 161, 1391 (2009).
  9. T. Madrakian, M.A. Zolfigol, M. Solgi, J. Hazard. Mater., 160, 468 (2008).
  10. Z. Karim, M. Mumtaz, A. Siddique, S. Anwer, U. Zafar, A. Karim, J. Saudi Chem. Soc., 12, 463 (2008).
  11. X. Chang, H. Luo, Y. Cui, X. Zhu, H. Zhai, Z. Hu, Q. He, J. Mol. Struct., 89, 145 (2008).
  12. Z. Sun, P. Liang, Microchim. Acta, 162, 121 (2008).
  13. C. Araneda, C. Fonseca, J. Sapag, C. Basualto, M. Yazdani-Pedram, K. Kondo, E. Kamio, F. Valenzuela, Sep. Purif. Technol., 63, 517 (2008).
  14. Y. Fazaeli, A.R. Jalilian, M.M. Amini, M.R. Aboudzadeh, S. Feizi, A. Rahiminezhad, K. Yousefi, J. Radioanal. Nucl. Chem., 295, 255 (2013).
  15. M. Ghaedi, A. Shokrollahi, F. Ahmadi, H.R. Rajabi, M. Soylak, J. Hazard. Mater. 150, 533 (2008).
  16. Organic sequestering agents. (1st Ed.), S.A. Chaberek, A.E. Martell, (1959), John Wiley & Sons Inc. New York.
  17. J. Vandegaer, S.A. Chaberck, A.E. Frost, J. Inorg. Nucl. Chem., 11, 210 (1959).
  18. J. Zhang, J. Shao, P. Guo, Y. Huang, Anal. Methods, 5, 2503 (2013).
  19. M. Nekoeinia, M. Kabiri Dehkordi, M. Kolahdoozan, S. Yousefinejad, Microchem. J., 129, 236 (2016).
  20. D. Huang, C. Deng, X. Zhang, Anal. Methods, 6, 7130 (2014).
  21. B.D. Cai, J.X. Zhu, Q. Gao, D. Luo, B.F. Yuan, Y.Q. Feng, J. Chromatogr. A, 1340, 146 (2014).
  22. J.K. Xu, F.F. Zhang, J.J. Sun, J. Sheng, F. Wang, M. Sun, Molecules, 19, 21506 (2014).
  23. M. Faraji, Y. Yamini, M. Rezaee, J. Iran. Chem. Soc., 7, 1 (2010).
  24. H. Veisi, S. Taheri, S. Hemmati, Green Chem., 18, 6337 (2016).
  25. Y. Wang, S. Wang, H. Niu, Y. Ma, T. Zeng, Y. Cai, Z. Meng, J. Chromatogr. A, 1283, 20 (2013).
  26. A.T. Nair, A.R. Makwana, M.M. Ahammed, Water Sci. Technol., 69, 464 (2014).
  27. M. Khajeh, A.R. Akbari Moghaddam, E. Sanchooli, Food Anal. Methods, 3, 133 (2010).
  28. S.L.C. Ferreira, R.E. Bruns, H.S. Ferreira, G.D. Matos, J.M. David, G.C. Brandão, E.G.P. da Silva, L.A. Portugal, P.S. dos Reis, A.S. Souza, W.N.L. dos Santos, Anal. Chim. Acta, 597, 179 (2007).
  29. G.D. Pierini, M.F. Pistonesi, M.S. Di Nezio, M.E. Centurion, Microchem. J., 125, 266 (2016).
  30. Analytical and instrumental techniques in agriculture, environmental and food engineering (2nd Ed.), A. Nag, 2015, PHI Learning Pvt. Ltd, P 31. ISBN: 9788120351189.
  31. H. Veisi, M. Pirhayati, A. Kakanejadifard, P. Mohammadi, M. R. Abdi, J. Gholami, S. Hemmati, ChemSelect, 3, 1820 (2018).
  32. M. Nekoeinia, S. Yousefinejad, M. R. Abdi, B. Ebrahimpour, Iranian journal of analytical chemistry, 4, 1 (2017).
  33. Ş. Tokalıoğlu, F. Gürbüz, Food Chem., 123, 183 (2010).
  34. M. Szybowicz, M. Koralewski, J. Karo, L. Melnikova, Acta Phys. Pol., A, 127, 534 (2015).
  35. W. Horwitz, G.W. Lafimer, Official Methods of Analysis of AOAC International, Sections 971.20, 990.05 (2002).
  36. K.S. Ferreira, J.C. Gomes, J.B.P. Chaves, Food Chem., 92, 29 (2005).
  37. M.A. Radwan, A.K. Salama, Food Chem Toxicol., 44, 1273 (2006).