Enhanced Naproxen Elimination in Water by Catalytic Ozonation Based on NiO Films

: This study evaluates naproxen (NP) degradation efﬁciency by ozonation using nickel oxide ﬁlms ( NiO ( F ) ) as a catalyst. The NiO ﬁlms were synthesized by chemical vapor deposition and characterized by X-ray diffraction, scanning electron microscopy, atomic force microscopy and X-ray photoelectron spectroscopy. NP degradation was conducted for 5 min using 10 films of NiO ( F ) comparing against ozonation using 100 mg/L NiO powder in suspension ( NiO ( S ) ) and conventional ozonation ( O 3 -conv). Total organic carbon analysis demonstrated a mineralization degree of 12% with O 3 -conv, 35% with NiO as powder and 22% with NiO ( F ) after 60 min of reaction. The films of NiO ( F ) were sequentially used 4 times in ozonation demonstrating the stability of the synthesized material, as well as its properties as a catalyst for ozonation. A proposed modeling strategy using robust parametric identification techniques allows the comparison of NP decomposition pseudo-monomolecular reaction rates.


Introduction
The growing interest in catalytic ozonation is a consequence of their remarkable advantages with respect to conventional ozonation whose benefits are (a) enhancing the recalcitrant pollutants removal in water, (b) reducing reaction period, and (c) increasing the mineralization degree [1,2]. The presence of suspended metallic oxides particles in a catalyst breaks down the ozone molecule leading to radical generation (·OH mainly). These radicals are more reactive than molecular ozone due to their higher oxidation potential [3]. Today, different catalysts have been used in the ozonation for decomposing a wide variety of pollutants in water, including suspended metal oxides [4], supported metal oxides [5,6], activated carbon [7], zeolites and clay minerals [8] are the most extended ones due to their fast and efficient radicals generations.
Nickel oxide has shown outstanding results degrading herbicides [9], pharmaceutical compounds [10] and polycyclic aromatic hydrocarbons [11] in combination with ozone. One additional problem associated with heterogeneous catalytic ozonation reactions is the lack of accurate models describing the complete process. The understanding of the process conditions and the expected products can contribute to the operation control and quality of the contaminant degradation results [32]. Classical modeling techniques of ozonation include the least square method to obtain an accurate curve that fits the experimental data better. However, high variability presented in such a complex chemical system (including serial and parallel reactions) can lead to inaccurate models, which may limit their applicability for real treatments. To overcome these problems, a robust differentiation technique was implemented.
Catalytic ozonation with NiO in suspension represents an effective alternative to eliminate NP from water [33], but the recovering and eventual reuse of suspended catalysts needs to be resolved. Therefore, it is feasible to propose catalytic ozonation based on NiO thin films aimed to decompose NP in aqueous solutions. Here, NiO films (NiO (F) ) were supported on glass substrates by the CVD method. A comparison was established between conventional and catalytic ozonation using NiO in powder and films.

