Abstract

Two coordination compounds, copper (II) fumarate (CuFum) and copper (II) tartrate (CuTart), synthesized from copper (II) with fumaric acid and tartaric acid as ligands and using the slow evaporation method have been applied to study the adsorption of phenacetin in aqueous solution. These compounds were characterized by elemental analysis, IR-FT spectroscopy, and X-ray powder diffraction. The melting points of the synthesized coordination compounds were found to be above 350°C. The influence of parameters such as the initial pH, the contact time, and the initial concentration on the adsorption of phenacetin in an aqueous solution has been studied. The studies showed that adsorption equilibrium was reached after 80 minutes for both coordination compounds; the adsorption capacity increased with increasing phenacetin concentration, and the maximum adsorption capacity was obtained in the acidic medium at pH 4. The adsorbed amount of phenacetin on copper (II) fumarate (CuFum) was 25.158 mg/g while that on copper (II) tartrate (CuTart) was 25.906 mg/g. Nonlinear regression analysis showed the best fit for the Freundlich model isotherm for CuTart with R2 of 0.963 and a Chi-square test (χ2) of 0.529 while for the CuFum material, it is the Redlich-Peterson model with R2 of 0.975 and Chi-square test (χ2) of 0.263. The kinetic study shows that the pseudo-second-order model better describes the adsorption of the two materials. The results show that physisorption and chemisorption participate in the adsorption of phenacetin and that these materials can be used for the elimination of phenacetin in solution.

1. Introduction

Phenacetin, also known by the designations pethoxyacetanilide, aceto-4-phenetidine acetophenetidin, N-acetyl-p-phenetidine, and acetophenetidine, is a non-opioid antipyretic and analgesic. It can be isolated from Bursera grandifolia, a herbal remedy with antipyretic properties [1]. The antipyretic effect was observed by its action on the brain via a decrease in the temperature set point [2]. Phenacetin is also used principally as an analgesic for the treatment of fever and related complications [1, 3]. The analgesic effects were observed through its activities on the sensory tracts of the spinal cord. However, the long-term and chronic consumption of phenacetin led to several toxicological complications ranging from nephrotoxicity to carcinogenicity [4]. The carcinogenicity was observed in the urinary tract and renal pelvis [5]. As a result of these severe complications, phenacetin and drugs containing phenacetin were withdrawn from the market by the order of the U.S. Food and Drug Administration in 1983 [6].

Due to its low cost, phenacetin is still being used for research into the physical and refractive properties of crystals [7]. Conventional methods are being used to remove pharmaceutical compounds from wastewater; some examples of which are chlorination, oxidation, biological treatment, electrochemical treatment, and adsorption. However, the adsorption method has many advantages over the others, such as easy operation, high efficiency, and no risk of highly toxic sub-products [8]. The most used adsorbents for the removal of pharmaceutical compounds from aqueous solution are activated carbons [911] clay minerals [12], polymeric resins [13], and other adsorbents, such as molecular polymers [14], mesoporous materials [15], and metal-organic frameworks [16]. Despite all these, it is still of interest to study modified adsorbents or new adsorbents that may be more efficient at removing pharmaceutical compounds from water than conventional ones.

