Alberto
Arce†
,
Francesco
Fornasiero
,
Oscar
Rodríguez†
,
Clayton J.
Radke
and
John M.
Prausnitz
*
Chemical Engineering Department, University of California, Berkeley, and Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720. E-mail: prausnit@cchem.berkeley.edu; Fax: +1 510 642 4778; Tel: +1 510 642 3592
First published on 27th November 2003
Using a gravimetric technique, experimental data at 35°C are reported for isothermal sorption equilibria and sorption kinetics of water vapor in poly(N,N-dimethyl methacrylamide)
(PDMAA), poly(2-dimethyl aminoethyl methacrylate)
(PDMAEMA), poly(acrylic acid)
(PAA), and in a typical membrane of a commercial soft-contact lens made of poly(2-hydroxyethyl methacrylate)
(PHEMA). The highest sorption of water vapor is in PDMAA. The Flory–Huggins model and the Zimm–Lundberg clustering theory are used to interpret the equilibrium data. Over the entire range of water activity, the four materials show clustering functions larger than −1 indicating that water molecules cluster together in the polymers. The least hydrophilic polymer, PDMAEMA, and the most hydrophilic one, PDMAA, show respectively, the highest and the lowest tendency of water molecules to form clusters. Water diffusion into the polymer matrix is faster in PDMAEMA (diffusion coefficient, D
=
10–20
×
10−8 cm2 s−1) and in PDMAA (D
=
4.9–8.7
×
10−8 cm2 s−1) than in PHEMA-lens material (D
=
0.55–3.4
×
10−8 cm2 s−1) and PAA (D
=
0.98–3.5
×
10−8 cm2 s−1). The measured diffusion coefficients increase with water activity in the PHEMA-lens, decrease with activity in PDMAEMA, and show a more complex concentration dependence in PDMAA and PAA. By writing the diffusion coefficient as a product of an intrinsic mobility, Đ, and a non-ideality thermodynamic factor, Γ, the concentration dependence of D is explained by the interplay of two factors: rising plasticization with water activity, which causes an increase in Đ, and the effect of composition on Γ, decreasing with water activity.
This work reports experimental, isothermal vapor-sorption equilibria and diffusion coefficients of water in three homopolymer films whose monomers may be included in the formulation of continuous-wear contact lenses, and in a typical membrane used in commercial soft contact lenses made of poly(2-hydroxyethyl methacrylate)
(PHEMA). The homopolymers are poly(N,N-dimethyl methacrylamide)
(PDMAA), poly(2-dimethyl aminoethyl methacrylate)
(PDMAEMA), and poly(acrylic acid)
(PAA). Hydrophilic polymers are used in contact-lens technology because adequate water content in a contact lens is essential for eye comfort.5 Experiments are carried out at 35°C, the average temperature of a contact lens on the eye, in a quartz-spring gravimetric apparatus.6 We use a well-defined, thin-film polymer-sample geometry and measure spring elongation as a function of time for different water-vapor activities. For each solvent-polymer system, the equilibrium data are used to calculate the binary Flory parameter.7 Toward obtaining some information on the molecular state of dissolved water, we use the clustering theory of Zimm and Lundberg.8,9 Water vapor diffusion coefficients are calculated from a diffusion model derived from Fick's second law. Although similar data have been reported for PAA at different temperatures,10 to our knowledge, no previous data have been reported for PDMAA and PDMAEMA.
Polymer | Supplier | CAS No./Serial number | MW | ρ/g cm−3 | T g/°C |
---|---|---|---|---|---|
a Estimated. b SP2: Scientific Polymer Products, Inc. c Reported by supplier. d Equilibrium water content in pure water of 38 weight % (reported by Bausch & Lomb). | |||||
PDMAA | Polymer Source | 26793-34-0 | 4000 | 1.15a | 105 |
PDMAEMA | SP2b | 25154-86-3 | 100![]() |
1.10a | 18 |
PAA | Aldrich | 9003-01-4 | 250![]() |
1.15c | 108 |
PHEMA-lens materiald | Bausch & Lomb | 2248-SH-140A | — | 1.15c | 118 |
LC![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() | (1) |
![]() | ||
Fig. 1 Schematic of the apparatus for water-vapor sorption and diffusion measurements. |
A U-tube manometer is used to measure the pressure. The manometer liquid is 1-octadecene, selected because of its low vapor pressure, immiscibility with water, suitable density for the pressure range of interest here, and for a suitable U-tube length. Its density is calculated with an equation given by Daubert and Danner.11
The spring elongation and the manometer-liquid-height difference are measured with a cathetometer (Wild, model KM274) with a precision of ±0.1 mm. The equipment is enclosed in an air bath with an Omega CN9000 PID controller (accuracy ±0.5°C) to keep the experimental temperature very close to 35
°C.
After vacuum testing is satisfied, water vapor is injected by opening the solvent-flask valve. Spring elongation and pressure are measured at several time steps until equilibrium is reached, i.e. when several successive measurements agree within experimental error. Then a new injection is made. The experiments correspond to a range of vapor activities, aw, from 15 to 85%. Water activity is calculated assuming an ideal gas phase of pure water vapor, neglecting the tiny vapor pressure of 1-octadecene:
![]() | (2) |
PDMAA | PDMAEMA | PAA | PHEMA-lens material | ||||
---|---|---|---|---|---|---|---|
a w | w w | a w | w w | a w | w w | a w | w w |
0.14 | 0.025 | 0.14 | 0.004 | 0.29 | 0.037 | 0.19 | 0.015 |
0.24 | 0.053 | 0.24 | 0.008 | 0.48 | 0.086 | 0.40 | 0.038 |
0.38 | 0.095 | 0.38 | 0.019 | 0.59 | 0.138 | 0.66 | 0.083 |
0.49 | 0.136 | 0.49 | 0.032 | 0.69 | 0.200 | 0.72 | 0.102 |
0.58 | 0.178 | 0.58 | 0.048 | 0.80 | 0.292 | 0.78 | 0.126 |
0.68 | 0.235 | 0.68 | 0.076 | 0.86 | 0.395 | 0.84 | 0.154 |
0.76 | 0.293 | 0.76 | 0.112 | ||||
0.83 | 0.372 | 0.83 | 0.174 |
The equilibrium isotherms show a large degree of upturn at high activities (Fig. 2). Two phenomena are suggested to interpret the equilibrium isotherm upturn: (1) plasticization of the polymer by water, and/or (2) clustering of water molecules. These phenomena provide opposite effects for the concentration dependence of the diffusion coefficient.14–16
The Flory–Huggins theory7 is widely used to correlate experimental phase equilibria for solvent-polymer systems:
![]() | (3) |
![]() | ||
Fig. 3 Flory interaction parameter χversus water-vapor activity for PDMAA, PDMAEMA, PAA, and PHEMA-lens material at 35![]() |
Water molecules have a tendency to form clusters when absorbed in a polymer.18 To provide a measure for clustering, Zimm and Lundberg8,9 defined the clustering integral Gww. The ratio of the clustering integral to the partial molecular volume Gww/vw, called the clustering function, is obtained from equilibrium data:
![]() | (4) |
![]() | ||
Fig. 4 Clustering function Gww/vw at various water-vapor activities for PDMAA, PDMAEMA, PAA, and PHEMA-lens material at 35![]() |
The quantity ΦwGww/vw has been interpreted as the mean number of solvent molecules in the neighborhood of a given solvent molecule in excess of those provided by the mean solvent concentration.8,9 The sum (1+
ΦwGww/vw), therefore, is the mean cluster size.22Fig. 5 shows ΦwGww/vw as a function of water vapor activity for four polymers. PDMAEMA shows the highest value of ΦwGww/vw that increases rapidly with water activity, reaching values in excess of 2 at high activities. In PDMAA, the mean cluster size is smaller with ΦwGww/vw always below 0.5. For PAA, ΦwGww/vw reaches a maximum near 1 at aw
∼
70% and is larger than ΦwGww/vw for PHEMA-lens material until near 80% activity. Then, ΦwGww/vw for PHEMA-lens material increases to 1.2 at the highest activity. These results for PHEMA-lens material are similar to those reported elsewhere13 for water in PHEMA homopolymer.
![]() | ||
Fig. 5 Mean excess of water molecules around a central one, ΦwGww/vw, at various water-vapor activities for PDMAA, PDMAEMA, PAA, and PHEMA-lens material at 35![]() |
![]() | ||
Fig. 6 Sorption kinetics of water in polymers at 35![]() |
Diffusion coefficients, D, are calculated using a model suggested by Crank23 to interpret kinetic absorption data “when the surface concentration changes rapidly but not instantaneously, a situation which usually arises when an instantaneous change is attempted in an experiment.”
We assume that at a fixed water activity, at the start, the water surface concentration approaches the equilibrium water concentration, C0, as an exponential function of time:
C(x![]() ![]() ![]() ![]() ![]() ![]() | (5) |
![]() | (6) |
A sigmoidal shape in experimental diffusion data can also indicate anomalous kinetics. We cannot exclude completely this possibility. However, for our systems, only the very beginning would be eventually affected. Moreover, eqn. (6) can be considered a simplified version of the model proposed by Long and Richman24 to interpret anomalous sorption kinetics. When the instantaneous surface concentration C0 in eqn. (2) of ref. 24 is set equal to zero, the Long and Richman model reduces to eqn. (6).
Therefore, whether this sigmoidal shape is due to some non-Fickian effect or to the impossibility to raise the pressure instantaneously (as we believe in absence of stronger evidence of anomalous kinetics), the adopted model captures the essential features of both phenomena.
Fig. 6 shows water-sorption kinetics in the four polymers for a similar water activity. Points show experimental data whereas lines represent fitting curves obtained with eqn. (6). Table 3 reports the fitted coefficients D and β at each water activity together with the film thickness for each polymer. We use the Sigmaplot© program with a least-squares method to calculate D and β with eqn. (6). In general, the diffusion coefficient is concentration dependent; thus, the diffusion coefficient reported here are averages for a concentration range of each sorption step, from an initial water concentration to the equilibrium concentration at the corresponding water activity.
Polymer | a w | 108D/cm2 s−1 | 108Đ/cm2 s−1 | β −1/s |
---|---|---|---|---|
PDMAA | 0.17 | 7.1 | 8.0 | 0 |
L![]() ![]() |
0.39 | 4.9 | 6.5 | 0 |
0.58 | 8.7 | 14 | 5 | |
0.77 | 8.1 | 17 | 122 | |
PDMAEMA | 0.21 | 20 | 35 | 80 |
L![]() ![]() |
0.40 | 14 | 29 | 81 |
0.61 | 13 | 36 | 1 | |
0.76 | 10 | 29 | 9 | |
PAA | 0.58 | 0.98 | 2.3 | 27 |
L![]() ![]() |
0.70 | 3.5 | 9.0 | 196 |
0.80 | 2.9 | 8.2 | 167 | |
PHEMA-lens material | 0.19 | 0.62 | 0.72 | 127 |
L![]() ![]() |
0.40 | 0.55 | 0.78 | 8 |
0.78 | 3.2 | 7.6 | 8 | |
0.84 | 3.4 | 8.9 | 208 |
The inverse of β measures the characteristic time required to reach the equilibrium water content at the polymer/water-vapor interface. In all cases, 1/β is less than 3.5 min, much smaller than the diffusion time scale, L2/4D, and consistent with the time required for increasing the pressure at the desired level during the experiment. Therefore, the water concentration on the film surface reaches equilibrium after at most the first few measurements of water uptake in the diffusion experiments and, except for a few initial points (2–3), the experimental sorption curves have Fickian characteristics.
In general, the water diffusion coefficient in these four materials is concentration dependent. D in PAA reaches a maximum of 3.5×
10−8 cm2 s−1 at 69% water activity. Chang et al.10 report diffusion coefficients for water in PAA at 40
°C reaching a maximum of 1
×
10−8 cm2 s−1 at 60% activity. For PDMAA, D tends to increase at high activities. However, it has a minimum of 4.9
×
10−8 cm2 s−1 at 58% activity. Diffusion coefficients in PDMAA are larger than reported by Ichikawa et al.25
(0.98
×
10−8 cm2 s−1) at aw
=
50%. The calculated diffusion coefficient for water in the PHEMA-lens material increases with water activity from a constant 0.6
×
10−8 cm2 s−1 at lower activities to 3.4
×
10−8 cm2 s−1 at higher activities. These values are similar to those obtained for PAA. They are in good agreement with those reported previously for PHEMA homopolymers.13,25 The water diffusion rate is faster in PDMAEMA than in the other three materials. D decreases with water activity from 20
×
10−8 cm2 s−1 to 10
×
10−8 cm2 s−1 within the measured activity range.
For binary systems, the Fickian diffusion coefficient can be written as the product of two terms, D=
ΓĐ;26 the first one is a thermodynamic non-ideality factor, Γ, defined by
![]() | (7) |
The non-ideality factor can be calculated from the clustering function:
![]() | (8) |
Fig. 7 shows that the thermodynamic factor is a decreasing function of water activity for our four systems. The same trend is observed for the diffusion coefficient in PDMAEMA. Therefore, the concentration dependence of the water diffusion coefficient reflects the concentration dependence of the thermodynamic factor for PDMAEMA.
![]() | ||
Fig. 7 Thermodynamic factor Γversus water-vapor activity for PDMAA, PDMAEMA, PAA, and PHEMA-lens material at 35![]() |
D has the opposite trend in PHEMA-lens material, increasing at larger water activity, whereas the concentration dependence of the diffusion coefficient is more complex for PDMAA and PAA. An increasing diffusion coefficient is often associated with polymer plasticization because it augments the intrinsic mobility Đ. Therefore, while plasticization appears to dominate the concentration dependence of D in PHEMA-lens, the complex behavior observed in PAA and PDMAA is probably a result of the competing concentration dependence of Đ and Γ.
The mobility Đ increases from 6.5–8×
10−8 at low water activity to 14–17
×
10−8 cm2 s−1 at high water activity for PDMAA, and from 0.55–0.62
×
10−8 to 3.2–3.4
×
10−8 cm2 s−1 for PHEMA-lens material. Despite some scatter in the reported values (easily accounted for by the possible source of error in the calculation), the mobility is nearly constant (29–36
×
10−8 cm2 s−1) with water activity in PDMAEMA, whereas for PAA, it rises sharply from 2.3 to 9.0
×
10−8 cm2 s−1 in a relatively narrow activity range.
The observed mobility trends are probably due to polymer plasticization. To support this conclusion, we have estimated the water uptake required for lowering the glass-transition temperature below 35°C for PDMAA and PAA. We used the Flory–Fox28 equation with a glass-transition temperature of water Tgw equal to −135
°C:
![]() | (9) |
The required water content is, respectively, 13.1 and 13.5 weight %. For PHEMA-lens material, we expect that the glassy-to-rubbery transition occurs at 35°C for a water content similar to that measured for HEMA homopolymers, i.e. around 8 weight %.29 Therefore, all three polymers undergo a glassy-to-rubbery transition during equilibrium sorption experiments, roughly at water activity equal to 50%, 60% and 65% respectively. Also, some of the kinetic water-uptake experiments are performed in the glassy region (the first and second activity for PDMAA; the first and second activity for PHEMA-lens) some other experiments span a water-content range embedding the glassy-to-rubbery transition (most likely, the third diffusion step for PDMAA, the first one for PAA and the third one for PHEMA-lens material), whereas the remaining measurements are in the rubbery region.
While the equilibrium isotherms do not present noticeable features related to the glassy-to-rubbery transition (for example, inversion of curvature as reported for sorption of gases in polymers), there is a clear correlation between the activity region where the kinetic sorption experiment is performed and the mobilities. In particular, the calculated mobilities for PDMAA, PAA and PHEMA-lens material are much larger in the rubbery region than those in the glassy region; mobilities increase with water content in the polymer probably because of plasticization induced by water sorption.
By contrast, PDMAEMA mobility is roughly constant with water content. PDMAEMA has a glass-transition temperature well below 35°C; the polymer is rubbery in the entire composition range at the experimental temperature and, therefore, there is no significant plasticization accompanying water sorption.
Footnote |
† Permanent address: Department of Chemical Engineering, University of Santiago de Compostela, 15782 Santiago de Compostela, Spain. |
This journal is © the Owner Societies 2004 |