About this Journal Submit a Manuscript Table of Contents
Journal of Thermodynamics
Volume 2012 (2012), Article ID 921693, 6 pages
http://dx.doi.org/10.1155/2012/921693
Research Article

Phase Behavior at High Pressure of the Ternary System: CO𝟐, Ionic Liquid and Disperse Dye

1Department of Chemical Engineering, State University of Maringá, Avenue Colombo 5790, 87020-900 Maringá PR, Brazil
2Institute for Research and Technology, ITP, PEP/UNIT, Avenue Murilo Dantas 300, 49032-490 Aracaju SE, Brazil
3Department of Chemistry, State University of Maringá, Avenue Colombo 5790, 87020-900 Maringá PR, Brazil
4Department of Chemical Process, School of Chemical Engineering, State University of Campinas, Avenue Albert Einstein 500, 13081-970 Campinas SP, Brazil

Received 27 April 2011; Accepted 22 June 2011

Academic Editor: Ramesh Gardas

Copyright © 2012 Helen R. Mazzer 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.

Abstract

High pressure phase behavior experimental data have been measured for the systems carbon dioxide (CO2) + 1-butyl-3-methylimidazolium hexafluorophosphate ([bmim] [PF6]) and carbon dioxide (CO2) + 1-butyl-3-methylimidazolium hexafluorophosphate ([bmim] [PF6]) + 1-amino-2-phenoxy-4-hydroxyanthraquinone (C.I. Disperse Red 60). Measurements were performed in the pressure up to 18 MPa and at the temperature (323 to 353 K). As reported in the literature, at higher concentrations of carbon dioxide the phase transition pressure increased very steeply. The experimental data for the binary and ternary systems were correlated with good agreement using the Peng-Robinson equation of state. The amount of water in phase behavior of the systems was evaluated.

1. Introduction

In the past years the development of new technologies to minimize water consumption and wastewater emissions received more attention. In this context, the textile industry has a significant stake in the generation of effluents with highly pollution potential, using large amounts of water with chemical additives [1]. The process of dyeing using supercritical fluids is an alternative method that does not use water and any additive is frequently not required. In addition, such process allows the recovering and the reusing of the dyes. According to van der Kraan et al. [2], in dyeing of textiles consisting of nonpolar materials such as polyester, supercritical CO2 acts as a blowing agent and increases the rate of diffusion of the dye molecules inside the fiber. Several studies have been reported on this textile-dyeing technology [36]. The literature describes applications of supercritical technology in combination with the ionic liquids in both the extraction and the separation processes and as an alternative solvent for organic and inorganic compounds in reaction mediums [79]. The ionic liquids are described as “designer solvents,” where their characteristics can be adjusted to make the change either of cations or anions, or both. In this way, dyeing process involving association of supercritical technology with ionic liquids (ILs) can be used to improve the solubility of the dyes in the dyeing process of textile fibers and minimize the amount of energy demand.

In this context, the objective of this work is to report new experimental data of CO2 (1) + 1-butyl-3-methylimidazolium hexafluorophosphate [bmim] [PF6] + (2) and CO2 (1) + [bmim] [PF6] (2) C.I. Disperse Red 60 systems in the temperature range of 323 to 353 K. C.I. Disperse Red 60 is indeed a popular dye and one of the widely used dyes in dyeing process of synthetic textile fibers [2, 10, 11]. [bmim] [PF6] is one of the most studied ionic liquids (ILs) in the literature with great availability of physical and chemical data as well as phase equilibrium data [1214]. Keskin et al. [7], Huddleston et al. [15], and Jaitely et al. [16] have reported that the ionic liquid consisting of imidazole as a cation and [PF6] as an anion presents hydrophobic characteristics compatible to that of disperse dyes with wide uses in dyeing of synthetic fibers. The synthetic static method with variable-volume cell was used in the experimental measurements of phase equilibrium. The amount of water present in ILs was determined by Karl Fischer analysis before and after equilibrium measures once that the ionic liquid [bmim] [PF6] absorbs the moisture during the synthesis process due to its hygroscopic property [17]. A small amount of water in the ionic liquid can have a dramatic effect on the phase equilibrium behavior [13]. A technique to recover [bmim] [PF6] from the [bmim] [PF6] + C.I. Disperse Red 60 mixture was also proposed here once the reuse of [bmim] [PF6] makes the dyeing process more environmentally benign and economically attractive. Peng-Robinson equation of state with quadratic mixing rules was used to correlate the experimental data.

