Quantification of Water-Ion Pair Interactions in Polyelectrolyte Multilayers Using a Quartz Crystal Microbalance Method.

Water existing within thin polyelectrolyte multilayer (PEM) films has significant influence on their physical, chemical, and thermal properties, having implications for applications including energy storage, smart coatings, and biomedical systems. Ionic strength, salt type, and terminating layer are known to influence PEM swelling. However, knowledge of water's microenvironment within a PEM, whether that water is affiliated with intrinsic or extrinsic ion pairs, remains lacking. Here, we examine the influence of both assembly and post-assembly conditions on the water-ion pair interactions of poly(styrene sulfonate)/poly(diallyldimethylammonium) (PSS/PDADMA) PEMs in NaCl and KBr. This is accomplished by developing a methodology in which quartz crystal microbalance with dissipation monitoring is applied to estimate the number of water molecules affiliated with an ion pair (i), as well as the hydration coefficient, πsalt H2O. PSS/PDADMA PEMs are assembled in varying ionic strengths of either NaCl and KBr and then exposed post-assembly to increasing ionic strengths of matching salt type. A linear relationship between the total amount of water per intrinsic ion pair and the post-assembly salt concentration was obtained at post-assembly salt concentrations >0.5 M, yielding estimates for both i and πsalt H2O. We observe higher values of i and πsalt H2O in KBr-assembled PEMs due to KBr being more effective in doping the assembly because of KBr's more chaotropic nature as compared to NaCl. Lastly, when PSS is the terminating layer, i decreases in value due to PSS's hydrophobic nature. Classical and ab initio molecular dynamics provide a microstructural view as to how NaCl and KBr interact with individual polyelectrolytes and the involved water shells. Put together, this study provides further insight into the understanding of existing water microenvironments in PEMs and the effects of both assembly and post-assembly conditions.

Introduction

Polyelectrolyte complexes (PECs) prepared into thin films are called polyelectrolyte multilayers (PEMs).[1−4] Assembly occurs when polycations and polyanions form intrinsic ion pairs, leaving other uncompensated charge sites available for extrinsic ion compensation by smaller counterions.[1,5,6] These PEMs may be applied in electrochemical systems, drug delivery platforms, smart coatings, and water and gas barriers.[7−11] They are highly tunable and responsive to changes in assembly conditions (e.g., polymer concentration, salt type and concentration, pH, and temperature) and also post-assembly environment (e.g., temperature, humidity, and ionic strength).[4,12−16] Thus, to ensure the integrity of a PEM in its final application, it is important to understand its stability from assembly to end use in response to external stimuli.

PECs and PEMs have been described as “saloplastics” due to the strong plasticizing effect of salt, which reduces the number of intrinsic ion pairs among polyelectrolyte chains.[17,18] This is often described as doping.[19,20] The effects of increasing salt content on a macroscopic level are evident in PEM growth profiles (linear to exponential growth), mechanical properties (rigid to soft), layer structure (stratified to interdigitated layers), surface roughness (rough to smooth), and glass transition temperature (Tg).[18,21−24] Interestingly, post-assembly exposure of the PEM to salt may either result in swelling, deswelling, or even disassembly of the PEM depending on the ionic strength.[25−27]

Similar to salt, water acts as a plasticizer in PEMs by increasing the free volume for polymer mobility and weakening electrostatic attraction between oppositely charged polyelectrolytes.[28,29] Our past studies have shown a universal dependence of key PEM properties, such as glass transition temperature (Tg) and relaxation time (τ), on water content for both strong and weak polyelectrolytes at varying pH and ionic strength.[16,22,30] Specifically, the inverse of Tg was proportional to the ratio of total water molecules to total number of intrinsic ion pairs,[22,30] pointing to the significance of water’s association with the polycation–polyanion intrinsic ion pair.

A number of studies classify the total water content in PEMs into smaller microenvironments.[22,31−33] In our recent work, by observing the hydrogen bonding strengths of the deconvoluted OD stretch peak obtained for poly(styrene sulfonate)/poly(diallyldimethylammonium) (PSS/PDADMA) PEMs using attenuated total reflectance Fourier transform infrared (ATR-FTIR) spectroscopy, we identified three water microenvironments around ion pairs: high frequency (tightly bound) water, low frequency (loosely bound) water, and bulk (free) water.[31] A prominent finding from the study is the lack of any bulk free water for both partially and fully immersed PEMs, indicating that all water molecules present in the PEMs are associated with and influenced by either intrinsic or extrinsic ion pairs.[31] Although the nature of the water microenvironment was examined, the study could not distinguish whether water was located at extrinsic or intrinsic ion pairs.

Experimentally, the Schlenoff group developed a method to distinguish between intrinsic and extrinsic water–ion pair interactions using ATR-FTIR spectroscopy.[27] Two polycations, PDADMA, and poly(4-vinylmethylpyridinium bromide) (P4VMP) were assembled with PSS at 1.0 M NaCl, annealed in 1.0 M NaCl for a week, and then exposed to varying concentrations of a range of 14 different salts. The purpose of the annealing step was to produce PEMs with a 1:1 polycation:polyanion stoichiometric ratio. The number of total water molecules per intrinsic ion pair, rH, was determined from the areas of the OH and SO3– peaks and plotted against the post-assembly salt concentration. By extrapolating the linear portion of the resulting curve to the y-axis, the number of water molecules associated with intrinsic ion pairs for an undoped PEM, “i”, was determined by the intercept and a hydration coefficient, πsaltH, determined by the slope. πsaltH represents the ability of a salt to influence an influx of water into a PEM. For PSS/PDADMA or PSS/P4VMP PEMs, the values of i were 6.9 ± 1.7 or 2.5 ± 0.6 water molecules per intrinsic ion pair, respectively, showing that P4VMP is a more hydrophobic polycation. It was observed that i remains constant irrespective of the post-assembly salt whereas πsaltH is dependent on the post-assembly salt type. This study serves as motivation to broaden the available methodology to determine i and the factors influencing i. In our case, we employ quartz crystal microbalance with dissipation (QCM-D) monitoring, which provides a method to directly measure mass changes with both assembly and post-assembly conditions, leading to an alternative method to assess water at the intrinsic ion pairs.