NiO (F) Characterization
The CVD technique used to deposit NiO allows the obtaining of a film with a uniform gray color. The structural analysis was complemented with the X-ray diffraction (XRD) study. Figure 2 shows the XRD pattern of NiO (F) . The analyzed spectrum presents five characteristic peaks at 2θ = 37.52, 43.58, 63.14, 75.7 and 79.84 • . These signals agree with the commercial NiO (S) , corresponding to the (111), (200), (220), (311) and (222) planes, respectively. These signals are indexed to the cubic phase of NiO (Joint Committee on Powder Diffraction Standards No. 78-0423, 73-1519) [34]. The presence of only well-defined five peaks indicates the high purity of the NiO (F) . According to the XRD pattern, the NiO (F) obtained by CVD has preferential orientation in 111 direction and the orientation does not change after ozonation. Furthermore, the films present compressive stress due to the slight slippage of the peaks. If we consider (111) planes, NiO (S) presents the peak at position 37.20, while O 3 -NiO (F) at 37.36 and NiO (F) at 37.52. The results indicate that NiO (F) has 0.86% compressive stress and O 3 -NiO (F) after the ozonation, there is a 0.43% relaxation compared to NiO (S) powders. Additionally, the X-ray diffraction analysis in the software Match3 enables determination of a density of 6.79 g/cm 3 for NiO (F) . The crystallinity stability was demonstrated in the NiO (F) XRD pattern after oxidation with ozone (O 3 -NiO (F) ). The well-defined five peaks detected in the sample without ozonation were observed after the treatment. This pattern confirmed the stability of the film after ozonation. The surface morphology and thickness of NiO (F) were evaluated by Scanning Electron Microscopy (SEM), Figure 3a. The image showed the aggregates of NiO on the glass substrates, revealing a dense and continuous coating. The optical image captured at 200,000× demonstrates the presence of NiO aggregates with an equivalent diameter smaller than 100 nm. These aggregates are conformed by little spherical particles (Figure 3a). A similar surface growth was reported in the synthesis of NiO films by electro-deposition on ITO coated [35]. The cross-sectional image of the film (insert Figure 3a) demonstrated the deposition of NiO on the glass slide was around 206 nm of thickness. This value was confirmed by spectroscopic ellipsometry. The average thickness was 240.97 ± 55.56 nm (n = 3, an average of five points by sample). After ozonation, NiO (F) (Figure 3b) showed apparent cracks. Lin et al. [35] in the synthesis of NiO (F) by galvanostatic electro-deposition method, observed cracks on the surface. This fact was attributed to the relaxation of the stress in the film. This behavior agrees with the XRD results, moreover, the film stability was demonstrated using the constant diminution of total organic carbon (TOC) when the reused films were used as the catalyst. Figure 4 shows the Atomic force microscopy (AFM) images of surface morphology of NiO (F) (Figure 4a) and O 3 -NiO (F) (Figure 4b). In these figures, no cracks are observed after ozonation. The values of average roughness (Ra) and the root mean square roughness (Rrms) for NiO (F) were 4.66 nm and 5.89 nm, respectively and Ra was 4.82 nm and Rrms was 6.14 nm for O 3 -NiO (F) . These results unlike that observed in SEM analysis, suggest that the surface is not affected during the ozonation.
X-ray photoelectron spectroscopy (XPS) results for NiO (F) before and after ozonation are shown in Figure 5. The Ni2p region in the XPS spectrum displays characteristic signals attributed to NiO, Figure 5a. Ozone did not change the spectrum in the Ni region which is similar to previous DRX results. Furthermore, the comparison NiO (F) and NiO (S) showed similar spectra. The O1s region (Figure 5b) is decomposed in two main signals for NiO (F) and NiO (S) : the first at 529.3 eV corresponding to oxygen atoms lattice of NiO and the second around 531 eV assigned to defective oxygen, for instance oxygen atoms due to nickel vacancies [9]. When the ozonation was carried out, NiO (F) displays two additional peaks (at 531.8 eV (-CO) and 532.3 eV (-COH)) which are attributed to oxygen coordinated to carbon atoms by adsorbed byproducts generating to NP decomposition.    Figure 6a shows the NP degradation profiles (obtained by high efficiency liquid chromatography) for conventional (O 3 -conv) and catalytic ozonation with NiO (F) and NiO (S) . The NP signal was only detected along the first five minutes of the reaction. In the presence of NiO (F) and NiO (S) , the catalyst did not modify the NP degradation profile in comparison with O 3 -conv. The determined NP degradation time was faster than the other advanced oxidation process [31,36,37]. The apparent absence of changes in the degradation profile of the parent compound agrees with previous results [38,39]. This fact was attributed to: (1) ozone fast reaction with the initial contaminant (NP) and (2) generation of oxidant species (ozone decomposition with catalyst) and reaction with NP ozonation byproducts.

