[Paper Review] Hydration of polyelectrolytes studied by molecular dynamics simulation

Oliver Biermann, E. HAEDICKE+3 more|arXiv (Cornell University)|Jan 9, 2001

Electrostatics and Colloid Interactions1 references8 citations

TL;DR

This study uses all-atom molecular dynamics simulations to compare the hydration behavior of sodium carboxymethyl cellulose (CMC) and sodium polyacrylate (PAA) in dilute aqueous solutions. Despite similar molecular weights, PAA forms significantly more hydrogen bonds with water (0.083 bonds/amu) than CMC (0.029–0.036 bonds/amu), leading to stronger solvation and higher counterion association, while CMC’s conformation depends critically on substitution pattern, yielding either extended or compact globular structures.

ABSTRACT

Molecular dynamics simulations of diluted (~2.5 weight percent) aqueous solutions of two polyelectrolytes, namely sodium carboxy methyl cellulose (CMC) and sodium poly(acrylate) (PAA) have been performed. Water and counterions were taken into account explicitly. For CMC the substitution pattern and starting conformation is all-important. Two simulations of CMC oligomers resulted in different structures: One molecule takes a stretched conformation, while the second one keeps a globule-like, toroidal one. PAA is stretched during the whole simulation, with an average characteristic ratio of 8.3. On a local atomistic scale CMC and PAA have different hydrogen-bond properties. The COO- groups of PAA can only act as hydrogen bond acceptors, but due to the high negative charge density there are still more water molecules assembled around PAA than around CMC. There are 0.036 bonds/amu respectively 0.029 bonds/amu to water for the two CMC oligomers, but more than twice as many for PAA: 0.083 bonds/amu. Beside intermolecular hydrogen bonding, there is a significant amount of intramolecular H-bonding for CMC, which is influenced by the COO- groups, which act as strong H-acceptor. In contrast to hydroxy- and carboxylic groups, ether oxygens are hardly involved into hydrogen bonding.

Motivation & Objective

• To understand the atomistic-level hydration structure and dynamics of polyelectrolytes in water.
• To compare the solvation behavior of natural polyelectrolyte CMC with synthetic PAA.
• To investigate how substitution patterns in CMC affect its conformation and hydrogen bonding.
• To quantify differences in hydrogen bonding, counterion association, and solvation shell structure between CMC and PAA.
• To generate atomistic data for future coarse-grained modeling of PAA.

Proposed method

• Perfomed all-atom molecular dynamics simulations of dilute aqueous solutions (≈2.5 wt% of CMC and PAA).
• Explicitly modeled water molecules and counterions (Na+) to capture ion-specific and solvent effects.
• Simulated two CMC oligomers with different random substitution patterns (DS ≈ 0.7) to assess structural diversity.
• Used radial distribution functions (RDFs) to analyze spatial distribution of water and Na+ ions around polymer chains.
• Calculated hydrogen bond statistics (number per amu) and monitored intramolecular and intermolecular H-bonding dynamics.
• Analyzed chain dimensions via characteristic ratio and conformational changes over time.

Experimental results

Research questions

• RQ1 How does the substitution pattern in CMC influence its conformation in aqueous solution?
• RQ2 What are the differences in hydrogen bonding between CMC and PAA with water at the atomistic level?
• RQ3 How does the charge density and ion association differ between CMC and PAA?
• RQ4 What is the impact of intramolecular hydrogen bonding on CMC’s solvation and dynamics?
• RQ5 How do the solvation shells and radial distribution functions of Na+ differ between CMC and PAA?

Key findings

• CMC oligomers adopted distinct conformations: one remained globular and toroidal, while the other stretched into an extended form, depending on substitution pattern.
• PAA maintained a stretched conformation throughout the simulation with a characteristic ratio of 8.3, indicating high chain flexibility.
• PAA formed 0.083 hydrogen bonds per amu with water, more than twice the 0.029–0.036 bonds/amu observed for the two CMC oligomers.
• The COO− groups in PAA acted only as hydrogen bond acceptors, but due to high charge density, they attracted more water molecules than CMC.
• CMC II, with its compact globular structure, formed more intramolecular hydrogen bonds and fewer intermolecular H-bonds with water than CMC I.
• Despite similar molecular weight, PAA attracted approximately 13 Na+ ions within 0.64 nm (screening length), while CMC bound only 0.41 counterions on average.