QCM-D over the years has become a relevant tool for studying the response of polymers and composites to changes in external environment. QCM-D, which operates on the converse-piezoelectric effect of quartz crystal, serves as an effective way to measure real-time changes in PEM thickness, mass, and viscoelasticity.[34,35] From our group, past works of Reid et al. involved the study of the reversible swelling behavior of PSS/PDADMA PEMs prepared at 0.5 M NaCl upon post-assembly exposure to different concentrations of varying divalent salt solutions.[36] Also, O’Neal et al. studied the swelling–deswelling behavior and the compositional changes of PSS/PDADMA PEMs prepared at 0.5 M NaCl upon exposure to different concentrations of monovalent salts (KBr, NaBr, NaCl, and KCl).[25] None of these studies quantitatively examined the water content or the water’s microenvironment with regard to intrinsic or extrinsic ion pairs.

In this study, we estimate the number of water molecules associated with intrinsic ion pairs in PSS/PDADMA PEMs prepared in NaCl and KBr using QCM-D. We monitor both the layer-by-layer (LbL) assembly and the post-assembly swelling or deswelling in solutions of different salt concentrations in real time. Molecular dynamics (MD) simulations provide a microscopic viewpoint of water at the intrinsic pair. The goal is to develop a method by which QCM-D can estimate the number of water molecules affiliated with an ion pair, as well as the hydration coefficient, i and πsaltH, respectively. Once obtained, we discuss how i and πsaltH are influenced by salt identity and the PEM’s terminating layer.

Experimental Section

Materials

Poly(diallyldimethylammonium) (PDADMA, Mw = 200,000–350,000 g/mol, 20 wt % solution), poly(styrene sulfonate) (PSS, Mw = 500,000 g/mol), and linear polyethylenimine (LPEI, Mw = 25,000 g/mol) were purchased from Sigma-Aldrich, Scientific Polymer Products, and Polysciences, Inc., respectively. Sodium chloride (NaCl) and potassium bromide (KBr) were purchased from Sigma-Aldrich and Alfa Aesar, respectively. Silicon dioxide-coated QSensors (QSX 303 SiO2) QCM-D substrates were purchased from Biolin Scientific.

Preparation of Freestanding Polyelectrolyte Multilayers (PEMs)

An automated Carl Zeiss HMS slide stainer was used to fabricate free-standing PSS/PDADMA PEMs on Teflon substrates. First, PDADMA and PSS solutions were prepared at 1 g/L concentration and at identical NaCl or KBr ionic strength. The ionic strength for each salt was varied from 0.25–1.0 M (0.25, 0.5, 0.75, and 1.0 M). Eight samples of PSS/PDADMA PEMs were prepared: four in the presence of NaCl and four in the presence of KBr. The Teflon substrates were sonicated in ethanol for 15 min followed by thorough rinsing with Milli-Q water. The first layer was fabricated by dipping the substrates in PDADMA solution of a certain ionic strength for 15 min, followed by three separate rinse steps for 2, 1, and 1 min in Milli-Q water at matching ionic strength. The same process was repeated for the second layer by replacing the PDADMA solution with a PSS solution to form a layer pair. A total of 140 layer pairs were prepared with a final rinse step in pure Milli-Q water for 5 s to remove any excess salt. The PEMs were dried under ambient conditions overnight and then dried under vacuum at 115 °C for 3 h. The completed PEMs are denoted as (PSS/PDADMA)140.

Quartz Crystal Microbalance with Dissipation (QCM-D) Monitoring

Quartz crystal microbalance with dissipation (QCM-D) monitoring was used to monitor both LbL assembly and post-assembly behavior using a QSense E4 instrument. The PEMs were assembled on SiO2-coated AT-cut quartz crystals with a resonant frequency of 4.95 MHz. The quartz crystals were cleaned by immersing them in 2 vol % sodium dodecyl sulfate (SDS) for 10 min, rinsing with Milli-Q water, drying with nitrogen, and one last step of 15 min O2-plasma treatment. PDADMA and PSS polyelectrolyte solutions were prepared at a concentration of 0.1 g/L at varying ionic strengths of NaCl and KBr (0.25–1.0 M) similar to that of the freestanding films. Rinse solutions were set at a matching ionic strength to that of the polyelectrolyte solutions. LPEI solution was adjusted to pH 5.5 to obtain a stable clear solution. All QCM-D experiments were carried out at room temperature. All polyelectrolyte and rinse solutions were flowed through a peristaltic pump at a constant approximate flowrate of 150 μL/min. Milli-Q water (pH 5.5) was first allowed to flow over the quartz crystal for 30 min as a baseline for each measurement. An anchor layer of LPEI was then deposited onto the crystal for 10 min before rinsing for 5 min with pure Milli-Q water. The rest of the LbL film followed an alternating deposition of PSS and PDADMA solutions for 10 min each separated by a 5 min rinse at a matching ionic strength.

