- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Recently Accepted Articles ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
Journal of Biomedicine and Biotechnology
Volume 2012 (2012), Article ID 693848, 11 pages
Production of Chemoenzymatic Catalyzed Monoepoxide Biolubricant: Optimization and Physicochemical Characteristics
School of Chemical Sciences and Food Technology, Faculty of Science and Technology, Universiti Kebangsaan Malaysia, Bangi, 43600 Selangor, Malaysia
Received 26 July 2011; Revised 22 August 2011; Accepted 24 August 2011
Academic Editor: Rumiana Koynova
Copyright © 2012 Jumat Salimon et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Linoleic acid (LA) is converted to per-carboxylic acid catalyzed by an immobilized lipase from Candida antarctica (Novozym 435). This per-carboxylic acid is only intermediate and epoxidized itself in good yields and almost without consecutive reactions. Monoepoxide linoleic acid 9(12)-10(13)-monoepoxy 12(9)-octadecanoic acid (MEOA) was optimized using D-optimal design. At optimum conditions, higher yield% (82.14) and medium oxirane oxygen content (OOC) (4.91%) of MEOA were predicted at 15 μL of H2O2, 120 mg of Novozym 435, and 7 h of reaction time. In order to develop better-quality biolubricants, pour point (PP), flash point (FP), viscosity index (VI), and oxidative stability (OT) were determined for LA and MEOA. The results showed that MEOA exhibited good low-temperature behavior with PP of . FP of MEOA increased to comparing with of LA. In a similar fashion, VI for LA was 224 generally several hundred centistokes (cSt) more viscous than MEOA 130.8. The ability of a substance to resist oxidative degradation is another important property for biolubricants. Therefore, LA and MEOA were screened to measure their OT which was observed at 189 and , respectively.
Nowadays, people are concerned about green technology that is more environmental friendly and can save our environment. Green technology concept is good for industry which can be realized with sustainable and renewable resources. Epoxidized oils are currently produced by epoxidation of unsaturated plant oils, soybean, linseed oil , and Jatropha curcas seed oil . Although there are several methods available to epoxidize the double bonds of unsaturated fatty acids, the only method applied on industrial scale is the Prileshajev epoxidation reaction. In this reaction, a peracid from a short-chain fatty acid and hydrogen peroxide under strong acidic conditions is used as the oxidizing agent. The presence of the strong acid in the reaction mixture, however, causes the formation of side products such as vicinal diols, estolides, and other dimmers. Although a carful choice of the peracid and the reaction conditions can help to minimize the epoxides loss, the selectivity of industrial epoxidation of plant oils rarely exceeds 80% . Monoepoxidized oil has gained widely application including functional fluids, fuel additives, polyol replacements, pharmaceutical molecules, reactive diluents in cationic and UV cure applications, surfactants, adhesives, sealants, coatings, and biolubricants .
The current used of epoxy biolubricant is costly and harmful to environment. This study will be one of the solutions to obtain cheaper biolubricant to replace the existing epoxy biolubricant. A milder and more selective epoxidation process has been suggested , wherein a lipase is used to catalyze the peracid formation. In the lipase-mediated epoxidation of unsaturation fatty acid, the acid itself reacts with hydrogen peroxide to form the peracid, which then epoxidized the double bond. The reaction is therefore often referred to as “self-epoxidation reaction,” in spite of the fact that the second step proceeds predominantly via an intermolecular process . Among the various lipases studied so far, Novozym 435, a commercial preparation of Candida Antarctica lipase B, has been shown to be the most effective. Most of the investigations have involved dilution of the substrate in organic solvent. Recently, lipase-mediated epoxidation in a solvent-free medium has also been reported .
This study focuses on the physicochemical characteristics such as pour point (PP), flash point (FP), viscosity index (VI), and oxidative stability (OT) of the monoepoxidation process, with the aim to determine optimal reaction conditions with regard to reaction efficiency and enzyme stability using D-optimal design. In this paper, producing of MEOA over diepoxide as biolubricants has been studied for improvement of the physicochemical characteristics (i.e., PP, FP, and VI). Linoleic acid (LA; 18 : 2) was used as the model substrate. Figure 1 demonstrates the scheme for the self-epoxidation reaction of LA. Monoepoxidation of LA results in the mixture of two monoepoxides (cis-9, 10-epoxy 12 c- 18 : 1 and cis-12, 13 epoxy 9 c- 18 : 1).
2.1. Experimental Procedure
The enzymatic monoepoxidation was carried out using Novozym 435, a commercial catalyst made up of lipase, from Candida antartica, immobilized on a polyacrylate resin. 3 Factors (variables) such as hydrogen peroxide (L, ), enzyme (mg, ), and ration time (h, ) were obtained according to . Table 1 shows different ratios of hydrogen peroxide, enzyme, and reaction time using D-optimal design. In a typical chemoenzymatic monoepoxidation of LA 9(12)-10(13)-monoepoxy 12(9)-octadecanoic acid (MEOA), the LA (1.4 g) was dissolved in 10 mL of toluene, and the lipase was added. After stirring for 15 min, 30% H2O2 were added, and every 15 min the addition was repeated. Afterwards, the lipase was removed by filtration, the mixture was washed with water to remove excess H2O2, and the organic phase was dried over anhydrous sodium sulfate, and solvent was evaporated in a vacuum rotary evaporator.
After the laboratory reaction was completed, the actual yield to the theoretical yield was compared to determine the percent yield
2.1.2. Oxirane Oxygen Content
OOC% was calculated according to .
2.1.3. Iodine Value
Iodine value was determined according to .
2.2. Experimental Design and Statistical Analysis
To explore the effect of the operation variables on the response in the region of investigation, a D-optimal design at three levels was performed. Hydrogen peroxide (L, X1), amount of enzyme (mg, X2), and reaction time (h, X3) were selected as independent variables. The range of values and coded levels of the variables are given in Table 1.
A polynomial equation was used to predict the response as a function of independent variables and their interactions. In this work, the number of independent variables was three, and therefore the response for the quadratic polynomials becomes where ; ; ; are constant, linear, square, and interaction regression coefficient terms, respectively, and and are independent variables. The Minitab software version 14 (Minitab Inc., USA) was used for multiple regression analysis, analysis of variance (ANOVA), and analysis of ridge maximum of data in the response surface regression (RSREG) procedure. The goodness of fit of the model was evaluated by the coefficient of determination and its statistical significance that was checked by the F-test.
2.3. Physicochemical Characteristics
2.3.1. Fourier Transforms Infrared Spectroscopy
Fourier transforms infrared spectroscopy (FTIR) has been carried out according to . FTIR of the products was recorded by using Perkin Elmer Spectrum GX spectrophotometer in the range of 400–4000 cm−1 A very thin film of products was covered on NaCl cells (4 mm thickness) and was used for analysis.
2.3.2. Nuclear Magnetic Resonance Spectroscopy
Nuclear magnetic resonance spectroscopy for proton 1H and 13C NMR has been carried out according to . 13C and 1H NMR spectra were recorded using a 400 MHz Bruker 300 NMR spectrometer. For these analysis, 20 mg of samples were dissolved in 1 mL of CDCl3 solvent and introduced into the NMR tube.
2.3.3. Low Temperature Operability
The pour point (PP) is defined as the lowest temperature at which a liquid remains pourable (meaning it still behaves as a fluid). This method is routinely used to determine the low-temperature flow properties of fluids. PP values were measured according to  using a phase technology analyzer, Model PSA-70S (Hammersmith Gate, Richmond, BC, Canada). Each sample was run in triplicate, and average values rounded to the nearest whole degree are reported. For a greater degree of accuracy, PP measurements were done with a resolution of 1°C instead of the specified 3°C increment. Generally, materials with lower PP exhibit improved fluidity at low temperatures than those with higher PP.
2.3.4. Flash Point Values
The flash point of a volatile liquid is the lowest temperature at which it can vaporize to form an ignitable mixture in air. Flash point determination was run according to the American National Standard Method using a Tag Closed Tester . Each sample was run in triplicate, and the average values rounded to the nearest whole degree are reported.
2.3.5. Viscosity Index Measurements
Viscosity index (VI) is an arbitrary measure for the change of kinematic viscosity with temperature. It is used to characterize the lubricating oil in the automotive industry. Automated multirange viscometer tubes HV M472 obtained from Walter Herzog (Germany) were used to measure viscosity. Measurements were run in a Temp-Trol (Precision Scientific, Chicago, Ill, USA) viscometer bath set at 40.0 and 100.0°C. The viscosity and viscosity index were calculated using ASTM method . Triplicate measurements were made, and the average values were reported.
2.3.6. Oxidative Stability
Pressurized DSC (PDSC) experiments were accomplished using a DSC 2910 thermal analyzer from TA Instruments (Newcastle, Del, USA). Typically, a 2 μL sample, resulting in a film thickness of <1 mm, was placed in an aluminum pan hermetically sealed with a pinhole lid and oxidized in the presence of dry air (Gateway Airgas, St. Louis, Mo, USA), which was pressurized in the module at a constant pressure of 1378.95 kPa (200 psi). A 10°C min−1 heating rate from 50 to 350°C was used during each experiment. The oxidation onset (OT, °C) was calculated from a plot of heat flow (W/g) versus temperature for each experiment. The sample was run in triplicate, and average values rounded to the nearest whole degree are reported.
3. Results and Discussion
In this study, following the reaction experiments, the response surface is approximated by D-optimal design. Hydrogen peroxide is an important reactant for the formation of peracids from fatty acids; hence, the influence of its amount on the monoepoxidation reaction was studied. The addition of H2O2 solution to the reaction medium (toluene with LA) leads to the formation of two distinct phases: an organic phase and an aqueous phase. The designs and the responses yield% of MEOA (Y1), OOC% (Y2), and IV mg/g (Y3) are given in Table 2.
The Novozym 435, being adsorbed on a hydrophobic carried, is mainly present in the organic phase, while H2O2 will be partitioned in both the aqueous and the organic phases, with the concentration being higher in the aqueous phase (Figure 2) . This was determined performing varying amounts of H2O2 (15, 17.5, and 20 L) which has been added every 15 min, and the addition was repeated 24 times, Novozym 435 (80, 100, and 120 mg), and different times (6, 7, and 8 h). A stoichiometric excess of the required amount of the peroxide was used to compensate for its possible decomposition by light and temperature.
Table 2 showed that the yield percentage of MEOA, Y1, has increased to 82.14%, while OOC%, Y2, 4.91 and iodine value, Y3, 66.65 which are considerably compared to the theoretical (OOCt) 9.02% and initial iodine value (IV°) 157.35 mg/g. Subsequent experiments were performed using different amounts of H2O2 15, 17.5, and 20 L for every 1.4 gm of LA in one single step. As seen in Table 2, there was a clear increase in the reaction rate (OOC%) and a decrease (IV mg/g) with increasing H2O2 amount. With 15 L, monoepoxidation was achieved at 7 h using 120 mg Novozym 435, while total epoxidation was observed within 10 h using 30 L. Increasing the peroxide amount used for the reaction results in increasing peracid formation. In the state of partial epoxidation, the amount of peracids accumulated is not significant (less than 2% of the total) , because the chemical reaction in which they are consumed is very fast. But once all the double bonds are epoxidized, the remaining peracid is not consumed.
The quadratic regression coefficient obtained by employing a least squares method technique to predict quadratic polynomial models for the yield% of MEOA (Y1), OOC% (Y2), and IV mg/g (Y3) is given in Tables 3, 4, and 5. For the yield% of MEOA (Y1), the linear term of Novozym 435 catalyst amount (X2), quadratic terms of H2O2 (X11) and Novozym 435 catalyst amount (X22) were significant (). The interaction between H2O2 (X11) and Novozym 435 catalyst amount (X12) and the interaction between H2O2 (X11), and reaction time (X13) were significant (), while its quadratic term of reaction time (X33) was highly significant ().
Highly significant () terms of OOC% (Y2) and IV mg/g (Y3) for the H2O2 (X1) were linear, while linear term of IV mg/g for the reaction time (X3) was significant (). The coefficients of independent variables (H2O2; X1, catalyst Novozym 435; X2, and reaction time; X3) were determined for the quadratic polynomial models Tables 3, 4 and 5.
The lack of fit F-value for all the responses showed that the lack of fit is not significant () relatively to the pure error. This indicates that all the models predicted for the responses were adequate. Regression models for data on responses Y1, Y2, and Y3 were highly significant () with satisfactory . However, for Y2 (0.69) was lower although the model was significant. Table 6 summarizes the analysis of variance (ANOVA) of all the responses of this study.
These results suggest that linear effect of hydrogen peroxide is the primary determining factors for MEOA. Reference  also concluded that this variable had a very large effect on the results of their monoepoxidation study. Final equations in terms of actual factors are:
RSM is one of the best ways of evaluating the relationships between responses, variables, and interactions that exist. Significant interaction variables in the fitted models (Tables 3 to 5) were chosen as the axes (amount of H2O2; X1, catalyst Novozym 435; X2 and reaction time; X3) for the response surface plots. The relationships between independent and dependent variables are shown in the three-dimensional representation as response surfaces. In a contour plot, curves of equal response values are drawn on a plane whose coordinates represent the levels of the independent factors. Each contour represents a specific value for the height of the surface above the plane defined for a combination of the levels of the factors. Therefore, different surface height values enable one to focus attention on the levels of the factors at which changes in the surface height occur .
Figures 3 to 5 are the design-expert plots for all the responses. In the MEOA, performing the technique using low amount of H2O2 would give the desired OOC% of MEOA as shown in Figure 4, while IV (Figure 5) was higher at this condition. As shown in Figures 4 and 5, the increasing amount of H2O2 led to an increase in the OOC% and reduction of IV mg/g. The relationships between the parameters and MEOA were linear or almost linear. High OOC% could be obtained by using high amount of H2O2 at high reaction time. Experimental variables should be carefully controlled in order to recover a medium percentage of MEOA of interest with reasonable yield .
Optimum conditions of the experiment to obtain higher yield% of MEOA and medium OOC% were predicted at amount of H2O2 L of 15, catalyst Novozym 435 of 120 mg, and 7 h of reaction time. At this condition, the yield% of MEOA was 82.14%, 4.91% of OOC, and 66.65 mg/g of IV.
In order to prove the present oxirane ring of MEOA, final product was tested by FTIR. The comparison between LA (a) and MEOA (b), FTIR spectra is shown in Figure 6. The main peaks and their assignment to functional groups are given in Table 7. Oxirane ring of MEOA can be detected at wavenumber 820 cm−1.
For the carboxylic acid carbonyl functional groups (C=O), FTIR spectrum showed absorption bands of LA and MEOA at 1719 and 1711 cm−1, respectively, while stretching vibration peak of C=C can be detected at wavenumber 3009 cm−1 , and peaks at 2927–2856 cm−1 indicated the CH2 and CH3 scissoring of LA and MEOA based on Figures 6(a) and 6(b). FTIR spectrum also showed absorption bands at 722, 723 for (C-H) group.
Figure 7(a) indicates the 13C NMR spectrum of LA. The 13C spectroscopy shows the main signals assignment of the LA as shown in Table 8(a). The signals at 180.49 ppm refer to the carbon atom of the carbonyl group (carboxylic acid). The signals at 128.27–130.38 ppm refer to the unsaturated carbon atoms (olefinic carbons); 24.86–29.79 ppm due to methylene carbon atoms in fatty acid moieties of LA .
Figures 7(b) and 7(c) can be confirmed by the oxirane ring of MEOA 54.59–57.29 ppm and diepoxide LA at about 54.61–57.32 ppm. Indeed, it appeared that the signals were present in the MEOA, as four peaks of roughly equal intensity (132.89, 132.72, 130.15, and 124.02 ppm) were observed in the alkenic carbon region in the 13C NMR spectrum Figure 7(b), while they disappeared in the diepoxide LA (Figure 7(c)) . The 13C-NMR spectra indicate the existence of carbonyl group (carboxylic acid) in their structure MEOA 179.32 ppm and diepoxide LA at about 178.79 ppm. The other distinctive signals were aliphatic carbons MEOA at about 22.69–29.38 ppm and di-epoxide LA at about 22.77–29.44 ppm, which are common for these types of compounds .
1H NMR spectroscopy shows the main signals assignments in LA, MEOA and di-epoxide LA as shown in Table 9. The 1H NMR spectra for the products show some of the key features for a typical (-CH-O-CH-) at about 2.92–3.12 ppm of MEOA, and about 2.99–3.13 ppm of di-epoxide LA (Figure 8(b)). The distinguishable groups are the protons of the terminal methyl of the fatty acid chain. The signals at 0.88–0.86 ppm referred to the methylene group (-CH3) of LA Figure 8(a) which also appear in MEOA 0.86–0.88 ppm and di-epoxide LA 0.88–0.92 ppm (Figures 8(b) and 8(c)) next to the terminal methyl (-CH2) at 1.30–2.77 ppm of LA, 1.29–2.33 of MEOA, and 1.34–2.36 ppm of di-epoxide LA.
However, the methane proton signals (-CH=CH-) were shifted upfield at about 5.35–5.36 ppm of LA and 5.38–5.49 ppm of MEOA , while they disappeared in di-epoxide LA. Another distinctive feature is the hydroxyl proton (-COOH) of the carboxylic acid at about 7.27 ppm.
This approach is used here to improve the low-temperature flow behavior of fatty acids by monoepoxide ring. MEOA improves the PP at −41°C significantly comparing with LA at −2°C. Monoepoxidation to produce MEOA was the most effective decreasing the PP to −41°C (Table 10). It can be assumed that the presence of oxirane ring on the fatty acid creates a steric barrier around the individual molecules and inhibits crystallization, resulting in lower PP .
Flash point is another important factor in determining how well oil will behave as a potential biolubricant. The flash point is often used as a descriptive characteristic of oil fuel, and it is also used to describe oils that are not normally used as fuels. Flash point refers to both flammable oils and combustible oils. There are various international standards for defining each, but most agree that oils with a flash point less than 43°C are flammable, while those having a flash point above this temperature are combustible . Table 10 has shown the improvement in flash point of MEOA which increased to 128°C comparing with 115°C of LA that means the result agrees with the various international standards.
The efficiency of the biolubricant in reducing friction and wear is greatly influenced by its viscosity. Generally, the least viscous biolubricant which still forces the two moving surfaces apart is desired. If the biolubricant is too viscous, it will require a large amount of energy to move; if it is too thin, the surfaces will rub and friction will increase. The viscosity index highlights how a biolubricants viscosity changes with variations in temperature. Many biolubricant applications require performing across a wide range of conditions, for example, in an engine. Automotive biolubricants must reduce friction between engine components when it is started from cold (relative to engine operating temperatures) as well as when it is running (up to 200°C). The best oils (with the highest VI) will not vary much in viscosity over such a temperature range and therefore will perform well throughout. In MEOA, decreased viscosity index (VI) to 130.8 of MEOA is the result of less double bond (Table 10).
The ability of a substance to resist oxidative degradation is another important characteristic of biolubricants. Therefore, MEOA was screened to measure their oxidation stability using PDSC through determination of OT. PDSC is an effective method for measuring oxidation stability of oleochemicals in an accelerated mode . The OT is the temperature at which a rapid increase in the rate of oxidation is observed at a constant, high pressure (200 psi). A high OT would suggest high oxidation stability of the material. OT was calculated from a plot of heat flow (W/g) versus temperature that was generated by the sample upon degradation and by definition. In this chapter, oxirane ring of MEOA did not show improvement in the oxidation stability of OT for MEOA at 168°C comparing with LA at 189°C (Table 10).
RSM and 3-factor D-optimal design were employed for optimization of MEOA. From the present study, it is evident that hydrogen peroxide is the most critical parameter influencing the chemoenzymatic monoepoxidation reaction. Increase in the hydrogen peroxide amount has a strong effect on the reaction kinetics; however, a large excess of hydrogen peroxide results in accumulation of peracid in the final product. Based on the results obtained, MEOA had a positive influence on the low-temperature properties PP and FP.
The authors thank UKM and the Ministry of Science and Technology for research Grants UKM-GUP-NBT-08-27-113 and UKM-OUP-NBT-29-150/2011.
- M. R. G. Klaas and S. Warwel, “Complete and partial epoxidation of plant oils by lipase-catalyzed perhydrolysis,” Industrial Crops and Products, vol. 9, no. 2, pp. 125–132, 1999.
- P.-P. Meyer, N. Techaphattana, S. Manundawee, S. Sangkeaw, W. Junlakan, and C. Tongurai, “Epoxidation of soybean oil and Jatropha oil,” Thammasat International Journal of Science Technology, vol. 13, pp. 1–5, 2008.
- A. Adhvaryu and S. Z. Erhan, “Epoxidized soybean oil as a potential source of high-temperature lubricants,” Industrial Crops and Products, vol. 15, no. 3, pp. 247–254, 2002.
- S. Warwel and M. R. G. Klaas, “Chemo-enzymatic epoxidation of unsaturated carboxylic acids,” Journal of Molecular Catalysis B, vol. 1, no. 1, pp. 29–35, 1995.
- M. R. G. Klaas and S. Warwel, “Lipase-catalyzed preparation of peroxy acids and their use for epoxidation,” Journal of Molecular Catalysis A, vol. 117, no. 1–3, pp. 311–319, 1997.
- C. Orellana-Coca, S. Camocho, D. Adlercreutz, B. Mattiasson, and R. Hatti-Kaul, “Chemo-enzymatic epoxidation of linoleic acid: parameters influencing the reaction,” European Journal of Lipid Science and Technology, vol. 107, no. 12, pp. 864–870, 2005.
- S. W. Lin, T. T. Sue, and T. Y. Ai, PORIM Test Methods, vol. 1, PORIM, Bandar Baru Bangi, Malaysia, 1995.
- J. Salimon and N. Salih, “Modification of epoxidized ricinoleic acid for biolubricant base oil with improved flash and pour points,” Asian Journal of Chemistry, vol. 22, no. 7, pp. 5468–5476, 2010.
- ASTM Standard D5949, Standard Test Method for Pour Point of Petroleum (Automatic Pressure Pulsing Method), ASTM, West Conshohocken, Pa, USA.
- ASTM Standard D 56-79, Standard Test Method for Flash Point of Liquids with a Viscosity Less than, 45 Saybolt Universal Seconds (SUS) at (that don’t contain suspended solids and don’t tend to form a surface film under test).
- ASTM D 2270-93, Standard Practice for Calculating Viscosity Index from Kinematic Viscosity at 40 and 100°C, ASTM, West Conshohocken, Pa, USA.
- M. Wu, H. Ding, S. Wang, and S. Xu, “Optimization conditions for the purification of linoleic acid from sunflower oil by urea complex fractionation,” Journal of the American Oil Chemist' Society, vol. 85, pp. 677–684, 2008.
- G. Socrates, Infrared and Raman Characteristic Group Frequencies: Tables and Charts, John Wily & Sons, Chichester, UK, 3rd edition, 2001.
- G. Du, A. Tekin, E. G. Hammond, and L. K. Woo, “Catalytic epoxidation of methyl linoleate,” Journal of American Oil Chemical Society, vol. 81, no. 5, pp. 477–480, 2004.
- K. M. Doll, B. K. Sharma, and S. Z. Erhan, “Synthesis of branched methyl hydroxy stearates including an ester from bio-based levulinic acid,” Industrial and Engineering Chemistry Research, vol. 46, no. 11, pp. 3513–3519, 2007.
- H. S. Hwang and S. Z. Erhan, “Synthetic lubricant basestocks from epoxidized soybean oil and Guerbet alcohols,” Industrial Crops and Products, vol. 23, no. 3, pp. 311–317, 2006.
- B. K. Sharma, K. M. Doll, and S. Z. Erhan, “Ester hydroxy derivatives of methyl oleate: tribological, oxidation and low temperature properties,” Bioresource Technology, vol. 99, no. 15, pp. 7333–7340, 2008.
- J. Salimon, N. Salih, and E. Yousif, “Chemically modified biolubricantbasestocks from epoxidized oleic acid: Improved low temperature properties and oxidative stability,” Journal of Saudi Chemical Society, vol. 15, pp. 195–201, 2011.
- Y. Y. Zhang, T. H. Ren, H. D. Wang, and M. R. Yi, “A comparative study of phenol-type antioxidants in methyl oleate with quantum calculations and experiments,” Lubrication Science, vol. 16, no. 4, pp. 385–392, 2004.