Coarse-grained model of water diffusion and proton conductivity in hydrated polyelectrolyte membrane

How fork-length asymmetry affects solvent connectivity and diffusion in grafted polymeric model membranes

The hydrophilic pore morphology and solvent diffusion within model (amphiphilic) polymer membranes are simulated by dissipative particle dynamics (DPD). The polymers are composed of a backbone of 18 covalently bonded A beads to which at regular intervals side chains are attached. The side chains are composed of linear Ap chains (i.e., -A1-A2…Ap) from which two branches, [AsC] and [ArC], split off (s ≤ r). C beads serve as functionalized hydrophilic pendent sites. The branch lengths (s + 1 and r + 1) are varied. Five repeat unit designs (with general formula A3[Ap[AsC][ArC]]) are considered: A2[A3C][A3C] (symmetric branching), A2[A2C][A4C], A2[AC][A5C], A2[C][A6C] (highly asymmetric branching), and A4[AC][A3C]. The distribution of water (W) and W diffusion through nanophase segregated hydrophilic pores is studied. For similar primary length p, an increase in side chain symmetry favors hydrophilic pore connectivity and long-range water transport. C beads located on the longer [ArC] branches reveal the highest C bead mobility and are more strongly associated with water than the C beads on the shorter [AsC] branches. The connectivity of hydrophilic (W and W + C) phases through mapped replica of selected snapshots obtained from Monte Carlo tracer diffusion simulations is in line with trends found from the W bead diffusivities during DPD simulations. The diffusive pathways for protons (H+) in proton exchange membranes and for hydronium (OH-) in anion exchange membranes are the same as for solvents. Therefore, control of the side chain architecture is an interesting design parameter for optimizing membrane conductivities.

Simulated and Experimental Trends Regarding Water Uptake in Polymeric Electrolyte Membranes

Polymeric membranes in an anion or a proton exchange membrane fuel cell need sufficient hydration in order to provide a high hydroxide ion or proton conductivity. The water uptake for six model ionomer membranes, all of the same ion exchange capacity, is modeled by dissipative particle dynamics. The architectures cover three types of families that are of potential interest in fuel cell membrane research. All architectures consist of connected hydrophobic backbone A beads, to which side chains are grafted. For the type I family, the hydrophilic (functional) C beads are pendent on (amphiphilic) [AxC] side chains. The type II architecture contains both hydrophobic [A4] and short hydrophilic [C] side chains. For type III, the C beads are embedded along various locations within the [AxCAy] side chains (x + y = constant). For similar equilibrium time, the membrane water volume fraction increases with side chain length x for type I, and for type III, it increases with the distance x that C beads are separated from the backbone. Among the architectures (types I and III) for which the number of covalent C-A bonds are the same, the water uptake increases with the average number of A-A and A-C bonds (dpd springs) between A beads and the nearest C bead. A picture emerges in which for similar ion exchange capacity model membranes water uptake increases as a function of ⟨Nbondphob-phyl⟩.

Extraction of Polyphenols from Slovenian Hop (Humulus lupulus L.) Aurora Variety Using Deep Eutectic Solvents: Choice of the Extraction Method vs. Structure of the Solvent

Polyphenols from Slovenian hops (Humulus lupulus L.) of the Aurora variety were extracted by different methods and using classical solvents and several deep eutectic solvents (DES) based on choline chloride as the hydrogen bond acceptor component. The obtained extract solutions were analyzed by HPLC for the content of extracted α- and β-acids and extracted xanthohumol. It was found that choline chloride:phenol DES concentrated aqueous solution had an extraction efficiency close to that of diethyl ether, which is considered one of the best classical extraction solvents for polyphenols from hops. The comparison of the extraction efficiency with other choline chloride-based DESs showed that the chemical similarity of the phenol ring in the solvent DES with the polyphenols in hops may be crucial for a highly efficient extraction with choline chloride:phenol DES. On the other hand, the choice of extraction method and the viscosity of the solvents tested seem to play only a minor role in this respect. As far as we know, this is the first study to attempt to relate extraction efficiency in the extraction of hydrophobic solutes to the compressibility of the DES extractants, the latter of which may be correlated with the extent of hydrophobic hydration around the DES components. In addition, using the heating and stirring method for the preparation of choline chloride-based DES concentrated aqueous solutions we found no support for the occurrence of water in two different roles (in the structural and in the dilution role) in these solvents.