The final number of layers is designated “x” in (PSS/PDADMA). For PEMs prepared in NaCl, x = 6, and for those prepared in KBr, x = 4 in order for the film thicknesses to stay within the QCM-D operational limits.[34] This is because the QCM-D guide suggests that this thickness limit exists at approximately 1 μm, beyond which overtone signals begin to become unstable, starting with the highest overtone (13th overtone).[34,35] In order to accurately model the resulting QCM-D data for a viscoelastic film, a minimum of 3 overtones are required. We, therefore, aimed to have all PEM thicknesses below 700 nm during both assembly and post-assembly experiments. Also, some PEMs were terminated with either PDADMA for (PSS/PDADMA) or PSS for (PSS/PDADMA) .

Post-assembly, PEMs were exposed to solutions of increasing salt concentrations (0–2.0 M for NaCl and 0–1.0 M for KBr) with an exposure time of 30 min for each salt. The low salt region (0–0.1 M) was in some cases eliminated if the corresponding swelling thickness went beyond the permissible thickness range. In such cases, the lowest salt concentration explored were 0.001 M NaCl and 0.1 M KBr. At the end of all QCM-D measurements, each sample was rinsed in pure Milli-Q water to remove any excess salt on the crystal.

Determining the Salt Hydration Coefficient, πsaltH, and the Number of Water Molecules Associated with Intrinsic Ion Pairs, i

From the raw QCM-D data, the hydrated thickness of the PEM was calculated using a viscoelastic model under each exposure condition. The PEMs were then dried, and the dry thickness of the film was measured using profilometry (see the SI). Using the reported density of dry PSS/PDADMA PEMs (ρdry = 1.27 kg/m3), the dry mass of the PEMs was calculated from its dry thickness.[20,37] The water content of the PEM at a given state was obtained by subtracting the dry PEM mass from the hydrated mass. Dry, freestanding PEMs were made separately, and a representative composition was obtained using NMR spectroscopy and XPS (see the SI). Assuming that all of the limiting polyelectrolyte in each sample engages in intrinsic charge compensation only, rH was estimated from the composition and the PEM’s water content. The combination of all the experimental procedures outlined above makes it possible to directly estimate the hydration coefficient, πsaltH, and the number of water molecules associated with intrinsic ion pairs in an undoped complex, i. Ref (27) provides a derivation toward these two parameters. πsaltH is defined by the equation below:where [MA] is the salt concentration and rH is the ratio of moles of water in the PEM to moles of PSS monomer units in the PEM. In this study, non-stoichiometric PEMs with PDADMA in excess were prepared; thus, PSS represents the limiting polyelectrolyte in the PEM. We therefore assume that all PSS chains in the PEM form intrinsic ion pairs. By rearranging eq , a linear relationship between rH and [MA] is obtained, in which the slope is equal to πsaltH and in which the y intercept is equal to i:

For most polyelectrolyte systems studied, a deviation from this linear relationship between rH and [MA] is observed at lower salt concentrations (<0.1 M) due to electrostatic repulsion within the PEMs causing increased swelling. Another explanation is that this may occur as a result of osmotic pressure, which brings in additional water into the multilayer films at low salt concentrations. Therefore, only the linear portion of the curve (existing at higher salt concentrations) was extrapolated to the y axis to determine both i and πsaltH.

Simulation Methodology

The simulation protocol and methodology for both MD and ab initio molecular dynamics (AIMD) simulations are detailed in the Supporting Information.

Results

Layer-by Layer-Assembly

QCM-D was applied to monitor the growth and the swelling responses of the (PSS/PDADMA) PEMs in real time. First, a layer of LPEI at pH 5.5 was deposited onto a bare SiO2-coated QCM-D sensor to promote subsequent layer growth. Following this, LbL assembly of PSS and PDADMA was conducted at varying NaCl or KBr concentrations (0, 0.25, 0.5, 0.75, 1.0 M). After assembly, the PEM was exposed to solutions of successively higher salt concentrations and the swelling response was examined. Figure shows examples of the raw data obtained from this procedure for (PSS/PDADMA) PEMs prepared at 0.5 M NaCl (Figure a) and 0.5 M KBr (Figure b) that were then exposed to solutions of NaCl (0–2 M) and KBr (0–1 M), respectively. Generally, a negative frequency change (Δf < 0) corresponds to an increase in hydrated mass and thickness, and a positive dissipation change (ΔD > 0) corresponds to a softening of the film. In Figure , we observed negative changes in frequency with each polyelectrolyte deposition step, confirming an increase in film thickness. For the PDADMA adsorption step, specifically, we observed a more prominent decrease in frequency relative to the PSS step. Rinsing showed small, but noticeable, positive deviations in frequency that we associate with the removal of loosely bound material. At the later stages of film growth, a larger drop in frequency with PDADMA deposition resulted, pointing to an exponential growth trend; for PSS deposition at later stages, frequency did not change much, also consistent with prior reports of exponential growth.[4,24,38−41] Other assembly salt concentrations explored herein generally followed the same responses shown in Figure but to varying magnitudes. The post-assembly swelling responses of the PEMs to solutions of varying salt concentrations and types will be discussed later below.

Figure 1