Evaluation of NiO (F) in the NP Removal Using Ozone
The pH changed during the ozonation from the initial NP solution which had a pH of 5.0 ± 0.5 to the final pH of 2.5 ± 0.5 after the treatment (60 min). Under these experimental conditions (acid pH), the direct ozonation mechanism was favored and ozone decomposition (formation of oxidant species) was not enforced. On the other hand, the oxalic acid accumulated in ozonation as the main final product which yields the decreasing of pH. This result coincides with the reported results by other researchers and our previous results [38,40].
Few articles have studied NP degradation by catalytic ozonation varying pH. The study proposed in [41] reported NP (50 mg/L) degradation with 0.48 mg/s of ozone at pH of 5, 6, 7, 8 and 9. The treatment at pH = 9 allowed near 100% of NP degradation after 4 min, while only 30% of NP was eliminated with pH = 5 during the same reaction time. This difference was explained considering that under basic pH, direct and indirect mechanisms occur simultaneously. In the case of catalytic ozonation of NP, [33] demonstrated that 1 g/L of TiO 2 yields a higher mineralization degree at pH of 5 and it decreased near 50% at pH 3. Moreover, the same authors proposed two periods of mineralization; the first period, named the fast, corresponding to the fast reaction of easy oxidable compounds; in this period pH plays a main role (while at greater pH = 7 increased the removal of TOC). The second mineralization stage includes the reaction of refractory byproducts and it was independent of pH. Notice that some of these results agree with the outcomes attained in this study.
The NiO effect can be observed in the oxalic acid profile over time ( Figure 6b). This organic acid is a recalcitrant compound to O 3 -conv [42,43]. Oxalic acid is a main final product formed as a consequence of the aromatic ring outbreak or loss of acid naproxen substituent in the reaction, (mechanism reported by Jallouli et al. [44]). The O 3 -NiO (F) duplicates the oxalic acid concentration in comparison with O 3 -conv during the first min of reaction. The oxalic acid concentration is similar after 10 min in catalytic and conventional ozonation. The differences in the oxalic acid profile in both ozonation systems justify the proposed inclusion of several oxidant species in the catalytic ozonation. According to [38], the combination of ozone and NiO produced hydroxyl radicals (·OH) which favors the elimination of refractory organic compounds due to their oxidation potential (2.8 eV) [45].
Oxalic acid as a NP final product can be generated by two main routes: (1) the breakout of the propionic acid substituent of NP and the release of a 3-carbon fragment, the subsequent 1-carbon addition form the malic acid [44] by the oxidant species in the reaction. This byproduct is susceptible to oxidation with ozone and hydroxyl radicals to produce oxalic acid (a recalcitrant final product) (2). The electrospray ionization mass spectrometry (ESI-Ms) analysis of NP degradation effluents at 5 min showed the formation of aromatic byproducts. The molecules proposed for these ions m/z of 237.11, 217.09, 177.66 and 149.06 only present one aromatic ring.
The fragments obtained in the breakout and degradation of the aromatic compounds can be degraded to oxalic acid. The increase of the oxalic acid concentration in the first two minutes of ozonation is possible due to the first route described above, where the break of aromatic rings is not needed. The subsequent decrease and slight increment of this compound concentration might be related to the aromatic byproduct's degradation after 5 min which was confirmed by the HPLC and ESI-Ms studies. As an example (Figure 7), a signal of a non-identified intermediate byproduct detected by HPLC (12.9 min of retention time) disappeared after 10 min. Its maximal concentration was determined at 2 min.   In a previous article regarding NP degradation with different proportions of water and ethanol as solvents [38], the studies of ESI-Ms and ESI-mass-mass (ESI-Ms-Ms) spectrometry were done in negative mode. Different aromatic byproducts such as

