Next Article in Journal
Highly Skin-Conformal Laser-Induced Graphene-Based Human Motion Monitoring Sensor
Next Article in Special Issue
Gas Porosimetry by Gas Adsorption as an Efficient Tool for the Assessment of the Shaping Effect in Commercial Zeolites
Previous Article in Journal
Graphene Oxide/Reduced Graphene Oxide Enhanced Noniridescent Structural Colors Based on Silica Photonic Spray Paints with Improved Mechanical Robustness
Previous Article in Special Issue
Preparation of Carbon-Covered Phosphorus-Modified Alumina with Large Pore Size and Adsorption of Rhodamine B
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Design of Hybrid PAH Nanoadsorbents by Surface Functionalization of ZrO2 Nanoparticles with Phosphonic Acids

1
College of Natural and Health Sciences, Zayed University, P. O. Box 144534 Abu Dhabi, United Arab Emirates
2
CP2M-ESCPE Lyon, CNRS-UMR 5128, Université Claude Bernard Lyon 1, 43 Bd du 11 Novembre 1918, 69616 Villeurbanne, France
3
ISA, CNRS-UMR 5280, Université Claude Bernard Lyon 1, 5 rue de la Doua, 69100 Villeurbanne, France
*
Authors to whom correspondence should be addressed.
Nanomaterials 2021, 11(4), 952; https://doi.org/10.3390/nano11040952
Submission received: 10 March 2021 / Revised: 1 April 2021 / Accepted: 2 April 2021 / Published: 8 April 2021
(This article belongs to the Special Issue Nanostructured Materials for Adsorption)

Abstract

:
This study focuses on the preparation of innovative nanocomposite materials based on surface modification of commercial nano-ZrO2 optimized from Brønsted acid–base surface reactions. This surface modification was carried out by direct grafting of suitable phosphonic acids bearing a vinylic or phenylic substituent in aqueous solution. Different loading quantities of the anchoring organophosphorus compounds were applied for each materials synthesis. The resulting nanohybrids were thoroughly characterized by infrared spectroscopy (DRIFT), solid-state nuclear magnetic resonance (NMR), nitrogen adsorption-desorption (BET), thermogravimetric analysis (TG), and X-ray photoelectron spectroscopy (XPS), demonstrating the reliability and efficient tunability of the surface functionalization based on the starting Zr/P ratio. Our nanocomposite materials exhibited a high specific surface area as well as complex porosity networks with well-defined meso-pore. The as-prepared materials were investigated for the adsorption of a mixture of 16 polycyclic aromatic hydrocarbons (PAHs) at 200 ng·mL−1 in an aqueous solution. Adsorption kinetics experiments of each individual material were carried out on the prepared PAHs standard solution for a contact time of up to 6 h. Pretreatments of the adsorption test samples were performed by solid-phase extraction (SPE), and the resulting samples were analyzed using an ultrasensitive GC-orbitrap-MS system. The pseudo-first-order and the pseudo-second-order models were used to determine the kinetic data. The adsorption kinetics were best described and fitted by the pseudo-second-order kinetic model. The correlation between the nature of the substituent (vinylic or phenylic) and the parameters characterizing the adsorption process were found. In addition, an increase of PAHs adsorption rates with phosphonic acid loading was observed.

1. Introduction

Polycyclic aromatic hydrocarbons (PAHs) are hydrophobic organic molecules that resulted from the condensation of at least two aromatic rings, being white or yellowish solids [1]. Their ubiquitous and persistent nature in the environment as well as their proven toxicity towards humans and ecosystems have raised general concern in recent years [2]. Hence, 16 of these chemicals were listed as priority pollutants by the US Environmental Protection Agency (US EPA) due to their common detection in the environment, their acute toxicity, and their ability to accurately represent the vast majority of the several hundred derivatives from their chemical class [3,4]. Their toxicity is mainly due to their carcinogenic properties, and it comes as well from their ability to suppress the immune system [5]. PAHs can occur from many anthropogenic or natural sources like forest fire, volcano eruption or vehicular emission, petroleum-related activity, storage, and transport [6]. The presence of these pollutants in surface water and soil mainly comes from the deposition of these airborne particles, which results in general pollution of PAHs in the environment [6]. Unfortunately, traditional wastewater treatments are not efficient enough to ensure their remediation in an aqueous environment [7]. Therefore, numerous methods have been proposed for the removal of PAHs from environmental matrices [8,9]. One of the most promising techniques is the adsorption method, which could have significant economic and environmental advantages over other proposed methods. PAHs removal by adsorption has been efficiently achieved with activated carbon, biochar, clay derivatives, but each of these proposed adsorbents implies significant drawbacks that need to be addressed prior to their democratization [10]. Among them, nanohybrid adsorbents have gained a lot of popularity during the last decades. These nano-adsorbents present interesting features that make them suitable candidates for PAHs remediation in the environment [11,12].
Different supports have been proposed, whether carbonaceous ones (outside the scope of this article) or natural [13] or manufactured inorganic supports functionalized on the surface or not [14,15]. ZrO2, because of its high chemical stability against corrosion and biocompatibility, has been shown to be highly suitable for applications involving the removal of organic pollutants from the environment [16,17,18]. Thus, unlike titania, which can undergo photocatalytic reaction leading to the degradation of the grafted organic part, zirconia is much less reactive under photoradiation [19]. In this regard, phosphonic acids (PAs) that cannot undergo homocondensation and that show a high affinity and stability for metal oxides have been widely studied for the surface modification of metal oxides such as titania, alumina, iron oxides, indium-tin-oxide (ITO), and zirconia [20,21,22,23,24]. In addition, PAs possess (like organosilicates) very stable P-C bonds to hydrolysis that provide a large variety of functional organic fragments available for functionalization. Finally, the high acid stability of the resulting hybrid has allowed to extend the field of application of these materials to complex oxidants and acidic matrices [24].
In this study, we report the synthesis and characterization of zirconia-based nanomaterials functionalized by phosphonic acids bearing a vinylic (VPA) or phenylic (PPA) substituent as nano-adsorbents suitable for PAH adsorption trough π–π interactions. The surface modification has been achieved by a simple approach of post-functionalization of commercial nanosuspensions to allow easy control of the final properties of the resulting nano-hybrids and take into account the economic prerequisites for possible integration in an industrial process. The grafting of the organic moieties has been achieved with different loading quantities to accurately assert the influence of the nature and the amount of the organic part on the synthesized materials’ resulting physicochemical properties. All materials were therefore fully characterized using a wide variety of characterization techniques (nitrogen adsorption-desorption (BET), infrared spectroscopy (DRIFT), solid state 1H 31P NMR, TGA,Differential Scanning Calorimetry (DSC) XPS, and elemental analysis) and then tested as an effective nanoadsorbent for the removal of 16 PAHs in aqueous solution (200 ppb); each sample being analyzed using a new ultra-sensitive gas chromatography (GC-orbitrap-MS) apparatus [25].

2. Materials and Methods

2.1. Chemicals and Reagents

All chemicals were of technical or analytical grade and used without further purification. Phenylphosphonic acid (PPA) (>98.0%, C6H7O3P, Mw = 158.09 g/mol) and Vinylphosphonic acid (VPA) (>95.0%, C2H5O3P, Mw = 108.03 g/mol) were purchased from TCI (Tokyo Chemical Industry). Zirconium (IV) oxide nanoparticles dispersion (99.9%; 10% by weight aqueous solution of a ZrO2 nanosuspension at pH 4–5) and ethyl acetate (C4H8O2, ACS reagent, >99.8%) were both of analytical grade and purchased from VWR. The 16 US EPA PAH mixture (500 µg/mL in Acetonitrile:Acetone:Toluene (6:3:1)) were purchased from LGC Standards. Deionized 18.2 MΩ·cm water was used throughout the whole experiment.

2.2. Syntheses of Surface Modified ZrO2 Nanoparticules by Phosphonic Acid

The general procedure is as follows: 12.33 g of the ZrO2 nanosuspension (10 mmol of ZrO2) were added to a 50 mL round-bottom flask. Different amounts of VPA solution diluted in distilled water or PPA powder (0.05, 0.1, 0.2, 0.4 mmol of VPA, respectively) were added slowly to the ZrO2 nanodispersion, and the mixture was stirred for 4 h at 100 °C. The molar ratio ZrO2/VPA and ZrO2/PPA ranged from 200:1 to 25:1. At the end of the reaction, the suspension obtained was centrifuged for 10 min at 12,000× g rpm. the solid phase was recovered then washed by adding 30 mL of distilled water, redispersing by sonication for 5 min, and then separating by centrifugation again. The washing procedure was repeated twice with water (2 × 30 mL). The white solid was then retrieved and dried in an oven at 95 °C for 24 h. The resulting pale-yellow solid was ground to a fine powder using a mortar and pestle and then stored in a labeled vial. The resulting materials were designated, respectively, by (ZrO2)x(VPA) and (ZrO2)x(PPA) where x = 200, 100, 50, 25 is the molar ratio related to the organic part grafted onto the inorganic support.