Δf and ΔD versus time for the 3rd, 5th, and 7th overtones (n) of the QCM-D response from the layer-by-layer deposition and post-assembly salt treatment of (a) (PSS/PDADMA)6 PEMs prepared at 0.5 M NaCl, post-assembly salt treatment in NaCl; (b) (PSS/PDADMA)4 PEMs prepared at 0.5 M KBr, post-assembly salt treatment in KBr. (c,d) Resulting hydrated thickness (e,f) and viscosity from viscoelastic modeling of the QCM-D response in (a) and (b), respectively.

Figure c–f shows the resulting modeled film thickness and viscosity. We observed that the film viscosity increased with film thickness as more polymer layers were adsorbed during assembly. Hence, during assembly, each deposition step caused an increase in the amount of adsorbed and hydrated polymer, leading to an increase in the viscous component of the PEM. However, during post-assembly monitoring, when exposed to pure water, an influx of water molecules causes swelling and reduced the viscosity. Likewise, with the introduction of salt ions to the PEMs post-assembly, as the PEMs either swelled or deswelled, the viscosity either decreased or increased, respectively.

Figure shows the growth profile of the PSS/PDADMA PEMs extracted from data similar to that shown in Figure . For all samples prepared, we observed a linear growth regime at early stages, after which a non-linear exponential (or parabolic) growth regime manifested. This growth trend mirrors the large frequency drops that were observed upon the deposition of PDADMA at later stages, Figure a,b. Comparing Figure a to 2b, it becomes clear that PEMs assembled in KBr have a thickness larger than those assembled in NaCl for an equivalent number of layers. Quantitatively, for x = 4 layer pairs (8 layers) and an assembly salt concentration of 1.0 M, the average hydrated film thickness was 130 nm for NaCl-assembled PEMs and 200 nm for KBr-assembled PEMs. For comparison, Figure S1 shows the final dry thickness of PEMs measured using profilometry for each of the assembly conditions, in which the final dry thickness increased as the assembly salt concentration increased.

Figure 2

Growth profiles of PSS/PDADMA PEMs assembled in (a) NaCl and (b) KBr at varying concentrations. Hydrated thicknesses were obtained from viscoelastic modeling of QCM-D data, such as those shown in Figure . All experiments were repeated three times, and the error bars represent the standard deviation. Odd layers represent PSS deposition, and even layers represent PDADMA deposition.

Taken together, these growth profiles show the influence of salt concentration and salt type on the growth mechanism and growth rate of a PEM. From Figure a,b, it is evident that increasing salt concentration leads to thicker layer pairs, consistent with previous studies.[4,42−44] This is because salt screens the charges on the polyelectrolytes, weakening their interactions and changing their conformation from an extended chain to a more coil-like conformation.[23,42,45,46] As for the salt type, NaCl-assembled PSS/PDADMA PEMs grew more slowly than the KBr-assembled ones, Figure . This is because NaCl is more kosmotropic than KBr, so NaCl is less effective at extrinsic charge compensation, leading to the lower growth rate.[4,47] This has previously been reported by experiments and simulations alike, in which the more hydrated ions (Na+, Cl–) generally bind more weakly to the polyelectrolyte charge sites.[45] Similar results have been identified for complex doping with NaCl and KBr.[19,48]

Composition of PSS/PDADMA PEMs Assembled at Various NaCl and KBr Concentrations

Because we desire to understand the swelling response of the PEM on a molecular level, it is important to identify the PEM’s composition. Thus, two measurement techniques were applied to calculate the PDADMA and PSS compositions: 1H NMR spectroscopy and X-ray photon spectroscopy (XPS).

First, 1H NMR spectroscopy was performed on freestanding (PSS/PDADMA)140 PEMs at each studied assembly salt concentration following the same procedure as past studies from our group,[22,25]Table and Figure S2. For PEMs assembled at all concentrations of NaCl, the molar percentage of PSS in the PEM ranged from 43.8 to 45.5%. Similarly, for PEMs assembled from KBr, the molar percentage of PSS was 35.7 to 45.8%. For the assembly salt concentrations studied here (0–1 M), we generally observed no significant influence on the film’s composition. Our past investigation of PEMs assembled in only 0.5 M NaCl or KBr yielded a similar percentage of PSS.[49]

Table 1

PDADMA and PSS Composition from 1H NMR Spectroscopy of Freestanding (PSS/PDADMA)140 PEMs Prepared at Varying Concentrations of NaCl and KBr

assembly salt	assembly salt concentration (M)	PSS/PDADMA	PSS mol %	PDADMA mol %
NaCl	0.25	0.8	45.5	54.5
	0.5	0.8	44.1	55.9
	0.75	0.8	44.1	55.9
	1.0	0.8	43.8	56.3
KBr	0.25	0.6	36.8	63.3
	0.5	0.6	39.2	60.9
	0.75	0.6	35.7	64.3
	1.0	0.8	45.8	54.2

The NMR results above were taken from PEMs comprising 140 layer pairs, which is quite different from the PEMs comprising 4–6 layer pairs used in QCM-D measurements. To examine PEM composition for the samples used in QCM-D, XPS was performed. The areas of the high-resolution peaks of sulfur (S) in PSS and nitrogen (N) in PDADMA were computed relative to the high-resolution peak of carbon (C), shown in Figure S3. XPS results presented in Table S2 give a PSS percentage of 39.9 to 58.4% for NaCl-assembled PEMs and 44.8 to 62.6% for KBr-assembled PEMs. Table S2 shows two data outliers, which we attribute to experimental error or the limited penetration depth of the XPS method.