acetic acid and (d) 4-Methylphenylacetic acid
were detected in the first 5 min of the reaction. However, these byproducts were not detected at 60 min of treatment (conventional and catalytic ozonation). Similar results were obtained with the NiO films for the NP elimination. These results confirm that all these intermediate compounds are close to be eliminated within the time window of the ozonation reaction (This fact is confirmed with the ESI-Ms-Ms analysis in samples taken after 60 min of ozonation). What is more important is that the decomposition of these compounds yields to the formation of organic acids (such as oxalic acid), aldehydes and carbon dioxide. All these compounds are considered to be more biodegradable.
The global efficiency of water treatments can be determined by the TOC variation. The role of NiO (F) was estimated also by TOC analysis (Figure 8). In the O 3 -conv, TOC removal after 60 min of the reaction was 12.3%, while introducing 10 films of NiO increases this removal to 22%. The presence of 100 mg/L of NiO (S) approximately achieved 35% of TOC removal, indicating that the solid catalyst was a better system for the NP elimination.
The low and incomplete NP mineralization has been described as a regular condition in catalytic ozonation. Rosal et al. reported near 40% of TOC removal in O 3 -conv and around 55% using 1000 mg/L of TiO 2 Degussa P25 after 60 min of reaction. Therefore, the use of TiO 2 as a catalyst increased the TOC removal about 15% [33]. This percentage is similar to the observed in this work with 100 mg/L of NiO in suspension.
The use of NiO (F) showed less TOC removal compared to NiO (S) . In this study, to compare correctly the catalytic effectiveness of both catalytic cases, it is worth taking into account what decomposed mass of NP corresponds to the mass of the catalyst in the reaction, it means a class of NP decomposition yield. A rough calculation of the catalyst on the glass slide showed a mass of 1.2 mg per film and about 12 mg over the 10 used films. The calculated amount of catalyst in the treatment is 10 times less than that used in suspension (40 mg). Corresponding, the catalytic activity of NiO (F) is 1.66 mg of the decomposed NP and only 0.5 mg for NiO (S) at 1 mg of catalyst. Therefore, the catalytic activity of NiO (F) is 3.3 times higher than NiO (S) . This result shows the possibility of using lower concentrations of the catalyst when depositing it on thin films, in addition to facilitate its recovery after ozonation.
The stability of NiO (F) was studied by TOC removal after 4 consecutive cycles of ozonation, Figure 8. In all cases, NiO (F) removed similar TOC percent after each cycle (22.2% ± 2.9% of TOC removal). Such a result indicates that ozone-NP and ozone-byproducts reactions did not inhibit the activity of NiO (F) after 4 cycles. The use of films as catalysts in the ozonation is scarcely described yet. Guzman et al. (2020) [19], in the 4-chlorophenol and 4-phenolsulfonic acid degradation by ozone (120 min) with six ceria oxide films (26.4 × 76.2 mm) demonstrated a similar TOC removal in five cycles of treatment (about 44% for both compounds). The results confirm the stability of the film in the catalytic ozonation.  Table 1 summarizes all parameters used in the Super-Twisting algorithm (STA) to estimate the kinetics of the NP degradation by ozonation. The robustness of STA produces similar values for the needed parameters for tuning the differentiator. Figure 9 shows the time variation of signals produce by the modeling process that yields the kinetics and the approximate model of NP degradation by conventional and catalytic ozonations. In the first step, data are collected from the HPLC. These data are marked with a blue circle in all subfigures. Then, the black continuous line is obtained by an interpolation process (spline functions as an approximation basis). This step is done to have enough data to feed STA applied as a differentiator. The accuracy of the STA is related to the quality of the collected data. The derivative is presented in the second graph. During the first 40 s of simulation, STA presents an overshoot during the adaptation phase, which means it is the time when the value x 1 (t) − C NP (t) is reaching the zero value. After this time, the derivative is exact, and one may calculate the maximum rate of NP degradation. This rate is marked with a red circle in each graph and it is reported in Table 1. This process was repeated three times to compare the behavior of the ozonation cases studied in this manuscript. Notice that all reaction rates constants are not varying significantly. However, their tendencies confirm what one may detect as a critical outcome of this and similar studies on supported catalysts over solid films: films of catalyst may be more efficient than conventional (no catalytic), but less than using the powder catalyst. This condition can be clearly justified by the mass transfer effect which is not limited in the O 3 -NiO (S) case, but it could be an issue for the O 3 -NiO (F) process. Nevertheless, the process efficiency in terms of the contaminants and byproducts removals concerning the catalyst mass is relevantly higher for the case of thin films. Moreover, the stability of such films regarding reusing experiments is a promising outcome attained here.

Materials and Methods
This section summarizes all materials and methods used in the catalytic ozonation of NP using supported catalysts forming thin films. Also, the kinetic characterization of catalytic ozonation of NP is proposed using a simplified model based on approximate pseudo-first-order dynamics.

Synthesis of Nickel Oxide Films
The synthesis of NiO (F) was carried out in a quartz tube within a CVD set-up reported elsewhere. Nickel (II) acetylacetonate Ni(acac) 2 (Sigma-Aldrich, Mexico City, Mexico 99.0%) was used as metalorganic precursor of the corresponding metallic oxide. The conditions for NiO deposition were: (a) 200 • C of Ni(acac) 2 sublimation temperature, (b) 550 • C of substrate temperature, (c) mixture of Ar-O 2 (60sccm-60sccm). With these conditions, the allowed pressure in the equipment was 10 mTorr. This parameter was kept constant along the reaction period. The deposition time was 15 min.
The NiO (F) grew up on glass plates of 2.5 × 3.0 cm as the substrate. The substrate was washed previously to the deposition with sonication cycles of 15 min using different solvents: acetone, ethanol, and distillate water. The films were used in ozonation without any additional techniques.