Metal-organic frameworks (MOFs) are a new class of adsorbent materials, with high porosity and high specific area that is formed by metal ions linked to organic ligand bridges, to generate coordination compounds that extend into space in all dimensions. These materials allow a more flexible design of chemical functionality and easy modification, which makes them more promising materials for the adsorption of pharmaceutical compounds. Among the MOFs, a porous Cu-based MOF based on a pentacarboxylate ligand 2,5-bis(3,5-dicarboxyphenyl)-benzoic acid (H5L), namely, [(CH3)2NH2]{[Cu2(L)·(H2O)2x solvent}n, has been used as a sorbent to remove diclofenac sodium and chlorpromazine hydrochloride from the aqueous solution [17]. These compounds exhibit excellent gas adsorption capacities taking advantage of large permanent porosity and excellent diclofenac adsorption in an aqueous solution.

In this work, we used copper (II) fumarate and copper (II) tartrate as novel adsorbents. These adsorbents have never been used for the adsorption of any pollutants. Hence, this work aimed to evaluate the adsorption properties of these complexes for the removal of phenacetin from an aqueous solution. Isotherms studies were carried out, and the experimental data were analyzed using the two-parameter non-linear regression kinetics, and two and three-parameter non-linear isotherm models.

2. Materials and Methods

2.1. Phenacetin

Phenacetin with a purity of 98% was purchased from Aldrich. The stock solution was prepared by dissolving 250 mg of phenacetin in 100 mL of hot distilled water (75°C). Working solutions of concentration varying between 20–50 mg/L were freshly prepared by diluting the stock solution. NaOH (0.10 M) and HCl (0.10 M) were used to adjust the pH value of the phenacetin solutions. Distilled water was used throughout this study.

2.2. Adsorbents

Copper (II) fumarate (CuFum) was synthesized hydrothermally as follows: In a 50 mL round-bottomed flask containing about 20 mL of distilled water, were introduced fumaric acid, C4H4O4 (196.07 mg; 1.23 mmol) and sodium hydroxide, NaOH (80 mg; 2 mmol) which dissolved instantly under magnetic stirring. To the colorless solution obtained was added copper (II) nitrate trihydrate, Cu(NO3)2.3H2O (384.7 mg, 1.23 mmol). The resulting solution was kept under continuous stirring for 24 hours. At the end of this time, a blue precipitate was formed which was collected by filtration. Subsequently, the blue powder obtained (396.4 mg) was washed with distilled water and dried at room temperature. The yield from this synthesis was 60%. The balanced equation of this reaction is as follows:

Copper (II) tartrate (CuTart) was prepared according to the method described by Al-Dajani et al. [18] as follows: In a round-bottomed flask containing about 20 mL of distilled water, tartaric acid was introduced, C4H6O6 (150 mg, 1 mmol) and sodium hydroxide, NaOH (80 mg, 2 mmol) which gradually dissolved under magnetic stirring. To the colorless solution obtained was added copper (II) chloride dihydrate, Cu(Cl)2.2H2O (170.5 mg, 1 mmol). A pale green precipitate was gradually formed upon magnetic stirring for over 24 hours. The mixture was filtered and the residue dried at room temperature. A green powder (174 mg) was obtained, with a yield of 44%. The balanced equation for this reaction is as follows:

2.3. Samples Characterization

Elemental analysis was performed using a Euro Vector CHNS-O element analyzer. FT-IR spectra were recorded using a Perkin-Elmer FT-IR 100 spectrophotometer. The melting point of the compound was obtained using an SMP3 Stuart Scientific melting point apparatus, while the UV-vis spectrophotometer model GENESYS 10S was employed to record the electronic spectra of the investigated compounds. Powder X-rays diffraction was performed using STOE Stadi-p X-ray powder diffractometer with Cu Kα1 radiation (λ = 1.54056 Å; Ge monochromator; flat samples) in transmission geometry.

2.4. Adsorption

Adsorption experiments were carried out by mechanical agitation at room temperature. For each run, 25 mL of phenacetin solution of a known initial concentration between 20 and 50 mg/L was treated with 35 mg of the two compounds. After agitation, the solution was filtered, and the filtrate was analyzed to obtain the concentration of the residual phenacetin by using the UV/Vis spectrophotometer (GENESYS 10S). Similar measurements were carried out by varying the pH of the solution and contact time. The amount (Qt) of phenacetin adsorbed was calculated using the following expression:where Co is the initial concentration of the phenacetin, Ct is the concentration at time t, V is the volume of the solution, and Qt is the quantity of phenacetin adsorbed at time t while m is the mass of the absorbent.

2.4.1. Effect of Initial pH

To determine the effect of pH, the adsorption of phenacetin by each adsorbent was investigated over a pH range of 2 to 10 at ambient temperature. For each adsorbent, 35 mg was treated with 25 mL of an aqueous solution of 40 mg/L of phenacetin. The initial pH of the solution was adjusted by adding HCl or NaOH solution at a concentration of 0.01 mol/L.

2.4.2. Effect of Contact Time

To determine the effect of agitation time on the adsorption process, 35 mg of adsorbent was agitated in a 25 mL solution of phenacetin with an initial concentration of 40 mg/L for different contact times varying between 5 and 120 minutes. After each contact time t, the solution was rapidly filtered and the residual concentration was determined by spectrophotometry. The amount (Qe) of phenacetin adsorbed was calculated using equation (4).

2.5. Equilibrium Isotherms

To study the equilibrium of the adsorption of phenacetin onto each adsorbent, data from the equilibrium experiments were analyzed using non-linear equilibrium models of two- to three parameters.

Langmuir isotherm: the Langmuir isotherm is often used for the equilibrium adsorption of solutes from solutions. It is expressed as [19]:where Qe is the adsorption capacity at the equilibrium concentration (mg/g); Ce is the equilibrium concentration of adsorbate in solution (mg/L), while Qm is the maximum adsorption capacity (mg/g) and K is the Langmuir constant (L/m).

Freundlich isotherm: the Freundlich isotherm is an empirical equation employed to describe multilayer adsorption. This model predicts that the pharmaceutical concentration on the adsorbent will increase with the increase in the adsorbate concentration in the solution. The model equation is given as [20]:where Kf (L/mg) is the Frendlich isotherm constant and 1/n is the heterogeneity factor which can vary between 0 and 1.

Elovich isotherm: the Elovich isotherm is an equation employed to describe monolayer adsorption. This model predicts that chemical bonds are formed between adsorbents and adsorbates. The model equation is given as [21]:

Here, KE is the Elovich constant.

Temkin isotherm: the Temkin isotherm assumes that the decrease in the heat of adsorption is linear and the adsorption is characterized by a uniform distribution of binding energies. It is expressed by the following equation [22]:where b = RT/B is related to the heat of adsorption (J/mol), R is the gas constant (8.314 J/mol·K), T is the absolute temperature (K), and A is the Temkin equilibrium constant (L/g) corresponding to the maximum binding energy.

Redlich-Peterson isotherm: the Redlich-Peterson is an empirical isotherm that incorporates three parameters. It may be used to represent adsorption equilibrium over a wide concentration range. It combines some elements from both the Langmuir and Frendlich equations, and consequently, it can be employed either in heterogeneous or homogenous systems [23]. It can be described as follow:where ARP (L/g) and BRP (L/mg) are Redlisch–Peterson isotherm constants, β is an exponent which lies between 0 and 1.

Langmuir-Freundlich isotherm: the Langmuir-Freundlich isotherm includes the knowledge of adsorption on heterogeneous surfaces. It describes the distribution of adsorption energy onto the heterogeneous surface of the adsorbent [24]. It can be described as follow:where Qe is the equilibrium adsorbed quantity (mg/g); Ce is the equilibrium adsorbate concentration (mg/L); β is the heterogeneity parameter and it has values between 0 and 1 and KLF is the equilibrium constant of the Langmuir-Freundlich equation (Lβ/mgβ) for a heterogeneous solid.

2.6. Kinetic Models

Four different kinetic models were considered for this study: pseudo-first order, pseudo-second order, Elovich, and intra-particle diffusion.

Pseudo-first-order model: the pseudo-first-order kinetic model describes an adsorption process base on multilayer adsorption. Its equation is generally expressed as follows:where Qe and Qt are the adsorption capacities at equilibrium and at time t, respectively (in mg/g) and K1 is the rate constant for the pseudo-first-order adsorption (L/mg.min)

Pseudo-second-order model: the pseudo-second-order kinetic model was initially proposed as a second order rate equation for the removal of heavy metals from water using natural zeolites and it was based on the strong bond between adsorbent and adsorbate [25].where K2 is the rate constant for the pseudo-second-order adsorption (L/mg.min).

Elovich model: the Elovich empirical equation model was firstly for the adsorption of carbon monoxide onto manganese dioxide [26]. However, this equation is now generally known as the Elovich equation and has been extensively applied to chemisorption data. This equation can be expressed mathematically as follows:where α (mg/g.min) represents the initial rate of adsorption and β (mg/g.min) the desorption rate constant.

Intra-particular diffusion model: the intraparticle diffusion model developed by Weber and Morris is presented as follows [19, 20]:where C (mg/g) is a constant associated with the thickness of the boundary layer and its higher value corresponds to a greater effect on the limiting boundary layer. Kid (mg/g.min−1/2) is the rate constant of the intra-particle diffusion model.

2.7. Error Functions

In this study, non-linear regression was applied using the Microsoft Excel Solver function from Excel Microsoft 2013 for fitting the curve. The best fit for experimental data was determined from the coefficient of determination (R2), residual root mean square error (RMSE), Chi-square test (χ2), Sum Square of Errors (SSE), and Average Relative Error (ARE). The expressions of the error functions are given as follows:where Qe,exp and Qe,cal (mg/g) are the equilibrium capacity of adsorption obtained from the experiment and by calculating from the model, respectively, and N is the number of data points.

3. Results and Discussions

3.1. Characterization of Copper (II) Carboxylate
3.1.1. Copper (II) Fumarate (CuFum)

The synthesized copper (II) fumarate was a blue powder that was found to be insoluble in toluene, distilled water, ethanol, acetone, and acetonitrile. It was soluble in DMSO and its boiling point was as high as 350°C.

(1) Infrared Spectroscopy. The IR spectrum of the copper (II) fumarate (Figure 1) shows the presence of organic ligands used in the synthesis through the typical vibrations of carboxylic groups.

In the analysis region, the complex shows a low-intensity peak around 3566.78 cm−1 corresponding to the O-H group, which could indicate that this compound contains water molecules. This frequency is slightly lower than 3600 cm−1, indicating that these water molecules may be linked to the central metal or engaged in interactions with its surroundings. This spectrum also exhibits a broad band with a maximum centered around 2895.62 cm−1; which is characteristic of aliphatic C-H vibrations of elongation. This band, expected around 3000 cm−1 is slightly shifted towards low frequencies. The bands observed at 1395.41 cm−1 and 1533.41 cm−1 could correspond respectively to asymmetric vibration of COO- and symmetric vibration of COO- in the fumarate, with a difference of 138.01 cm−1 which reflects a bridging coordination mode of this grouping [27, 28]. The absence of any band around ≈1710 cm−1 indicates the total deprotonation of all carboxylate groups in the complex [29]. The band at 1214.78 cm−1 can be assigned to the C-O group of the fumarate. The band observed in the interval 998.97–989.84 cm−1 represents the valence vibrations, that are characteristic of the C=C deformation group. Finally, the absorption bands between 691.92 and 579.38 cm−1 are the characteristic frequencies of the Cu-O bond. From the attributions from the IR spectrum, it can be said that the product formed is a copper complex based on the fumarate ligand.

(2) Elemental Analysis. The experimental results obtained from the elemental analysis of the synthesized compound and the relative theoretical values are presented in Table 1 below.

The comparison of the percentages of the different elements (C and H) obtained experimentally and theoretically, as well as infrared spectroscopy permitted us to propose the empirical formula [Cu(O4C4H2)2(H2O)4] for the compound. The analysis we were able to have was about (C,H,N), but the result shows us that the compound contained only C and H with is in accordance since the reactant used does not contain Nitrogen.

Figure 2 below shows the XRD of the adsorbent copper (II) fumarate obtained between 5 to80°.

(3) Powder X-ray Diffraction. This diffractogram shows sharp Bragg peaks of good intensity indicating that the material has high crystallinity.

From the empirical formula proposed by elemental analysis, and the functional group shown by Infrared spectroscopy, the following structure was proposed for the compound.

3.1.2. Copper.(II) Tartrate (CuTart)

The results of the UV-Vis spectroscopy, IR spectroscopy, and powder XRD of the synthesized copper(II) tartrate are similar to those obtained from the literature [30]. The structure proposed from the analyses is the following:

3.2. Adsorption Study
3.2.1. Effect of Initial pH

The effect of initial pH on the adsorption of phenacetin by the two porous coordination compounds was studied by varying the initial pH of the phenacetin solution. The result is shown in Figure 3. According to this figure, the removal of phenacetin from an aqueous solution is highly dependent on the pH of the solution. Figure 4 shows that the adsorption is maximum at pH 4 with 17.928 mg/g and 19.631 mg/g adsorbed for copper (II) fumarate and copper (II) tartrate, respectively. The decrease in the amount of phenacetin adsorbed at pH greater than 4 can be explained by repulsive electrostatic interactions between an adsorbent and the adsorbate. The amount of phenacetin adsorbed by copper (II) tartrate is greater than that for copper (II) fumarate. This could be explained by the presence in copper (II) tartrate of hydroxide groups which are absent in copper (II) fumarate.

3.2.2. Effect of Contact Time

The time of agitation is a very significant parameter in the adsorption process because it determines the time necessary to reach equilibrium adsorption. The adsorption of phenacetin from an aqueous solution of an initial concentration of 40 mg/L and pH 4 was carried out on 35 mg of each carboxylate at the rotational speed of 150 rpm. The result obtained for the increasing stirring times is shown in Figure 5 which gives the amount of phenacetin adsorbed (Qt) as a function of time.

Figure 5 shows that the adsorption takes place in three phases. The first rapid phase occurs during the first 20 minutes and the second phase takes place a little more slowly, stopping at equilibrium which is 80 min for the metal carboxylates. The rapidity of the first phase can be explained by a large number of sites available at the start of the adsorption process. This availability favors their occupation by the molecules of phenacetin. The same trend was observed by Gadipelly et al. [16] during the adsorption of ciprofloxacin onto MOF-5 (zinc (II) benzene dicarboxylate). As the adsorption continues, the number of available sites decreases and the adsorbent becomes saturated. The amount of phenacetin adsorbed during this phase can be attributed to the π-π interactions between the organic compounds and adsorbents. The adsorption capacity of copper (II) tartrate is 21.381 mg/g while that of copper (II) fumarate is 16.225 mg/g.

3.2.3. Effect of Initial Concentration

The effect of initial concentration on the adsorption of phenacetin was investigated and the results are shown in Figure 6 below.

The increasing adsorption capacities of the porous coordination compounds with the increase in the concentration of phenacetin may be due to π-π interactions between the organic compounds and adsorbents. The π-π interactions in most cases are responsible for the mechanism of adsorption of aromatic compounds [31]. This result can also be explained by saying that, increasing the initial concentration of phenacetin leads to an increase in the driving force of mass transfer and hence an increase in the rate at which phenacetin molecules pass from solution to the particle surface [32]. The quantity of phenacetin adsorption onto all samples of porous coordination compounds did not show a plateau, suggesting that there was no monolayer formation on the surfaces of the adsorbents. The increase in adsorption capacity with an increase in initial concentration indicates that the porous coordination compounds have a high potential for the removal of phenacetin from their solutions.

3.3. Isotherm Study

In general, the adsorption isotherm indicates how the molecules are distributed between the liquid and solid phase when the adsorption processes attain equilibrium. This part aims at finding the models that can accurately describe the data obtained from the experiments of the adsorption of phenacetin onto copper (II) fumarate and copper (II) tartrate. Three models of two-parameter isotherms and two models of three-parameter isotherms were studied. The non-linear plots of these isotherms are given in Figures 7 and 8.

The analysis in Table 2 shows that among the two-parameter isotherms studied, the Freundlich model is the one that best describes the metal carboxylate adsorption phenomenon. This is due to its coefficient of determination which tends toward 1 and the values of the error functions which are lower than those obtained by the other models. The validation of the Freundlich model allows us to say that it is a multilayer adsorption on a heterogeneous surface for the two adsorbents, and the interactions between the adsorbent’s surface and the adsorbate play a significant role in the adsorption process [33]. The value of 1/n less than 1 in the Freundlich model shows that there is an affinity between the surface of the materials and the adsorbate solution. These values ​​support the hypothesis that the surfaces of the different materials are heterogeneous [33]. The change in adsorption energy ∆Q obtained from the Temkin model is positive regardless of the material, which means that the adsorption process is exothermic.

By still using Table 3, and considering the three-parameter isotherms, it can be deduced that based on the low error values ​​and the high values ​​of the coefficient of determination, the best representation of the experimental adsorption data is given by the Redlich-Peterson model. This implies that the adsorption of the phenacetin by the different materials takes place on heterogeneous surfaces. Its heterogeneity parameter between 0 and 1 shows that the model for the adsorption of phenacetin by the adsorbents cannot be reduced to the Langmuir isotherm.

This confirms the hypothesis of the two-parameter isotherm, according to which the Freundlich model is that which best describes the adsorption process. The distribution of residues between [−2; 2], around the x-axis shows that the Freundlich and Redlich-Peterson models are more suitable for describing the adsorption process. The asymmetry at the level of the distribution of residuals around the x-axis suggests that the basic model hypothesis is not followed during the adsorption process.where ARP and BRP are models constants of Redlich-Peterson; β is heterogeneity parameter; A is Temkin model constant; KL, Langmuir constant, KLF, Langmuir-Freundlich constant and KF, Freundlich constant.

3.4. Batch Studies

Pseudo-first order, pseudo-second order, intra-particle diffusion, and Elovich kinetic models were used to investigate and describe the adsorption mechanism. These kinetic models have been presented in non-linear forms (Figure 9) and their calculated parameters are shown in Table 3 below.

The pseudo-first-order model is based on multilayer adsorption on the surface of each adsorbent material. This type of adsorption is based on van der Waal-type interactions between the phenacetin and the copper MOFs. These interactions thus highlight physical adsorption between the adsorbate and the adsorbent. The low values of the error functions of the pseudo-first-order model as well as the coefficient of determination R2 which tends toward 1 suggest that this model best describes the adsorption of phenacetin by copper (II) fumarate. Also, the pseudo-second-order model highlights π-π interaction between the phenacetin and copper (II) fumarate. These interactions, which are of a chemical nature imply that there is monolayer adsorption. The low values of the error functions of the pseudo-second-order model as well as the coefficient of determination R2 which tends towards unity plead in favor of this model for the adsorption of phenacetin by copper (II) tartrate. The equilibrium adsorbed amounts, which are close to the experimental amounts in the case of this model suggest that the adsorption of phenacetin by the adsorbents depends on π-π interactions. The Elovich model assumes that the solid surfaces are energetically heterogeneous.

The comparison of the adsorbed quantity of phenacetin from this present work with those obtained from literature are shown in Table 4 below.

By comparing the adsorbed quantities of phenacetin by copper (II) fumarate (25.158 mg/g) and by copper (II) tartrate (25.906 mg/g) and that adsorbed by activated carbons BP1M (13.18 mg/g), P1M (7.40 mg/g) and B1M (4.11 mg/g) [34], the adsorbed quantities obtained from this work is greater than those obtained in the literature. These coordination compounds are potentials adsorbent of phenacetin.

4. Conclusion

In this work, two porous coordination compounds, CuFum and CuTart have been used for the adsorption of phenacetin in an aqueous solution. These compounds were characterized using elemental analysis, IR spectroscopy, thermogravimetric analysis, and powder XRD. The different results obtained enabled us to propose the chemical formulas and then to deduce the structures of the two coordination compounds. Batch adsorption studies showed that increasing the initial concentration of phenacetin enhances the interaction between it and the coordination compounds, resulting in an increase in adsorption capacity. According to the different error functions and coefficients of determination, the kinetic studies show that the adsorption data best fitted by the pseudo-second-order model is more appropriate to describe the adsorption of phenacetin by CuTart, while for CuFum the pseudo-first-order model is the most appropriate. This indicates that the adsorption of phenacetin on CuTart may be due to chemical interactions while for CuFum it may be due to physical interactions. Isotherm models confirmed that the system is heterogeneous and the adsorption is physical in nature. This study shows that these coordination polymers are suitable for the adsorption of pharmaceutical compounds.

Data Availability

The data used to support the findings of the study can be obtained from the corresponding author upon request.

Conflicts of Interest

The authors have declared no conflicts of interest regarding the publication of this study.

Acknowledgments

This work was the result of the thesis work of the student (now Dr. Alvine Djoumbissie).