2.3. PAHs Batch Adsorption Test Procedure

Diluted solutions of the PAH mixture have been precisely conducted to assert the reproducibility of each test. An amount of 400 µL of the 16 US EPA PAHs mixture (500 µg·mL−1) contained in the commercial amber sealed bottles was collected using a micropipette and introduced into a 1 L volumetric flask containing 40 mL of absolute anhydrous ethanol. The resulting solution was manually shaken to obtain a homogeneous mixture and completed with distilled water to obtain a solution of 200 ng·mL−1 of PAH in water: ethanol (24:1) mixture. Adsorption tests for each material were achieved by introducing 10 mg of our synthesized nanoadsorbent (weighed on a precision balance) in a series of 125 mL flasks disposed on a multi-stirring platform with 100 mL of the 200 ng·mL−1 standard PAH solution with a stirring rate of 1000 rpm. For the kinetic study, each flask was taken out at a predefined time interval of 5, 15, 30 min and 1, 2, 4, 6 h and 20 mL samples of the resulting mixture were collected into centrifugation tubes. For the isotherm study, different PAH solutions were prepared in the same aqueous medium. The concentrations were varied from 50 to 400 ppb. An amount of 10 mg of ZrO2 or (ZrO2)50(PPA) was suspended in 100 mL of each PAH solution. After stirring the suspension at room temperature for 30 min at 300 rpm, the resulting mixture was collected into centrifugation tubes. The spent adsorbent particles were separated by centrifuging the collected samples at 12,000× g rpm for 10 min. An amount of 10 mL of the supernatant was introduced into a glass tube to be treated by solid-phase extraction (SPE) on commercial Strata-X 33u polymeric sorbent cartridges to extract non-adsorbed PAHs remaining in the sample from the aqueous matrix. An amount of 1 mL of isopropanol (10 vol% of the sample) was added to the tube containing the sample prior to SPE analysis, as our previous tests revealed that Heavy Molecular Weight (HMW) PAHs tend to be sorbed on the walls of the tube without any addition of organic solvents. The SPE procedure included three steps: conditioning the cartridge with successively 4 mL of ethyl acetate and 4 mL of 2.0 vol% isopropanol in water for removing substantial impurities present in the column and also activating the immobile phase; loading the column with the water sample at a controlled flow rate of 1–2 droplets per seconds (10 mL·min−1) followed by a drying step to remove moisture in the cartridge; and then collecting the analytes by eluting the cartridge with 8 mL of ethyl acetate. No preconcentration treatment was needed prior to Gas Chromatography-Mass Spectrometry (GC-MS) analysis since the analyzer system is sensitive enough to detect trace levels (ppb) [26]. Standard solutions of 5, 20, 50, 100, 150, and 200 ng·mL−1 were further elaborated to precisely calibrate the GC-orbitrap-MS system by diluting the starting 200 ng·mL−1 solution with water and extracting by SPE.

2.4. Equipment and Characterization

DRIFT experiments of our surface-modified nano-ZrO2 materials were carried out on a Nicolet 6700 FT-IR spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA) equipped with a Mercury Cadmium Telluride (MCT) detector (number of scans: 64, resolution: 4 cm−1). Nitrogen gas adsorption-desorption isotherms were carried out using a ASAP 2020 instrument (Micromeritics, Merignac, France) at 77 K based on the BET method for specific surface area measurement. Prior to analysis, solids were pretreated at 120 °C for 3 h with a thermal ramp of 10 °C·min−1 allowing their degasification. Determination of the organic grafting moieties of our nanohybrid materials were carried out by thermogravimetric analysis using TGA 2 STARe System (Mettler Toledo, Greifensee, Switzerland) with medium pressure alumina crucibles of 70µL (Ø 7 mm). A heating rate of 10 °C·min−1 was used from 25 °C to 900 °C under air atmosphere. Elemental analyses were carried out by ICP-OES Activa (Horiba Jobin Yvon, Longjumeau, France) in IRCELYON (Villeurbanne). 31P solid-state NMR spectra of our ZrO2 nanomaterials functionalized with phosphonic acid were obtained on a AVANCE III 300 MHz NMR spectrometer (Bruker, Billerica, MA, USA) with Magic Angle Spin (MAS) technique. Solid state MAS 31P NMR spectra were recorded using frequency of 121.4 MHz, at rotation speed of 10 KHz with a 4 mm 15N−31P/1H. The centrifugation of our suspension resulting from our synthesis and batch adsorption test were carried out using an 2.5-Liter Multilab centrifuge LISA (AFI centrifuge, Château-Gontier, France) equipped with an Angular Rotor 8x50–25° at rotation speed of 12,000× g rpm. Samples resulting from our batch adsorption tests were extracted using Rapid Trace SPE Workstation (Zymark)with commercial Strata X 33u polymeric sorbent cartridges (Phenomenex, Torrance, CA, USA). Analyses were carried out on a Trace 1310 GC equipped with a Triplus RSH autosampler and coupled to an Orbitrap-MS analyzer (Thermo Fisher Scientific, Waltham, MA, USA). All devices were from Thermo Fisher Scientific. A TG-5SilMS (Thermo Fisher Scientific, Waltham, MA, USA) column was used (30 m × 0.25 mm I.D, film thickness 0.25 µm) with helium as carrier gas (Helium flow: 1.2 mL/min). The injection was performed in splitless mode (1 µL) at a temperature of 325 °C with a split flow rate of 50 mL/min and splitless time of 0.8 min. The heating gradient was 25.0 °C/min to 240 °C; 5.0 °C/min to 280 °C; 25.0 °C/min to 325 °C (hold time: 8.40 min); and 25.0 °C/min to 350 °C for total running time of 27 min. Orbitrap-MS analyzer was employed using a 70 eV electronic impact set in full scan mode with a solvent delay of 3 min. The ion source and transfer line temperature were set to 250 °C. The scanned mass range was from 50 to 600 m/z.

2.5. Adsorption Kinetics

2.5.1. Pseudo-First-Order Kinetic Equation

This well-known kinetic equation may be depicted as follows [27]:
dq t dt =   k 1   ( q e q t ) ,
where qe and qt are the adsorption capacity of the adsorbent at equilibrium and at time t, respectively [mg/g]. k1 is the pseudo-first-order rate constant [min−1]. When integrating Equation (1), the linear form of the equation becomes:
  ln (   q e q t ) = lnq e   k 1 t ,
The adsorption rate constant, k1 [min−1], can be obtained from the slope of the linear plot of ln (qe − qt) versus t.

2.5.2. Pseudo-Second-Order Equation

The pseudo-second-order model can be expressed as follows [27]:
  dq t dt =   k 2   ( q e   q t ) 2 ,
where k2 is the pseudo-second-order rate constant [g/(mg min)]. On integration, this equation can be re-arranged to give the linear form:
  t q t =   1 k 2 q 2 e   +   t q e   ,
The plot of t/qt versus t gives a linear relationship, which allows the values of qe and k2 to be computed.

2.6. Adsorption Isotherm

Adsorption isotherms measurements were carried out by varying the PAH molecules concentration between 50 and 400 ppb. All the experiments are carried out at room temperature for 30 min according to the kinetic study described in the previous paragraph. A reliable prediction of adsorption parameters can be obtained from equilibrium adsorption isotherm equations. In this study, the evaluation of the experimental results was based on Langmuir and Freundlich isotherms models. Langmuir [28] isotherm supposes monolayer adsorption onto a surface containing a finite number of identical and uniform adsorption sites with no transmigration of adsorbate in the plane of the surface and no interaction between the adsorbed molecules on the adsorbent surface. The linearized Langmuir isotherm equation is given as:
1 q e =   1 q m   +   1 q m K L C e
where qe and Ce are the adsorption capacity (mg g−1) and the PAH concentration (mg L−1) at equilibrium, respectively, qm and KL represent the maximum adsorption capacity of adsorbents (mg g−1), and the Langmuir adsorption constant (L mg−1), respectively. The Freundlich [29] isotherm is an empirical equation. It is another approach that is used for the description of the multilayer and heterogeneous adsorption of molecules to the adsorbent surface. It is expressed by the following linear form:
lnq e =   lnK F +   1 n F   lnC e
where KF and nF are Freundlich isotherm constants. KF (mg·g−1) indicates the adsorption capacity of the adsorbent toward the adsorbate and n is an indicator for the degree of the surface heterogeneity and describes the distribution of the adsorbed molecules on the adsorbent surface. A value of higher n reflects the higher intensity of adsorption.