NiO (F) Characterization
Characterization of NiO (F) was made with the following techniques: XRD, XPS, AFM (Brunker Zeiss Bioscope catalyst), spectroscopic ellipsometry and SEM (JEOL JSM 6701F, Peabody, Massachusetts, USA, 5 kV). These techniques allowed to study the morphology and crystal structure of the obtained NiO (F) . The SEM determined the thickness of the NiO film. This parameter was compared with the result obtained by spectroscopic ellipsometry in the HORIBA's UVISEL equipment where the visible-UV spectra were acquired from a range of 1.5-5.5 eV.
XRD analysis was carried out in the Bruker D8 advance equipment with copper radiation operating at 35 kV and 25 mA (CuKα, λ = 1.54 Å) in 2θ between 20-90 • . The XRD pattern obtained from NiO (F) was compared with NiO (S) with the Bruker D8 AXS equipment using the same operating conditions. XPS (Thermo Fisher Scientific, Waltham, Massachusetts, USA) was applied to determine changes in chemical binding types (Al Kα X-ray monochromatic source, 1487 eV; charge correction C1s: 284.6 eV).
AFM and NP-10 probes tips (Bruker, www.brukerafmprobes.com) with V-shaped cantilever and 20 nm radius pyramidal-geometry were used. The measurements were done in Tapping mode. The software NanoScope Analysis v1.4 (Bruker, Billerica, Massachusetts, USA, 2018) was used to calculate the average roughness which determines the deviation in height. Root mean square roughness represents the standard deviation of surface heights from AFM images. The thickness of NiO (F) was obtained by spectroscopic ellipsometry.

Catalytic Ozonation
The stability condition of NiO (F) was studied in catalytic ozonation of NP solutions (Sigma-Aldrich, 99.8%; 20 mg/L dissolved in distillate water). The initial pH of the solution was 5.0 ± 0.5, this parameter was not controlled during the ozonation. In O 3 -NiO (F) treatments, ten coated slides were assembled on a Teflon sheet. They were put in a glass bubbled column reactor. This configuration allowed the contact of NiO (F) with NP solution (400 mL) and ozone. The mixture O 3 -O 2 (flow: 0.5 L/min and [O 3 ] = 5.5 ± 0.5 mg/L, HTU500 G ozone generator corona discharge type-AZCO Industries) bubbled through a ceramic porous filter located at the bottom of the self-designed reactor (500 mL). NP degradation by O 3 -NiO (F) was compared with O 3 -NiO (S) using 100 mg/L of NiO nanopowder <50 nm Sigma-Aldrich, 99.0%.

Analytical Methods
The monitoring of NP and oxalic acid (the main byproduct formed by the NP decomposition by conventional ozonation) by HPLC was carried out in the Perkin Elmer Flexar equipment. The column for the NP analysis was a Platinum C18 Alltech (250 mm × 4.6 mm) and 0.3 mL/min of water:acetonitrile 50:50 with a pH of 2.5 (phosphoric acid) as mobile phase. The separation followed at 240 nm. The byproducts detection was made using a Prevail TM Alltech-Grace column (150 × 4.6 mm and 5 µm). The mobile phase was a 0.25 mM of KH 2 PO 4 Buffer pH of 2.5 (phosphoric acid), the determination followed with 1.0 mL/min at 210 nm. The mineralization degree was calculated using TOC data obtained by the Torch IQOQ Teledyne Tekmar equipment by direct injection. All samples were previously filtered with membranes of 0.2 µm.