Overall, despite the large difference in the number of layer pairs and the method of preparation, both NMR spectroscopy and XPS techniques resulted in relatively similar compositions for the PEMs. Early studies on the growth regimes of the PEMs suggest that PEMs at later growth stages such as these at 140-layer pairs (for NMR spectroscopy) exist in the exponential growth regime in which polyelectrolyte chains exist in a more intertwined state. However, the PEMs made at x = 4 or 6 (for XPS and QCM-D) are formed in earlier growth stages such that the PEMs are still influenced by the individual polyelectrolyte chains and interactions with the substrate. Given these considerations, we used the results obtained from 1H NMR spectroscopy for our further analysis to remain consistent with our previous works.

It is notable that PDADMA is the majority component in the PEMs. This indicates that most (or all) of the PSS is occupied in intrinsic ion pairing and that PDADMA participates in both intrinsic and extrinsic ion pairs. The presence of excess PDADMA is further confirmed through examination of the raw QCM-D data shown in Figure a,b, where each PDADMA deposition step led to a greater frequency decrease than the following PSS deposition step. While complexes and multilayers with a 1:1 stoichiometric have been obtained using an extra annealing step,[50,51] many reports show a non-stoichiometric composition of PSS/PDADMA PEMs similar to observations herein.[4,17,52]

Post-Assembly Swelling Response

As earlier described, PSS/PDADMA PEMs were exposed to increasing concentrations of salt matching the assembly salt type, either NaCl (0–2.0 M) or KBr (0–1.0 M). These concentration ranges were chosen taking into consideration the critical salt concentration for PSS/PDADMA complexes in NaCl (∼2.0 M) and KBr (∼1.6 M) as well as the QCM-D thickness limit.[26,53,54]Figure presents the hydrated thickness at each post-assembly salt concentration obtained from an extended viscoelastic model applied to the QCM-D data. In general, as the assembly salt concentration increases, the hydrated thickness increases. Specifically, PEMs assembled at 1.0 M NaCl and KBr (green triangles) both exhibited hydrated thicknesses greater than those assembled at 0.25 M (black squares), regardless of post-assembly salt exposure. This is consistent with the growth profiles shown in Figure . By observation, three major swelling regions are observed for the range of post-assembly salt conditions, indicated by the different shaded regions. These swelling regions have been observed elsewhere in experimental and simulation studies of doping, swelling, and ion-pairing properties.[25,27,31,45,55]

Figure 3

Hydrated thickness of (a) (PSS/PDADMA)6 PEMs assembled in NaCl and (b) (PSS/PDADMA)4 PEMs assembled in KBr at each post-assembly salt concentration. Blue shading represents the zero to low salt concentration region; orange represents the intermediate salt concentration region; and green represents the mid to high salt concentration region for both NaCl and KBr post-assembly exposures.

In the first region (blue), at low post-assembly salt concentrations (0–0.1 M NaCl and 0–0.01 M KBr), the PEM experiences a large influx of water and counterions – so much so that several samples (assembled at>0.75 M NaCl and >0.5 M KBr) swelled beyond the range of our QCM-D instrument’s measurement range. Recalling that PDADMA exists in excess within the PEM, the film swells to allow for charge compensation of uncompensated extrinsic PDADMA sites at low salt concentrations of the contacting solution.[22,25] However, charge compensation is restricted due to the larger Debye screening length, leading to a small fraction of PDADMA units that are not ion-paired.[25] This results in electrostatic self-repulsion and significant swelling.

In the second region (orange), at intermediate salt concentrations (0.1–0.25 M NaCl and 0.01–0.1 M KBr), a slight deswelling of the PEM is observed. This can be attributed to the sufficient charge screening provided by the salt counterions, causing the PEMs to become more closely packed with less free volume for water. As more salt is added to the PEM from the contacting solution, the previously uncompensated chains become compensated, and electrostatic self-repulsion is minimized.

In the last region (green), at high salt concentrations (0.25–2.0 M NaCl and 0.1–1.0 M KBr), the PEMs begin to swell again. In this case, the swelling results from excess counterions, leading to the formation of additional extrinsic sites. This region exhibits a linear trend of the hydrated thickness with the post-assembly salt concentration. This process by which intrinsic ion pairs are broken to form extrinsic ion pairs due to the introduction of salt ions into the PEMs is also known as doping. The amount of doping and the overall amount of salt in a complex follows a Hofmeister series, in which more hydrophobic salt ions are better dopants.[20] The tendency of a salt to break intrinsic ion pairs is given by the doping constant Kunpair:where y is the fraction of intrinsic ion pairs broken or the doping level, [PE]PEC is the concentration of polyelectrolyte in complex, and [salt] is the NaCl or KBr salt concentration.