3. Results

3.1. Characterization

Commercial ZrO2 nanosuspension (10 wt% in H2O, pH 4–5, dp < 100 nm) purchased from Sigma Aldrich (Wassergasse, Switzerland) was previously characterized with XRD powder diffractograms in literature and this analysis confirmed the presence of a mixture of tetragonal (~58%) and monoclinic (~42%) ZrO2 [24]. The commercial ZrO2 nanosuspension was also thoroughly described in the same study with TEM imaging, revealing that the sample was rather polydispersed in size (20 ± 10 nm). The XRD powder diffractograms and TEM images obtained during this study revealed that the grafting process did not modify the morphology and structure of the ZrO2 starting support.
Nitrogen adsorption-desorption isotherms spectra for our synthesized materials and bare ZrO2 are shown in Figure A1 and recorded after the degasification of our compounds at 120 °C for 3 h. All of the resulting nanomaterials exhibited almost identical isotherms to the starting material. All compounds without discrimination displayed type IV(a) isotherms based on the IUPAC terminology [30], representing mesoporous adsorbents where pore condensation occurs just after the initial monolayer-multilayer adsorption on the mesopore walls. Our isotherms exhibit a final inflection point instead of a simple saturation plateau that reveals the complex pore network of our materials. Furthermore, the nanohybrids isotherm presents an H2(b)-type hysteresis loop, which further highlights the complex pore structure in which the network effects are significant. Hysteresis loops generally occur during multilayer physisorption and are usually associated with capillary condensation resulting from either adsorption metastability and/or pore blocking effects. In our case, the H2(b)-type hysteresis loop seems to support the hypothesis of very complex pore networks in our nanomaterials. The isotherm study agrees with the pore size measurement based on the BET method, revealing that all our materials are mesoporous, which is beneficial for our PAH adsorption application. However, the major pore’s width distribution seems to be shifted from the starting material, revealing the pore structure’s interparticle origin that is affected by the surface functionalization of the starting ZrO2 support. As shown on Table A1, the specific surface areas of the synthesized materials are in the same range as the starting material (49–57 m2·g−1), suggesting that our grafting process did not affect the morphology of the inorganic support.
Furthermore, Barrett, Joyner, and Halenda (BJH) pore distribution exhibited in Figure 1 revealed that pore structure was slightly changed in the case of ZrO2-VPA materials but kept in the polydispersed range of 20 ± 10 nm, being consistent with the result of the previous study [24]. However, for PPA derivatives, the major pore size distribution shifts as a function of the Zr/P ratio with a linear increase in pore width for x = 25, a trend that was not revealed by the BET analysis of the VPA derivatives compounds. Since pore structure is defined as intergranular porosity for our nanomaterials, pore width seems to be directly related to the nature of grafting moieties. This is also attested by the isotherm spectra of our PPA materials showing desorption branches of the hysteresis tending towards a rather steep desorption branch, a characteristic of H2(a) loops for highly loaded materials. This observation results from a pore-blocking phenomenon in a narrow range of pore necks, which can be explained by the size of phenyl substituents being wider than vinyl ones [31].
Characterization of the effectiveness of the surface functionalization of our nanohybrid materials has first been carried out using DRIFT analysis, since this technique allows the direct analysis of our powders without prior preparation and result in brief analysis where the data are collected on the overall sample. Both ZrO2-VPA and ZrO2-PPA spectra, shown in Figure 2, were recorded in absorbance to clearly emphasize on representative peak of the chemical transformation between bare ZrO2 and the resulting ZrO2-VPA and ZrO2-PPA nanohybrids.
All DRIFT spectra exhibited the same bands regardless of the Zr/P ratio of the materials. Strong absorption bands around 800 and 1464 cm−1 can be easily ascribed to the Zr-O bond since these bands are exhibited in both the bare ZrO2 and ZrO2-VPA or ZrO2-PPA spectra. This assignment was further supported on these ZrO2-PPA spectra, since we observed a clear linear disappearance of the Zr–O band at 1464 cm−1 as a function of the increase in the Zr/P ratio. In addition, the broad band centered at 3300 cm−1 is characteristic of the O-H bond and can be attributed here to surface hydroxyl for bare zirconia and to synthetic solvents. The P–O stretching peak, at 986 cm−1 for ZrO2-VPA, was slightly shifted to 998 cm−1 for ZrO2-PPA materials and was even more covered by the characteristic aromatic absorption occurring in the fingerprint between ~1000 and ~1100 cm−1 [31]. The usual characteristic weak to medium absorption band of P–O–H located near 2600 and 2200 cm−1 was not displayed on our spectra, which may be due to the low concentration of phosphonic acid in our samples. Same observations were made for TiO2-VPA materials [24]. In the case of ZrO2-VPA spectra, the 3087 cm−1 band results in the adsorption of the C–H stretching vibration of the vinylic C sp2, represented by the adsorption peak of the C=C stretching vibration at 1616 cm−1. In the case of ZrO2-PPA spectra, the medium absorption peak at 3062 cm−1 results from the C–H stretching vibration of the aromatic function, further evidenced by the weak adsorption band of the C=C stretching vibration occurring at 1597 cm−1. The medium absorption peak at 1438 cm−1 was characteristic and resulted from the P–C bond’s vibration linked to an aromatic substituent. Moreover, the band intensity linearly increases with the Zr/P ratio of our composite materials, supporting the evidence that our grafting was efficiently achieved. It is worth noting that this band’s absorption frequency was shifted based on the substituent bearing the carbon linked to the phosphorus; 1403 cm−1 for the vinyl and 1438 cm−1 for the phenyl. As with ZrO2-VPA, ZrO2-PPA compounds exhibited the absorption of the vibrational stretching band of O–P–O found at 1081 cm−1, which is still moderately covered by the zirconia absorption area. Furthermore, characteristic P=O stretching vibration was also retrieved in the case of ZrO2-PPA materials but significantly shifted at 1153 cm−1, whereas it was found at 1273 cm−1 for ZrO2-VPA samples, which can induce a different binding conformation between the two organophosphorus compounds in regard to zirconia surface modification since this bond has been reported to be very sensitive to association effects [32]. The remaining band at 1273 cm−1 is highly characteristic of the P=O stretching vibration, which is almost independent of the type of compound or the size of the substituents and can easily help us address the grafting efficiency. In addition, the presence of a P=O function for our material induces that the grafting occurs preferentially by the reaction of the P-OH function on surface hydroxyls, whether they are monocoordinated or bicoordinated. Finally, the intensity of peaks ascribed to phosphonic acid grafting was significantly and linearly shifted based on the Zr/P ratio, which could quantitatively address the grafting quantity present in the final material.
Solid state 31P NMR has been used to better apprehend the coordination of the anchoring group onto inorganic support thanks to a wide variety of observed nuclei in organic grafting moieties. The solid-state 31P MAS NMR spectra of our ZrO2-VPA and ZrO2-PPA samples are given in Figure 3. In the case of ZrO2-VPA, all our materials spectra displayed a distinct resonance downfield shift at 7.4 ppm and broadened in regard of the starting vinylphsophonic acid reagent, which, according to the literature, gave a sharp peak at 18 ppm [33]. This downfield shift underlines the fact that the grafting of the vinylphosphonic acid was correctly achieved on the ZrO2 surface. Moreover, the signal’s broad characteristic results in the sensitivity of the analysis to the variation of bond angles; therefore, coordination sites with slightly distorted geometries induce many different peaks that result in a broad signal covering each of them [19,20,21,22,23,24]. Studies with 31P NMR have shown that the predominant binding state for monophosphonic acids on zirconia surface was generally bidentate and resulted in a peak at around 7 ppm, which correlates with our 31P NMR study and in good agreement to our DRIFT data that revealed P=O stretching vibration. Furthermore, a well-defined peak appeared on the (ZrO2)200(VPA) spectrum at around 12 ppm, which seems to be induced by a tridentate coordination mode based on an identical study conducted on surface-modified nano-TiO2 by VPA [31]. In the case of ZrO2-PPA, the spectra displayed a distinct downfield resonance shift at 8.7 ppm and broadened in comparison to the bare phenylphsophonic acid reagent that was reported in the literature to give a precise, sharp peak at 20.8 ppm [33]. The observed shift between our VPA and PPA materials is in the same range of those reported in the literature [24]. For the most PPA-loaded samples, a distinct downfield peak at −5.60 ppm reveals molecular zirconium phosphonate presence and the corrosion of ZrO2, which could be explained by the addition of pure non-solubilized PPA [34].
TGA experiments were carried out on all synthesized ZrO2-VPA and ZrO2-PPA materials to quantitatively estimate the VPA and the PPA content. Figure 4 shows the thermogravimetric analysis recorded under air from 25 °C to 900 °C with a thermal ramp of 10 °C/min. In the case of ZrO2-VPA, all curves displayed two defined major mass losses, one from 25 °C to 280 °C and the other from 280 °C to 580 °C. In the case of ZrO2-PPA, TGA curves reveal three major mass losses step ranging from 25 °C to 280 °C, 280 °C to 430°C, and 430 °C to 580 °C, respectively. The first major mass loss for both ZrO2-VPA and ZrO2-PPA was ascribed to the evaporation of remaining water and solvents in the samples, further evidenced by the endothermic step appearing on the DSC analysis conducted afterward for this range of temperature (Figure A2). For ZrO2-VPA, the second mass loss was due to the combustion of grafted surface vinyl phosphonate as well as the dehydration of the surface hydroxyl groups of the zirconia. The mass loss evolution was in agreement with the starting Zr/P ratio of the materials since the most concentrated organic material shows the highest mass loss in this temperature range. DSC experiments were carried out to try to precisely attribute the different temperature ranges where the surface hydroxyl and PAs were degraded. Based on our study conducted for TiO2 [23], DSC should have revealed two distinct exothermic peaks that could have helped us determine each mass loss and efficiently measure our nanohybrids’ organic part. However, in our case, the organic moiety quantity was too low to induce a well-defined DSC curve for our materials, even for the most concentrated one (ZrO2)25(VPA). Unlike VPA derivatives, our ZrO2-PPA compounds have presented two distinct mass losses in the range of 280 °C to 580 °C. The mass loss occurring from 280 °C to 430 °C has been attributed to the combustion of the remaining surface hydroxyl function after surface modification, whereas the one from 430 °C to 580 °C was clearly assigned to the organic grafting degradation since the mass loss quantity is linearly proportional to the initial Zr/P ratio of the considered materials. By plotting the mass loss between 25 and 280 °C and that between 280 and 580 °C as a function of the theoretical Zr/P ratios (Figure A3), we can efficiently demonstrate that the mass losses occurring followed a linear trend with a high regression coefficient of 0.90 and 0.94, respectively, for the first and combined second and third mass losses. Nevertheless, these mass losses cannot be used to accurately calculate the number of organic moieties present on the surface of our materials since phosphate residues are still present after calcination at 900 °C as shown by previous TGA studies on similar materials [24]. Thus, trying to quantify the grafting part with given TGA mass losses efficiently would lead to an underestimated value of the organic moieties present in the nanohybrid.
In order to accurately determine the Zr/P ratio of our nanomaterials, elemental analyses were carried out on the most representative samples of the VPA and PPA-based derivatives. The elemental analyses confirmed the grafting of a quantity of P (Table 1) in adequacy with the initial ratios in the case of ZrO2-VPA. However, in the case of ZrO2-PPA, they showed lower than theoretical amounts of P, which could confirm that some of the PPA were eliminated as molecular species such as phosphonates observed in 31P NMR.
The XPS study was carried out to match the observation of our solid NMR analysis and to accurately describe the surface coordination mode of our phosphonic acids on our nanozirconia, based on a study that examined the use of XPS to efficiently discriminate the binding modes of PPA on ZrO2 [35]. XPS studies have been carried out on our (ZrO2)xPPA materials (Figure A4). They have allowed us to confirm that the total amount of phosphorus detected in XPS increases as the Zr/P ratio of our materials increases, proving once again that our approach is effective (Table 1). Finally, it would seem that the majority of the PPA was grafted according to a tridentate mode and no monodentate mode (Figure A5).

