Carbamazepine (CBZ), 5H-dibenzo[b,f]azepine-5-carboxamide (Figure 1), is an
Figure 1. Chemical structure and emission spectra of CBZ: (a) CBZ aqueous solution (5.0 µg∙mL−1) without irradiation; (b) CBZ aqueous solution (5.0 µg∙mL−1) + NaOH 0.1 M without irradiation; (c) solution a irradiated and (d) solution b irradiated.
anticonvulsant and mood-stabilizing drug used primarily in epilepsy treatment and bipolar disorder, as well as trigeminal neuralgia in old age patients . It is also used for a variety of indications, including attention-deficit hyperactivity disorder, schizophrenia, phantom limb syndrome, paroxysmal extreme pain and post-traumatic stress disorders  . CBZ is a first-choice anticonvulsant because of its relatively low psychological toxicity and the rarity of serious adverse effects .
Several analytical methods have been applied to CBZ determination or its photodegradation products in different matrices; i.e. chromatographic, spectrophotometric, electrochemical and capillary electrophoresis methods    . CBZ is not fluorescent in solution but can be converted into a fluorescent compound using UV irradiation, thus fluorimetric methods have been developed in a Nylon surface or through the formation of fluorescent derivatives by transformation assisted with different oxidant agents  - .
Flow injection analysis (FIA) systems are recognized as excellent automated methods since they show advantages in terms of precision along with significant reduction in reagents and samples consumption, thus minimizing wastes and sampling times; unlike conventional batch mode. Furthermore, FIA schemes are easy to couple with different optical detectors  .
Optimization of the experimental conditions of an FIA system using the ‘‘one variable at a time’’ methodology require various experiments to modify each factor versus the response studied disregarding the influence that each variable may exert on another and even then, the optimum may not be achieved. The multivariate statistical technique has arisen as an important tool for optimizing analytical processes. One of the most important within them is the response surface methodology (RSM). This method involves various statistical and mathematical techniques based on fitting a polynomial equation and symmetric models to the experimental data to describe the behavior of the independent variables. The central composite design (CCD) matrix is a type of RSM characterized by an embedded factorial or fractional factorial design with center points that is augmented with a group of “star points” that allow estimation of curvature. If the distance from the center of the design space to a factorial point is ±1 unit for each factor, the distance from the center of the design space to a star point is |α| > 1. The precise value of α depends on certain properties desired for the design and on the number of factors involved   .
The aim of the present work was to develop a method for CBZ determination using photoinduced fluorescence which is based on the photochemical transformation by UV irradiation. In order to make the proposed methodology suitable for the pharmaceutical analysis, experimental parameters affecting the sensitivity and sampling frequency were studied and optimized by a multivariate experimental scheme. The results were compared with those given by the official methods.
2.1. Instrumentation and Software
A Shimadzu spectrofluorometer model RF-5301PC (Kyoto, Japan), equipped with a Xenon discharge lamp was used for the fluorescent measurements and a 120 µL volume flow-cell unit for the continuous measurements were used to record transient signals. The CBZ fluorescence measurements were carried out operating the spectrofluorometer in the time-course mode (transient signals; λex = 251 nm; λem = 434 nm, slits widths = 3/3 nm). In batch studies, quartz cells with an optical path length of 1 cm were used.
The propulsion system consisted of two peristaltic pumps (Gilson Minipuls 3) with Polytetrafluoroethylene (PTFE) tubing. A Rheodyne model 5041 six-port two-way rotary valve (Rohnert Park, CA) was employed.
A home-made continuous photochemical reactor was made as follows: a UV-C source that emits radiation of 100 to 280 nm was used. The lamp (15 W high pressure Hg, PHILIPS) ignited with a suitable starter and chock, and was surrounded by a 2.5 m PTFE tubing (internal diameter = 0.75 mm and outer diameter = 1.0 mm), and covered by an aluminum foil, to permit the maximum reflectance of UV light.
In order to minimize interferences by matrix effect, a mini-column (mC) of multiwalled carbon nanotubes (MWCNTs) as adsorbent material was used. It was previously activated following the procedure described by Acosta et al. .
A combined glass electrode and a pH meter Orion Expandable Ion Analyzer (Orion Research, Cambridge, MA, USA, Model EA940) was used for pH adjustments.
In order to validate the methodology, a chromatographic analysis was performed using an HPLC system equipped with a Gilson 322 pump controller, a Rheodyne injector with 20 µL loop, a Gilson 170 UV-vis diode array detector (Buenos Aires, Argentina). The separative column used was a Phenomenex Luna 5µ CN (250 × 4.6 mm) (Buenos Aires, Argentina). The mobile phase used was a mixture (85:12:3) of water (milli Q): methanol (Sintorgan Bs. As., Argentina): tetrahydrofuran (JB Baker, New Jersey, USA), to which 0.22 mL of formic acid (Anedra, Buenos Aires, Argentina) and 0.5 mL of triethylamine (Merck, Darmstadt, Germany) were subsequently added.
The Stat-Ease Design-Expert Trial Version 11.0 was used for design and analysis of experimental runs. The predicted and adjusted R2 coefficient values were the criteria for the quality of the polynomial model equation was determined by an F-test.
CBZ (C15H12N2O, MW 236.27 g∙mol−1) was supplied by Parafarm Droguería Saporiti S.A.C.I.F.I.A. rev. 07, Bs. As., Argentina. In the optimization study, the pH was adjusted using 1 mol∙L−1∙NaOH (Mallinckrodt Chemical Works, New York, Los Angeles, St. Louis, USA) and 0.1 mol∙L−1 HCl (Merck, Darmstadt, Germany) solutions until the desired pH value. Multiwalled carbon nanotubes (MWCNTs) (diam = 110 - 170 nm, length = 5.9 micron, 90 + %) were purchased from Sigma Chemical Co. (St. Louis, MO, USA). All other reagents employed in this study were of analytical grade.
2.3. Stock Solutions
CBZ stock solutions containing 1 mg∙mL−1 were weekly prepared by dissolution of suitable amounts of the drug in 50% v/v ethanol water mixture (Merck Darmstadt, Germany). The standard working solutions of 0.1 mg∙mL−1 were prepared daily by stepwise dilutions with doubly distilled water.
2.4. Sample Preparation
A syrup sample volume of 0.25 mL (Tegretol® Novartis S.A., Bs. As., Argentina, batch H5595, labeled as 2%) was transferred to 50 mL volumetric flask, and diluted with 50% v/v ethanol water mixture to volume. The CBZ final concentration was about 0.1 mg∙mL−1. This solution was further diluted with doubly distilled water and was adjusted to pH 11 with NaOH (1 mol∙L−1) to obtain suitable experimental conditions for CBZ quantification.
Ten tablets of Tegretol® (Novartis S.A., Bs. As., Argentina, batch 405731) labeled to contain 200 mg CBZ per tablet, were weighted and finely powdered. A portion of the powder, equivalent to one tablet (0.29 g), was accurately weighted, dissolved in 100 mL of 50% v/v ethanol water mixture and filtered to remove insoluble material. The filtered solution was transferred to a 200 mL volumetric flask and taken to volume with doubly distilled water. The CBZ final concentration was about 1 mg∙mL−1. This solution was further diluted and was adjusted to pH 11 with NaOH (1 mol∙L−1) to obtain suitable experimental conditions for CBZ quantification.
2.5. General Procedure
Standard and samples were injected into a flowing buffer stream conformed by HCl 0.003 mol∙L−1. The continuous flow diagram is shown in Figure 2. In the “loading” position, the carrier was impulse by the first peristaltic pump (PP) to the fluorescence detector generating the baseline at λem = 434 nm and λex = 251; while the second peristaltic pump aspirate the sample plug, which is pushed towards the photoreactor using NaOH pH 11.0. Then, when leaving the photoreactor, the sample was collected into the MWCNTs mC located among 1 and 4 positions at the injection valve. In the “injection” position carrier flowed through the mC and carried the photoproduct to the detector. The second PP is turned off.
2.6. Multivariate Optimization
2.6.1. Screening Phase: Full Factorial Design (FFD)
Experimental factorial design is a factorial analysis involving 2k experiments (with k = number of studied variables). This type of analysis allows defining important variables affecting the analytical response(s), and its possible interactions. Moreover, a good experimental design provides a simple, efficient, and systematic approach to optimize designs for performance, quality, and cost. Consequently, an FFD was built for the evaluation of the main factors affecting the sensitivity and sampling frequency that may be quantified through Equations (1) and (2), respectively:
A two-level full factorial design with 2(4) = 16 experiments (Table 1) is described here for the variables: A) NaOH concentration, mol∙L−1 (standard/sample carrier concentration); B) Photoreactor flow rate, mL∙min−1; C) HCl concentration, mol∙L−1 (eluent concentration); D) eluent flow rate, mL∙min−1; while sensitivity, cps∙L∙mol−1 (S) and sampling frequency, samples∙h−1 (f) were regarded as the dependent variables. Maximum and minimum levels of the four factors were determined based on preliminary studies. All experiments were performed in random order to minimize the effects of uncontrolled factors that may introduce bias on the measurements.
2.6.2. Optimization Phase: Central Composite Design (CCD)
Once influencing factors were established with the factorial analysis, an RSM was used to find the best conditions for Sensitivity and Sampling frequency at
Figure 2. Schematic diagram of the on-line system for CBZ photoinduced determination. PhR: Photoreactor; mC: microcolumn; PP: peristaltic pumps; V: load/injection valve ((a) load position and (b) injection position); (d) spectrofluorimeter detector.
Table 1. Full factorial design (FFD) built for the responses Sensitivity (S) and Sampling frequency (f).
once. A spherical CCD was used here consisting of 27 experiments; i.e. combinations of 4 factorial points, 4 axial points, and replicates of a central point. In the common way, a simple response is analyzed, and the model analysis indicates areas in the design region where the process is likely to give desirable results. To do this, the Derringer’s desirability function (D) was used . The desirability approach is a widespread method that gives scores (Y) to a set of responses and chooses factor settings that maximize that score (Equation (3)) with di, as the individual desirability and k denoting the number of responses.
The response surface was built with the matrix of data acquired according to a Central composite design that consisted of 27 experiments. Experiments were combinations of the independent variables in the following ranges: NaOH concentration 0.01 - 0.1 mol∙L−1; Photoreactor flow rate 0.6 - 1.4 mL∙min−1; HCl concentration, 0.01 - 0.1 mol∙L−1; eluent flow-rate, 0.6 - 1.4 mL∙min−1 (Table 2). The levels studied were selected considering the results of the FFD. Sensitivity and Sampling frequency were thus fitted to polynomial models and the maximum in the experimental domain was then assigned according to the desirability function maximum obtained through simplex optimization. All experiments were performed in random order to minimize the effects of uncontrolled factors that may introduce bias in the measurements.
2.6.3. Outliers Removed after Cook’s Distance Analysis
Cook’s distance, Di, is used in regression analysis to find outliers in a set of predictor variables. It is a way to identify points that negatively affect a regression model. The measurement is a combination individual leverage and residual values; the higher the leverage and residuals, the higher the Cook’s distance. In this study, the models used as stated in Sections 2.6.1 and 2.6.2 were analyzed by Cook’s distance test and outliers were removed.
3. Results and Discussion
3.1. Fluorescence Characteristics of CBZ in Aqueous Medium
Figure 1 shows CBZ emission spectra and chemical structure of CBZ. Figure 1(a) and Figure 1(b) show the emission spectra of the aqueous and basic solution of CBZ without irradiation. As expected, upon exposure to UV-radiation CBZ in a basic medium (Figure 1(d)) was quickly derived to fluorescent compound (photochemical product) with a maximum emission wavelength of 434 nm when it is excited at 251 nm. This phenomenon was less observed in the irradiated aqueous solution (Figure 1(c)). The fluorescent intensity has been improved 10 times by inducing hydrolysis in basic medium, which would favor the formation of the photoproducts described by Pan et al. . These wavelengths were selected to measure the fluorescence intensity for the following assays. Based on these observations, the optimum medium for the on-line photochemical
Table 2. Central Composite Design with seven replicates of the central point built for Sensitivity (S) and Sampling frequency (f).
reaction of CBZ was NaOH solution which was used to be on-line merged with the sample zone. The NaOH concentration (standard/sample carrier) and the photoreactor flow rate were optimized.
3.2. Optimization of FIA System and Flow Conditions
Several solvents (such as water, ethanol, acid and basic buffers and mixtures of them) were tested as eluents to achieve complete release of the analyte from the minicolumn. Weak fluorescence signals (evaluated as peak area) were obtained when water, ethanol, alkaline buffers, or mixtures of them were employed. However, adequate elution was achieved when an HCl solution was employed. Therefore, its concentration and flow rate were optimized through multivariate experimental design. Furthermore, considering the results of some earlier works , the influence of the elution flow direction on the sensitivity and peak shape was studied. The best results were obtained when the elution was performed in a countercurrent manner (see Figure 2).
3.3. Experimental Design for CBZ Determination
Factorial designs are models that allow the size of an experiment being reduced by obtaining critical information about the effects of controllable factors on a response of interest (i.e., sensitivity and sampling frequency). In this research, the chosen design was the two-level factorial, which consists of 2(k) experiments, where k is the number of factors under study. Based on preliminary studies and experiments, four factors were considered including: NaOH concentration, mol∙L−1 (standard/sample carrier concentration, factor A); photoreactor flow rate, mL∙min−1 (factor B); HCl concentration, mol∙L−1 (eluent concentration, factor C); eluent flow rate, mL∙min−1 (factor D) (Table 1). In order to reduce the effect of uncontrolled variables, the experiment runs were performed in a random manner.
The evaluation consisted of the analysis of stock standard solutions of CBZ in ways to miming the analytical sample concentrations, and in all cited conditions. ANOVA tests (data not shown) were applied to experimental data in order to establish the significance of the factors studied (p < 0.05) upon the analytical responses. In order to get deeper insights, percentage normal probability plots for analyzed responses were built (Figure 3). Results indicated that all factors studied influenced significantly upon the two responses, but with different signs; i.e., Factors A and D, had positive effects on sensitivity; while B and C had negative. On the other hand, B, C and D effects had positive effects on sampling frequency whilst A had negative. Those observations were consistent with qualitative observations of the experimental setup, and with previous experiments within the studied experimental domain.
Once the level for main factors was determined, a central composite design was built, in which each design variable was assayed in the following ranges: NaOH concentration 0.01 - 0.1 mol∙L−1; Photoreactor flow rate 0.6 - 1.4 mL∙min−1; HCl concentration, 0.01 - 0.1 mol∙L−1; eluent flow-rate, 0.6 - 1.4 mL∙min−1. Polynomial models were used to fit the responses for all the 27 experiments (see Table 2). The model coefficients were calculated by multiple regressions and validated by the analysis of variance (ANOVA). In all cases, irrelevant terms were eliminated (alpha > 0.1), although some of them were maintained to ensure hierarchy. The optimization procedure was thus carried out and response surfaces were obtained. The CCD models yield polynomial equations in which the model terms are preceded by the adjusted coefficients, including the intercept. As can be observed the lack of fit was not significant (p < 0.05).
Figure 3. Normal Plots for the independent variables included in the FFD.
The responses in terms of studied factors effects were obtained (Equations (4) and (5)):
Desirability Function (DF) is one of the most popular chemometrics tools used to find the experimental conditions to reach simultaneously the optimal value for all the evaluated factors . Figure 4 shows the DF for each individual factor under study, while maintaining the others at their optimal values. The experimental conditions corresponding to one maximum (DF = 0.707); [NaOH] = 0.01 mol∙L−1; Photoreactor FR = 0.6 mL∙min−1; [HCl] = 0.03 mol∙L−1 and EFR = 1.4 mL∙min−1. These conditions were determined through simplex optimization using design points as starting conditions. Twenty-nine solutions were suggested, predicting values between 19 and 20 samples∙h−1 for sampling frequency, and 1.4 × 106 - 6.1 × 106 cps∙L∙mol−1 for sensitivity.
Figure 4. Response surface plots for the Desirability function in terms of (a) NaOH concentration and photoreactor flow rate; (b) NaOH and HCl concentrations (c) NaOH concentration and eluent flow rate (d) HCl concentration and eluent flow rate. Desirability function is showed in the four plots (a) (b) (c) and (d) in terms of two variables while the remaining were left at their optimal values; i.e. (a): NaOH concentration = 0.01; (b): Photoreactor FR = 0.6; (c): HCl Concentration = 0.03; and (d): Eluent FR = 1.4.
3.4. Analytical Figures of Merit
Figure 5 shows a flow diagram obtained for the injection of standard solutions of CBZ in the optimized conditions, showing transient signals with minimum dispersion. The triplicate signals put in evidence good repeatability. The equation for the calibration graph was obtained by the least-squares regression analysis, employing the areas of the analyte standard fluorescent signals: F = 19.138c + 20.982 where F is the relative fluorescence intensity and c the concentration of CBZ. The method was linear in the range 0.13 - 40 μg∙mL–1 CBZ. The correlation coefficient was 0.9956. The detection (0.04 μg∙mL–1) and quantification limits (0.13 μg∙mL–1) were calculated from the residual standard deviation for the calibration curve and its slope . The main advantages of the proposed method are its ease of sampling, which enables a throughput of 20 samples per hour achieving analytical characteristics comparable with other available methodologies. In contrast, the present method does not require prohibitive apparatus price nor pretreatment samples involve tedious time-consuming steps. Additionally, it shows a low consumption of reagents and the samples generating a low volume of waste. Therefore, the methodology represents an eco-friendly alternative without using toxic additives.
3.5. Application to Real Samples and Validation
The developed on-line photoinduced fluorescence approach was applied to determinate CBZ concentration in tablets and syrup.
Figure 5. Transient signal (FIAgram) obtained under optimized conditions. Standard solutions of CBZ (µg∙mL−1) = (a) 0.1; (b) 0.5; (c) 1.0; (d) 2.0; (e) 4.0; (f); 5.0; (g) real sample of CBZ syrup without standard solution addition. λex = 251 nm and λem = 434 nm; slits = 3/3 nm.
In order to study the accuracy, recovery studies were carried out by the standard addition method. Known amounts of analyte were added on sample solutions of different commercial formulations. Results showed quantitative recoveries for all real samples (Table 3), indicating good accuracy of the proposed procedure. The intra-day precision of the method, based on repeatability, were determined by replicating injections (n = 6) of the sample solutions prepared by the standard addition method. A relative standard deviation (RSD) < 4.36% was obtained in all cases. The method was also validated against the USP HPLC official method  ; Table 4 shows the results of the CBZ determination in pharmaceuticals by the developed spectrofluorimetric method and HPLC official methods. All values were in good statistical agreement with the nominal values, indicating good tolerance to common excipients. In comparison with the methodology proposed, the official method (USP) describes an HPLC method involving hight-cost instrument, using organic solvents, and is laborious and time-consuming. The proposed method is simple, sensitive and offers higher throug put sample than the official method. In addition, it requires small amounts of samples and reagents, so it is safe and environmentally friendly.
An on-line photoinduced fluorescence method for CBZ determination in pharmaceuticals coupled to multivariate calibration has been developed and validated. In alkaline aqueous medium, photo-degradation of CBZ is rapid and leads to the formation of a fluorescent photoproduct. The automation of the methodology was achieved by employing a very simple scheme of flow injection analysis, obtaining a sample throughput of 20 samples per hour (considering the required time for CBZ photoreaction). A proposed strategy was optimized in this
Table 3. Determination of CBZ in pharmaceutical samples by the developed method.
aMean value, n = 6. bRecovery = [(found−base)/added] × 100. RSD: Relative standard deviation.
Table 4. Analysis of CBZ in pharmaceutical samples by the developed spectrofluorimetric method and HPLC official methods.
aSpectrofluorimetric method developed in this work. bUSP official method using HPLC. cPercent relative error (calculated considering that the tablets contain the amount reported by the manufacturing laboratories).
study using the response surface methods. This method advantages other available as it is fast, economical and produces little waste.
Authors gratefully acknowledge at Instituto de Química de San Luis– Consejo Nacional de Investigaciones Científicas y Tecnológicas (INQUISAL–CONICET, Project PIP-CONICET: 11220130100605 C0) and Universidad Nacional de San Luis (Project 22/Q228) for financial support.
 Maiti, R., Mishra, B.R., Sanyal, S., Mohapatra, D., Parida, S. and Mishra, A. (2017) Effect of Carbamazepine and Oxcarbazepine on Serum Neuron-Specific Enolase in Focal Seizures: A Randomized Controlled Trial. Epilepsy Research, 138, 5-10.
 Ghoraba, Z., Aibaghi, B. and Soleymanpour, A. (2017) Application of Cation-Modified Sulfur Nanoparticles as an Efficient Sorbent for Separation and Preconcentration of Carbamazepine in Biological and Pharmaceutical Samples Prior to Its Determination by High-Performance Liquid Chromatography. Journal of Chromatography B, 1063, 245-252.
 Frag, E.Y.Z., Zayed, M.A., Omar, M.M., Elashery, S.E.A. and Mohamed, G.G. (2012) Spectrophotometric Determination of Carbamazepine and Mosapride Citrate in Pure and in Pharmaceutical Preparations. Arabian Journal of Chemistry, 5, 375-382.
 Taibon, J., Schmid, R., Lucha, S., Pongratz, S., Tarasov, K., Seger, C., Timm, C., Thiele, R., Mark Herlan, J. and Kobol, U. (2017) An LC-MS/MS Based Candidate Reference Method for the Quantification of Carbamazepine in Human Serum. Clinica Chimica Acta, 472, 35-40.
 Datar, P.A. (2015) Quantitative Bioanalytical and Analytical Method Development of Dibenzazepine Derivative, Carbamazepine: A Review. Journal of Pharmaceutical Analysis, 5, 213-222.
 Lozano, V.A. and Escandar, G.M. (2013) Second-Order Advantage with Excitation-Emission Photoinduced Fluorimetry for the Determination of the Antiepileptic carbamazepine in Environmental Waters. Analytica Chimica Acta, 782, 37-45.
 Escandar, G.M., GonzÁlez Gómez, D., Espinosa Mansilla, A., Muñoz de la Peña, A. and Goicoechea, H.C. (2004) Determination of Carbamazepine in Serum and Pharmaceutical Preparations Using Immobilization on a Nylon Support and Fluorescence Detection. Analytica Chimica Acta, 506, 161-170.
 Zhang, Z.-Q., Liang, G.-X., Ma, J., Lei, Y. and Lu, Y.-M. (2006) A Sensitive Flow Injection Fluorimetry for the Determination of Carbamazepine in Human Plasma. Analytical Letters, 39, 2417-2428.
 Walash, M.I., El-Enany, N. and Askar, H. (2015) Validated Spectrofluorimetric Method for the Determination of Carbamazepine in Pharmaceutical Dosage Forms after Reaction with 4-chloro-7-nitrobenzo-2-oxa-1,3-diazole (NBD-Cl). Luminescence, 30, 1119-1124.
 Hurtado-SÁnchez, M.C., Lozano, V.A., Rodríguez CÁceres, M.I., DurÁn-MerÁs, I. and Escandar, G.M. (2015) Green Analytical Determination of Emerging Pollutants in Environmental Waters Using Excitation-Emission Photo Induced Fluorescence Data and Multivariate Calibration. Talanta, 134, 215-223.
 Peralta, C.M., Henestrosa, C., Gil, R.A., FernÁndez, L.P. and Acosta, G. (2017) Novel Spectrofluorimetric Method for Boldine Alkaloid Determination in Herbal Drugs and Phytopharmaceuticals. Spectrochimica Acta Part A, 184, 101-108.
 amini, Y., Faraji, M., Shariati, S., Hassani, R. and Ghambarian, M. (2008) On-Line Metals Preconcentration and Simultaneous Determination Using Cloud Point Extraction and Inductively Coupled Plasma Optical Emission Spectrometry in Water Simples. Analytica Chimica Acta, 612, 144-151.
 Bezerra, M.A., Santelli, R.E., Oliveira, E.P., Villar, L.S. and Escaleira, L.A. (2008) Response Surface Methodology (RSM) as a Tool for Optimization in Analytical Chemistry. Talanta, 76, 965-977.
 Martins, A.C., Bukman, L., Vargas, A.M.M., Barizão, É.O., Moraes, J.C.G., Visentainer, J.V. and Almeida, V.C. (2013) The Antioxidant Activity of Teas Measured by the FRAP Method Adapted to the FIA System: Optimising the Conditions Using the Response Surface Methodology. Food Chemistry, 138, 574-580.
 Acosta, G., Silva, R., Gil, R.A., Gómez, R. and FernÁndez, L.P. (2013) On-Line Enantioseparation of Chlorpheniramine Using β-Cyclodextrin and Carbon Nanotubes after Multivariate Optimization. Talanta, 105, 167-172.
 Pan, Y., Cheng, S.S., Yang, X., Ren, J., Fang, J., Shang, C., Song, W., Lian, L. and Zhang, X. (2017) UV/Chlorine Treatment of Carbamazepine: Transformation Products and Their Formation Kinetics. Water Research, 116, 254-265.