In Figure , we present the post-assembly hydrated thicknesses of the PEMs normalized against the hydrated thicknesses of those freshly prepared at the assembly salt concentration indicated. This data represents a normalization of the hydrated thickness in Figure to that of their original hydrated states. In the first swelling region (blue), the normalized thickness increased with the assembly salt concentration. Specifically for post-assembly exposure to 0 M salt concentration (Figure a), the normalized thickness for PEMs assembled in 0.25 M NaCl was 1.25 ± 0.39, and for those assembled in 0.5 M NaCl, the normalized thickness was 1.77 ± 0.05. In line with our expectations, a 4-layer pair PSS/PDADMA PEM prepared in 1.0 M NaCl doubled in thickness when exposed to water at 0 M salt concentration.[55] Taken together, this shows an increase in the swelling of the PEM due to an influx of water. In the second swelling region (orange), the normalized thickness drops below 1.0, suggesting that the film shrinks beyond its initial hydrated thickness. Here, we observe that all the NaCl assembly concentrations overlap to a similar degree. Last, in the third swelling region (green), we observe a steady increase in normalized thickness with increasing post-assembly salt. The behavior for PEMs assembled from and exposed to KBr (Figure b) followed a similar trend to that for NaCl but with some slight differences. Data for the first swelling region is harder to interpret due to a high degree of swelling beyond the instrument’s limits. Also, the large error bars allow for only a qualitatively interpretation of the effects of the assembly salt on the post-assembly swelling of KBr.

Figure 4

Post-assembly swelling thickness normalized to the final assembly thickness of (a) (PSS/PDADMA)6 assembled in NaCl and (b) (PSS/PDADMA)4 assembled in KBr. Specifically, the final assembly thickness is that of the PEM freshly prepared at the assembly salt concentration indicated but before the PEM’s exposure to media of varying salt concentrations.

Similar observations have been made from our past QCM-D studies of PSS/PDADMA PEMs prepared at only at a singular assembly condition of 0.5 M NaCl.[25,36] From the data displayed in Figures and , we now discuss the effect of varying assembly conditions (salt type: NaCl and KBr, salt concentration: 0.25–1.0 M) on the post-assembly swelling behavior. Whereas Figure shows an increase in hydrated thickness with increasing assembly salt concentration, Figure shows that actual swelling is independent of the assembly salt concentration. The overlapping of the normalized thickness plots across all salt concentrations points to a general swelling behavior of PSS/PDADMA PEMs in NaCl and KBr. Therefore, we conclude that the swelling in PSS/PDADMA PEMs appear to be “general”, regardless of the ionic strength of the assembly salt.

Water Content

We next isolated the percentage of water in the PEM by mass, wH. This calculation requires knowledge of the total PEM mass and thickness in dry and hydrated states (obtained from Figure ). We obtained the dry PEM thickness from profilometry of multilayers dried after the completion of QCM-D experiments, Figure S1. The water content, wH, was calculated using the density values of PSS/PDADMA PEMs in the dry and hydrated states (ρdry = 1270 kg/m3 and ρhydrated = 1050 kg/m3, from refs (37) and (25), respectively) to yield Figure .

Figure 5

Percentage of water by mass (wH) in (a) (PSS/PDADMA)6 PEMs assembled in NaCl and (b) (PSS/PDADMA)4 assembled in KBr.

In Figure , we observe that wH follows a trend consistent with the three swelling regions described above. This clarifies that the changes in thickness observed are not solely driven by changes in polymer conformation but also the influx of water into the PEMs. Prior to post-assembly salt exposure, wH for all PEMs assembled in NaCl was 53 ± 3% and for those assembled in KBr, it was 59 ± 6%. When exposed to pure water, wH was 58.4 ± 1.3% and 72.3 ± 1.8% for PEMs assembled in 0.25 M and 0.5 M NaCl, respectively. With the increase in ionic strength of the exposure solution, wH dropped to its lowest value of 45% at 0.25 M NaCl post-assembly concentration. Above this, wH increased with increasing post-assembly salt concentration. For the PEMs assembled in 0.25 M KBr, a similar trend was observed, but for lower KBr concentrations, the PEMs swelled outside of the instrument’s acceptable parameters (data not shown).

The amount of water in a PEM is influenced by the molecular weight, charge density, and packing density of the polyelectrolyte, as well as the contacting solution’s salt type and concentration.[56] Using neutron reflectometry to measure changes in thickness, Hariri et al. observed a water content of 80% in PSS/PDADMA complexes prepared in 2.5 M NaCl after exposure to pure water.[57] Elsewhere for PSS/PDADMA PEMs, FTIR spectroscopy was used to observe the ratio of the OH stretch peak to the SO3– peak against post-assembly salt concentration. Specifically, the trend in the OH/SO3– ratio followed a trend with post-assembly salt similar to that displayed in Figure a.[27] This supports the conclusion that the PEM swelling observed herein is a consequence of water and salt.

The PEMs herein are swollen with water, but the exact microenvironment of each water molecule remains undescribed. The molecular distribution of total water content in a PEC or PEM has been captured by various methods based on porosity, water diffusion, neutron reflectometry scattering length density, spectroscopic properties, and thermal characteristics.[22,31,33,58] Large pores have been identified in PSS/PDADMA PEMs using pulse field gradient (PFG) NMR spectroscopy.[58] Many studies have categorized the total water content in a PEM as either pore water or free water.[58,59] In other works, similarly described water states were named void and swelling water using neutron reflectometry.[14,33,60] Through the deconvolution of the OD stretch peak obtained using FTIR spectroscopy of fully immersed PSS/PDADMA PEMs, three water states based on binding energy have been identified.[31] Also, using differential scanning calorimetry (DSC), three water states were identified based on the freezing temperature of PSS/PDADMA and PAH/PAA complexes at varying hydrations.[22,32,61] While these works are aimed at separating bulk water into various water microenvironments within the PEMs, these works do not directly identify whether water molecules are associated with intrinsic or extrinsic sites.

Water Associated with Intrinsic Ion Pairs, i, and Hydration Coefficient, πsaltH