3.2. Adsorption Efficiency and Kinetic Study

Adsorption kinetics studies have been conducted to determine the extraction efficiency and the selectivity of our (ZrO2)x(VPA) and (ZrO2)x(PPA) nanomaterials towards the targeted PAHs. These tests have been carried out with standard aqueous solutions of the 16 US EPA PAHs at 200 ng·mL−1 (ppb) for different contact times of 5, 15, 30, 60, 120, 240, and 360 min with 1000 rpm stirring. Adsorption tests have shown an efficiency higher than 90% for all PAHs. In general, the adsorption equilibrium was reached within 30 min of contact without saturation of the nanoadsorbent, regardless of the Zr/P starting ratio, since extraction percentages were found in the same range for each material after 6 h of contact.
The kinetic parameters of PAHs adsorption by our adsorbents were determined by applying the pseudo first-order model (Equation (2)) and the pseudo second-order model (Equation (4)), respectively, to the experimental data. It was found that the values of the qe calculated from the pseudo first-order model did not fit the corresponding experimental qe values. In addition, the correlation coefficient (R2) values determined for this model were less than 0.85. On the other hand, the pseudo second-order model’s application led to similar experimental and calculated qe values, as well as R2 values that were within the range 0.95–0.99. This result indicates that the PAH adsorption mechanism is better explained in terms of the pseudo second-order reaction model. The main assumptions of this kinetic model are the following:
Adsorption occurs only at specific binding sites located on the surface of the adsorbent.
The adsorption energy does not depend on the formation of a layer on the adsorbent surface.
The second-order equation governs the rate of the adsorption process.
No interaction occurs between the molecules adsorbed on the surface of the adsorbent.
Figure A6 shows an example of applying the pseudo-second-order model to the data in the form of t/qt versus t plots. These were linear in all cases, allowing the values of k2 and qe to be calculated from the corresponding slopes and intercepts. All kinetic parameters obtained using pseudo-second-order models are listed in Table 2.
Although the surface loading of ZrO2-PPA is low compared to VPA (Table 2), all (ZrO2)xPPA nanomaterials displayed better affinities to PAHs than (ZrO2)xVPA with adsorption rates between 50 and 230 g·mg−1·min−1 % for Heavy Molecular Weight (HMW) PAHs (Figure 5), whereas the adsorption rate constant of ZrO2-VPA was in the range of 10–70 g·mg−1·min−1. Figure A6 shows the adsorption of benz(a)anthracene using (ZrO2)100(PPA) or (ZrO2)100(VPA). (ZrO2)100(PPA) showed 97% of adsorption during the first 5 min of the test compared to (ZrO2)100(VPA), which showed 85% of adsorption during the same short period. In addition, (ZrO2)50(PPA) and (ZrO2)100(PPA) showed higher adsorption rates than (ZrO2)200(VPA), revealing some influence of the organic charge on the kinetic part of the adsorption. Figure A7 shows the adsorption curves of benzo(a)pyrene using (ZrO2)50(PPA), (ZrO2)100(PPA), and (ZrO2)200(PPA). Adsorption curves showed that a higher PPA content resulted in a more stable system with a higher adsorption capacity.
Furthermore, adsorption rate constants presented in Figure 5 showed competitive adsorption between Light Molecular Weight (LMW) and Heavy Molecular Weight (HMW) PAHs, which is more significant in the case of ZrO2-PPA nano adsorbents. HMW PAHs have more significant affinities to the adsorbent surface than LMW PAHs which have an adsorption rate constant less than 30 g·mg−1·min−1 in both ZrO2-PPA and ZrO2-VPA. However, for adsorption capacities, both showed similar values in the range of 14.5–19.5 µg·g−1.
In general, there is little difference in the case of VPA, both in terms of capacity and affinity. However, in the case of PPA, capacity remains almost unchanged, but affinity is increased by increasing the graft content to reach a pseudo-plateau from a ratio greater than 1/100. This behavior can be correlated with phosphorus species’ content on the surface of ZrO2 determined by XPS. As shown in Figure A5, the amount of PO3- and PO2- on the ZrO2 surface increased with the amount of PPA added until a plateau was reached. It started between the ratio 1/100 and the ratio 1/50.
Isotherm parameters studied for ZrO2 and (ZrO2)50(PPA) were obtained by fitting the adsorption equilibrium data to the isotherm models, and are listed in Table A3. It can be noticed that the adsorption isotherms are fitted better by the Freundlich model since the Langmuir parameters had negative values, suggesting that PAH adsorption onto ZrO2 and (ZrO2)50(PPA) is a heterogeneous coverage. The correlation coefficient (R2) values for the Freundlich model are higher than 0.95. Figure 6 presents the variation of Freundlich parameters as a function of molecular weight of the 16 PAH molecules.

4. Discussion