Changes in Ion Concentrations upon the Binding of Short Polyelectrolytes on Phospholipid Bilayers: Computer Study Addressing Interesting Physiological Consequences

This computer study was inspired by the experimental observation of Y. Qian et al. published in ACS Applied Materials and Interfaces, 2018 that the short positively charged β-peptide chains and their oligomeric analogues efficiently suppress severe medical problems caused by antimicrobial drug-resistant bacteria despite them not penetrating the bacterial membrane. Our coarse-grained molecular dynamics (dissipative particle dynamics) simulations confirm the tentative explanation of the authors of the experimental study that the potent antimicrobial activity is a result of the entropically driven release of divalent ions (mainly magnesium ions essential for the proper biological function of bacteria) into bulk solution upon the electrostatic binding of β-peptides to the bacterial membrane. The study shows that in solutions containing cations Na⁺, Ca²⁺ and Mg²⁺, and anions Cl⁻, the divalent cations preferentially concentrate close to the membrane and neutralize the negative charge. Upon the addition of positively charged oligomer chains (models of β-peptides and their analogues), the oligomers electrostatically bind to the membrane replacing divalent ions, which are released into bulk solvent. Our simulations indicate that the entropy of small ions (which controls the behavior of synthetic polyelectrolyte solutions) plays an important role in this and also in other similar biologically important systems.

Perspective: Morphology and ion transport in ion-containing polymers from multiscale modeling and simulations

Ion-containing polymers are soft materials composed of polymeric chains and mobile ions. Over the past several decades they have been the focus of considerable research and development for their use as the electrolyte in energy conversion and storage devices. Recent and significant results obtained from multiscale simulations and modeling for proton exchange membranes (PEMs), anion exchange membranes (AEMs), and polymerized ionic liquids (polyILs) are reviewed. The interplay of morphology and ion transport is emphasized. We discuss the influences of polymer architecture, tethered ionic groups, rigidity of the backbone, solvents, and additives on both morphology and ion transport in terms of specific interactions. Novel design strategies are highlighted including precisely controlling molecular conformations to design highly ordered morphologies; tuning the solvation structure of hydronium or hydroxide ions in hydrated ion exchange membranes; turning negative ion-ion correlations to positive correlations to improve ionic conductivity in polyILs; and balancing the strength of noncovalent interactions. The design of single-ion conductors, well-defined supramolecular architectures with enhanced one-dimensional ion transport, and the understanding of the hierarchy of the specific interactions continue as challenges but promising goals for future research.

Anion Exchange Membranes for Fuel Cells Based on Quaternized Polystyrene-b-poly(ethylene-co-butylene)-b-polystyrene Triblock Copolymers with Spacer-Sidechain Design

This work studied the polystyrene-b-poly(ethylene-co-butylene)-b-polystyrene (SEBS) triblock copolymers functionalized by butyl quaternary ammonium (C₄Q) groups and alkyl side chains of different chain lengths (C_n, n = 0 to 24). The hydrated membrane morphology was modeled by dissipative particle dynamics simulation at hydration levels from 10 to 30. A hydroxide model was devised to characterize the diffusivity of anions under the coarse-grained framework. In general, the ionomers with alkyl side chains provided ion conductivity of a similar level at a lower ion exchange capacity. All hydrated SEBS–C₄Q–C_n ionomers showed clear phase separation of the hydrophobic and hydrophilic domains, featuring 18.6 mS/cm to 36.8 mS/cm ion conductivity. The hydrophilic channels expanded as the water content increased, forming more effective ion conductive pathways. Introducing excess alkyl side chains enhanced the nano-segregation, leading to more ordered structures and longer correlation lengths of the aqueous phase. The membrane morphology was controlled by the length of alkyl side-chains as well as their tethering positions. Ionomers with functionalized side chains tethered on the same block resulted in well-connective water networks and higher conductivities. The detailed structural analysis provides synthesis guidelines to fabricate anion exchange membranes with improved performances.

Coarse-Grained Modeling of Ion-Containing Polymers