Figure shows a plot of the ratio of the total amount of water present in the PEMs at each post-assembly salt concentration to the amount of PSS repeat units, nH/nPSS or rH. PSS is used here as a representation of the amount of intrinsic ion sites present in the PEM because PSS is the limiting component. By plotting this ratio against the post-assembly salt concentration, we obtain a curve with a high (nH/nPSS) ratio (∼60–25) at low salt concentrations and a linear portion at mid to high salt concentrations. The slope of the linear portion is of particular interest because it can be used to estimate the πsaltH, which is defined by eq and describes the doping power of a salt for a specific PEM, as well as how much water a particular salt ion would bring into a PEM upon exposure.

Figure 6

Ratio of total moles of water to moles of PSS repeat units plotted against post-assembly salt concentrations for (a) (PSS/PDADMA)6 assembled in NaCl and (b) (PSS/PDADMA)4 assembled in KBr. The dashed lines show linear fits to the higher-salt regions; the slope is taken as πsaltH, and the y intercept is taken as i.

From the hydration coefficients calculated here, we observed that regardless of the assembly salt conditions of the PEMs, there is little influence on the hydration coefficient for KBr and NaCl. PEMs assembled in 0.5–1.0 M NaCl exhibited hydration coefficients of 8.5 ± 0.7 L/mol, and PEMs assembled in 0.25–1.0 M KBr exhibited hydration coefficients of 20.3 ± 1.4 L/mol. Notable exceptions include PEMs assembled at 0.25 M NaCl (13.9 ± 0.9 L/mol) and 1.0 M KBr (10.3 ± 0.9 L/mol). For comparison, Schlenoff et al.’s FTIR spectroscopy study yielded a hydration coefficient of 3.8 for a PSS/PDADMA PEM assembled in 1.0 M NaCl with a 1:1 stoichiometry.[27] The difference between the two methods may be attributed to FTIR spectroscopy measuring more strongly bound water and QCM-D measuring both strongly and loosely bound water, or else the different PEM compositions.

Using eq also, the y intercept of the linear portion of the curve in Figure is equal to i, the number of water molecules affiliated with an intrinsic ion pair. i ranges from 15 to 24 water molecules per intrinsic ion pair for NaCl-assembled PEMs, and i ranges from 20 to 29 water molecules per intrinsic ion pair, for KBr-assembled PEMs (Figure S4). Other studies indicate lower values of i. From FTIR spectroscopy, i was 6.9 ± 1.7 for PSS/PDADMA PEMs assembled in 1.0 M NaCl for all post-assembly salts studied.[27] Using neutron reflectivity, the amount of tightly bound, “immobile” water was 3.1 water molecules per PAH/PSS intrinsic ion pair for protonated layers and 1.5 water molecules per monomer pair for deuterated layers.[62] Using NMR spectroscopy, the number of water molecules per ion pair for PSS/PDADMA PEMs exposed to ambient conditions was estimated to be between 5.5 and 7.5 for 3–5 layer pairs.[63] The authors also show that for water-saturated environments, 20–25 water molecules per ion pair are expected.[63] Taken together, i might vary depending on the observation method as well as whether the PEM is immersed or not.

For another viewpoint, we conducted molecular simulations to examine the number of water molecules near PSS and PDADMA charge groups. We observed in MD simulations (non-polarizable, fixed point charge model) that the counterion position of Na+ differs from that of K+, with Na+ preferring a location close to a single PSS charge group and K+ readily shared between two PSS charge groups, at a bridging configuration. The AIMD simulations with polarizable, explicit quantum chemical solvation and ionic group descriptions show a similar tendency for a closer location of Na+ than that of K+ for the studied configurations, see Figure . In the radial distribution function g(r) plot of Figure S5, this manifests as the Na+ peak maximum being at shorter distance than that of the K+ counterion. The difference in the counterion location translates to also a binding affinity difference and affects additionally the average distance between the PSS charge groups. These in turn influence PSS conformations and hydration, e.g., coordination number. Additionally, the higher g(r) peak for MD results (Figure S5) suggests stronger binding of Na+ than K+ ions to PSS, resulting in a larger number of Na+ ions condensed around PSS chains, which is in line with the previous studies.[64,65] Considering that at moderate salt concentration, the entropy gain related to the release of counterions is a driving force for complexation, the larger fraction of condensed Na+ ions can lead to formation of a larger number of intrinsic ion pairs with PDADMA, which was observed experimentally.[49]

Figure 7

Representative snapshots of the counterion positions around PSS and PDADMA charge groups obtained from the ab initio molecular dynamics (AIMD) simulations. A single Na+, K+, Cl–, or Br– counterion, highlighted in blue, pink, orange, and red, respectively, is present in each system. The coordination numbers (CNs) of the charge group closest to the counterion were calculated between the central atom of the PE charge group (S in PSS and N in PDADMA) and oxygen of water. Similarly, the CNs of the counterions were calculated between the ion and oxygen of water.

The AIMD simulations give also more insight into the solvation of the PE charge groups and counterions. The simulations indicate a notably larger number of water molecules in the first solvation shells of the charge groups of PDADMA (ca. 20–23 water molecules) compared to the charge groups of PSS (ca. 6–7 water molecules). Although the small system size and short duration of the AIMD simulations cannot reveal this, we would expect counterion specific differences in solvation to influence also the binding distances and PE dynamics. Indeed, when comparing the solvation of the counterions in the PE–ion systems, slightly more water resides around K+ (ca. 6 water molecules) than around Na+, Cl–, and Br– (ca. 4–5 water molecules). The K+ ion also brings a slightly increased amount of water in comparison to Na+ when binding with the charge groups of PSS. This reflects indirectly in the g(r) peak positions corresponding to PE charge group – counter ion positions of both the MD and AIMD simulations and the peak heights in the MD results (Figure S5).