Mathematical Modeling and Parameter Estimation
To compare the NP removing efficiency enforced by the three studied methods, the NP ozonation rate constants were determined. For this study, it is used as a simplified approximation of the kinetic model based on pseudo-first-order dynamics. Then, with the data obtained from the decomposition of NP by ozonation, one can propose a set of finite parameters that correspond to the reaction rate constants. The following kinetic model considers the interaction between the measurable components in the reaction (contaminant, byproducts and ozone).
It is known that the reaction of either molecular ozone or hydroxyl radicals with contaminants satisfies a second-order kinetics, i.e.: where k NP,O 3 , min −1 is the reaction rate constant between ozone in the liquid phase O 3,l , mole/L and NP while k NP,OH is the reaction rate constant between hydroxide radicals OH, mole/L and NP. Nevertheless, due to the interaction of ozone with the catalyst leading to the corresponding formation of hydroxyl radicals, the explicit quantification of O 3 and OH independently is a complex task in the presence of NiO films. Therefore, in this case, we propose an alternative method to characterize the decomposition of NP. Consider that k(t) = k O 3 O 3 (t) + k OH OH(t) is the pseudo-mono molecular time dependent reaction rate parameter. Now, if one considers that k NP = max t≥0 k(t) it is feasible to propose that the following mathematical structure is given by: where c NP represents the concentration of the contaminant in time during the ozonation and k NP is the kinetic constant associated with the decomposition rate of NP. Using the differential equations representations of chemical processes and considering that advanced oxidation processes can be sufficiently approximated by pseudo-first-order dynamics. The inherent problem is the available information along with heterogeneous catalytic ozonation reaction. Usually, available measurements along the reaction are concentrations of just a few chemical species including the initial contaminant as well as its byproducts [46]. However, this information is not enough to get an accurate estimation of the first-order reaction rates. Commonly, it is necessary to get estimates of the contaminants' decomposition velocity. With the aim of getting a precise estimation of such decomposition velocity, a robust differentiation technique can be implemented.
Numerical differentiation is a well-studied topic in systems theory. The most common differentiation technique is based on the Euler approximation [47]. However, this method is highly affected by noises in the measurement. This manuscript proposes the application of the STA, which is a robust exact differentiator [48] based on the sliding mode theoretical results. The advantages of this algorithm are its robustness against parametric uncertainties that are common in the chemical process, well-posed discretization equivalents and finite-time convergence. Therefore, by means of the STA, the reaction kinetics are obtained. With the concentration and the information provided by the STA, a least square regression is performed to obtain a more accurate model representation of the heterogeneous catalytic ozonation.
STA [48] can determine the time derivative of a signal even in the presence of certain class of noises and model uncertainties. The algorithm is given by the following set of differential equations In Equation (3), d c NP is the actual time derivative of the concentration c NP (t), x 1 and x 2 are internal variables, k 1 and k 2 are positive constants that define the rate of convergence and the quality of estimation of the algorithm. The function sign(·) is defined as Notice that d dt c NP (t) = d c NP (t) when the algorithm has converged to the real concentration c NP (t), the estimation error defined as x 1 (t) − c NP (t) is equal zero. Notice the necessity of applying the STA differentiator as many times as the number of compounds identified along with the heterogeneous catalytic ozonation characterization: one for the degradation of NP in the three experimental setups: The complete procedure to estimate the reaction rates is summarized in the Algorithm 1. An interpolation is made to enhance the data available for analysis −→ c NP ; 3.
The STA is explicitly discretized using Euler method (See the work in [47]); 4.
The value of d NP (t) is the time derivative of c NP ; 5.
The maximum value of d NP (t) corresponds to the maximum reaction rate; 6.
A LS method approximates the reaction model. 7.
The approximated model is compared with the original data end /* The for cycle runs three times, one for each procedure (O 3 -conv, O 3 -NiO (S) and O 3 -NiO (F) . */ /* The LS method can be applied using the software Matlab R */ For further references, the theoretical background of the STA working as a differentiator can be found in [47,48].

Conclusions
CVD method allowed the synthesis of thin NiO films (240 nm thick) with relevant catalytic characteristics in NP degradation by ozone. NP elimination was carried out in five minutes of reaction with and without NiO (powder and film), due to its fast reaction with molecular ozone. The byproduct analysis and TOC removal showed the catalytic activity of the NiO in the ozonation. The treatments with 100 mg/L of NiO (S) generated a TOC removal of 35% while in the presence of 10 films of NiO, TOC removal increased 10% in comparison with O 3 -conv. However, NiO mass supported on the films was close to a tenth of the NiO used in suspension. The reuse of films (4 cycles) attained similar TOC removal in each cycle (22 ± 2.9%). This fact showed the stability of the catalyst films. Moreover, the film characterization after ozonation demonstrated some changes in the surface morphology but XRD and XPS patterns confirmed the crystalline stability. The results obtained suggest the possibility of using the deposition of catalyst in the film form in the catalytic ozonation of recalcitrant pollutants with the reduction of several detected disadvantages of the suspended catalyst (elimination and reuse). This study also presented a simplified model characterizing the ozonation kinetics of NP. A novel method yields to estimate a pseudo-monomolecular reaction rate constant which established a comparison between the proposed ozonation of NP with (suspension and films) and without catalysts. This model seems to offer a primary method to explain the reaction kinetics of ozone and NP in the presence of NiO (F) as a catalyst.  Acknowledgments: The authors kindly acknowledge to Instituto Politécnico Nacional, Centro de Nanociencias y Micro y Nanotecnologías-IPN, Consejo Nacional de Ciencia y Tecnología (Posdoctoral studies of Claudia M. Aguilar Melo) and PhD. Denisse Fabiola González Ramírez for economical supporting and contribution in this work.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript: NP Naproxen CVD Chemical vapor deposition O 3 -conv Conventional ozonation NiO (F) NiO film O 3 -NiO (S) Catalytic ozonation with suspension NiO O 3 -NiO (F) Catalytic ozonation with film NiO