Ion-containing polymers have continued to be an important research focus for several decades due to their use as an electrolyte in energy storage and conversion devices. Elucidation of connections between the mesoscopic structure and multiscale dynamics of the ions and solvent remains incompletely understood. Coarse-grained modeling provides an efficient approach for exploring the structural and dynamical properties of these soft materials. The unique physicochemical properties of such polymers are of broad interest. In this review, we summarize the current development and understanding of the structure-property relationship of ion-containing polymers and provide insights into the design of such materials determined from coarse-grained modeling and simulations accompanying significant advances in experimental strategies. We specifically concentrate on three types of ion-containing polymers: proton exchange membranes (PEMs), anion exchange membranes (AEMs), and polymerized ionic liquids (polyILs). We posit that insight into the similarities and differences in these materials will lead to guidance in the rational design of high-performance novel materials with improved properties for various power source technologies.

Rhamnolipid Biosurfactants for Oil Recovery: Salt Effects on the Structural Properties Investigated by Mesoscale Simulations

Rhamnolipids (RLs) are biosurfactants produced by Pseudomonas. The biodegradability and the variety of their functionality make them suitable for environmental remediation and oil recovery. We use dissipative particle dynamics simulations to investigate the aggregation behaviors of ionic RL congeners with nonane in various operating conditions. Under zero-salinity conditions, all RL congeners studied here form small ellipsoidal clusters with detectable free surfactants. When salt ions are present, the electrostatic repulsion between the ionized heads is overcome, resulting in the formation of larger aggregates of unique structures. RLs with C10-alkyl tails tend to form elongated wormlike micelles, while RLs with C16-alkyl tails tend to form clusters in spherical symmetry, including vesicles. Di-rhamnolipids (dRLs) require stronger solvation than monorhamnolipids (mRLs) to form clusters, and the resulting size of micelles is decreased. The morphology of the mixed dRL/mRL/oil systems is controlled based on the type of the congeners and the oil contents. In addition, the divalent calcium ions are found to be influential to the structure of the micelles through different mechanisms. For 5 wt % salinity, the ionic RLs can form oil-swollen micelles up to a 1:1 surfactant-to-oil ratio, suggesting that ionic RLs are superb to act as cleaning agents for petroleum hydrocarbons in the marine area. These key findings may guide the design for RL-based washing techniques in enhanced oil recovery.

DPD Modelling of the Self- and Co-Assembly of Polymers and Polyelectrolytes in Aqueous Media: Impact on Polymer Science

This review article is addressed to a broad community of polymer scientists. We outline and analyse the fundamentals of the dissipative particle dynamics (DPD) simulation method from the point of view of polymer physics and review the articles on polymer systems published in approximately the last two decades, focusing on their impact on macromolecular science. Special attention is devoted to polymer and polyelectrolyte self- and co-assembly and self-organisation and to the problems connected with the implementation of explicit electrostatics in DPD numerical machinery. Critical analysis of the results of a number of successful DPD studies of complex polymer systems published recently documents the importance and suitability of this coarse-grained method for studying polymer systems.

Interactions of Crosslinked Polyacrylic Acid Polyelectrolyte Gels with Nonionic and Ionic Surfactants

The morphology and stability of surfactant-loaded polyelectrolyte gels are of great interest for a variety of personal care, cosmetic, and pharmaceutical products. However, the mechanisms of surfactant interactions with gel-forming polymers are poorly understood and experimentally challenging. The aim of this work is to explore in silico the specifics of surfactant absorption within polyelectrolyte gels drawing on the examples of typical non-ionic octaethylene glycol monooctyl ether (C₈E₈) and anionic sodium dodecyl sulfate (SDS) surfactants and polyacrylic acid modified with hydrophobic sidechains mimicking the practically important Carbopol polymer. Using the systematically parameterized coarse-grained dissipative particle dynamics models, we generate and characterize the equilibrium conformations and swelling of the polymer films in aqueous solutions with the surfactant concentrations varied up to the critical micelle concentration (cmc). We discover the striking difference in interactions of Carbopol-like polymers with nonionic and ionic surfactants under mildly acidic conditions. The sorption of C₈E₈ within the polymer film is found substantial. As the surfactant concentration increases, the polymer film swells and, close to cmc, becomes unstable due to the formation and growth of water pockets filled with surfactant micelles. Sorption of SDS at the same bulk concentrations is found much lower, with only about 1% of surfactant mass fraction achieved at cmc. As the SDS concentration increases further, a lamellae structure is formed within the film, which remains stable. Reduced swelling and higher stability indicate better prospects of using SDS-type surfactants with Carbopol-based gels in formulations for detergents and personal care products.