2. Experimental

2.1. Materials

Commercial CO2 and nitrogen (99%) were purchased from White Martins (Rio de Janeiro, Brazil). 1-Methylimidazole (99.0% GC) and 1-Bromobutane (96.0% GC) were purchased from Fluka (Sigma-Aldrich Ltda, Brazil). Ethyl acetate (99.5% P.A.), acetone (99.5% P.A.), activated charcoal (P.A.), anhydrous methanol (99.8% HPLC), dichloromethane (99.5% P.A.), and ethyl ether were purchased from Nuclear (CAQ-Chemical of Home Ind. Com. Ltda, Brazil). Potassium hexafluorophosphate (98% GC) was obtained from Aldrich (Sigma-Aldrich Ltda, Brazil). Silica gel mesh 0.035–0.070, 60 Å 99.5% was obtained from Acros (Sigma-Aldrich Ltda, Brazil). The Karl Fisher reagent (99.99% HPLC) was purchased from Merck (Merck S.A-Brazil). All water concentrations were measured through Karl Fischer (Analyzer, KF-1000) analysis. The 1-amino-2phenoxy-4-hydroxyanthaquinone (disperse Red 60) was supplied by DyStar (Brazil).

The ionic liquid [bmim] [PF6] was synthesized from 1-methylimidazol, bromobutane and potassium hexafluorophosphate (KPF6) by the process described in the literature [18]. [bmim] [PF6] was dried under vacuum at 333 K and stored under dried nitrogen. It was purified in a chromatography column with silica gel Acros (Sigma-Aldrich Ltda, Brazil). (silica gel 60; mesh 70–230) and dichloromethane. Hydrogen and Carbon Nuclear Magnetic Resonance (1H and 13C NMR) analyses of the [bmim] [PF6] obtained were δH (300 MHz; CDCl3) 8.34 (sl, 3 N–CH3–N), 7.31 (t, CH2=CH), 7.28 (t, CH=CH2), 4.09 (t, N–CH2), 3.83 (s, N–CH3), 1.78 (quintet, CH2–CH2), 1.27 (sextet, CH2–CH3), 0.85 (t, CH2–CH3) for δC (75 MHz, CDCl3) 135.8 (N–CH2–N), 123.8 (CH2=CH), 122.5 (CH=CH2), 49.8 (N–CH2), 36.1 (N–CH3), 31.8 (CH2–CH2), 19.3 (CH2–CH3), 13.3 (CH2–CH3). The analyses using NMR were performed in a Varian spectrometer, model Mercury Plus BB using CDCl3 as a solvent.

2.2. Equipment and Procedure

Phase equilibrium experiments (cloud points) were performed through the synthetic static method [19] in a high-pressure variable-volume view cell. The experimental apparatus and experimental procedure are similar to those used in previous studies of phase equilibrium carried out by our group [2022]. Briefly, the experiment begins with the introduction of a predetermined amount of compounds (CO2, IL, and dye) in the equilibrium cell so that the global composition of the mixture is known. In the sequence, the system was stirred for 12 h at room temperature for complete solubilization of CO2 in [bmim] [PF6]. The procedure used here follows the work of Aki et al. [23]. After the initial procedure, pressure and temperature are adjusted to form a homogeneous solution. After reaching the thermal equilibrium of the system, pressure is slowly decreased until the appearance of a second phase is visually detected.

After the measurements of the phase transitions, the IL + Dye mixture was submitted to a process to recover the IL. The process of IL recovery consisted of the addition of a dichloromethane solution and active coal in the mass ration of 1 : 3 into the mixture of IL + Dye. Afterwards, the solution was stirred at room temperature for 48 h. Then, the mixture was filtrated using qualitative filter paper in a chromatographic column with silica gel (silica gel 60 Å; mesh 0.035–0.70 mm). The purified solution was evaporated under vacuum at 343 K for 10 h.

