Machine learning enabled quantification of the hydrogen bonds inside the polyelectrolyte brush layer probed using all-atom molecular dynamics simulations

All-atom molecular dynamics simulations of polymer and polyelectrolyte brushes

Densely grafted polymer and polyelectrolyte (PE) brushes, owing to their significant abilities to functionalize surfaces for a plethora of applications in sensing, diagnostics, current rectification, surface wettability modification, drug delivery, and oil recovery, have attracted significant attention over the past several decades. Unfortunately, most of the attention has primarily focused on understanding the properties of the grafted polymer and the PE chains with little attention devoted to studying the behavior of the brush-supported ions (counterions needed to screen the PE chains) and water molecules. Over the past few years, our group has been at the forefront of addressing this gap: we have employed all-atom molecular dynamics (MD) simulations for studying a wide variety of polymer and PE brush systems with specific attention to unraveling the properties and behavior of the brush-supported water molecules and ions. Our findings have revealed some of the most fascinating properties of such brush-supported ions and water molecules, including the most remarkable control of nanofluidic transport afforded by the specific ion and water responses induced by the PE brushes grafted on the inner walls of the nanochannel. This feature article aims to summarize some of our key contributions associated with such atomistic simulations of polymer and PE brushes and brush-supported water molecules and counterions.

Hydrophilic and Apolar Hydration in Densely Grafted Cationic Brushes and Counterions with Large Mobilities

We employ an all-atom molecular dynamics (MD) simulation framework to unravel water microstructure and ion properties for cationic [poly(2-(methacryloyloxy)ethyl) trimethylammonium chloride] (PMETAC) brushes with chloride ions as counterions. First, we identify locally separate water domains (or first hydration shells) each around {N(CH3)3}+ and the C═O functional groups of the PMETAC chain and one around the Cl- ion. These first hydration shells around the respective moieties overlap, and the extent of the overlap depends on the nature of the species triggering it. Second, despite the overlap, the water molecules in these domains demonstrate disparate properties dictated by the properties of the atoms and groups around which they are located. For example, the presence of the methyl groups makes the {N(CH3)3}+ group trigger apolar hydration as evidenced by the corresponding orientation of the dipole of the water molecules around the {N(CH3)3}+ moiety. These water molecules around the {N(CH3)3}+ group also have enhanced tetrahedrality compared to the water molecules constituting the hydration layer around the C═O group and the Cl- counterion. Our simulations also identify that there is an intervening water layer between the Cl- ion and {N(CH3)3}+ group: this layer prevents the Cl- ion from coming very close to the {N(CH3)3}+ group. As a consequence, there is a significantly large mobility of the Cl- ions inside the PMETAC brush layer. Furthermore, the C═O group of the polyelectrolyte (PE) chain, due to the partial negative charge on the oxygen atom and the specific structure of the PMETAC brush system, demonstrates strongly hydrophilic behavior and enforces a specific dipole response of water molecules analogous to that experienced by water around anionic species of high charge density. In summary, our findings confirm that PMETAC brushes undergo hydrophilic hydration at one site and apolar hydration at another site and ensure large mobility of the supported Cl- counterions.