Multiscale simulations of polyzwitterions in aqueous bulk solutions and brush array configurations

Zwitterionic polymers are very promising candidates for antifouling materials that exhibit high chemical stability as compared to polyethylene glycol-based systems. A number of simulation and experimental studies have emerged over recent years for the investigation of sulfobetaine-based zwitterionic polymers. Investigating the structural and thermodynamic properties of such polymers requires access to broad time and length regimes, thus necessitating the development of multiscale simulation strategies. The present article advocates a mesoscopic dissipative particle dynamics (DPD) model capable of addressing a wide range of time and length scales. The mesoscopic force field was developed hand-in-hand with atomistic simulations based on the OPLS force field through a bottom-up parameterization procedure that matches the atomistically calculated strand-length, strand-angle and pair distribution functions. The DPD model is validated against atomistic simulations conducted in this work, and against relevant atomistic simulation studies, theoretical predictions and experimental correlations from the literature. Properties examined include the conformations of SPE polymers in dilute bulk aqueous solution, the density profile and thickness of brush arrays as functions of the grafting density and chain length. In addition, we compute the potential of mean force of an approaching hydrophilic or hydrophobic foulant via umbrella sampling as a function of its position relative to the poly-zwitterion-covered surface. The aforementioned observables lead to important insights regarding the conformational tendencies of grafted polyzwitterions and their antifouling properties.

Dissipative particle dynamics simulations in colloid and Interface science: a review

Dissipative particle dynamics (DPD) is one of the most efficient mesoscale coarse-grained methodologies for modeling soft matter systems. Here, we comprehensively review the progress in theoretical formulations, parametrization strategies, and applications of DPD over the last two decades. DPD bridges the gap between the microscopic atomistic and macroscopic continuum length and time scales. Numerous efforts have been performed to improve the computational efficiency and to develop advanced versions and modifications of the original DPD framework. The progress in the parametrization techniques that can reproduce the engineering properties of experimental systems attracted a lot of interest from the industrial community longing to use DPD to characterize, help design and optimize the practical products. While there are still areas for improvements, DPD has been efficiently applied to numerous colloidal and interfacial phenomena involving phase separations, self-assembly, and transport in polymeric, surfactant, nanoparticle, and biomolecules systems.

Effects of metal-polymer complexation on structure and transport properties of metal-substituted polyelectrolyte membranes

Morphological and transport properties of hydrated metal-substituted Nafion membranes doped with metal ions of different valency and coordination strength are explored using coarse-grained dissipative particle dynamics simulations. To incorporate the effects of metal-polymer complexation, we introduce a novel metal ion complexation model, in which the charged central metal ion is surrounded by dummy sites that coordinate with ligands. The model parameters are determined by matching the metal-ligand running coordination numbers and the diffusion coefficients obtained from atomistic simulations and/or experiments. The increase of valency and coordination strength is found to strongly influence both the morphology and transport characteristics of the membrane at all hydration levels. The membrane segregation into hydrophobic and hydrophilic sub-phases is affected by metal-sulphonate coordination induced crosslinking at the hydrophilic/hydrophobic interface. The simulation results indicate that the interfacial crosslinking influences the interfacial tension and thereby affect the growth and coalescence of water clusters upon the increase of hydration. Multivalent complexation hinders water and ion mobility and causes anomalous sub-diffusion and dramatic decrease of the water permeability and ionic conductivity. Our DPD model is found efficient in elucidating the mechanisms of coordination-induced cross-linking and complexation and predicting on a semi-quantitative level the morphological and transport properties of metal-substituted Nafion membranes depending on the ion valency and coordination strength. The proposed model can be further advanced and adopted for other polyelectrolyte systems, such as sulfonated block-copolymers, polysaccharide solutions and composites, and biopolymer assemblies.

Micellization of Rhamnolipid Biosurfactants and Their Applications in Oil Recovery: Insights from Mesoscale Simulations