Based on the water coordination numbers of PSS and PDADMA charge groups, see Figure , the number of water molecules around a PSS/PDADMA intrinsic ion pair in a PE assembly can be expected to be slightly smaller than the sum of the water coordination numbers corresponding to PDADMA and PSS separately. This is because of the overlap of their hydration shells when forming an intrinsic ion pair. Considering both the simulations and experimental results here, the number of water molecules around an intrinsic ion pair is expected to be significantly larger than that found, e.g., by Schlenoff et al. via FTIR spectroscopy.[27,63] This is because the FTIR and neutron reflectivity method measures only strongly bound water, whereas the QCM-D method likely captures those water molecules both loosely and strongly associated with intrinsic ion pairs.

We next discuss sources of error in our measurements and analysis. First, eq was originally applied to a 1:1 stoichiometric PSS/PDADMA PEM in which the y intercept was taken for a purely undoped PEM with no existing extrinsic ion pairs.[27] However, in our case, the PEMs produced from QCM-D deposition are non-stoichiometric (Table and Table S2). This suggests that our values could be influenced by the presence of extrinsic ion pairs at the condition of no added salt. Second, the assumption that all PSS is in intrinsic ion pairing used in the determination of the amount of intrinsic ion pairs in Figure could also influence the validity of the results as some of the PSS might pair with Cl–. However, XPS did not reveal the presence of any Cl–. Third, any errors in the experimental measurement of thickness would also influence the determination of the number of intrinsic ion pairs used in calculation. Last, the density used in QCM-D modeling and calculations would influence the modeled hydrated thickness of the PEM.

Effect of Terminating Layer on πsaltH and i

PEMs for the preceding results were terminated with PDADMA (x = 4 or 6). In Figure , we explore the influence of the terminating layer by assembling PEMs in 0.5 M NaCl and KBr with one extra PSS layer (x = 4.5 or 6.5). Overall, a notable decrease in the water content of the PSS-terminated PEMs is produced. Figure a shows a decrease in πsaltHfrom 9.25 ± 0.55 to 4.11 ± 0.48 L/mol with an added PSS layer for NaCl. Meanwhile, Figure b shows hydration coefficients within error from 19.6 ± 1.9 to 17.3 ± 1.1 L/mol with an added PSS layer for KBr. Also, for PSS-terminated PEMs was systematically lower than PDADMA-terminated PEMs. Thus, i decreases in value with an added PSS layer, from 19.7 ± 0.6 to 0.63 ± 0.19 and from 21.7 ± 1.2 to 11.5 ± 0.7 moles of water per intrinsic ion pairs with an added PSS layer for NaCl and KBr, respectively.

Figure 8

Ratio of total moles of water to intrinsic ion pairs plotted against post-assembly concentrations for (a) (PSS/PDADMA)6 and (PSS/PDADMA)6.5 assembled in NaCl and (b) (PSS/PDADMA)4 and (PSS/PDADMA)4.5 assembled in KBr.

Many studies have shown the effect of the terminating layer of PEMs on PEM physical and chemical properties. Bruening and Miller demonstrated that PDADMA-terminated or PSS-terminated PEMs in 0.5 M NaCl swelled by 380 ± 60% or 106 ± 9% when immersed in pure water, respectively.[66] This can first be attributed to having hydrophobic PSS existing as the outermost layer of the PEM. The hydrophobicity leads to repulsion of water, evident from the high static water contact angles of ∼70° for PSS-terminated PEMs and ∼20° for PDADMA-terminated PEMs.[67] Also, PDADMA-terminated PEMs exhibit a higher surface anion concentration than PSS-terminated PEMs, which have a net charge of zero.[68,69]

Conclusions

The response of the water content of PEMs to changes in salt concentration has been extensively studied because of the great impact these changes have on the properties of these multilayer films. In this work, we show the water content responses of PDADMA/PSS PEMs to changes in salt type, salt concentration, and terminating layer. Having understood that water exists in PEMs in various microenvironments, we also further determined both by experiments and simulations the amount of water molecules associated with intrinsic ion pairs, i, in the PEMs and the hydration coefficient, πsaltH, of the salts. Using QCM-D, we identified a fast growth rate and higher post-assembly water content of PEMs assembled in KBr compared to those assembled in NaCl. This is because KBr hydrates PDADMA/PSS PEMs more effectively than NaCl, evident in the significantly higher KBr hydration coefficient. Similarly, NaCl-assembled PEMs exhibited a lower i value when compared to KBr-assembled PEMs at the same ionic strength. MD and AIMD simulations provided support of the closer binding of Na+ to PSS and pointed toward a difference in counterion position with respect to the PSS charge groups. Both observations affect the solvation of the resulting extrinsic ion pair, and the difference may lead to the formation of less extrinsic sites for the PSS–Na system than for the PSS–K system. This allows for higher influx of water into the PEMs when assembled in KBr. When terminated with the more hydrophobic polyelectrolyte, PSS, both πsaltH and i were lower than when terminated at PDADMA, proving that the hydrating ability of the salt can also be limited by the exposed polyelectrolyte. We believe that the knowledge provided by this study helps further the understanding of the water–ion pair interactions in complexes.