Subscriber access provided by GAZI UNIV
Article
Conformational Properties of Comb-Like Polyelectrolytes: A Coarse-Grained MD Study Mahdi Ghelichi, and Michael H. Eikerling J. Phys. Chem. B, Just Accepted Manuscript • DOI: 10.1021/acs.jpcb.6b00568 • Publication Date (Web): 24 Feb 2016 Downloaded from http://pubs.acs.org on March 1, 2016
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
The Journal of Physical Chemistry B is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
Conformational Properties of Comb-like Polyelectrolytes: a Coarse-grained MD Study Mahdi Ghelichi†, Michael H. Eikerling†* †
Department of Chemistry, Simon Fraser University, 8888 University Drive, Burnaby, BC V5A
1S6, Canada KEYWORDS: Ionomer molecules, conformational properties, coarse-grained molecular dynamics, backbone radius of gyration, persistence length
ACS Paragon Plus Environment
1
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 2 of 33
Abstract
This article presents a coarse-grained molecular dynamics study of single comb-like polyelectrolyte or ionomer chains in aqueous solution. The model polymer is comprised of a hydrophobic backbone chain with grafted side chains that terminate in anionic headgroups. The combpolymer is modeled at a coarse-grained level with implicit treatment of the solvent. The computational study rationalizes conformational properties of the backbone chain and localization of counterions as functions of side chain length, grafting density of side chains, backbone stiffness, and counterion valence. The main interplay that determines the ionomer properties unfolds between electrostatic interactions among charged groups, hydrophobic backbone interactions, and steric effects induced by the pendant side chains. Depending on branching density, we have found two opposing effects of side chain length on the backbone gyration radius and local persistent length. Variation in comb-polyelectrolyte architecture is shown to have nontrivial effects on the localization of mobile counterions. Changes in Bjerrum length and counterion valence are also shown to alter the strength of Coulomb interactions and emphasize the role of excludedvolume effects on controlling the backbone conformational behavior. The simulation results are in qualitative agreement with existing experimental and theoretical studies. The comprehensive conformational picture provides a fundamental framework for future studies of combpolyelectrolyte systems.
ACS Paragon Plus Environment
2
Page 3 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
Introduction Polyelectrolytes are a class of polymers with a high fraction of periodically attached ionizable groups.1,2 Polyelectrolyte solutions exhibit intriguing physical properties and they find a wide range of applications as rheology modifiers and colloidal stabilizers, adsorbent agents, biocompatible coatings, and biosensors.1–5 The structure of polyelectrolyte chains in aqueous solution is dictated by a subtle interplay of electrostatic interactions between ionic species and hydrophobic interactions of the polymer backbone with the solvent.4,6–10 The ionizable groups can be located within the backbone chain or grafted to it. The first category, linear polyelectrolytes, forms the largest class of polyelectrolytes with paramount importance in biology and various commercial uses.11,12 Examples of linear polyelectrolytes include: poly(acrylic acid) and poly(methacrylic acid), poly(styrenesulfonate), DNA, RNA, and many other polyacids and polybases.1–3,11–14 Structural properties of linear polyelectrolytes have been explored extensively in theoretical studies4,6,8,9 and computer simulations.3,4,7,10,15–20 Graft polyelectrolytes, on the other hand, have received generally less attention in studies of basic ionomer properties. In graft polyelectrolytes, ionizable groups are either located along pendant side chains or at the termini of the side chains. Many synthetic and biological polymers belong to the latter group. Polypeptides containing amino acid monomers are prominent examples.21–24 The amino acid repeating unit comprises pendant side chains terminated with either positively (Lysine and Arginine) or negatively (Aspartic and Glutamic) charged units.21,22 Conformational properties of the polypeptide backbone control numerous biological functions of the protein. Polypeptides with longer side chains exhibit improved conformational backbone flexibility due to weakened electrostatic repulsion of ionic groups.21–24
ACS Paragon Plus Environment
3
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 4 of 33
The other notable class of polyelectrolytes with comb-like architecture corresponds to perfluorinated ionomers.25–28 Ionomers contain a smaller fraction of ionizable groups compared to polyelectrolytes. Ionomer applications span from paint suspensions27,29 to electrochemical energy applications25, 26,30,31,32–35. Ionomers are synthetized both in acid and basic forms.36–40 Acid-type ionomers are the key material in polymer electrolyte fuel cells (PEFCs).25,28,30,41 The basic forms are receiving escalating recent attention as anion exchange membranes (AEM) in alkaline membrane fuel cells.40,42 Excellent proton conductivity as well as mechanical and chemical stability render PFSA ionomers the material of choice in PEFC.43 PFSA-type ionomers consist of a PTFE backbone, to which the perfluorosulfonic acid side chains are grafted. Hydrocarbon-based ionomers, both in acidic and basic forms, have also been synthesized with controllable comb-like architectures and backbone rigidities, achieved mostly through incorporation of aromatic rings into the backbone chain.28,31,39,42,44,45–48 Variation in length and density of grafted side chains has been employed in both PFSA and hydrocarbon ionomers to tune membrane properties such as water sorption, mechanical and chemical durability, and ionic conductivity.49–51,28,39,44,47 Valency of counterion exerts a pronounced effect on mechanical and water sorption properties of the ionomer membrane.52,53 A recently developed poroelectroelastic model of water sorption in ionomer membranes relates these changes to the alteration of the elastic properties of the polymer walls formed by self-aggregated ionomer backbones.54 A coherent understanding of structure, properties and performance of ionomer membranes is complicated by the multiple scales involved. The three main scales are (i) molecular scale comprising ionomer chains, ionomer bundles, anionic charges and surrounding cloud of protons, (ii) pore level at which elastic properties of pore walls as well as proton distri-
ACS Paragon Plus Environment
4
Page 5 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
bution and mobility in the pore can be defined, and (iii) hydrated membrane scale at which sorption characteristics, proton conductivity and elasto-mechanical properties can be defined.55,56. Understanding basic structure and properties at the single chain level is of paramount importance. From a practical standpoint, it determines the ionomer solution behavior and elucidates how the structure and properties will evolve at all three levels. However, single chain studies are scarce; mostly because of the practically motivated focus on membrane-level investigations; but also due to spontaneous aggregation of hydrophobic chains that complicates single chain investigations. Recently, Minko and coworkers reported atomic force microscopy (AFM) imaging of singlechain adsorbed Nafion chains on mica and highly ordered pyrolitic graphite (HOPG) substrates.57 A compact globular structure was reported for adsorbed ionomer chains.57 Theoretically, conformational features of single chains of Nafion-type ionomers have been studied mainly through atomistic simulations. In the early 2000’s, Vishnyakov et al.58,59 studied the effect of aqueous and alcoholic solvents on the conformation of Nafion oligomers. Recently, density functional theory (DFT) combined with energy-loss spectroscopy (EELS) was employed to study the effect of different side chain chemistries on the conformational features of single PFSA backbone chains.60 Computational limitations associated with the incorporation of atomistic details along with illdefined simulations of inadequately short oligomeric chains have hindered drawing definite conclusions regarding the interplay of side chain length and density on the conformational features of comb-like ionomer molecules. Coarse-grained MD simulations, in which chemical details are foregone for the sake of increasing the time and length scale of simulations, have been able to capture structural changes in charged polymers. Changes in the radius of gyration of the backbone and in the persistence
ACS Paragon Plus Environment
5
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 6 of 33
length of single linear polyelectrolyte chains have been explored in relation to variations in chain rigidity61,62, charge fraction10,65,47,49,50, and dielectric properties of the solvent.10,17,19,20,49 For branched polyelectrolytes with charged side chains, conformational properties of the hydrophobic backbone have been studied as a function of the length and density of the side chains.63,66 In neutral polymers with pendant side chains (bottle-brush or comb-polymers) it is known that steric hindrances due to side chain congestion increase the gyration radius and stiffness of the backbone.67–69 For comb-polyelectrolytes with charged moieties at side chain terminals more complex conformational responses are expected due to the presence of both electrostatic and excludedvolume interactions. To the best of our knowledge, for hydrophobic comb-like polyelectrolytes with charge-terminated side chains, there have not been any attempts in developing a coherent picture of their conformational properties. This article strives to rationalize the backbone conformation and the local rigidity of single chains of comb-polyelectrolytes as a function of changes in architectural and environmental parameters. The Section “Chain model” introduces the coarse-grained chain model, the type of interactions considered, and the simulation protocol. Section “Results and discussion” provides a systematic evaluation of backbone conformational changes as a function of variation in pendant side chain length and density, strength of electrostatic interactions, and counterion valence. Chain model In a coarse-grained model, the interacting species are groups of atoms or beads that interact via effective forces. We model a bead-spring polyelectrolyte chain consisting of a hydrophobic backbone and pendant side chains in implicit solvent with explicit counterions. Figure 1 is a schematic illustration of the branched polyelectrolyte chain considered in the model.
ACS Paragon Plus Environment
6
Page 7 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
The model is devised to emulate PFSA-type ionomers. Each backbone bead resembles one (CF2-CF2) repeating unit with diameter of approximately 2.4 Å. The simulated backbone chain consists of 104 beads with regularly spaced pendant side chains, representing a real backbone contour length of ~ 25 nm. Anionic headgroups carry a -1e charge and they are represented by a single bead (red) with the same size as backbone and side chain beads. Ls denotes the number of beads in each side chain, and Nb is the number of backbone beads per ionomer repeat unit (containing 1 side chain), as illustrated in Figure 1. We explore values Ls = 1, 3, 5, and 7 and Nb = 1, 2, 3, 5. Since all side chains contain a charged unit, we define a charge fraction f = 1/Nb, for which we will consider values f = 1, 1/2, 1/3, 1/5, with the lower range resembling ionomers and the upper range corresponding to bottle-brush (f = 1) and polyelectrolyte systems (f = 1/2, 1/3). Electroneutrality is guaranteed by adding Nc = Ns/Zc counterions (yellow beads in Figure 1) where Zc is the counterion valence and Ns is the total number of side chains. No other electrolyte is added to the solution system. As a reference, we simulated all comb-like architectures in their uncharged form in electrolyte-free solution, for all considered combinations of f and LS. As a simplification, all beads are assumed to have the same diameter σ. Adjacent monomer beads in the backbone and side chains are connected by the finitely extensible nonlinear elastic (FENE) potential with spring constant k = 30 ε/σ2 and maximum extension R0 =1.5σ,70–72 U FENE
⎛ ⎛ ⎞2 ⎞ r 1 rij = − kR0 ln ⎜1− ⎜⎜ ij ⎟⎟ ⎟ ⎜ 2 R ⎟, ⎝ ⎝ 0⎠ ⎠
( )
(1)
The flexibility of the backbone is modeled with the bond angle potential
( )
(
(
Uθ rij = kθ 1− cos θ − θ 0
)) ,
(2)
ACS Paragon Plus Environment
7
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 8 of 33
where kθ is the bending rigidity and θ0 the equilibrium angle that is set to 180º for linear segments along backbone. We considered values kθ = 0, 15, 30, 45 kBT for the backbone chain. Pendant side chains are freely rotating. Van der Waals interactions between monomer beads are modeled through a standard shifted and truncated Lennard-Jones (LJ), ⎧ ⎡⎛ ⎞12 ⎛ ⎞6 ⎛ ⎞12 ⎛ ⎞6 ⎤ σ σ σ ⎥ ⎪ ⎢ σ ⎪ 4ε LJ ⎢⎜⎜ ⎟⎟ − ⎜⎜ ⎟⎟ − ⎜⎜ ⎟⎟ + ⎜⎜ ⎟⎟ ⎥, r r r r U LJ (rij ) = ⎨ ⎣⎝ ij ⎠ ⎝ ij ⎠ ⎝ c ⎠ ⎝ c ⎠ ⎦ ⎪ ⎪ , ⎩0
rij < rc rij > rc
(3)
where rij is the distance between beads i and j and σ is the bead diameter. The cut-off distance was set to rc = 2.5σ, for interactions between uncharged polymer beads (backbone and side chain beads), and rc = 21/6σ for all other pairwise interactions (neutral-anionic, neutral-cationic, anionic-anionic, anionic-cationic) resulting in the purely repulsive part of the Lennard-Jones (LJ) potential. The reference value of the Lennard-Jones interaction parameter is set to εLJ = 1 kBT. The units of energy, diameter and mass of all polymer beads are set to unity, i.e. kBT = 1, σ = 1, and m = 1. The θ solvent condition for the polymer backbone corresponds to εLJ = 0.33 kBT.19,73 The reference value of εLJ for interactions between uncharged beads was set to 1.0 kBT, representing an “intermediate” level of hydrophobicity. For interactions between charged beads, we used a Coulomb potential of the form,
( )
U Coul rij = k BT
λB qi q j rij
,
(4)
where λB = e02/4πε εskBT is the Bjerrum length, e0 is the unit charge and ε0 and εs are the permit0
tivity of vacuum and the relative dielectric constant of the solvent. The baseline value of the Bjerrum length is λB = 3.0 σ, which corresponds to about 0.7 nm in aqueous electrolyte solution
ACS Paragon Plus Environment
8
Page 9 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
(σ ≈ 0.24 nm). In our parametric analysis, we explore effects of LS, f, λB, ZC, and kθ . Table 1 summarizes the baseline values and parameter ranges explored. The particle–particle particle–mesh (PPPM)74 was used for the calculation of electrostatic interactions. Simulations were carried out in the NVT ensemble with periodic boundary conditions in all three directions for a cubic simulation box with length L = 110 σ. The temperature was maintained by a Langevin thermostat,75 in which the motion of beads is governed by
mi
!" dV i t
()
dt
! !" ! ! d ri t !" = −∇U i − Γ + Fi t dt ,
()
()
(5)
where mi is the mass of bead i and Ui is the potential energy experienced by the bead as the result
! of bonded and nonbonded interactions. Fi t is a random force with zero average value. The
()
friction coefficient of Γ = τLJ-1 was calibrated to maintain a reduced temperature of T* = kBT/εLJ = 1 where τLJ is the standard LJ time, τLJ = σ (m/εLJ)1/2. The equations of motion for beads in the system were integrated using the velocity-Verlet algorithm with a time step of Δt = 0.01 τLJ. The following procedures were employed in simulations, which were performed with LAMMPS code.76 A single polyelectrolyte chain was placed in the middle of the simulation box and simulated for 7×106 time steps. The first 5.0×106 timesteps we used for equilibration, tracked through convergence in the total energy and the radius of gyration. For equilibrium state analysis, data collections were performed for the last 2×106 timesteps and data were saved every 2×104 timesteps, resulting in 100 regularly time-spaced snapshots. Before presenting simulation results, a remark on the implicit solvent treatment is in order. The effects of solvent molecules are implemented through the friction term of the Langevin thermostat. This specific treatment, which is a common practice in coarse-grained MD simulations on
ACS Paragon Plus Environment
9
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 10 of 33
the polyelectrolyte systems, implies the absence of explicit solvent molecules. Absence of these molecules lowers the computational costs and thus longer simulation runs become practically possible, allowing for a proper equilibration of the system. However, as argued in Refs. [17,77,78], the absence of many-body effects imposed by solvent molecules can significantly influence the collapse and conformational dynamics of both neutral and polyelectrolyte solutions under poor solvent conditions. Since we aim at capturing the equilibrium conformational features of long comb-polyelectrolyte chains and not the relaxation dynamics of our model system, the use of implicit solvent molecules does not alter any of the reported conformational behaviors. We have also assumed a uniform polarizability of the implicit solvent by fixing the dielectric constant value. This assumption implies that the electrostatic interaction among any pair of charged groups depends only on their inter/intra separation distances. In reality, local charging effects result in polarization of solvent molecules and alteration of local solvent dipole moments. This implies subtle variations in dielectric constant and spatially dependent electrostatic interactions. Solvent polarization and local dielectric effects can be introduced to the studied system either through incorporation of explicit solvent molecules or by developing polarizable coarsegrained solvent models.65,79 Even though the consideration of polarization effects yields more realistic results for local conformational properties and needs to be the focus of future studies, the imposed computational complexity shortens the equilibration time period. For the sake of simulating long comb-like chains and capturing the global conformational trends, it seems that consideration of uniform dielectric media is a legitimate assumption. Results and Discussion Effects of side chain length and density. As discussed above, conformational behavior of a branched polyelectrolyte chain originates from both steric effects due to branching and the elec-
ACS Paragon Plus Environment
10
Page 11 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
trostatic interactions of charged units. To be able to discriminate both of these effects, we simulated grafted polymers in both their charged and neutral states (UCoul = 0). We used the backbone gyration radius, Rg (σ), as the main conformational variable, which is defined as
1 N R = ∑ ri − rCOM N i=1 2 g
2
, rCOM
1 N = ∑ ri N i=1
,
(6)
where i is the backbone bead index, N is the number of backbone beads, ri is the position vector of the ith bead, and rCOM the center-of-mass (COM) of the chain. Angular brackets indicate averaging over chain configurations. Figure 2 shows the results for Rg obtained in the case of freely rotating chains. Zero bending rigidity implies that all conformational changes stem from excluded-volume and electrostatic effects. Figure 2a shows that increasing LS and f increases Rg. These effects become more prominent as f increases. No significant change in Rg with LS is observed for f = 1/5. The increase in backbone Rg upon increasing LS and f for the neutral polymer is due to enhanced excludedvolume interactions of side chain beads. This trend reflects the classical behavior that was previously captured through extensive Monte Carlo (MC)67–69,80 and MD63,81 simulations. Figure 2b shows the changes in backbone Rg for comb-polyelectrolyte chains in the presence of monovalent counterions at λB = 3 σ. For comb-polyelectrolytes, higher f increases the backbone Rg for all LS values explored. Higher density of side chains results in closer proximity of charged units and stronger electrostatic repulsion that causes greater stretching of the backbone. Charged polymers maintain larger Rg values compared to the uncharged equivalent polymer, clearly demonstrating that eliminating the charges leads to smaller and more compact conformations.
ACS Paragon Plus Environment
11
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 12 of 33
Regarding the effect of side chain length, two different trends are observed depending on the fraction of side chains. At f = 1, increasing LS increases Rg, while f ≤ 1/2 results in the opposite trend. The change in LS induces a transition in the interplay of excluded-volume and Coulombic interactions. Increasing LS leads to a larger distance between charged headgroups and correspondingly weakened electrostatic repulsion, which favors a reduction in Rg. On the other hand, higher LS increases the steric congestion of side chains, which favors a stretching of the backbone chain, to increase Rg. At f = 1, the steric effect of side chain congestion outweighs the electrostatic effect and the increase in LS leads to enhanced Rg. The milder increase in Rg at f = 1 for the charged chain as compared to an uncharged chain reflects the opposite influence of increasing LS on the chain size. For f ≤ 1/2 excluded-volume effects are weakened due to the greater separation of pendant side chains and the electrostatic interactions among anionic headgroups exert the dominant effect on Rg. Weakened Coulombic repulsion at higher LS lowers the electrostatically induced stretching of the backbone and Rg decreases. The decrease in backbone size upon increasing the length of side chains was also experimentally reported for polypeptide chains bearing longer pendant side chains.21,22,24 Figure 2c shows the pure electrostatic contributions to the backbone Rg once the excluded-volume effects, Rg values in Figure 2a, are subtracted from Rg values presented in Figure 2b. There is a noticeable difference in the pure electrostatic contribution to Rg between short and long regimes of LS for all f values. This contribution is more pronounced for regimes of small LS values. The trend seen in Fig 2c shows that for short side chain polymers Coulombic effects prevail. This effect is owned to the closer proximity of anionic moieties at smaller LS values that strengthen the electrostatic
ACS Paragon Plus Environment
12
Page 13 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
repulsion among charged groups. The electrostatic contribution continuously diminishes as LS increases and the excluded-volume effect gains weight in determining Rg. Figure 3 shows the snapshots of the comb-polyelectrolyte chains at different LS and f values. The conformational behavior observed in Figure 2b is in line with the conformational trends in Figure 3. For f = 1/5 we can see the dominant effect of LS on the size of the chain: while chains with LS ≥ 3 exhibit a globular conformation, chains with LS = 1 assume an expanded coil-like conformation. For f = 1, increasing LS results in the transformation of an expanded chain at LS = 1 to a rod-like chain at LS = 7. Increase in the fraction of pendant side chains results in monotonic expansion of collapsed chains. All-atomistic15 and coarse-grained65 MD simulation studies of sodium polystyrenesulfonate (NaPSS) in water bearing different fractions of anionic-terminated branches demonstrated a similar conformational trend. Conformational structures obtained for f < 1/5 are indistinguishable from f = 1/5 and are not discussed in this work. The globular structure observed for f = 1/5 and LS ≥ 3 is in line with the globular structure of a single chain Nafion ionomer reported in Ref [57]. Distally located branches weaken the steric congestion of side chains and lower the Coulombic self-charge repulsions. Minimization of the surface energy of hydrophobic polymer triggers the collapse of the ionomer chain to form a globular structure. Figure 4 shows the effect of different bending rigidities, kθ , on Rg of both charged and uncharged comb-polymer for the explored range of f and LS. The value f = 1/2 closely resembles the charge fraction of the coarse-grained Nafion models employed by Malek et al82,83. Imposing bending rigidity to the backbone chain continuously increases Rg for both the charged and uncharged polymers.
ACS Paragon Plus Environment
13
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 14 of 33
More rigid chains exhibit a smaller sensitivity of Rg to the structural changes for both charged and uncharged polymers. For very stiff chains, Rg is nearly independent of f and LS. Increased backbone rigidity introduces an energetic barrier for the conformational changes introduced either through topological effects (excluded-volume effect) or through electrostatic repulsions. In what follows, the effect of changes in length and density of side chains on the local orientational correlation or local persistence length, Lkp , of the backbone chain is explored. The local !" persistence length is the scalar projection of the end-to end vector, Rb , onto each bond in the
!"
backbone, rb , which is defined by63,67
Lkp
!" ! rb !" ! • Rb , = !" rb
(7)
where k is the backbone monomer or bead index and the angular brackets, , denote averaging over chain configurations. Figure 5 shows the results for Lkp or the persistent length (LP) as a function of monomer index, k, along the backbone chain for neutral and charged comb-polymers. For both charged and uncharged polymers, greater f results in higher Lkp . Backbones of a combpolyelectrolyte maintain a greater local rigidity than the uncharged polymer due to Coulombicinduced chain stiffening. Increasing LS results in a monotonic increase in Lkp at all f values for neutral polymers and at f = 1 for the comb-like polyelectrolytes. Increasing LS increases the backbone flexibility at f ≤ 1/2 for charged chains. No special features are observed in Lkp for f = 1/5 at LS > 1 in both uncharged and charged chains. In the regime of low and distant side chain
ACS Paragon Plus Environment
14
Page 15 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
contents, Lkp values are close to zero, which suggests that the entire backbone is in a compact globular conformation. The Lkp values of f = 1, 1/2, 1/3 for the charged polymer, shown in Figure 5 reveal a curving of Lkp upon increase in LS. This curvature arises from lower rigidity at both ends of the backbone chains. The monomers at the ends of a chain experience reduced steric repulsion from the adjacent grafts. Shorter branching, i.e., smaller LS, yields more homogenous backbone rigidity and diminished end effects as the long-range Coulomb forces gain in importance Overall, smaller LS and larger f result in greater rigidity of the backbone chain for f ≤ 1/2. These two main features, smaller LS and higher f, resemble the key characteristics of short side chain (SSC) ionomers compared to longer side chain ionomers such as Nafion. At the single-chain level, SSC ionomers are expected to maintain greater rigidity than their longer side chain counterparts. Localization of counterions We explore how the length and fraction of pendant side chains as well as the local rigidity of the backbone influences the localization of monovalent counterions. Figure 6a shows the radial distribution function (RDF) of anionic headgroups and counterions, gh-c(r) for LS = 3 at various f. For an isotropic system, an RDF gA-B(r) describes the probability of finding a particle A in a specific radial distance, r, of a reference particle B. It is defined as84
N
()
g A−B r =
N
A B V ∑∑ δ r − rij 4πr 2 N A N B i=1 j=1
(
)
,
(8)
ACS Paragon Plus Environment
15
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 16 of 33
where NA and NB are the numbers of A and B particles in the simulated system, respectively, and the triangular brackets around the δ-function indicate averaging over the equilibrium system trajectory. Greater f increases the surface charge of the chain and creates a stronger surface potential, thereby causing greater Coulombic attraction of the counterion, a trend that was found in other simulation studies as well [15,65,79]. Figure 6b demonstrates the effect of LS on the localization of counterions for f = 1/2. Counterions have a stronger affinity to anionic headgroups for ionomers with shorter side chains. The enhanced surface charge density at lower LS is responsible for the creation of stronger surface potential and greater localization of counterions. The changes in the localization of counterions in the vicinity of anionic headgroups as a function of f and LS can also be seen in the snapshots depicted in Figure 3. Figure 6c shows how the intrinsic rigidity of the backbone chain influences the localization of counterions to anionic groups for f = 1/2 and LS = 3. Affinity of counterions to the headgroups decreases as the bending rigidity increases. Greater kθ results in more stretched backbone chains with farther separated side chains. This greater distance among charged moieties results in lower effective surface charge of the chains and thus, weaker electrostatic attraction to the counterions. Counterion localization exhibits a smaller sensitivity to kθ than to f and LS. These observations highlight the nontrivial effects of architectural variation in comb-polyelectrolyte chains on the distribution and possibly the dynamics of mobile counterions. Effect of electrostatic interaction strength In this section, we explore how the change in strength of electrostatic interactions introduced through the variation in λB influences Rg. Figure 7a shows Rg plotted over λB for f = 1/2 and for different LS. Several trends are noticeable in this plot. First, a nonmonotonic dependence of Rg on
ACS Paragon Plus Environment
16
Page 17 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
the λB is evident for all LS. Secondly, the sensitivity of the nonmonotonic dependence decreases as LS increases. A nonmonotonic dependence of Rg on λB was reported in Refs.[19,20,61,62] for linear polyelectrolyte chains. It was attributed to the initial expansion of the chain due to enhanced electrostatic repulsion upon increasing λB. Further increase in λB results in strong localization of counterions in the vicinity of the chain that screens anionic charges and lowers their effective electrostatic repulsion. In the regime of large and small λB values, i.e., λB = 0.5 σ and λB ≥ 9 σ, side chains with larger LS assume a greater radius of gyration. This behavior can be rationalized through the dominance of excluded-volume effects in the case of weak or screened electrostatic interactions. In the midrange of λB = 1 to 6 σ, enhanced Coulombic repulsions along with shorter distance between the charged units of chains with small LS results in higher Rg. Thus, we can see that increasing the side chain length does not always result in enhanced chain size and this behavior is sensitive to the strength of electrostatic interactions. Increasing LS results in weaker sensitivity of Rg to λB, which can be attributed to the reduced impact of Coulomb interactions. Figure 7b shows how comb-polyelectrolytes (f = 1/2, LS = 3) with different backbone rigidity behave upon variation of λB. We can see the nonmonotonic dependence for chains with different kθ values. However, the sensitivity to the change in Bjerrum length, and the extent of chain ex-
pansion and shrinkage, diminishes for more rigid backbones. Ou et al.62 also observed similar behavior and reported that increasing the backbone rigidity diminishes the dependence of Rg on λB. Effect of counterion valency The results for the change in Rg upon change in ZC for chains of different f and LS are shown in Figure 8. Figure 8a shows the change in Rg of comb-polyelectrolytes with different fraction of
ACS Paragon Plus Environment
17
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 18 of 33
side chains at fixed LS = 3. Increasing the valency of counterions, from ZC = 1 to ZC = 2 and further to ZC = 3 leads to a decrease in Rg. Increasing the counterion valence yields stronger correlation of counterions with charged headgroups. Enhanced counterion localization triggers an effective charge screening that results in decreased electrostatic repulsion. Weakened Coulomb interactions along with attractive short-range hydrophobic effects causes the comb-polyelectrolyte chain to collapse to a neutral-like structure. The reduction in Rg is more pronounced for chains with higher f, owed to the prevailing electrostatic contribution to Rg at high f (Figure 2c). Figure 8b shows the effect of an increase in ZC for f = 1/2 and different LS. The extent of the decrease in Rg is more pronounced with ZC upon decreasing LS. At ZC = 1, shorter branching yields greater Rg due to enhanced Coulombic repulsion. This trend is opposite of that observed at ZC = 2, 3 where excluded-volume effects primarily control the backbone Rg and higher LS yields larger backbone Rg. Conclusions The presented coarse-grained molecular dynamics study scrutinizes the effects of different parameters on the conformational properties of hydrophobic comb-like polyelectrolytes. Observed trends can be rationalized through the fine interplay of hydrophobic interactions of individual backbone units, excluded-volume effects of side chains, and Coulombic interactions among charged moieties. It is observed that increasing the fraction of pendant side chains increases the radius of gyration and rigidity of comb-polyelectrolyte chains, a trend similar to the neutral counterpart. For high density of side chains, increasing the length of side chains leads to backbone stretching. For lower density of side chains, increasing side chain length yields a smaller backbone radius of gyration, resulting in the formation of collapsed globular structures. Counterions possess stronger affinity to anionic moieties for shorter side chains and higher grafting den-
ACS Paragon Plus Environment
18
Page 19 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
sities. Chain radius of gyration is shown to demonstrate a nonmonotonic dependence on the strength of electrostatic interaction. In the regime of weak Coulomb interaction steric effects of side chains control the radius of gyration and backbone conformation. Multivalent counterions are shown to decrease the radius of gyration through the screening of anionic charges. Qualitative agreement in the conformational trends with both experimental and similar theoretical studies highlights the ability of the generic coarse-grained model in capturing the subtle conformational behavior of comb-polyelectrolytes and ionomer chains. Conformational understandings at the level of single chain are particularly essential for characterization of ionomer solutions and the underlying process of self-assembly. Finally, we hope that this study will trigger further fundamental and single chain studies of fluorinated and hydrocarbon ionomers. Author information Corresponding Author *E-mail:
[email protected]. Notes The authors declare to competing financial interest. Acknowledgment The work reported in this manuscript was performed with the financial support from an NSERC Discovery Grant, entitled Materials for Electrochemical Energy Conversion: From Fundamental Physics to Advanced Design. Simulations were performed using the Grex computational facility at the University of Manitoba, provided through WestGrid and Compute Canada.
ACS Paragon Plus Environment
19
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 20 of 33
References (1)
Hara, M. Polyelectrolytes: Science and Technology; Marcel Dekker: New York, 1993.
(2)
Holm, C.; Joanny, J. F.; Kremer, K.; Netz, R. R.; Reineker, P.; Seidel, C.; Vilgis, T. A.; Winkler, R. G. Polyelectrolyte Theory. Adv. Polym. Sci. 2004, 166, 67–111.
(3)
Khokhlov, A. R.; Khalatur, P. G. Solution Properties of Charged Hydrophobic/hydrophilic Copolymers. Curr. Opin. Colloid Interface Sci. 2005, 10, 22–29.
(4)
Dobrynin, A. V. Theory and Simulations of Charged Polymers: From Solution Properties to Polymeric Nanomaterials. Curr. Opin. Colloid Interface Sci. 2008, 13, 376–388.
(5)
Burak, Y.; Ariel, G.; Andelman, D. Onset of DNA Aggregation in Presence of Monovalent and Multivalent Counterions. Biophys. J. 2003, 85, 2100–2110.
(6)
Dobrynin, A. V; Colby, R. H.; Rubinstein, M. Scaling Theory of Polyelectrolyte Solutions. Macromolecules 1995, 28, 1859–1871.
(7)
Dobrynin, A. V. Effect of Counterion Condensation on Rigidity of Semiflexible Polyelectrolytes. Macromolecules 2006, 39, 9519–9527.
(8)
Dobrynin, A. V.; Rubinstein, M. Theory of Polyelectrolytes in Solutions and at Surfaces. Prog. Polym. Sci. 2005, 30, 1049–1118.
(9)
Dobrynin, A. V.; Rubinstein, M. Hydrophobic Polyelectrolytes. Macromolecules 1999, 32, 915–922.
(10)
Limbach, H.-J.; Holm, C. Single Chain Properties of Polyelectrolytes in Poor Solvent. J. Phys. Chem. B 2003, 107, 8041–8055.
(11)
Visakh, P. M. Polyelectrolyte: Thermodynamics and Rheology; Springer International Publishing, 2014.
(12)
Radeva, T. Physical Chemistry of Polyelectrolytes; CRC Press, 2001.
(13)
Tanford, C. Physical Chemistry of Macromolecules; Wiley: New York, 1961.
(14)
Katchalsky, A. Polyelectrolytes and Their Biological Interactions. Biophys. J. 1964, 4, 9–41.
(15)
Carrillo, J. M. Y.; Dobrynin, A. V. Detailed Molecular Dynamics Simulations of a Model NaPPS in Water. J. Phys. Chem. B 2010, 114, 9391–9399.
(16)
Carrillo, J.-M. Y.; Dobrynin, a V. Polyelectrolytes in Salt Solutions: Molecular Dynamics Simulations. Macromolecules 2011, 44, 5798–5816.
(17)
Chang, R.; Yethiraj, A. Dilute Solutions of Strongly Charged Flexible Polyelectrolytes in Poor Solvents: Molecular Dynamics Simulations with Explicit Solvent. Macromolecules 2006, 39, 821–828.
(18)
Liao, Q.; Dobrynin, A. V.; Rubinstein, M. Counterion-Correlation-Induced Attraction and Necklace Formation in Polyelectrolyte Solutions: Theory and Simulations. Macromolecules 2006, 39, 1920–1938.
(19)
Micka, U.; Holm, C.; Kremer, K. Strongly Charged, Flexible Polyelectrolytes in Poor Solvents: Molecular Dynamics Simulations. Langmuir 1999, 15, 4033–4044.
(20)
Stevens, M. J.; Kremer, K. The Nature of Flexible Linear Polyelectrolytes in Salt Free Solution: A Molecular Dynamics Study. J. Chem. Phys. 1995, 103, 1669.
(21)
Deming, T. J. Synthesis of Side-Chain Modified Polypeptides. Chem. Rev. 2016, 116, 786–808.
(22)
Zhang, Y.; Lu, H.; Lin, Y.; Cheng, J. Water-Soluble Polypeptides with Elongated, Charged Side Chains Adopt Ultrastable Helical Conformations. Macromolecules 2011, 44, 6641–6644.
(23)
Wang, Y.; Zheng, M.; Meng, F.; Zhang, J.; Peng, R.; Zhong, Z. Branched Polyethylenimine Derivatives with Reductively Cleavable Periphery for Safe and Efficient in Vitro Gene Transfer. Biomacromolecules 2011, 12, 1032–1040.
(24)
Lu, H.; Wang, J.; Bai, Y.; Lang, J. W.; Liu, S.; Lin, Y.; Cheng, J. Ionic Polypeptides with Unusual Helical Stability. Nat. Commun. 2011, 2, 206.
ACS Paragon Plus Environment
20
Page 21 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
(25)
Eikerling, M.; Kulikovsky, A. Polymer Electrolyte Fuel Cells: Physical Principles of Materials and Operation; CRC Press: Boca Raton, FL, 2014.
(26)
Eisenberg, A.; Kim, J.-S. Introduction to Ionomers; John Wiley & Sons: New York, 1998.
(27)
Nomula, S.; Cooper, S. L. Ionomer Solution Structure and Solution Diagram. Macromolecules 2001, 34, 2653–2659.
(28)
Kreuer, K. D. Ion Conducting Membranes for Fuel Cells and Other Electrochemical Devices. Chem. Mater. 2014, 26, 361–380.
(29)
Nomula, S.; Cooper, S. L. Influence of Ionic Content in Polyurethane Ionomer Solutions. J. Phys. Chem. B 2000, 104, 6963–6972.
(30)
Pu, H. Polymers for PEM Fuel Cells; John Wiley & Sons: Hoboken, NJ, 2014.
(31)
Zhang, L.; Brostowitz, N. R.; Cavicchi, K. A.; Weiss, R. A. Perspective: Ionomer Research and Applications. Macromol. React. Eng. 2014, 8, 81–99.
(32)
Hall, L. M.; Stevens, M. J.; Frischknecht, A. L. Dynamics of Model Ionomer Melts of Various Architectures. Macromolecules 2012, 45, 8097–8108.
(33)
Hall, L. M.; Seitz, M. E.; Winey, K. I.; Opper, K. L.; Wagener, K. B.; Stevens, M. J.; Frischknecht, A. L. Ionic Aggregate Structure in Ionomer Melts: Effect of Molecular Architecture on Aggregates and the Ionomer Peak. J. Am. Chem. Soc. 2012, 134, 574–587.
(34)
Hall, L. M.; Stevens, M. J.; Frischknecht, A. L. Effect of Polymer Architecture and Ionic Aggregation on the Scattering Peak in Model Ionomers. Phys. Rev. Lett. 2011, 106, 1–4.
(35)
Ting, C. L.; Stevens, M. J.; Frischknecht, A. L. Structure and Dynamics of Coarse-Grained Ionomer Melts in an External Electric Field. Macromolecules 2015, 48, 809–818.
(36)
Skalski, T. J. G.; Britton, B.; Peckham, T. J.; Holdcroft, S. Structurally-Defined, Sulfo-Phenylated, Oligophenylenes and Polyphenylenes. J. Am. Chem. Soc. 2015, 137, 12223–12226.
(37)
Wright, A. G.; Holdcroft, S. Hydroxide-Stable Ionenes. ACS Macro Lett. 2014, 3, 444–447.
(38)
Thomas, O. D.; Soo, K. J. W. Y.; Peckham, T. J.; Kulkarni, M. P.; Holdcroft, S. A Stable HydroxideConducting Polymer. J. Am. Chem. Soc. 2012, 134, 10753–10756.
(39)
Tsang, E. M. W.; Zhang, Z.; Shi, Z.; Soboleva, T.; Holdcroft, S. Considerations of Macromolecular Structure in the Design of Proton Conducting Polymer Membranes : Graft versus Diblock Polyelectrolytes. J. Am. Chem. Soc. 2007, 129, 15106–15107.
(40)
Merle, G.; Wessling, M.; Nijmeijer, K. Anion Exchange Membranes for Alkaline Fuel Cells: A Review. J. Memb. Sci. 2011, 377, 1–35.
(41)
Schlick, S. Ionomers: Characterization, Theory, and Applications; CRC Press: Boca Raton, FL, 1996.
(42)
Arges, C. G.; Ramani, V.; Pintauro, P. N. Anion Exchange Membrane Fuel Cells. Elctrochemical Soc. Interface 2010, 31–35.
(43)
Mauritz, K. a; Moore, R. B. State of Understanding of Nafion. Chem. Rev. 2004, 104, 4535–4585.
(44)
Elabd, Y. A.; Hickner, M. A. Block Copolymers for Fuel Cells. Macromolecules 2011, 44, 1–11.
(45)
Li, N.; Guiver, M. D. Ion Transport by Nanochannels in Ion-Containing Aromatic Copolymers. Macromolecules 2014, 47, 2175–2198.
(46)
Titvinidze, G.; Kreuer, K.-D.; Schuster, M.; de Araujo, C. C.; Melchior, J. P.; Meyer, W. H. Proton Conducting Phase-Separated Multiblock Copolymers with Sulfonated Poly(phenylene Sulfone) Blocks for Electrochemical Applications: Preparation, Morphology, Hydration Behavior, and Transport. Adv. Funct. Mater. 2012, 22, 4456–4470.
(47)
Chang, Y.; Mohanty, A. D.; Smedley, S. B.; Abu-Hakmeh, K.; Lee, Y. H.; Morgan, J. E.; Hickner, M. A.; Jang, S. S.; Ryu, C. Y.; Bae, C. Effect of Superacidic Side Chain Structures on High Conductivity Aromatic Polymer Fuel Cell Membranes. Macromolecules 2015, 48, 7117–7126.
(48)
Choi, U. H.; Ye, Y.; Cruz, D. S. de la; Liu, W.; Winey, K. I.; Elabd, Y. A.; Runt, J.; Colby, R. H. Dielectric
ACS Paragon Plus Environment
21
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 22 of 33
and Viscoelastic Responses of Imidazolium-Based Ionomers with Different Counterions and Side Chain Lengths. Macromolecules 2014, 47, 777−790. (49)
Luo, X.; Holdcroft, S.; Mani, A.; Zhang, Y.; Shi, Z. Water, Proton, and Oxygen Transport in High IEC, Short Side Chain PFSA Ionomer Membranes: Consequences of a Frustrated Network. PCCP 2011, 13, 18055–18062.
(50)
Merlo, L.; Ghielmi, A.; Cirillo, L.; Gebert, M.; Arcella, V. Resistance to Peroxide Degradation of Hyflon® Ion Membranes. J. Power Sources 2007, 171, 140–147.
(51)
Ghassemzadeh, L.; Marrony, M.; Barrera, R.; Kreuer, K. D.; Maier, J.; Müller, K. Chemical Degradation of Proton Conducting Perflurosulfonic Acid Ionomer Membranes Studied by Solid-State Nuclear Magnetic Resonance Spectroscopy. J. Power Sources 2009, 186, 334–338.
(52)
Fan, Y.; Tongren, D.; Cornelius, C. J. The Role of a Metal Ion within Nafion upon Its Physical and Gas Transport Properties. Eur. Polym. J. 2014, 50, 271–278.
(53)
Daly, K. B.; Panagiotopoulos, A. Z.; Debenedetti, P. G.; Benziger, J. B. Viscosity of Nafion Oligomers as a Function of Hydration and Counterion Type : A Molecular Dynamics Study. J. Phys. Chem. B 2014, 118, 13981–13991.
(54)
Eikerling, M. H.; Berg, P. Poroelectroelastic Theory of Water Sorption and Swelling in Polymer Electrolyte Membranes. Soft Matter 2011, 7, 5976–5990.
(55)
Gebel, G. Structural Evolution of Water Swollen Perfluorosulfonated Ionomers from Dry Membrane to Solution. Polymer (Guildf). 2000, 41, 5829–5838.
(56)
Eikerling, M.; Kulikovsky, A. Chapter 2. In Polymer Electrolyte Fuel Cells: Physical Principles of Materials and Operation; 2014, 59–154.
(57)
Koestner, R.; Roiter, Y.; Kozhinova, I.; Minko, S. AFM Imaging of Adsorbed Nafion Polymer on Mica and Graphite at Molecular Level. Langmuir 2011, 27, 10157–10166.
(58)
Vishnyakov, A.; Neimark, A. V. Molecular Dynamics Simulation of Nafion Olygomer Solvation in Equimolar Methanol – Water Mixture. J. Phys. Chem. B 2001, 105, 7830–7834.
(59)
Vishnyakov, A.; Neimark, A. V. Molecular Simulation Study of Nafion Membrane Solvation in Water and Methanol. J. Phys. Chem. B 2000, 104, 4471–4478.
(60)
Wang, C.; Duscher, G.; Paddison, S. J. Characterization of Chain Conformations in Perfluorosulfonic Acid Membranes Using Electron Energy Loss Spectroscopy. RSC Adv. 2015, 5, 2368–2373.
(61)
Liu, S.; Muthukumar, M. Langevin Dynamics Simulation of Counterion Distribution around Isolated Flexible Polyelectrolyte Chains. J. Chem. Phys. 2002, 116, 9975.
(62)
Ou, Z.; Muthukumar, M. Langevin Dynamics of Semiflexible Polyelectrolytes: Rod-Toroid-Globule-Coil Structures and Counterion Distribution. J. Chem. Phys. 2005, 123, 074905.
(63)
Elder, R. M.; Jayaraman, A. Coarse-Grained Simulation Studies of Effects of Polycation Architecture on Structure of the Polycation and Polycation − Polyanion Complexes. Macromolecules 2012, 45, 8083–8096.
(64)
Košovan, P.; Limpouchová, Z.; Procházka, K. Molecular Dynamics Simulation of Time-Resolved Fluorescence Anisotropy Decays from Labeled Polyelectrolyte Chains. Macromolecules 2006, 39, 3458– 3465.
(65)
Mantha, S.; Yethiraj, A. Conformational Properties of Sodium Polystyrenesulfonate in Water: Insights from a Coarse-Grained Model with Explicit Solvent. J. Phys. Chem. B 2015, 119, 11010–11018.
(66)
Košovan, P.; Limpouchová, Z.; Procházka, K. Conformational Behavior of Comb-like Polyelectrolytes in Selective Solvent: Computer Simulation Study. J. Phys. Chem. B 2007, 111, 8605–8611.
(67)
Connolly, R.; Bellesia, G.; Timoshenko, E. G.; Kuznetsov, Y. A.; Elli, S.; Ganazzoli, F. “Intrinsic” and “Topological” Stiffness in Branched Polymers. Macromolecules 2005, 38, 5288–5299.
(68)
Hsu, H. P.; Paul, W.; Binder, K. Standard Definitions of Persistence Length Do Not Describe the Local “Intrinsic” Stiffness of Real Polymer Chains. Macromolecules 2010, 43, 3094–3102.
ACS Paragon Plus Environment
22
Page 23 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
(69)
Yethiraj, A. A Monte Carlo Simulation Study of Branched Polymers. J. Chem. Phys. 2006, 125, 204901.
(70)
Kremer, K.; Grest, G. S. Dynamics of Entangled Linear Polymer Melts: A Molecular-Dynamics Simulation. J. Chem. Phys. 1990, 92, 5057.
(71)
Binder, K. Monte Carlo and Molecular Dynamics Simulations in Polymer Science; Oxford University Press: New York, 1995.
(72)
Dunweg, B.; Kremer, K. Molecular Dynamics Simulation of a Polymer Chain in Solution. Phys. Rev. Lett. 1993, 99, 6983–6997.
(73)
Grest, G. S.; Murat, M. Structure of Grafted Polymeric Brushes in Solvents of Varying Quality : A Molecular Dynamics Study. Macromolecules 1993, 26, 3108–3117.
(74)
Darden, T.; York, D.; Pedersen, L. Particle Mesh Ewald: An NŊ Log (N) Method for Ewald Sums in Large Systems. J. Chem. Phys. 1993, 98, 10089.
(75)
Grest, G. S.; Kremer, K. Molecular Dynamics Simulation for Polymers in the Presence of a Heat Bath. Phys. Rev. A 1986, 33, 3628–3631.
(76)
Plimpton, S. Fast Parallel Algorithms for Short-Range Molecular Dynamics. J. Comput. Phys. 1995, 117, 1– 19. http://lammps.sandia.gov.
(77)
Reddy, G.; Yethiraj, A. Implicit and Explicit Solvent Models for the Simulation of Dilute Polymer Solutions. Macromolecules 2006, 39, 8536–8542.
(78)
Reddy, G.; Yethiraj, A. Solvent Effects in Polyelectrolyte Adsorption: Computer Simulations with Explicit and Implicit Solvent. J. Chem. Phys. 2010, 132, 074903.
(79)
Park, S.; Yethiraj, A. Atomistic Simulations of Polyelectrolyte Solutions. J. Phys. Chem. B 2012, 116, 1–30.
(80)
Kramarenko, E. Y.; Pevnaya, O. S.; Khokhlov, A. R. Stoichiometric Polyelectrolyte Complexes as Comb Copolymers. J. Chem. Phys. 2005, 122, 084902.
(81)
E. Theodorakis, P.; G. Fytas, N. Molecular Dynamics Simulations of Bottle-Brush Polymers with a Flexible Backbone under Theta and Good Solvent Conditions. Am. J. Condens. Matter Phys. 2012, 2, 101–108.
(82)
Malek, K.; Eikerling, M.; Wang, Q.; Navessin, T.; Liu, Z. Self-Organization in Catalyst Layers of Polymer Electrolyte Fuel Cells. J. Phys. Chem. C 2007, 111, 13627–13634.
(83)
Malek, K.; Eikerling, M.; Wang, Q.; Liu, Z.; Otsuka, S.; Akizuki, K.; Abe, M. Nanophase Segregation and Water Dynamics in Hydrated Nafion: Molecular Modeling and Experimental Validation. J. Chem. Phys. 2008, 129, 204702.
(84)
Allen, M. P.; Tildesley, D. J. Computer Simulation of Liquids; Clarendon Press: New York, NY, USA, 1989.
ACS Paragon Plus Environment
23
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 24 of 33
Insert Table of Contents Graphic and Synopsis Here.
ACS Paragon Plus Environment
24
Page 25 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
Figure 1. Schematic illustration of comb-polyelectrolyte architecture. Backbone and side chain hydrophobic beads are shown in blue and green, respectively. Red beads represent anionic headgroups and yellow beads represent the counterions.
ACS Paragon Plus Environment
25
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 26 of 33
Table 1. List of system paramters and their baseline values along with the explored range.
Parameter
Baseline value
Parameter range explored in this work
LS
3
1-7
f
1/2
1 - 1/5
λB
3.0 σ
0.5 - 12 σ
kθ
0 kBT
0 - 45 kBT
ZC
1 e-
1 - 3 e-
ACS Paragon Plus Environment
26
Page 27 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
Figure 2. Backbone Rg for (a) neutral comb-like polymers, (b) comb-polyelectrolyte chains, and (c) pure electrostatic contribution to the backbone Rg.
ACS Paragon Plus Environment
27
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 28 of 33
Figure 3. Equilibrated comb-polyelectrolyte chains at different values of f and LS. Color-coding is the same as Figure 1.
ACS Paragon Plus Environment
28
Page 29 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
Figure 4. Effect of
kθ on backbone Rg. Comb-polyelectrolyte: (a) LS = 3 (b) f = 1/2. Neutral counterpart: (c) LS = 3 (d) f = 1/2. All
the plots share the legend of (a).
ACS Paragon Plus Environment
29
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 30 of 33
Figure 5. Local persistent length for the comb-polyelectrolyte and the neutral branched polymer at different side chain fractions and lengths.
ACS Paragon Plus Environment
30
Page 31 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
Figure 6. Radial distribution functions between anionic headgroups and counterions for (a) f = 1-1/5 and LS = 3 (b) f = 1/2 and LS = 1-7 (c) f = 1/2, LS = 3, and
kθ = 0 - 45 kBT.
ACS Paragon Plus Environment
31
The Journal of Physical Chemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Figure 7. Rg values as a function of Bjerrum length (a) f = 1/2 and LS = 1-7 (b) f = 1/2, LS = 3, and
Page 32 of 33
kθ = 0 - 45 kBT.
ACS Paragon Plus Environment
32
Page 33 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
Figure 8. Values of square root of backbone radius of gyration at different ZC (a) f = 1-1/5 and LS = 3 (b) f = 1/2 and LS = 1-7.
ACS Paragon Plus Environment
33