The dissipative particle dynamics (DPD) mesoscopic method is used to investigate the self-assembly of rhamnolipid congeners and their aggregation behaviors with paraffins including nonane and pentadecane. The coarse-grained force field is parameterized by combining molecular dynamics (MD) simulations, COSMOtherm calculations, and available experimental data. This model reproduces the vesicular formation of α-l-rhamnopyranosyl-β-hydroxydecanoyl-β-hydroxydecanoate (Rha-C10-C10) reported by all-atom MD simulations. The vesicle composed of Rha-C10-C10 is found to be most stable at a surfactant concentration of 100-146 mM based on asphericity analysis. The architecture of rhamnolipid congeners affects the morphology of their aggregates. Di-rhamno-di-lipidic dRha-C16-C16 forms vesicles with a thicker unilamellar layer of 3.2 nm. Rha-C16-C16 forms vesicles at a lower concentration of 70 mM, but the enclosed water space collapses when the surfactant concentration increases. dRha-C10-C10 forms wormlike micelles, which agglomerate into a torus and interconnected network at higher concentrations. In the presence of alkane molecules, dRha-C10-C10 maintains its wormlike micellar morphology with alkane molecules wrapped inside the aggregates. For Rha-C10-C10, Rha-C16-C16, and dRha-C16-C16, nonane molecules are distributed in the hydrophobic subdomain formed by rhamnolipid molecules. Spherical vesicles are formed at a surfactant concentration of 50 mM and then develop into ellipsoidal vesicles when the concentration increases to 125 mM. When mixed with pentadecane, the alkane molecules are aggregated and surrounded by surfactants forming a core-shell structure at a low surfactant concentration of 20 mM. At higher alkane and surfactant concentrations, the morphologies develop into disk micelles, wormlike micelles, and vesicles, with pentadecane molecules being distributed and packed with rhamnolipids. The obtained simulation results suggest that these biosurfactants have potential as environmental remediation agents.

Temperature-responsive morphology formation of a PS-b-PI copolymer: a dissipative particle dynamics simulation study

Self-assembly responsiveness to stimuli of polystyrene-block-polyisoprene (PS-b-PI) diblock copolymer materials is explored by means of classical molecular dynamics (MD) and dissipative particle dynamics (DPD) simulations. A concerted relationship between the parameters achieved from atomistic and DPD simulations is obtained for this molecular recognition as clearly pronounced in a phase transition. Effects of temperature, model size and composition on the morphological formation were systematically investigated for the diblock copolymeric system. Structural changes resulting in the evolution of rheology as well as an equilibrium ordered structure were analyzed in terms of order parameters and radial distribution functions. From our models, various morphologies were observed including discrete clusters (sphere-liked morphology), connected clusters (gyroid-liked morphology), hexagonally packed cylinders (HEX), connected cylinders, irregular cylinders, perfect lamellae, perforated lamellae and defected lamellae. Based on this finding, a bottom-up multi-scale simulation of the PS-b-PI diblock copolymer provides a link between equilibrium copolymeric morphologies and the crucial parameters.

Designing Highly Conductive Block Copolymer-Based Anion Exchange Membranes by Mesoscale Simulations

Hydroxide ion conductivity is a key aspect of anion exchange membranes and is mainly determined by the nanoscale membrane morphologies. Fundamental understanding of the structural and transport properties of membranes in terms of polymer architectures is crucial for future development of membrane-based applications. Using mesoscale simulations, this work predicts the mesostructure of the hydrated triblock copolymers; the designed polymers are composed of aromatic (polyphenylene oxide, PPO) or aliphatic (polystyrene-ethylene-butylene-styrene, SEBS) backbones, with cationic side chains being modified by hydrophobic or hydrophilic spacers. For PPO-based polymers, using octyl spacers creates a meshlike water network, yielding ion conductivity equal to 30.6 mS/cm at room temperature. For SEBS-based polymers, the nonmodified form is sufficient to produce ion-conducting pathways. Adding hydrophobic spacers further enhances the nanosegregation, and the membranes provide similar conductivity at a lower ion exchange capacity and water content. Adding hydrophilic spacers, however, has negative impacts on the ion transport. The side chains are in the stretched configurations, which sterically hinder the mobility of water and hydroxide ions. Such a resistance can be overcome by adapting multication side-chain designs, where large water channels are formed, yielding ion conductivity as high as 32.8 mS/cm.

A new equation of state for homo-polymers in dissipative particle dynamics