Phosphonic acids (PA) R-PO(OH)2 with R = vinyl (VPA), phenyl (PPA) were chosen in this study because they condense efficiently and selectively with the surface hydroxyl groups of the metal oxide, leading to the formation of a monolayer of unsaturated ligand and very stable ionocovalent P–O–M bonds [34]. VPA or PPA grafting was based on acid–base Brønsted surface reactions. This reaction was carried out at 100 °C for 4 h in aqueous solution, allowing the surface modifications of ZrO2 nanoparticles and thus allowing to obtain inorganic–organic hybrid nanomaterials (as a white powders). The commercial ZrO2 nanosuspension’s pH was set around 4–5 to avoid the support aggregation if it gets close to the isoelectric point (IEP) of zirconia (around 9) [24].
Depending on the yield of PPA added to this suspension, the pH can vary, as shown in Figure 7. By considering the PPA species present in the suspension under different pH conditions, the mechanism by which PPA acts can be understood.
The pKa values for PPA are 2.3 and 7.8 [36,37]. The deprotonation equilibria for PPA are as follows:
Ph–PO(OH)2 (aq) ↔ Ph–PO2 (OH) (aq) + H+(aq) (pKa1 = 2.3) [36]
Ph–PO2 (OH) (aq) ↔ Ph–PO3 (aq) + H+(aq) (pKa2 = 7.8) [36,37]
In the pH range of the resulting suspension, the dominant species of PPA are Ph-PO2 (OH). In this case, it can easily react with the ZrO2 surface, forming a bidentate coordination. As the isoelectric point of zirconia is around 9, the coordination mode can evolve from bidentate to tridentate coordination as found by XPS and NMR analysis. It is important to note that at a relatively high PPA amount (ratio 1/25), the dominant species is neutral PPA, which can lead to corrosion of ZrO2 and can form zirconium phosphonate as shown by NMR analysis.
Furthermore, it is worth noting that the substitution of the hydrophilic surface hydroxyls of the metal oxide with PA bearing an unsaturated hydrophobic function resulted in rapid precipitation of the synthesized nanohybrid for the materials with the highest organic charge. Additionally, this phenomenon was observed more for phenyl PA derivatives than vinyl PA because the aromatic function has a lower aqueous solubility than simple unsaturation.
The hydrophobic aspect of the aromatic function in the case of PPA allows having more affinity to the PAHs molecules compared to VPA as determined by the kinetic study carried out for the two series of nanomaterials containing PPA or VPA. Indeed, the kinetic study has shown that, unlike VPA, in the PPA grafting and thanks to the aromatic ring’s hydrophobicity, the resulting nanomaterials have more affinity to the PAH molecules, especially for the heavier ones. This grafting did not change the nanomaterials’ textural properties; therefore, no difference was observed in terms of adsorption capacity in both cases of VPA and PPA.
The properties of organic contaminants also affect the adsorption process. The hydrophobicity and hydrophilicity, the surface charge, and the functional groups of the pollutants can affect their adsorption behaviour. In the case of PAHs molecules, only the hydrophobicity and hydrophilicity could be a factor affecting the adsorption process. The octanol–water partition coefficient (Kow) is a useful physiochemical property to model the hydrophobicity and hydrophilicity of organic pollutants generally and PAH specifically. Figure 8 presents the variation of adsorption constant rate (k2) as a function of octanol–water partition coefficient (Kow) for the 16 PAH molecules in the case of ZrO2 or (ZrO2)50(PPA) tested at room temperature with a stirring rate of 300 rpm and an initial PAH concentration of 200 ppb. As shown in Figure 8 (and Table A2), the adsorption rate constant k2 is increasing in the same order as their octanol–water partition coefficients, which demonstrate that hydrophobic interaction can be a main mechanism for the adsorption process [38].
According to the variation of Freundlich parameters presented in Figure 6, the constant KF related to the adsorption capacity does not present notable variation between ZrO2 and (ZrO2)50(PPA), which is correlated with the kinetic study revealing a constant qe almost unchanged between grafted and ungrafted ZrO2. However, the values of nF are quite different between ZrO2 and (ZrO2)50(PPA), showing a different adsorption mechanism. Ungrafted ZrO2 is a polar species while PAHs are rather nonpolar molecules. In this case, the interaction that can form is probably the inductive force, also called Debye interaction. The values of nF increase as a function of the size of the PAH molecules. Indeed, the interaction is stronger since it increases with the polarizability, knowing that the polarizability rises with the size of the molecules, which clearly shows that the types of interaction between the ungrafted ZrO2 and the PAHs are indeed Debye interactions.
After surface functionalization by PPA, the polarity of the material surface decreases since phenyl groups are apolar. Thus, the interaction between the PAH molecules and the phenyl groups grafted on the surface of ZrO2 is an instantaneous dipole interaction, also called dispersion forces or London interaction. The dissociation energy of London interaction is less than 1 kcal.mol−1, which is coherent with nF having values less than 1 in the case of (ZrO2)50(PPA). This low energy interaction would facilitate the regeneration of the material for reuse.

5. Conclusions

In conclusion, a straightforward and eco-friendly approach for surface modification of commercial ZrO2 nanosuspension was developed to precisely tune the inorganic support’s surface functionalization without altering its structural and textural properties. Simple phosphonic acids bearing vinylic or phenylic substituents were used as a suitable anchoring group with a high affinity toward zirconia surface modification. The innovativedesigned nanocomposite materials were elaborated bearing different amounts of the organic grafting part to address the influence of their loading amount on the resulting nanohybrids materials’ structural and physicochemical properties. Further on, each material was thoroughly described with a wide variety of characterization: DRIFT, 1H 31P solid NMR, TGA, DSC, BET, XPS, and elemental analysis, to precisely define their different specific properties (high specific surface area, unique mesoporous networks, and various binding mode of the anchoring group). The synthesized nanomaterials were then used as nanoadsorbent to remediate the alarming PAHs pollution occurring worldwide and in numerous environmental matrices. Kinetics adsorption studies of each material towards a prepared standard aqueous solution containing trace levels of the 16 US EPA PAHs were carried out and showed that our systems are highly efficient. The application of various kinetic models to the adsorption data indicated that PAHs’ adsorption was best explained in terms of a pseudo-second-order kinetic mechanism. The adsorption kinetics parameters for the adsorption process have been investigated, indicating competitive adsorption between LMW and HMW PAHs, the higher affinity of ZrO2-PPA compared to ZrO2-VPA towards, and some influence of the organic charge towards the kinetic part of the adsorption.
The samples’ pretreatments were ensured by solid-phase extraction using commercial cartridges prior to their analyses by a novel and ultrasensitive GC-orbitrap-MS system that has recently shown high detection features toward PAHs trace levels in aqueous solution. This study lay the foundation of zirconia-based nanocomposite materials as highly efficient and eco-friendly adsorbents to address PAHs ubiquitous pollution.

Author Contributions

Conceptualization, N.B.O. and S.D.; methodology, N.B.O., S.D. and E.V.; validation, N.B.O., S.D. and E.V.; formal analysis, T.G., M.H., A.L. and F.L.; investigation, T.G., M.H., A.L. and F.L.; resources, N.B.O. and S.D.; writing—original draft preparation, M.H. and T.G.; writing—review and editing, M.H., N.B.O. and S.D.; supervision, F.L., N.B.O. and S.D.; project administration, N.B.O. and S.D.; funding acquisition, N.B.O. and S.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research and APC were funded by Zayed University (contract R19051).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in the study can be requested from the authors.

Acknowledgments