The values of water content refer to values obtained before and after the phase equilibrium measurements. For the binary system CO2 + [bmim] [PF6] the water content was 0.10 ± 0.02% (after phase equilibrium) and 0.23 ± 0.02 (before measurement). For the ternary system CO2 + [bmim] [PF6] + C.I. disperse Red 60 the water content corresponded to 0.43 ± 0.04 and 0.26 ± 0.02, after and before measurements, respectively. For the binary system with IL recovered the measurement were 0.19 ± 0.01 and 0.66 ± 0.02, after and before, respectively.

2.3. Thermodynamic Modeling

Peng-Robinson equation [24] of state was used to model the vapor-liquid equilibrium data: 𝑃=𝑅𝑇𝑎𝑣𝑏.𝑣(𝑣+𝑏)+𝑏(𝑣𝑏)(1) The parameters a and by were calculated through the following mixing rules: 𝑎=𝑛𝑐𝑖=1𝑛𝑐𝑗=1𝑦𝑖𝑦𝑗𝑎𝑖𝑎𝑗1𝑘𝑖𝑗𝑏=𝑛𝑐𝑗=1𝑦𝑗𝑏𝑗,(2) where 𝑦𝑖 is the molar fraction of the component 𝑖 and 𝑎𝑖 and 𝑏𝑖 are the Peng-Robinson parameters for the pure component 𝑖. The binary interaction parameters, 𝑘𝑖𝑗, were determined through the correlation of experimental data obtained in this work. The optimization procedure employed for fitting 𝑘𝑖𝑗 was based on the Nelder and Mead simplex method [25].

3. Results and Discussion

Experimental data of phase transitions at high pressures of the CO2 (1) + [bmim] [PF6] (2) binary system at 333 K are presented in Figure 1 and Table 1. All phase transitions observed were vapor-liquid type. The CO2 (1) solubility in [bmim] [PF6] (2) presented in Figure 1 is compared with experimental data reported by Blanchard and Brennecke [26] and Fu et al. [17]. Solubility obtained here is in agreement with Fu et al. [17] data. However, experimental data reported by Blanchard and Brennecke [26] show a significant divergence compared to the results obtained in this work. Fu et al. [17] state that the discrepancy in the solubility measurements can be attributed to different experimental techniques, uncertainties in the experimental procedure, impurity in the IL, among other factors. From own our point of view, the differences verified in this figure are mainly due to the water content in the ionic liquids used in the experiments. Blanchard and Brennecke [26] reported 0.15% water content in IL prior to experimental measures. However, water content data after the experiments have not yet been reported.

tab1
Table 1: Experimental data of (L→ LV) of binary system CO2 (1) + [bmim] [PF6] (2).
921693.fig.001
Figure 1: Experimental data obtained in this work (▪) and reported by Blanchard et al. [13] (), Fu et al. [17] () for the CO2 (1) + [bmim] [PF6] (2) binary system at 333 K isotherm. Thermodynamic modeling (—).

It is possible to notice that CO2 solubility in the IL increases continually with pressure, reaching 0.50 at 333 K and 9.83 MPa. Such behavior is also reported by Bermejo et al. [27], Kamps et al. [14], and Shariati and Peters [12].

Figure 2 and Table 2 show the experimental phase behavior for the ternary system CO2 (1) + [bmim] [PF6] (2) + disperse dye Red 60 (3). In this case, it is possible to notice the increase of CO2 solubility with increasing pressure, an analogue behavior to that of the binary system. Moreover, CO2 solubility decreases as temperature increases. The IL addition in supercritical environment induces a significant solubilization of the dispersed dye Red 60, which is almost insoluble in supercritical CO2 in the temperature and pressure range investigated. The solubility of the dye Red 60 in CO2 is around 10−6 in the investigated range [28]. On the other hand, in this range the solubility of the dye Red 60 in CO2 and IL has an average value of 10−4. For all experimental conditions (Figures 1 and 2), a good agreement of the Peng-Robinson equation with the experimental data is observed.

tab2
Table 2: Experimental data of liquid-vapor ternary system CO2 (1) + [bmim] [PF6] (2) + Red Dye 60 (3).
921693.fig.002
Figure 2: Experimental data of bubble-type phase transition at (♦) 323 K, (▪) 333 K, () 343 K, and (•) 353 K for the ternary system CO2 (1) + [bmim] [PF6] (2) + disperse dye Red 60 (3). Thermodynamic modeling using the Peng-Robinson equation of state (—).