A chain-revised Groot-Warren equation of state (crGW-EOS) was developed and tested to describe systems of homo-oligomeric chains in the framework of dissipative particle dynamics (DPD). First, thermodynamic perturbation theory is applied to introduce correction terms that account for the reduction in pressure with an increasing number of bonds at constant bead number density. Then, this EOS is modified by introducing a set of switching functions that yields an accurate second virial coefficient in the low-density limit. The crGW-EOS offers several improvements over the revised Groot-Warren equation of state and Groot-Warren equation of state for chain molecules. We tested the crGW-EOS by using it to predict the pressure of oligomeric systems and the B₂ virial coefficient of chain DPD particles for a range of bond lengths. Additionally, a method is developed for determining the strength of cross-interaction parameters between chains of different compositions and sizes and for thermal and athermal mixtures. We explored how different levels of coarse-graining affect the upper-critical solution temperature.

Coarse-grained model of nanoscale segregation, water diffusion, and proton transport in Nafion membranes

We present a coarse-grained model of the acid form of Nafion membrane that explicitly includes proton transport. This model is based on a soft-core bead representation of the polymer implemented into the dissipative particle dynamics (DPD) simulation framework. The proton is introduced as a separate charged bead that forms dissociable Morse bonds with water beads. Morse bond formation and breakup artificially mimics the Grotthuss hopping mechanism of proton transport. The proposed DPD model is parameterized to account for the specifics of the conformations and flexibility of the Nafion backbone and sidechains; it treats electrostatic interactions in the smeared charge approximation. The simulation results qualitatively, and in many respects quantitatively, predict the specifics of nanoscale segregation in the hydrated Nafion membrane into hydrophobic and hydrophilic subphases, water diffusion, and proton mobility. As the hydration level increases, the hydrophilic subphase exhibits a percolation transition from a collection of isolated water clusters to a 3D network of pores filled with water embedded in the hydrophobic matrix. The segregated morphology is characterized in terms of the pore size distribution with the average size growing with hydration from ∼1 to ∼4 nm. Comparison of the predicted water diffusivity with the experimental data taken from different sources shows good agreement at high and moderate hydration and substantial deviation at low hydration, around and below the percolation threshold. This discrepancy is attributed to the dynamic percolation effects of formation and rupture of merging bridges between the water clusters, which become progressively important at low hydration, when the coarse-grained model is unable to mimic the fine structure of water network that includes singe molecule bridges. Selected simulations of water diffusion are performed for the alkali metal substituted membrane which demonstrate the effects of the counter-ions on membrane self-assembly and transport. The hydration dependence of the proton diffusivity reproduces semi-qualitatively the trend of the diverse experimental data, showing a sharp decrease around the percolation threshold. Overall, the proposed model opens up an opportunity to study self-assembly and water and proton transport in polyelectrolytes using computationally efficient DPD simulations, and, with further refinement, it may become a practical tool for theory informed design and optimization of perm-selective and ion-conducting membranes with improved properties.

Theoretical analyses on water cluster structures in polymer electrolyte membrane by using dissipative particle dynamics simulations with fragment molecular orbital based effective parameters

The mesoscopic structures of polymer electrolyte membrane (PEM) affect the performances of fuel cells. Nafion® with the Teflon® backbone has been the most widely used of all PEMs, but sulfonated poly-ether ether-ketone (SPEEK) having an aromatic backbone has drawn interest as an alternative to Nafion. In the present study, a series of dissipative particle dynamics (DPD) simulations were performed to compare Nafion and SPEEK. These PEM polymers were modeled by connected particles corresponding to the hydrophobic backbone and the hydrophilic moiety of sulfonic acid group. The water particle interacting with Nafion particles was prepared as well. The crucial interaction parameters among DPD particles were evaluated by a series of calculations based on the fragment molecular orbital (FMO) method in a non-empirical way (Okuwaki et al., J. Phys. Chem. B, 2018, 122, 338–347). Through the DPD simulations, the water and hydrophilic particles aggregated, forming cluster networks surrounded by the hydrophobic phase. The structural features of formed water clusters were investigated in detail. Furthermore, the differences in percolation behaviors between Nafion and SPEEK revealed much better connectivity among water clusters by Nafion. The present FMO-DPD simulation results were in good agreement with available experimental data.