Authors thank Zayed University (R19051) for fundings and Christine LUCAS and Nesrine OUSLEATI (CP2M, Villeurbanne, France), M. Olivier BOYRON (CP2M, Villeurbanne, France), and Anass BENAYAD (CEA LETI, Grenoble, France) for solid-state NMR, TGA, and XPS experiments, respectively.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. Isotherms of (a) (ZrO2)x(VPA) and (b) (ZrO2)x(VPA) Samples (Zr/P ratio (x) = 25, 50, 100, 200) and bare ZrO2.
Figure A1. Isotherms of (a) (ZrO2)x(VPA) and (b) (ZrO2)x(VPA) Samples (Zr/P ratio (x) = 25, 50, 100, 200) and bare ZrO2.
Nanomaterials 11 00952 g0a1
Figure A2. Differential Scanning Calorimetry analysis of (a) (ZrO2)x(VPA) and (b) (ZrO2)x(VPA) Samples (Zr/P ratio (x) = 50, 75, 100, 200) in the range of 25–600 °C.
Figure A2. Differential Scanning Calorimetry analysis of (a) (ZrO2)x(VPA) and (b) (ZrO2)x(VPA) Samples (Zr/P ratio (x) = 50, 75, 100, 200) in the range of 25–600 °C.
Nanomaterials 11 00952 g0a2
Figure A3. Variation of organic weight loss in the range of 0–280 °C and 280–580 °C in thermogravimetric analysis as a function of the starting Zr/P ratio for (ZrO2)x(PPA) Samples.
Figure A3. Variation of organic weight loss in the range of 0–280 °C and 280–580 °C in thermogravimetric analysis as a function of the starting Zr/P ratio for (ZrO2)x(PPA) Samples.
Nanomaterials 11 00952 g0a3
Figure A4. X-Ray Photoelectron Spectroscopy study of (ZrO2)x(PAA) Samples (Zr/P ratio (x) = 25, 50, 100, 200).
Figure A4. X-Ray Photoelectron Spectroscopy study of (ZrO2)x(PAA) Samples (Zr/P ratio (x) = 25, 50, 100, 200).
Nanomaterials 11 00952 g0a4
Figure A5. Variation of the surface species content of ZrO2 as a function of the ratio used, according to X-ray Photoelectron Spectroscopy study of (ZrO2)x(PAA) Samples (Zr/P ratio (x) = 25, 50, 100, 200).
Figure A5. Variation of the surface species content of ZrO2 as a function of the ratio used, according to X-ray Photoelectron Spectroscopy study of (ZrO2)x(PAA) Samples (Zr/P ratio (x) = 25, 50, 100, 200).
Nanomaterials 11 00952 g0a5
Figure A6. (a) Evolution with time of benz(a)anthracene adsorption using (ZrO2)100(VPA) and (ZrO2)100(PPA) and (b) Pseudo-second-order kinetic plots for the adsorption of benz(a)anthracene using (ZrO2)100(VPA) and (ZrO2)100(PPA).
Figure A6. (a) Evolution with time of benz(a)anthracene adsorption using (ZrO2)100(VPA) and (ZrO2)100(PPA) and (b) Pseudo-second-order kinetic plots for the adsorption of benz(a)anthracene using (ZrO2)100(VPA) and (ZrO2)100(PPA).
Nanomaterials 11 00952 g0a6
Figure A7. (a) Evolution with time of benzo(a)pyrene adsorption using (ZrO2)x(PPA) Samples (Zr/P ratio (x) = 50, 100, 200) and (b) Pseudo-second-order kinetic plots for the adsorption of benzo(a)pyrene using (ZrO2)x(PPA) Samples (Zr/P ratio (x) = 50, 100, 200).
Figure A7. (a) Evolution with time of benzo(a)pyrene adsorption using (ZrO2)x(PPA) Samples (Zr/P ratio (x) = 50, 100, 200) and (b) Pseudo-second-order kinetic plots for the adsorption of benzo(a)pyrene using (ZrO2)x(PPA) Samples (Zr/P ratio (x) = 50, 100, 200).
Nanomaterials 11 00952 g0a7
Table A1. Textural properties for (ZrO2)x(VPA) and (ZrO2)x(PPA) nanohybrids.
Table A1. Textural properties for (ZrO2)x(VPA) and (ZrO2)x(PPA) nanohybrids.
AdsorbentSBET (m2 g−1)
ZrO258
(ZrO2)200(VPA)57
(ZrO2)100(VPA)51
(ZrO2)50(VPA)50
(ZrO2)25(VPA)56
(ZrO2)200(PPA)58
(ZrO2)100(PPA)57
(ZrO2)50(PPA)56
(ZrO2)25(PPA)53
Table A2. Adsorption capacities (qe) (µg.g−1), pseudo-second order rate constant (k2) (g.mg−1.min−1), and fit accuracy using a Pseudo-second-order Kinetic Plots for the Adsorption of 16 PAHs onto (ZrO2) and (ZrO2)50(PPA) nanocomposites.
Table A2. Adsorption capacities (qe) (µg.g−1), pseudo-second order rate constant (k2) (g.mg−1.min−1), and fit accuracy using a Pseudo-second-order Kinetic Plots for the Adsorption of 16 PAHs onto (ZrO2) and (ZrO2)50(PPA) nanocomposites.
PAHLog Kow *ZrO2(ZrO2)50(PPA)
qe
(mg.g−1)
k2
(g.mg−1 min−1)
R²qe
(mg.g−1)
k2
(g.mg−1 min−1)
R²
Naphtalene3.370.008735.80.99610.007738.50.9991
Acenaphtylene40.007355.10.99030.006687.40.9996
Acenaphtene3.920.008259.90.99540.0075102.80.9997
Fluorene4.180.008183.20.99570.0076179.00.9997
Anthracene4.540.0083139.20.99670.0079362.60.9998
Phenantrene4.570.0088264.40.99910.0087450.90.9999
Fluoranthene5.220.0088258.80.99910.0087330.70.9997
Pyrene5.180.0088255.00.99920.0086258.60.9995
Benz(a)anthracene5.910.0095192.40.99990.0096697.40.9999
Chrysene5.910.0094375.20.99990.0095576.10.9999
Benzo(a)pyrene5.910.0093253.10.99990.0100661.50.9999
Benzo(b)fluoranthene5.80.0094449.41.00000.0098869.90.9999
Benzo(k)fluoranthene60.0095279.11.00000.0095860.10.9999
Benzo(ghi)perylene6.50.0095161.10.99970.00981008.10.9999
Indeno(1.2.3-cd)pyrene6.50.0093396.90.99990.00981008.10.9999
Dibenz(a.h)anthracene6.750.0094493.40.99990.00991543.90.9999
* Kow values are obtained from reference [37].
Table A3. Freundlich parameters, and fit accuracy using the slope of the linear plot of lnqe versus lnCe for the adsorption of 16 PAHs onto (ZrO2) and (ZrO2)50(PPA) nanocomposites.
Table A3. Freundlich parameters, and fit accuracy using the slope of the linear plot of lnqe versus lnCe for the adsorption of 16 PAHs onto (ZrO2) and (ZrO2)50(PPA) nanocomposites.
PAHMolecular weight (g/mol)ZrO2(ZrO2)50(PPA)
R2KFnFR2KFnF
Naphtalene1280.9962.940.920.9614.680.33
Acenaphtylene1520.9862.390.870.9763.010.44
Acenaphtene1540.9893.000.880.9772.990.57
Fluorene1660.9652.980.630.9652.980.63
Anthracene1780.9955.400.710.9965.570.53
Phenantrene1780.9974.720.690.9714.210.59
Fluoranthene2020.9994.750.760.9465.420.52
Pyrene2020.9995.370.720.9866.250.53
Benz(a)anthracene2280.9992.061.280.9951.471.02
Chrysene2280.9911.461.610.9741.790.51
Benzo(a)pyrene2520.9961.781.350.9801.660.56
Benzo(b)fluoranthene2520.9813.041.530.9933.000.93
Benzo(k)fluoranthene2520.9991.122.190.9872.030.67
Benzo(ghi)perylene2760.9551.311.550.9851.520.75
Indeno(1.2.3-cd)pyrene2760.9991.541.930.9781.580.99
Dibenz(a.h)anthracene2780.9411.841.000.9412.290.48