Figure 3 shows a comparison of the phase diagram for the systems CO2 (1) + [bmim] [PF6] (2) and CO2 (1) + [bmim] [PF6] (2) + disperse dye Red 60 (3) at 333 K. It is possible to observe that CO2 shows lower solubility in ionic liquid for the ternary system when compared to the binary system.

921693.fig.003
Figure 3: Comparison of Px diagram for the binary system CO2 (1) + [bmim] [PF6] (2) (♦) and ternary system CO2 (1) + [bmim] [PF6] (2) + disperse dye Red 60 (3) (▪) at 333 K.

Figure 4 shows the phase behavior for the binary systems CO2 (1) + [bmim] [PF6] (2) and CO2 (1) + [bmim] [PF6]recovered (2). Such results are summarized in Table 3. It can be observed that the phase behavior of both systems (CO2 (1) + [bmim] [PF6] (2) and CO2 (1) + [bmim] [PF6]recovered (2)) are similar. Solubility differences between these two systems can be attributed to the water content present in each situation. Furthermore, the RMN spectra of IL after removal of the disperse dye Red 60 were consistent with values found for pure LI produced in our laboratory and with those reported in literature [29].

tab3
Table 3: Experimental data of liquid-vapor binary system carbon dioxide + [bmim] [PF6]recovered.
921693.fig.004
Figure 4: Phase transition Px diagram for the binary system CO2 (1) + [bmim] [PF6] (2) (▪) and CO2 (1) + [bmim] [PF6]recovered (2) () at 333 K.

The PR-EoS with the quadratic mixing rules was used to modeling the phase transition experimental data for all systems. The binary interaction parameters in the mixing rules were optimized using the P, x data for each system. The critical properties and acentric factors of the IL were obtained of the Valderrama and Robles [30], and the critical properties, and acentric factors of the disperse dye Red 60 were obtained of the Mishima et al. [31]. The critical properties are presented in Table 4 and Table 5 shows the interaction parameters for each system.

tab4
Table 4: Critical properties and acentric factor of carbon dioxide, ionic liquid and disperse Red 60.
tab5
Table 5: PR-EoS binary interaction parameters for the binary and ternary systems.

4. Conclusions

It was verified that CO2 solubility in the ionic liquid decreased by the presence of dye. Such phenomenon results in higher vapor-liquid pressure transitions for the ternary system when compared to the binary system. Water content (impurities) in the systems studied considerably change in the pressure values of phase transition, corroborating literature observations [19, 31, 32]. Therefore, water content in IL synthesis should be carefully controlled. IL’s purification protocol for dye removal presented good results without IL physical chemical changes, which enables IL to be reused. Peng-Robinson equation of state presented a good performance in describing the experimental data observed. The results obtained might be useful in the development of dyeing processes using supercritical fluids in association with ionic liquids. Finally, the results of partition of the undisclosed dye composition between SCCO2 and an ionic liquid demonstrated that it is possible to increase significantly the disperse dye solubility in SCCO2 using ionic liquids as an addictive. This effect can be used to improve the color aspects in the textiles, one of the most desirable characteristics in the dyeing of textiles.

Acknowledgments

The authors would like to thank CNPq for the financial support. This paper is in memory of Proffessor D. S. Octávio Augusto Ceva Antunes with thanks for the orientation, friendship and dedication during the investigation.