The mesoscopic structures of polymer electrolyte membrane (PEM) affect the performances of fuel cells.

Improving proton conduction pathways in di- and triblock copolymer membranes: Branched versus linear side chains

Phase separation within a series of polymer membranes in the presence of water is studied by dissipative particle dynamics. Each polymer contains hydrophobic A beads and hydrophilic C beads. Three parent architectures are constructed from a backbone composed of connected hydrophobic A beads to which short ([C]), long ([A₃C]), or symmetrically branched A₅[AC][AC] side chains spring off. Three di-block copolymer derivatives are constructed by covalently bonding an A₃₀ block to each parent architecture. Also three tri-blocks with A₁₅ blocks attached to both ends of each parent architecture are modeled. Monte Carlo tracer diffusion calculations through the water containing pores for 1226 morphologies reveal that water diffusion for parent architectures is slowest and diffusion through the di-blocks is fastest. Furthermore, diffusion increases with side chain length and is highest for branched side chains. This is explained by the increase of water pore size with 〈N_bond〉, which is the average number of bonds that A beads are separated from a nearest C bead. Optimization of 〈N_bond〉 within the amphiphilic parent architecture is expected to be essential in improving proton conduction in polymer electrolyte membranes.

Dependence of Solvent Diffusion on Hydrophobic Block Length within Amphiphilic-Hydrophobic Block Copolymer Membranes

Pore networks and water diffusion within model (amphiphilic-hydrophobic) diblock copolymer membranes in the presence of 16 vol % water is studied by dissipative particle dynamics in combination with Monte Carlo tracer diffusion calculations. The amphiphilic block (parent architecture (A[A₃C])₁₀) is composed of a backbone that contains 10 consecutively connected hydrophobic A beads; to each A bead, a side chain is grafted composed of three connected A beads and a pendant hydrophilic C bead. Hydrophobic blocks are constructed from x covalently bonded A beads, with x = 20, 30, or 50. Water diffusion through the pores is modeled by Monte Carlo tracer diffusion within more than 500 mapped morphologies. Long range water diffusion within the amphiphilic-hydrophobic ((A[A₃C])₁₀-A_x) diblock architectures increases with hydrophobic block length. Diffusion increases with Q = ⟨N_bond⟩|C||1 - C|^-1, where C is the hydrophilic C bead fraction and ⟨N_bond⟩ the average number of bonds that A beads are separated from the nearest C bead. These trends are also anticipated for amphiphilic parent architectures (ACA₃)₁₀, (A₂[C]A₂)₁₀, and (A₂[AC]A)₁₀. This is explained by the squeezing of water from the hydrophobic phase into the amphiphilic phase. Two characteristic distances are observed: The shorter distance corresponds to the interpore (or intercluster) separation within the "parent architecture-water" phase and obeys the earlier obtained linear relation between intercluster distance and ⟨N_bond⟩_amphi of the amphiphilic parent architecture. The longer distance is governed by the phase separation between the amphiphilic-water phase and hydrophobic blocks.

Water Diffusion Dependence on Amphiphilic Block Design in (Amphiphilic-Hydrophobic) Diblock Copolymer Membranes

Polyelectrolyte membranes (PEMs) are applied in polyelectrolyte fuel cells (PEFC). The proton conductive pathways within PEMs are provided by nanometer-sized water containing pores. Large-scale application of PEFC requires the production of low-cost membranes with high proton conductivity and therefore good connected pore networks. Pore network formation within four alternative model diblock (hydrophobic_amphiphilic) copolymers in the presence of water is studied by dissipative particle dynamics. Each hydrophobic block contains 50 consecutively connected hydrophobic (A) fragments, and amphiphilic blocks contain 40 hydrophobic A beads and 10 hydrophilic C beads. For one amphiphilic block the C beads are distributed uniformly along the backbone. For the other architectures C beads are located at the end of the side chains attached at regular intervals along the backbone. Water diffusion through the pores is modeled by Monte Carlo tracer diffusion through mapped morphologies. Diffusion is highest for the grafted architectures and increases with increase of length of the side chains. A consistent picture emerges in which diffusion strongly increases with the value of ⟨Nbond⟩ within the amphiphilic block, where ⟨Nbond⟩ is the average number of bonds between hydrophobic A beads and the nearest C bead.