1. Introduction
Tissue engineering is a rapidly expanding multidisciplinary scientific field with promising applications in biomedicine. Non-cancerous cells do not grow three-dimensionally; therefore, they need appropriate scaffolds to form three-dimensional tissues. The ideal scaffolds should allow the seeded cells to adhere and proliferate, enable the mobility of nutrients, and be made from biodegradable materials. Commonly used materials that satisfy the above requirements include polylactic acid (PLA), polyglycolic acid (PGA) and polycaprolactone (PCL) [
1,
2]. Scaffolds made from these materials often give satisfactory results; however, researchers worldwide are still seeking other materials suitable for particular applications. Polyurea copolymers may have optimal biocompatibility and the appropriate physical and chemical properties. Block polyurethanes could be promising candidates for nerve repair scaffolds [
3,
4]. Recently published results on their
in vivo and
in vitro applications have demonstrated satisfactory nerve regeneration with better performance than using autografts [
3,
4]. Nevertheless, the surface properties of such materials and their ability to adsorb and distribute liquids may not be optimal.
Cell adhesion and proliferation on scaffolds depends on numerous parameters, including the surface energy of the material used. Several authors have addressed the adhesion and multiplication of different biological cells on materials as a function of the material’s composition and structure [
5,
6,
7,
8,
9]. The reported results vary significantly; however, one robust conclusion is that the best results are obtained by using moderately hydrophilic substrates [
10,
11]. The majority of biocompatible polymers are moderately hydrophobic, often with water contact angles between 70° and 100°. PLA exhibits a contact angle in the range of 61° to 87° [
12,
13,
14,
15,
16], with a typical value of approximately 80° [
13,
14]; for PGA, a contact angle of 53° was reported [
17]; and for PCL fiber and PCL mesh, contact angles of 86° [
18] and 130° [
19,
20] were reported, respectively.
The surface properties of porous and/or fibrous polymers suitable for synthesizing scaffolds for tissue engineering can be tailored with different treatments; however, one of the best techniques is the application of gaseous plasma [
21]. Gaseous molecules partially ionize, dissociate and excite under plasma conditions [
22,
23], and the resulting reactive gaseous species interact chemically with the polymer surface to form new functional groups. The types of functional groups formed depends on the polymer’s properties, the type of gas and the plasma parameters. Polymers become more hydrophobic if they are functionalized with non-polar functional groups, usually by using fluorine-containing gases [
24]. Hydrophilization, on the other hand, is achieved by functionalization with polar groups, which is usually achieved by treatment with plasma using oxygen-containing gases [
25,
26,
27]. Because the required surface finish of materials used for scaffolds is moderately hydrophilic, oxygen plasma is often applied to polymeric scaffolds [
5,
6]. The plasma treatment may or may not lead to the preferred surface finish. Some authors report rapid functionalization as well as a change in surface wettability [
28], whereas others obtain rather unsatisfactory results. Poor functionalization with polar groups (or absence of these groups) upon treatment of fluorinated polymers by weakly ionized oxygen plasma is explained by extensive etching [
29,
30]: reactive oxygen species do interact chemically with fluorinated polymer but the molecular fragments containing oxygen desorb from the surface immediately. The functionalization of some other types of polymers often does not lead to satisfactory surface functionalization upon treatment with reactive oxygen species due to preferential etching [
31,
32]. The discrepancies in reported results may arise because most authors do not report the details of treatment conditions, but rather represent plasma as a black box. Hydrophobic recovery, which is typical in oxygen plasma treated polymers, is frequently neglected.
Although numerous authors have addressed the modification of the surface properties of materials that are suitable for tissue engineering using gaseous plasma treatment, there have only been a handful of reports on the functionalization of polyurethane and polyethylene glycol. Furthermore, to our knowledge, there is only one published paper in which polyurethane/urea copolymers based on polyethylene glycol were modified with plasma treatment [
33]. In that case, the polymer was treated with ammonia plasma; whereas in this paper, we report modifications of this type of polymer using oxygen plasma treatment while adjusting different treatment parameters.
3. Results and Discussion
The plasma in E-mode was characterized using optical emission spectroscopy, a catalytic probe and an electrical probe prior to experiments with polymer samples. Optical spectroscopy is a qualitative technique, whereas the probes give quantitative results. An optical spectrum of the plasma created at 33 Pa and a forward power of 150 W without a sample is shown in
Figure 4.
The spectrum reveals extensive emission arising from neutral oxygen and hydrogen atoms as well as a weak band in the UV region at 309 nm. This band corresponds to radiation from the excited OH radicals formed in the plasma as a result of the dissociation of water molecules caused by electron impact. The absence of nitrogen bands confirms that there is negligible leakage of the plasma reactor. The origin of the water is a slow desorption from internal reactor surfaces. The experimental setup does not allow for baking, therefore, the residual atmosphere in both the plasma and RGA chambers is predominantly water vapor. The concentration of water in the plasma reactor cannot be determined by mass spectrometry because the water molecules entering the RGA originate from both the plasma and RGA chambers. The water partial pressure can only be estimated by considering the base pressure in the plasma chamber, which is approximately 1 Pa. The partial pressure of water is therefore much smaller than that of oxygen. The optical spectrum presented in
Figure 4 indicates the high dissociation of both H
2O and O
2 molecules under plasma conditions. At a pressure of 33 Pa and forward power of 150 W, the absolute value of the O-atom density was determined with the catalytic probe to be approximately 1.2 × 10
21 m
−3. Such a high density is explained by extensive dissociation upon electron impact and the poor association of O atoms due to heterogeneous surface recombination. The association in the gas phase is negligible at 33 Pa because the reaction requires a three-body collision to preserve energy and momentum [
21,
34]. A forward power of 150 W does not allow for a high ionization fraction because the RF generator is mismatched. The reflected power meter on the generator showed a reflected power of 115 W; thus, the real discharge power was only approximately 35 W. The density of charged particles at these conditions was estimated using the double electrical probe to be only 4 × 10
16 m
−3.
PURPEG polymer samples were placed in the middle of the plasma reactor, directly under the collimating lens used for OES. The chamber was pumped down to a pressure of approximately 1 Pa and then filled with commercially available oxygen. The pressure was set to 33 Pa. Samples were treated with plasma generated with various forward RF power levels.
We attempted to treat samples in the powerful plasma created in H-mode at a forward power of approximately 900 W. The RF generator was matched well, and the reflected power was only 85 W. Such a powerful discharge created plasma of high luminosity. The integration time of the optical spectrometer was only 1 ms.
Figure 5 shows the optical spectrum acquired through treatment of the polymer sample at 900 W. The absence of oxygen atomic lines indicate a rather low electron temperature, which is typical for RF plasma in H-mode. The radiation from O atoms arises from excited states of excitation energy approximately 14 eV, thus it is highly unlikely that electrons of low temperature excite the radiative states. The optical emission was much larger for other plasma particles. Apart from the hydrogen Balmer series (H
α was left to saturation), there were emissions from CH, CN, OH and NH radicals as well as C
2 dimers. The presence of the dimers, CN radicals, and NH radicals indicates the thermal degradation of the polymer sample. If samples were oxidized, the predominant radiation would be from CO radicals. No mass spectrum was recorded during the experiment in H-mode because the sample visually degraded before the spectrum could be acquired. Therefore, sample treatment in H-mode at a real power close to 800 W caused thermal degradation rather than functionalization of the polymer sample. A similar thermal degradation of the PURPEG samples was also observed in H-mode ammonia plasma [
33].
Numerous samples were treated with oxygen plasma in E-mode. Typical OES spectra of oxygen plasma during the sample treatment in E-mode at a forward power of 150 W are shown in
Figure 6 and
Figure 7. The optical spectrum just after turning on the discharge was identical to the spectrum acquired in an empty chamber (
Figure 4). The spectrum in
Figure 4 comprises only a few spectral features: hydrogen Balmer lines, an OH band (A
2Σ–X
2Π transition) at 309 nm and oxygen atom lines at 777.4 and 844.6 nm corresponding to transitions 3p5P–3s5S° and 3p3P–3s3S°, respectively. The presence of oxygen emission lines is obvious because the chamber was filled with oxygen. However, the OH band and Balmer lines are observed because of the water vapor in the system, which is present due to the water desorption from reactor walls and the sample.
As the plasma treatment continued, some new spectral features emerged.
Figure 6 reveals the spectrum after 30 s of plasma treatment at a forward power of 150 W. Between 300 and 400 nm, a barely visible N
2 second positive band (C
3Π
u–B
3Π
g) indicates that the sample also contains nitrogen. All other subsequent spectral features are attributed to CO radicals, which can be explained by etching of the sample [
35,
36]. In the UV region of the spectra, there is a 3
rd positive band, which belongs to the transitions between b
3Σ
+ and a
3Π states. The central part of the spectrum from 350 to 800 nm is dominated by a broad continuum, which is known to be composed of overlapped bands of three CO systems: Angstrom (B
1Σ
+–A
1Π), Triplet (d
3Δ–a
3Π) and Asundi (a’
3Σ
+–a
3Π) bands [
37,
38]. The intensity of the CO bands in
Figure 6 is a rather weak, qualitative indication that the amount of reaction products (CO
2 and CO) is small in the initial stage of polymer processing.
Figure 7 represents the OES spectrum of plasma after a prolonged treatment of 600 s. The spectrum is similar to the one acquired at 30 s; however, the intensity of CO bands compared to O lines is much stronger. This observation is explained by the consumption of oxygen and accumulation of reaction products in the plasma chamber. Optical spectra were acquired continuously during sample treatment, and the intensity of the major features
versus treatment time is presented in
Figure 8. The intensity of the oxygen line at 777 nm and the H
α signal decreased with treatment time for the first few minutes and then stabilized, whereas the intensity of the CO bandhead at 519.5 nm (CO Angstrom line (0,2)) continued to increase. The behavior of the selected lines presented in
Figure 8 indicates the etching of the polymer during plasma treatment. However, no firm conclusions about the gas composition can be drawn because the technique is qualitative. Furthermore, the intensity of emission lines depends on the electron temperature, which in turn depends on the gas composition. Unfortunately, the home-made double electrical probe does not allow for reliable determination of the electron temperature.
The optical spectra were also recorded at different forward power levels of 100, 150 and 300 W. Considering the measured reflected power levels, these values correspond to real power levels of 25, 35 and 70 W, respectively.
Figure 9 presents the time evolution of the CO Angstrom line (0,2 vibrational transitions) normalized by the time evolution of the oxygen atom line (777.4 nm) for each RF power level. We conclude that the etching rate is larger when the forward RF power is higher. At higher power levels, the densities of both neutral atoms and ions are larger [
39], which explains the observed result.
The results of the OES measurements as summarized in
Figure 8 and
Figure 9 do not reveal the absolute concentration of stable gaseous molecules in the plasma chamber during polymer sample treatment. To quantify the results, we also performed measurements by mass spectrometry.
Figure 10 shows a mass spectrum that was measured in an empty chamber filled with oxygen at a pressure of 33 Pa, and the spectra shown in
Figure 11 and
Figure 12 were obtained in a chamber containing samples after 30 and 600 s of plasma treatment, respectively. The forward RF power was 150 W.
Peaks attributed to H
2, H
2O, O
2, CO
2 and CO can be observed in the mass spectra. The RGA spectrum recorded in an empty chamber contains only H
2, H
2O and O
2 peaks. Because of the cracking pattern of the RGA, molecules are not detected as only one mass peak. All of the detected peaks for a certain molecule are marked in
Figure 10,
Figure 11 and
Figure 12. For example, the peak at mass 16 is the sum of the contributions from the following stable molecules: H
2O, O
2, CO
2 and CO. Fortunately, the molecules H
2O, O
2 and CO
2 can be monitored by directly observing masses 18, 32 and 44, respectively. To observe CO molecules, which are produced during the plasma treatment of a sample, one must subtract 11.4% of the mass 44 peak intensity from mass 28. Thus, it is possible to quantify the evolution of the molecules’ partial pressures in the discharge chamber. The result is presented in
Figure 13. The results indicate behavior that is consistent with the observations from OES (
Figure 8). The O
2 mass peak decreases with time because oxygen is consumed during the oxidation of the sample. The products of sample etching are CO and CO
2, whose concentration increases with time. Although CO
2 reaches a maximum at approximately 250 s and then starts to decrease, CO continues to increase, thus indicating etching of the sample.
The etching of the polymer by plasma was also monitored by following the time evolution of the CO concentration.
Figure 14 represents the partial pressure of the CO molecule
versus treatment time at the power levels 100, 150, 200 and 300 W. The etching rate is higher if the sample is treated with plasma that is generated with a higher forward RF power. At 300 W, the concentration of CO reaches a maximum after a certain treatment time (approx. 350 s) and remains nearly constant with prolonged treatment. A feasible explanation for this observation is that all of the oxygen inside the discharge chamber has been consumed for polymer oxidation; therefore, CO can no longer be produced. After consuming the oxygen, the CO concentration slowly decreases as it is pumped out through the RGA chamber. This conclusion is supported by the time evolution of partial pressures at 300 W (
Figure 15). One can observe practically no remaining oxygen after approximately 350 s of treatment. Approximately 5% of the initial partial pressure of oxygen can easily be explained by the dissociation of CO in the discharge and recombination of O-atoms into O
2 on the surface of the discharge chamber.
The behavior of the oxygen partial pressure in
Figure 13 is rather unexpected and thus worth discussing. The oxygen partial pressure first increases with time, reaches a maximum after approximately 50 s of plasma treatment, and then decreases with further treatment. The increase in O
2 partial pressure during the first 50 s is probably an artifact of thermal effects. The discharge power initially heats the gas inside the plasma chamber. The real power (35 W) dissipated in the plasma reactor causes slow heating of the gas until such a temperature is achieved that heating is balanced by cooling with ambient room-temperature gas. The pressure in the plasma chamber depends on the density of molecules and gas temperature:
p =
nkBT. The oxygen partial pressure in
Figure 13 increased by approximately 10% from the initial value over a period of 50 s; therefore, one could estimate the gas temperature to be approximately 330 K at that treatment time. This gas temperature is an underestimate because oxygen is continuously being lost due to chemical reactions with the polymer sample (etching). A better estimate would account for this effect. The simplest method is to extrapolate the time derivative of the oxygen partial pressure towards zero, as shown in
Figure 16. The difference in the measured oxygen partial pressure and the cross-section of the extrapolation line with the y-axis is approximately 30%; thus, a better estimate of the gas temperature after a 50 s treatment time is 390 K. This value is typical for cold, non-equilibrium gas plasma.
Figure 16 can also be used to estimate the etching rate. As a rough approximation, one could assume that each oxygen molecule that is lost through polymer oxidation captures 3 polymer atoms (reaction –CH
2– + O
2→CO + H
2O). The true polymer composition is much more complex; therefore, the assumption is an oversimplification and underestimates the etching rate. However, it does allow us to make a rough estimate of the etching rate. The rate is calculated as follows:
where d
NO2/d
t is the slope of the line in
Figure 15,
xmon is the thickness of an atomic monolayer in the polymer,
S is the surface area of the sample and
Nmon is the surface density of atoms in the solid material. Considering the measured value (d
NO2/d
t = 2.4 × 10
16 s
−1 in
Figure 15) of
S = 2 × 10
−4 m
2, and typical values for
xmon = 3 × 10
−10 m and
Nmon = 1 × 10
19 m
−2, the etching rate is estimated to be 10 nm/s. This etching rate is typical for polymers and has been measured by much more sophisticated techniques, such as quartz crystal microbalance [
40].
It might be possible to estimate the etching rate in a simplified manner, i.e., by weighing a sample before and after plasma treatment. This procedure is common for smooth, inorganic materials that do not adsorb significant amounts of water. Our samples, however, are highly porous polymers; therefore, they adsorb significant amounts of water. The water molecules slowly desorb from the sample due to both vacuum conditions (equilibrium partial pressure of water vapor at 300 K is over 3000 Pa) as well as heating of the sample.
The behavior of the curve in
Figure 15 indicates linearity up to a treatment time of approximately 220 s (given a correction for thermal effects, see explanation above) and then a lower absolute value of the time derivative. The curve reflects the etching kinetics: up to approximately 220 s, the etching rate
ζ does not depend on time, and thereafter, it slowly decreases with increasing treatment time. The deviation from linearity is explained by a shortage of reactive oxygen species, which is due to the consumption of O
2 molecules and the formation of more stable CO
2 and CO molecules. The etching rate is therefore initially limited by the intensity of the surface reactions, and after approximately 220 s, by the arrival of reactive species onto the polymer surface.
The major reactants in the plasma created by these particular conditions are neutral oxygen atoms. They abound in the plasma at a density on the order of 10
21 m
−3. The corresponding flux of O-atoms on the surface is over 10
23 m
−2·s
−1; therefore, they should be capable of removing thousands of monolayers in a second. The positively charged ions are minor reactants, but they may contribute to etching acceleration due to synergistic effects with atomic oxygen and UV radiation from the plasma. The theoretical limit of the etching rate by O
2+ ions can be estimated from the measured density of charged particles. Because the Debye length in our plasma is orders of magnitude smaller than the lateral dimension of our samples, an approximation of infinitely large samples can be adopted. The ion flux onto the surface in this approximation is
j = ¼
n+<
v+>, where
n+ is the ion density in plasma and <
v+> is the ions’ average random velocity (almost 450 m/s). If all ions react chemically with the polymer, the etching rate would be:
The measured value is an order of magnitude larger; therefore, it is clear that the ions alone are not responsible for all of the polymer etching. The measured etching rate is thus explained by the interaction of PURPEG polymer with neutral reactive species (i.e., O-atoms in the ground state) and perhaps by synergistic effects with ions and/or UV radiation.
Whatever the mechanism involved, the oxygen plasma treatment not only causes etching of the polymer but also results in functionalization [
41,
42] that modifies the wettability of the material. The functionalization in our case was investigated by XPS. Survey spectra were measured for each sample soon after plasma treatment to minimize the influence of any ageing effects (hydrophobic recovery). All characterized samples were treated at a pressure of 33 Pa and a forward power of 150 W. The composition of the surface film as calculated from the XPS survey spectra is shown in
Figure 17. The results are not impressive. The original concentration of carbon is approximately 80 atomic %. It decreases monotonously, and the concentration after a prolonged treatment of 5 minutes is 65 atomic %. The opposite behavior is observed for oxygen and nitrogen. The nitrogen concentration is only up to 2 atomic %; therefore, considering the accuracy of this technique, the relative variations are not good representations of surface changes. Consequently, we limit our discussion to the evolution of the oxygen concentration. The oxygen concentration is approximately 17 atomic % for untreated samples, increasing to approximately 23 atomic % in the first 10 s of plasma treatment. The concentration continues increasing slowly until it reaches approximately 27 atomic % for the sample treated for 300 s. The survey spectra were measured at 3 spots on each sample, and the error bars in
Figure 17 indicate the statistical error only. The results presented in
Figure 17 therefore indicate that functionalization is not the predominant mechanism of interaction between reactive species from the oxygen plasma and the polymer surface. This type of polymer PURPEG does not exhibit much tendency to functionalize, which is explained by the instability of the oxidized surface material, which tends to desorb from the surface rather than accumulate much oxygen.
The types of newly formed functional groups can be deduced from XPS high resolution C1s peaks.
Figure 18 represents the evolution of this peak for different plasma treatment times. Three sub-peaks are observed. The largest peak at 284.8 eV corresponds to bonds between carbon atoms, the second peak at 286.4 eV to the single bonds of either C–N or C–O, and the small peak at 289 eV to –NH–CO–O– and –NH–CO–NH– groups in the original polymer structure. Because the nitrogen concentration is small and does not change with plasma treatment, it is clear that the modifications in the C1s peak are due to interactions with oxygen. Both peaks at 286.4 and 289 eV increase with plasma treatment time. There is no reason for the increased concentration of the peak at 289 eV; therefore, one can conclude that the increased peak at this binding energy corresponds to the formation of a O=C–O carboxylic or ester group. Unfortunately, the technique does not allow us to distinguish between these groups because the difference in binding energy between them is too small [
43,
44,
45]. The increasing intensity of the peak at 286.4 eV indicates that hydroxyl groups are also formed during plasma treatment.
Irrespective of the functional groups formed and despite the fact that functionalization is not reasonable, the plasma treatment caused significant modifications to the surface energy and the sorption kinetics. The water contact angle
versus plasma treatment time for the case of an oxygen pressure of 33 Pa and a forward power of 150 W is plotted in
Figure 19, with the corresponding water droplet absorption time in
Figure 20. The error bars in both figures indicate statistical error resulting from 5 repetitions.
The water contact angle (
Figure 19) roughly follows the concentration of oxygen on the film surface (
Figure 17). The contact angle for the untreated sample is approximately 85°, indicating the moderately hydrophobic character of as-synthesized PURPEG. The water contact angle drops to approximately 75° even with only 3 s of plasma treatment. Such a rapid effect is explained by a huge flux of O-atoms onto the sample surface. The segments of the polymer surface which were directly exposed to oxygen plasma quickly chemisorbed oxygen atoms and became saturated with polar functional groups. Further exposure causes etching, as revealed in
Figure 5,
Figure 6,
Figure 7,
Figure 8,
Figure 9,
Figure 10,
Figure 11,
Figure 12,
Figure 13,
Figure 14 and
Figure 15. As the treatment proceeds, the saturation of the polymer surface with oxygen functional groups slowly propagates into the pores. Once the entire surface in contact with the water droplet saturates, the contact angle assumes a constant value of approximately 30°, which is typical for this type of material. This value is much higher than for most aromatic polymers and comparable to the values for aliphatic polymers. Specifically, polymers with a high concentration of aromatic rings, such as polyethylene terephthalate or polystyrene, usually exhibit contact angles of approximately 10° when smooth and much smaller when nanostructured [
46,
47,
48]. PURPEG does not contain aromatic rings but just normal cyclohexane rings, which could explain the rather poor activation of the surface as revealed in
Figure 19. From application point of view, such a poor wettability may be beneficial because the best cell proliferation has been observed for moderately hydrophilic polymers [
10,
11]. Furthermore, a rather slow decrease in the contact angle with respect to treatment time allows us to select the desired wettability over a broad range of contact angles between 70° and 30°, corresponding to treatment times of 10 to 100 s.
Improved sorption of water droplets does not require the saturation of the entire material with functional groups.
Figure 20 reveals that the water droplet absorption time quickly drops from almost 9 s (for the untreated material) to 2 s (for samples exposed to plasma for 10 s). Further treatment does not have a significant influence on the adsorption time because even after 5 minutes of plasma treatment, the adsorption time remains at approximately 0.5 s.