References

  1. S. Papić, N. Koprivanac, A. Lončarić Božić, and A. Meteš, “Removal of some reactive dyes from synthetic wastewater by combined Al(III) coagulation/carbon adsorption process,” Dyes and Pigments, vol. 62, no. 3, pp. 291–298, 2004. View at Publisher · View at Google Scholar · View at Scopus
  2. M. van der Kraan, M. V. F. Cid, G. F. Woerlee, W. J. T. Veugelers, and G. J. Witkamp, “Equilibrium study on the disperse dyeing of polyester textile in supercritical carbon dioxide,” Textile Research Journal, vol. 77, no. 8, pp. 550–558, 2007. View at Publisher · View at Google Scholar · View at Scopus
  3. E. Bach, E. Cleve, and E. Schollmeyer, “Past, present and future of supercritical fluid dyeing technology—an overview,” Review of Progress in Coloration and Related Topics, vol. 32, pp. 88–102, 2002. View at Scopus
  4. A. Schmidt, E. Bach, and E. Schollmeyer, “Damage to natural and synthetic fibers treated in supercritical carbon dioxide at 300 bar and temperatures up to 160°C,” Textile Research Journal, vol. 72, no. 11, pp. 1023–1032, 2002. View at Scopus
  5. A. S. Özcan, A. A. Clifford, and K. D. Bartle, “Solubility of disperse dyes in supercritical carbon dioxide,” Journal of Chemical and Engineering Data, vol. 42, no. 3, pp. 590–592, 1997. View at Scopus
  6. T. Shinoda and K. Tamura, “Solubilities of C.I. Disperse orange 25 and C.I. Disperse blue 354 in supercritical carbon dioxide,” Journal of Chemical and Engineering Data, vol. 48, no. 4, pp. 869–873, 2003. View at Publisher · View at Google Scholar · View at Scopus
  7. S. Keskin, D. Kayrak-Talay, U. Akman, and Ö. Hortaçsu, “A review of ionic liquids towards supercritical fluid applications,” Journal of Supercritical Fluids, vol. 43, no. 1, pp. 150–180, 2007. View at Publisher · View at Google Scholar · View at Scopus
  8. L. A. Blanchard, D. Hancu, E. J. Beckman, and J. F. Brennecke, “Green processing using ionic liquids and CO2,” Nature, vol. 398, no. 6731, pp. 28–29, 1999. View at Scopus
  9. J. Dupont, R. F. De Souza, and P. A. Z. Suarez, “Ionic liquid (molten salt) phase organometallic catalysis,” Chemical Reviews, vol. 102, no. 10, pp. 3667–3692, 2002. View at Publisher · View at Google Scholar · View at Scopus
  10. A. Ferri, M. Banchero, L. Manna, and S. Sicardi, “Dye uptake and partition ratio of disperse dyes between a PET yarn and supercritical carbon dioxide,” Journal of Supercritical Fluids, vol. 37, no. 1, pp. 107–114, 2006. View at Publisher · View at Google Scholar · View at Scopus
  11. P. Bao and J. Dai, “Relationships between the solubility of C. I. disperse red 60 and uptake on PET in supercritical CO2,” Journal of Chemical and Engineering Data, vol. 50, no. 3, pp. 838–842, 2005. View at Publisher · View at Google Scholar · View at Scopus
  12. A. Shariati and C. J. Peters, “High-pressure phase equilibria of systems with ionic liquids,” Journal of Supercritical Fluids, vol. 34, no. 2, pp. 171–176, 2005. View at Publisher · View at Google Scholar · View at Scopus
  13. L. A. Blanchard, Z. Gu, and J. F. Brennecke, “High-pressure phase behavior of ionic liquid/CO2 systems,” Journal of Physical Chemistry B, vol. 105, no. 12, pp. 2437–2444, 2001. View at Scopus
  14. A. P. S. Kamps, D. Tuma, J. Xia, and G. Maurer, “Solubility of CO2 in the ionic liquid [bmim][PF6],” Journal of Chemical and Engineering Data, vol. 48, no. 3, pp. 746–749, 2003. View at Publisher · View at Google Scholar · View at Scopus
  15. J. G. Huddleston, A. E. Visser, W. M. Reichert, H. D. Willauer, G. A. Broker, and R. D. Rogers, “Characterization and comparison of hydrophilic and hydrophobic room temperature ionic liquids incorporating the imidazolium cation,” Green Chemistry, vol. 3, no. 4, pp. 156–164, 2001. View at Publisher · View at Google Scholar · View at Scopus
  16. V. Jaitely, A. Karatas, and A. T. Florence, “Water-immiscible room temperature ionic liquids (RTILs) as drug reservoirs for controlled release,” International Journal of Pharmaceutics, vol. 354, no. 1-2, pp. 168–173, 2008. View at Publisher · View at Google Scholar · View at Scopus
  17. D. Fu, X. Sun, J. Pu, and S. Zhao, “Effect of water content on the solubility of CO2 in the ionic liquid [bmim][PF6],” Journal of Chemical and Engineering Data, vol. 51, no. 2, pp. 371–375, 2006. View at Publisher · View at Google Scholar · View at Scopus
  18. J. Dupont, P. A. Z. Suarez, C. S. Consorti et al., “Preparation of 1-butyl-3-methyl imidazolium-based room temperature ionic liquids,” Organic Syntheses, vol. 79, pp. 236–243, 2002. View at Scopus
  19. M. A. McHugh and V. J. Krukonis, Supercritical Fluid Extraction: Principles and Practice, Butterworth-Heinemann, Boston, Mass, USA, 1994.
  20. R. Favareto, P. H. Fregadolli, V. F. Cabral, O. A. Antunes, and L. Cardozo-Filho, “Phase equilibria of acrylonitrile and p-bromobenzaldehyde in carbon dioxide,” Journal of Chemical and Engineering Data, vol. 53, no. 5, pp. 1080–1084, 2008. View at Publisher · View at Google Scholar · View at Scopus
  21. A. T. Souza, M. L. Corazza, L. Cardozo-Filho, R. Guirardello, and M. A. A. Meireles, “Phase equilibrium measurements for the system clove (Eugenia caryophyllus) Oil + CO2,” Journal of Chemical and Engineering Data, vol. 49, no. 2, pp. 352–356, 2004. View at Publisher · View at Google Scholar · View at Scopus
  22. M. L. Corazza, L. C. Filho, O. A. C. Antunes, and C. Dariva, “High pressure phase equilibria of the related substances in the limonene oxidation in supercritical CO2,” Journal of Chemical and Engineering Data, vol. 48, no. 2, pp. 354–358, 2003. View at Publisher · View at Google Scholar · View at Scopus
  23. S. N. V. K. Aki, A. M. Scurto, and J. F. Brennecke, “Ternary phase behavior of ionic liquid (IL)-organic-CO2 systems,” Industrial and Engineering Chemistry Research, vol. 45, no. 16, pp. 5574–5585, 2006. View at Publisher · View at Google Scholar · View at Scopus
  24. D. Y. Peng and D. B. Robinson, “A new two-constant equation of state,” Industrial and Engineering Chemistry Fundamentals, vol. 15, no. 1, pp. 59–64, 1976. View at Scopus
  25. J. A. Nelder and R. Mead, “A Simplex-method for function minimization,” Computer Journal, vol. 7, pp. 308–313, 1965.
  26. L. A. Blanchard and J. F. Brennecke, “Recovery of organic products from ionic liquids using supercritical carbon dioxide,” Industrial and Engineering Chemistry Research, vol. 40, no. 1, pp. 287–292, 2001. View at Scopus
  27. M. D. Bermejo, M. Montero, E. Saez et al., “Liquid-vapor equilibrium of the systems butylmethylimidazolium nitrate -CO2 and hydroxypropylmethylimidazolium nitrate -CO2 at high pressure: influence of water on the phase behavior,” Journal of Physical Chemistry B, vol. 112, no. 43, pp. 13532–13541, 2008. View at Publisher · View at Google Scholar
  28. S. N. V. K. Aki, B. R. Mellein, E. M. Saurer, and J. F. Brennecke, “High-pressure phase behavior of carbon dioxide with imidazolium-based ionic liquids,” Journal of Physical Chemistry B, vol. 108, no. 52, pp. 20355–20365, 2004. View at Publisher · View at Google Scholar · View at Scopus
  29. R. B. Gupta and J. J. Shin, Solubility in Supercritical Carbon Dioxide, CRC Press, New York, NY, USA, 2007.
  30. J. O. Valderrama and P. A. Robles, “Critical properties, normal boiling temperatures, and acentric factors of fifty ionic liquids,” Industrial and Engineering Chemistry Research, vol. 46, no. 4, pp. 1338–1344, 2007. View at Publisher · View at Google Scholar · View at Scopus
  31. K. Mishima, K. Matsuyama, H. Ishikawa, K. I. Hayashi, and S. Maeda, “Measurement and correlation of solubilities of azo dyes and anthraquinone in supercritical carbon dioxide,” Fluid Phase Equilibria, vol. 194–197, pp. 895–904, 2002. View at Publisher · View at Google Scholar · View at Scopus
  32. J. L. Anthony, E. J. Maginn, and J. F. Brennecke, “Solution thermodynamics of imidazolium-based ionic liquids and water,” Journal of Physical Chemistry B, vol. 105, no. 44, pp. 10942–10949, 2001. View at Publisher · View at Google Scholar · View at Scopus