References

  1. Grimmer, G. Environmental Carcinogens: Polycyclic Aromatic Hydrocarbons, 1st ed.; CRC Press: Boca Raton, FL, USA, 2018. [Google Scholar]
  2. Menzie, C.A.; Potocki, B.B.; Santodonato, J. Exposure to carcinogenic PAHs in the environment. Environ. Sci. Technol. 1992, 26, 1278–1284. [Google Scholar] [CrossRef]
  3. U.S. EPA (US Environmental Protection Agency). Polycyclic Aromatic Hydrocarbons (PAHs)-EPA Fact Sheet; US Environmental Protection Agency: Washington, DC, USA, 2008.
  4. Keith, L.H. The Source of U.S. EPA’s Sixteen PAH Priority Pollutants. Polycycl. Aromat. Compd. 2015, 35, 147–160. [Google Scholar] [CrossRef]
  5. U.S. EPA (US Environmental Protection Agency). Guidelines for Carcinogen Risk Assessment; US Environmental Protection Agency: Washington, DC, USA, 2005; 630/P-03/001F.
  6. Abdel-Shafy, H.I.; Mansour, M.S.M. A review on polycyclic aromatic hydrocarbons: Source, environmental impact, effect on human health and remediation. Egypt. J. Pet. 2016, 25, 107–123. [Google Scholar] [CrossRef] [Green Version]
  7. Rubio-Clemente, A.; Torres-Palma, R.A.; Peñuela, G.A. Removal of polycyclic aromatic hydrocarbons in aqueous environment by chemical treatments: A review. Sci. Total Environ. 2014, 478, 201–225. [Google Scholar] [CrossRef]
  8. Adeola, A.O.; Forbes, P.B.C. Advances in water treatment technologies for removal of polycyclic aromatic hydrocarbons: Existing concepts, emerging trends, and future prospects. Water Environ. Res. 2020. [Google Scholar] [CrossRef]
  9. Manousi, N.; Zachariadis, G.A. Recent advances in the extraction of polycyclic aromatic hydrocarbons from environmental samples. Molecules 2020, 25, 2182. [Google Scholar] [CrossRef]
  10. Lamichhane, S.; Bal Krishna, K.C.; Sarukkalige, R. Polycyclic aromatic hydrocarbons (PAHs) removal by sorption: A review. Chemosphere 2016, 148, 336–353. [Google Scholar] [CrossRef]
  11. Mahgoub, H.A. Nanoparticles used for extraction of polycyclic aromatic hydrocarbons. J. Chem. 2019, 2019. [Google Scholar] [CrossRef]
  12. Akbari, A.; Arsalani, N. Organic inaorganic condensed polyhedral oligomeric silsesquioxane-based nanohybrid: Synthesis, characterization and dye removal properties. Polym. Plast. Technol. Eng. 2016, 55, 1586–1594. [Google Scholar] [CrossRef]
  13. Farrokhzadeh, S.; Razmi, H.; Jannat, B. Application of marble powder as a potential green adsorbent for miniaturized solid phase extraction of polycyclic aromatic hydrocarbons from water samples. Sep. Sci. Technol. 2020, 55, 2737–2745. [Google Scholar] [CrossRef]
  14. Hedayati, M.S.; Li, L.Y. Removal of polycyclic aromatic hydrocarbons from aqueous media using modified clinoptilolite. J. Environ. Manag. 2020, 273, 111113. [Google Scholar] [CrossRef]
  15. Li, C.; Sun, M.; Ji, X.; Han, S.; Feng, J.; Guo, W.; Feng, J. Triazine-based organic polymers@SiO2 nanospheres for sensitive solid-phase microextraction of polycyclic aromatic hydrocarbons. J. Sep. Sci. 2020, 43, 622–630. [Google Scholar] [CrossRef]
  16. Gupta, P.K.; Khan, Z.H.; Solanski, P.R. Prospects of nanostructured ZrO2 as a point-of-care-diagnostic. In Recent Trends in Nanomaterials: Synthesis and Properties; Springer: Singapore, 2017; Volume 83, pp. 285–305. ISBN 978-981-10-3842-6. [Google Scholar]
  17. Petit, M.; Monot, J. Functionalization of Zirconium Oxide Surfaces. Chemistry of Organo-Hybrids; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2015; pp. 168–199. [Google Scholar] [CrossRef]
  18. Catauro, M.; Barrino, F.; Dal Poggetto, G.; Milazzo, M.; Blanco, I.; Vecchio Ciprioti, S. Structure, drug absorption, bioactive and anti bacterial properties of sol-gel SiO2/ZrO2 materials. Ceram. Int. 2020, 46, 29459–29465. [Google Scholar] [CrossRef]
  19. Mingalyov, P.G.; Lisichkin, G.V. Chemical modification of oxide surfaces with organophosphorus(V) acids and their esters. Russ. Chem. Rev. 2006, 75, 541. [Google Scholar] [CrossRef]
  20. Mutin, P.H.; Guerrero, G.; Vioux, A. Hybrid materials from organophosphorus coupling molecules. J. Mater. Chem. 2005, 15, 3761–3768. [Google Scholar] [CrossRef]
  21. Guerrero, G.; Alauzun, J.G.; Granier, M.; Laurencin, D.; Mutin, P.H. Phosphonate coupling molecules for the control of surface/interface properties and the synthesis of nanomaterials. Dalton Trans. 2013, 42, 12569–12585. [Google Scholar] [CrossRef]
  22. Lomoschitz, C.J.; Feichtenschlager, B.; Moszner, N.; Puchberger, M.; Müller, K.; Abele, M.; Kickelbick, G. Directing Alkyl Chain Ordering of Functional Phosphorus Coupling Agents on ZrO2. Langmuir 2011, 27, 3534–3540. [Google Scholar] [CrossRef] [PubMed]
  23. Bou Orm, N.; Trieu, Q.A.; Daniele, S. TiO2-Based Hybrid Nanocomposites Modified by Phosphonate Molecules as Selective PAH Adsorbents. Molecules 2018, 23, 3046. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Trieu, Q.A.; Pellet-Rostaing, S.; Arrachart, G.; Traore, Y.; Kimbel, S.; Daniele, S. Interfacial study of surface modified ZrO2 nanoparticles with thioctic acid for the selective recovery of palladium and gold from electronic industrial wastewater. Sep. Purif. Technol. 2020, 237, 116353. [Google Scholar] [CrossRef]
  25. Domínguez, I.; Javier Arrebola, F.; Martínez Vidal, J.L.; Garrido Frenich, A. Assement of wastewater pollution by gas chromatography and high resolution Orbitrap mass spectrometry. J. Chromatogr. A 2020, 1619, 460964. [Google Scholar] [CrossRef]
  26. Yang, L.; Wang, S.; Peng, X.; Zheng, M.; Yang, Y.; Xiao, K.; Liu, G. Gas chromatography-Orbitrap mass spectrometry screening of organic chemicals in fly ash samples from industrial sources and implications for understanding the formation mechanisms of unintentional persistent organic pollutants. Sci. Total Environ. 2019, 664, 107–115. [Google Scholar] [CrossRef]
  27. Ezzati, R. Derivation of Pseudo-First-Order, Pseudo-Second-Order and Modified Pseudo-First-Order rate equations from Langmuir and Freundlich isotherms for adsorption. Chem. Eng. J. 2020, 392, 123705. [Google Scholar] [CrossRef]
  28. Langmuir, I. Modelisation of adsorption. Phys. Rev. 1915, 6, 79–80. [Google Scholar]
  29. Freundlich, H.M.F. Over the adsorption in solution. J. Phys. Chem. 1906, 57, 385–471. [Google Scholar]
  30. Thommes, M.; Kaneko, K.; Neimark, A.V.; Olivier, J.P.; Rodriguez-Reinoso, F.; Rouquerol, J.; Sing, K.S.W. Physisorption of gases, with special reference to the evaluation of surface area and pore size distribution (IUPAC Technical Report. Pure Appl. Chem. 2015, 87, 1051–1069. [Google Scholar] [CrossRef] [Green Version]
  31. Al-Ghouti, M.A.; Da’ana, D.A. Guidelines for the use and interpretation of adsorption isotherm models: A review. J. Hazard. Mater. 2020, 393, 122383. [Google Scholar] [CrossRef] [PubMed]
  32. Lin-Vien, D.; Colthup, N.B.; Fateley, W.G.; Grasselli, J.G. Organophosphorus coumponds. In The Handbook of Infrared and Raman Characteristic Frequencies of Organic Molecules; Elsevier: Amsterdam, The Netherlands, 1991; pp. 263–276. ISBN 978-0-12-451160-6. [Google Scholar]
  33. Bou Orm, N. Nano-Materiaux Hybrides pour L’elimination de Micro-Polluants Organiques (HAP) dans les Effluents Aqueux. Ph.D. Thesis, University of Lyon, Lyon, France, 2012. [Google Scholar]
  34. Sheikh, J.A.; Bakhmutov, V.I.; Clearfield, A. Layered metal(IV) phosphonate materials: Solid-state 1H, 13C, 31P NMR spectra and NMR relaxation. Magn. Reson. Chem. 2018, 56, 276–284. [Google Scholar] [CrossRef]
  35. Wagstaffe, M.; Thomas, A.G.; Jackman, M.J.; Torres-Molina, M.; Syres, K.L.; Handrup, K. An Experimental Investigation of the Adsorption of a Phosphonic Acid on the Anatase TiO2 (101) Surface. J. Phys. Chem. C 2016, 120, 1693–1700. [Google Scholar] [CrossRef]
  36. Glover, C.F.; Williams, G. Inhibition of corrosion-driven organic coating delamination and filiform corrosion on iron by phenyl phosphonic acid. Prog. Org. Coat. 2017, 102, 44–52. [Google Scholar] [CrossRef]
  37. Glover, C.F.; Williams, G. Inhibition of Localized Corrosion of Hot Dip Galvanized Steel by Phenylphosphonic Acid. J. Electrochem. Soc. 2017, 164, C407–C417. [Google Scholar] [CrossRef] [Green Version]
  38. Sangster, J. Octanol-water partition coefficients of simple organic compounds. J. Phys. Chem. Ref. Data. 1989, 18, 1111–1229. [Google Scholar] [CrossRef]
Figure 1. BJH pore size distributions for (a) (ZrO2)x(vinylic (VPA)) and (b) (ZrO2)x(phenylic (PPA)) samples (Zr/P ratio (x) = 25, 50, 100, 200) and bare ZrO2.
Figure 1. BJH pore size distributions for (a) (ZrO2)x(vinylic (VPA)) and (b) (ZrO2)x(phenylic (PPA)) samples (Zr/P ratio (x) = 25, 50, 100, 200) and bare ZrO2.
Nanomaterials 11 00952 g001
Figure 2. Infrared spectroscopy (DRIFT) spectra of (a) (ZrO2)x(VPA) and (b) (ZrO2)x(PPA) samples (Zr/P ratio (x) = 25, 50, 100, 200) and bare ZrO2.
Figure 2. Infrared spectroscopy (DRIFT) spectra of (a) (ZrO2)x(VPA) and (b) (ZrO2)x(PPA) samples (Zr/P ratio (x) = 25, 50, 100, 200) and bare ZrO2.
Nanomaterials 11 00952 g002
Figure 3. 31P CP MAS NMR spectra of (a) (ZrO2)x(VPA) and (b) (ZrO2)x(PPA) samples (Zr/P ratio (x) = 25–200).
Figure 3. 31P CP MAS NMR spectra of (a) (ZrO2)x(VPA) and (b) (ZrO2)x(PPA) samples (Zr/P ratio (x) = 25–200).
Nanomaterials 11 00952 g003
Figure 4. Thermogravimetric analysis of (a) (ZrO2)x(VPA) and (b) (ZrO2)x(PPA) samples (Zr/P ratio (x) = 25, 50, 100, 200) and bare ZrO2 in the range of 25–900 °C.
Figure 4. Thermogravimetric analysis of (a) (ZrO2)x(VPA) and (b) (ZrO2)x(PPA) samples (Zr/P ratio (x) = 25, 50, 100, 200) and bare ZrO2 in the range of 25–900 °C.
Nanomaterials 11 00952 g004
Figure 5. Pseudo-second order rate constant of PAHs adsorption by (a) (ZrO2)x(VPA) and (b) (ZrO2)x(PPA) samples (Zr/P ratio (x) = 50, 100, 200).
Figure 5. Pseudo-second order rate constant of PAHs adsorption by (a) (ZrO2)x(VPA) and (b) (ZrO2)x(PPA) samples (Zr/P ratio (x) = 50, 100, 200).
Nanomaterials 11 00952 g005
Figure 6. Variation of Freundlich parameters: (a) nF and (b) KF as a function of molecular weight of the 16 PAH molecules.
Figure 6. Variation of Freundlich parameters: (a) nF and (b) KF as a function of molecular weight of the 16 PAH molecules.
Nanomaterials 11 00952 g006
Figure 7. pH variation as a function of the concentration of PPA in solution and predominance of phenylphosphonate species as a function of pH.
Figure 7. pH variation as a function of the concentration of PPA in solution and predominance of phenylphosphonate species as a function of pH.
Nanomaterials 11 00952 g007
Figure 8. Variation of adsorption constant rate (k2) as a function of octanol–water partition coefficient (Kow) for the 16 PAH molecules in the case of ZrO2 or (ZrO2)50(PPA).
Figure 8. Variation of adsorption constant rate (k2) as a function of octanol–water partition coefficient (Kow) for the 16 PAH molecules in the case of ZrO2 or (ZrO2)50(PPA).
Nanomaterials 11 00952 g008
Table 1. Quantitative analysis of (ZrO2)x(VPA) and (ZrO2)x(PPA) nanocomposites using XPS and elemental chemical analysis. (*** = not done)
Table 1. Quantitative analysis of (ZrO2)x(VPA) and (ZrO2)x(PPA) nanocomposites using XPS and elemental chemical analysis. (*** = not done)
Elemental Chemical AnalysisXPS
AdsorbentAtomic Ratio Zr/PMeasured P (wt%)PO3-PO2(OH)∑P
TheoreticalMeasured
(ZrO2)200(PPA)200******1.720.071.79
(ZrO2)100(PPA)1003000.082.390.723.11
(ZrO2)50(PPA)501220.193.631.415.04
(ZrO2)25(PPA)25******4.101.485.58
(ZrO2)100(VPA)100770.31*********
(ZrO2)50(VPA)50470.5*********
Table 2. Adsorption capacities (qe) (µg·g−1), pseudo-second order rate constant (k2) (g·mg−1·min−1), and fit accuracy using a pseudo-second-order kinetic plot for the adsorption of 16 polycyclic aromatic hydrocarbons (PAHs) onto (ZrO2)x(VPA) and (ZrO2)x(PPA) nanocomposites.
Table 2. Adsorption capacities (qe) (µg·g−1), pseudo-second order rate constant (k2) (g·mg−1·min−1), and fit accuracy using a pseudo-second-order kinetic plot for the adsorption of 16 polycyclic aromatic hydrocarbons (PAHs) onto (ZrO2)x(VPA) and (ZrO2)x(PPA) nanocomposites.
Adsorbent(ZrO2)xVPA(ZrO2)xPPA
Zr/P RatioPAH Pollutantqek2R2qek2R2
50Naphtalene15.550.997117.560.9976
10016.450.990216.8150.9941
20017.260.990218.520.9977
50Acenaphthylene14.7100.999217.280.9968
10015.780.994916.4210.9961
20017.270.997917.250.9903
50Acenaphthene17.5110.999618.8100.9952
10017.9110.997818.2260.9981
20018.770.997119.230.9954
50Fluorene17.4160.999618.8120.9966
10018.160.996518.3340.9991
20018.790.99791950.9973
50Phenanthrene17.2290.999518.4200.9979
10017.6130.999118570.9994
20018.4120.998618.670.9982
50Anthracene18.3330.999618.3130.9971
10017.2120.998218.4600.9997
20018.6120.997518.6120.9992
50Pyrene17.9670.999918.3460.9993
10018270.999818.51540.9999
20018.5250.999518.4200.9996
50Fluoranthene17.7510.999718.1420.9992
10017.8250.999818.41290.9998
20018.4210.999318.3170.9993
50Benz(a)anthracene19.170.997118.9590.9996
10018.970.9966191360.9995
20019.5200.999915.590.9905
50Chrysene19.460.996119.11480.9996
10019.290.994119810.9989
20019.7160.999913.5100.9976
50Benzo(b)fluoranthene18.7150.999218.81180.9999
10018.6260.999618.62280.9998
20018.8730.995316.9350.9971
50Benzo(k)fluoranthene19.150.994918.8980.9997
10019.180.990118.8540.9991
20019.3160.999914.3100.9944
50Benzo(a)pyrene19.650.994419.31390.9969
10019.490.993519.2780.9991
20019.9160.999914.3100.9944
50Dibenz(a,h)anththracene19.840.997519.5720.9983
10018.9100.998719.6530.9965
20020.390.999810.1100.9157
50Indeno(1,2,3-cd)pyrene18.180.997817.81400.9996
10017.790.999417.9890.9994
20019.7100.997915210.9946
50Benzo(ghi)perylene18.650.993117.91890.9997
10018.4150.992718.4380.9987
20018.9180.999811.5230.9618
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bou Orm, N.; Gréa, T.; Hamandi, M.; Lambert, A.; Lafay, F.; Vulliet, E.; Daniele, S. Design of Hybrid PAH Nanoadsorbents by Surface Functionalization of ZrO2 Nanoparticles with Phosphonic Acids. Nanomaterials 2021, 11, 952. https://doi.org/10.3390/nano11040952

AMA Style

Bou Orm N, Gréa T, Hamandi M, Lambert A, Lafay F, Vulliet E, Daniele S. Design of Hybrid PAH Nanoadsorbents by Surface Functionalization of ZrO2 Nanoparticles with Phosphonic Acids. Nanomaterials. 2021; 11(4):952. https://doi.org/10.3390/nano11040952

Chicago/Turabian Style

Bou Orm, Nadine, Thomas Gréa, Marwa Hamandi, Alexandre Lambert, Florent Lafay, Emmanuelle Vulliet, and Stéphane Daniele. 2021. "Design of Hybrid PAH Nanoadsorbents by Surface Functionalization of ZrO2 Nanoparticles with Phosphonic Acids" Nanomaterials 11, no. 4: 952. https://doi.org/10.3390/nano11040952

APA Style

Bou Orm, N., Gréa, T., Hamandi, M., Lambert, A., Lafay, F., Vulliet, E., & Daniele, S. (2021). Design of Hybrid PAH Nanoadsorbents by Surface Functionalization of ZrO2 Nanoparticles with Phosphonic Acids. Nanomaterials, 11(4), 952. https://doi.org/10.3390/nano11040952

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop