Structure and dynamics of water surrounding the poly(methacrylic acid): A molecular dynamics study

Conformations and Solvation of Synthetic Polymers in Water by Generalized Born Implicit-Solvent Molecular Dynamics Simulations: Stereoisomers of Poly(acrylic acid) and Poly(methacrylic acid)

We present the GB-OBC model as an approach for implicit-solvent MD simulations of a synthetic macromolecule in water. The model is tested and found to be successful in reproducing the chain dimensions and predicting the free energy of solvation of carboxylic acid vinyl polymers. The influence of stereochemistry and the hydrophobic nature of the polymer was investigated as a function of chain length (20 < N < 600) for poly(acrylic acid) (PAA) and poly(methacrylic acid) (PMA). The dimensionless parameters of the GB-OBC model were parameterized to be applicable to PAA and PMA. Scaling relations for chain dimensions obtained using implicit-solvent MD simulations in this study are in good agreement with those from experiments, theory of solvated chains in good solvents, and all-atom MD simulations in explicit water. Results show that ⟨R_g²⟩/NL² is greater for the atactic chain as compared to the isotactic chain, for PAA as well as PMA. ⟨R_g²⟩/NL² values of chains attain constancy in water for N = 200, with the values being greater for PMA. The PMA chain is conformationally more perturbed than the PAA chain, for both isotactic and atactic stereochemistry. The solvation free energy ΔG_hyd of PAA and PMA in water is negative for all chain lengths (N = 20-600) and becomes more favorable with an increase in molecular weight. The ΔG_hyd values for isotactic and atactic chains are identical at lower values of N but differ slightly for N > 300. Irrespective of the hydrophobic nature of the polymer, the atactic chain is thermodynamically more soluble in water as compared to the isotactic chain. The isotactic chain is less hydrophilic as compared to the atactic chain due to the closer proximity of the COOH groups along the backbone. This implicit solvent method is an effective way to accurately simulate the configurational properties and solvation of synthetic polymers in water.

The Role of Sulphonic and Phosphoric Pendant Groups on the Diffusion of Monovalent Ions in Polyelectrolyte Membranes: A Molecular Dynamics Study

Lithium-ion consumption has risen significantly in recent years due to its use in portable devices. Alternative sources of lithium, which include the recovery from brine using the sustainable and eco-friendly electrodialysis technology, has been explored. This technology, however, requires effective cation-exchange membranes that allow the selective permeation of lithium ions. In this study, we have investigated, via molecular dynamics simulations, the role of the two common charged groups, the sulfonic and the phosphoric groups, in promoting the adsorption of monovalent ions from brine comprising Li+, Na+, Mg2+, and Ca2+ ions. The analysis of the mean square displacement of the ions revealed that Li+ and Na+ ions exhibit superior diffusion behaviors within the polyelectrolyte system. The O-atoms of the charged groups bind strongly with the divalent ions (Mg2+ and Ca2+), which raises their diffusion energy barrier and consequently lowers their rate of permeation. In contrast, the monovalent ions exhibit weaker interactions, with Na+ being slightly above Li+, enabling the permeation of Li+ ions. The present study demonstrates the role of both charged groups in cation-exchange membranes in promoting the diffusion of Li+ and Na+ ions, and could serve as a guide for the design of effective membranes for the recovery of these ions from brine.

Adsorption Conformation of Comb-Shaped Polycarboxylate Ethers on Ettringite (100) Surface: An Atomic Scales Simulation

The adsorption of five polycarboxylate ethers (PCEs) oligomers with different side-chain number and side-chain length on the ettringite (100) surface in explicit solution is studied by all-atom molecular dynamics (MD) simulations. The adsorption conformations, adsorption energies, the radial distribution functions (RDF) between PCEs and ettringite surface and density distribution profile of water perpendicular to the substrate are analyzed. After dynamic equilibrium, negatively charged carboxylate groups are absorbed on the surface of the disordered ettringite crystal and the side-chains are extended to solution. The influence of the number of side-chain on the adsorption strength of PCEs on the ettringite (100) surface is more significant than that of the length of side-chain. The less number of grafted side-chains, the weaker electrostatic shielding, the electrostatic interaction between PCEs and ettringite (100) surface is stronger, which is favorable for the adsorption. The conformation of adsorbed PCEs is closely related to the length of side chain. The greater the length of side-chains is, the larger will be the coverage area of PCE on surface and the degree of water reduction. The mechanism of adsorption of PCEs on ettringite surface is concluded and the major contribution to the adsorption includes ion pairs, hydrogen bonds and an entropic compensation.

Effect of sodium gluconate on molecular conformation of polycarboxylate superplasticizer studied by the molecular dynamics simulation

Sodium gluconate (SG) has been accepted as one of the main additional components in polycarboxylate superplasticizer (PCE) system, due to its excellent retarding effect. While the negative effect on dispersion of PCE was reported in the literature, the reason was not completely revealed. In this study, molecular dynamics simulation was used to investigate the mutual influence between SG and PCE in calcium hydroxide (CH) solution. Radial distribution function (RDF) was used to analyze the effects of SG on the complexation of PCE with Ca²⁺. Radius of gyration (Rg) was adopted to characterize the conformations of the backbone and side chains of PCE in CH solution. Finally, several adsorption and dispersion models were proposed. The results showed that the presence of SG would perturb adsorption of PCE, which was one of the main reasons that affected the dispersion ability of PCE. SG could preferentially combine with Ca²⁺ so that less amount of Ca²⁺ is available for combination of PCE, and this could extend the main chain of PCE and show advantage for PCE adsorption. Besides, adding SG could squeeze the side chains of PCE, which would put a negative effect on the dispersion. These findings gave deeper insight into understanding the dispersion mechanism of PCE-SG system.

Molecular dynamics study of coil-to-globule transition in a thermo-responsive oligomer bound to various surfaces: hydrophilic surfaces stabilize the coil form

The structural and dynamical properties of 40-mer of thermo-responsive polymer PNIPAM covalently bound to different surfaces have been studied, at different temperatures, by means of molecular dynamics simulations. Evolution of the radius of gyration, Rg, of the polymer chain and radial distribution functions (RDFs) calculated for the carbon atoms of the PNIPAM backbone with water oxygens and for the hydrogen atom of the amide groups with water oxygens indicate that functionalized surfaces affect the coil-to-globule transition of PNIPAM, by means of electrostatic interactions, increasing the lower critical solution temperature (LCST) of the polymer. Such interactions, mainly represented by a H-bond, hinder the transition in the globular form while hydrophobic groups on the surface, such as -OCH3, contribute to the globular collapse. A significant alteration in the arrangement of water molecules around the polymer is testified by: (i) the absence of the second peak in the RDF between the C atoms of the PNIPAM backbone and the O atoms of water at the same temperature at which the radius of gyration decreases; (ii) the height of both the first and the second peak of the RDF between the H atom of the amide groups and water O atoms decreases when the temperature increases above the LCST. Finally, the H-bond autocorrelation function indicates that: (i) hydrogen bonds between the bound-to-surface PNIPAM acceptor groups (O[double bond, length as m-dash]C[double bond splayed right]) and the H atoms of water molecules are less persistent than H-bonds formed between the free PNIPAM acceptor groups and water; (ii) H-bonds between the PNIPAM acceptor groups and hydroxyl groups on the quartz surface are longer lived than those formed on graphene oxide.

Binding of calcium cations with three different types of oxygen-based functional groups of superplasticizers studied by atomistic simulations

This work investigated interactions between calcium cations (Ca²⁺) and three common types of oxygen-based functional groups of concrete superplasticizers using density functional theory (DFT) calculations and all-atom molecular dynamics (MD) simulations. The three common types of oxygen-based functional groups were modeled as three hypothetical, low-molecular-weight organic molecules, each containing a methyl-terminated oxyethylene dimer and an adsorbing head of two oxygen-based functional groups, and are referred to as carboxylate, sulfonate, and phosphate groups, respectively, following the usual terminology in the field of concrete admixtures. Our DFT results show that the binding strength of the three groups with calcium cations follows (from high to low) phosphate>carboxylate>sulfonate, and both the electrophilic attack and the chemical reactivity of the three groups contribute significantly to the binding strength. The MD simulation results indicate that the adsorption of the three small molecules on the calcite (1 0 4) surface in aqueous solution shares a similar pattern in the sense that just two oxygen atoms of two adjacent anchor groups adsorb on the calcium atoms on the top layer of the crystal. The adsorption strength among the three types of functional groups follows the same order as the binding strength obtained from DFT calculations; both results corroborate a similar rule-of-thumb established by experiments. Furthermore, interactions of the three types of groups with water molecules suggest that strong hydrogen-bonding interactions exist in those systems. Graphical abstract Binding of calcium cations with three different types of oxygen-based functional groups of superplasticizersᅟ.

Adsorption of organic molecules on mineral surfaces studied by first-principle calculations: A review

First-principle calculations, especially by the density functional theory (DFT) methods, are becoming a power technique to study molecular structure and properties of organic/inorganic interfaces. This review introduces some recent examples on the study of adsorption models of organic molecules or oligomers on mineral surfaces and interfacial properties obtained from first-principles calculations. The aim of this contribution is to inspire scientists to benefit from first-principle calculations and to apply the similar strategies when studying and tailoring interfacial properties at the atomistic scale, especially for those interested in the design and development of new molecules and new products.

Hydrogen bond directed surface dynamics at tactic poly(methyl methacrylate)/water interface

The complexity of induced ordering for tactic poly(methyl methacrylate) (PMMA) thin films in contact with water is examined through all-atom molecular dynamics with validated potentials. We observe that for the water molecules that are hydrogen bonded to the PMMA surface, the isotactic and atactic PMMA show a 33% longer relaxation time compared to syndiotactic PMMA. Almost 94% of hydrogen bonds are with the carbonyl groups of PMMA, irrespective of temperature and tacticity. The stability in re-orientation and nature of hydrogen bond participation for the carbonyl groups as well as about 20% higher interaction energies of carbonyl group hydrogen bonded with water for atactic form indicates existence of cooperative effects. Quantifying the dynamics of hydrogen bond at the tactic interface is important in understanding the role tacticity plays in controlling adhesion and biocompatibility, a design choice that has been gaining ground in the soft material science community.

Effect of methyl groups on conformational properties of small ionized comb-like polyelectrolytes at the atomic level

Comb-like polycarboxylate ether (PCE) molecules with different content of methyl groups substituted on backbone and different location of methyl groups substituted on the side chains, respectively, were designed and were studied in explicit salt solutions by all-atom molecular dynamics simulations. Methyl groups substituted on the backbone of PCE have a great effect on the conformation of PCE. Stiffness of charged backbone was not only affected by the rotational freedom but also the electrostatic repulsion between the charged COO^- groups. The interaction of counterions (Na⁺) with COO^- groups for PCE3 (with part of AA substituted by MAA on the backbone) was stronger and the screen effect was great, which decided the smaller size of PCE3. The interaction between water and COO^- groups was strong regardless of the content of AA substituted by MAA on the backbone. The effect of methyl groups substituted on the different location of side chains on the conformation of PCE was less than that of methyl groups substituted on the backbone. The equilibrium sizes of the four PCE molecules with methyl groups substituted on the side chains were similar. Graphical Abstract Effect of methyl groups on conformational properties of small ionized comb-like polyelectrolytes at the atomic level.

Investigation of pH-induced conformational change and hydration of poly(methacrylic acid) by analytical ultracentrifugation

Analytical ultracentrifugation was performed on poly(methacrylic acid) (PMAA) with a series of weight average molar masses (Mw) in aqueous solutions as a function of pH. The scales of the sedimentation coefficient (s) and the diffusion coefficient (D) to Mw at infinite dilutions were obtained at different pH values, indicating that PMAA chains adopt a collapsed structure at low pH values, and stretch at pH higher than 5.2. Our results show that the sedimentation coefficient exhibits a minimum at pH ∼ 6.0, presumably due to the effect of the conformational change and the hydration state of PMAA chains. When pH increases from 6.0 to 8.5, PMAA chains with high molar mass shrink a little bit, presumably because the sodium ions act as a bridging agent between nonadjacent carboxylate groups. Furthermore, the weight average molar mass of PMAA at pH 8.5 increases by one fold than that at pH 4.0, indicating the condensation of sodium ions and the increase in the number of hydration water molecules around carboxylate groups at high pH values.

The x-ray absorption spectroscopy model of solvation about sulfur in aqueous L-cysteine

The environment of sulfur in dissolved aqueous L-cysteine has been examined using K-edge x-ray absorption spectroscopy (XAS), extended continuum multiple scattering (ECMS) theory, and density functional theory (DFT). For the first time, bound-state and continuum transitions representing the entire XAS spectrum of L-cysteine sulfur are accurately reproduced by theory. Sulfur K-edge absorption features at 2473.3 eV and 2474.2 eV represent transitions to LUMOs that are mixtures of S-C and S-H σ∗ orbitals significantly delocalized over the entire L-cysteine molecule. Continuum features at 2479, 2489, and 2530 eV were successfully reproduced using extended continuum theory. The full L-cysteine sulfur K-edge XAS spectrum could not be reproduced without addition of a water-sulfur hydrogen bond. Density functional theory analysis shows that although the Cys(H)S⋯H-OH hydrogen bond is weak (∼2 kcal) the atomic charge on sulfur is significantly affected by this water. MXAN analysis of hydrogen-bonding structures for L-cysteine and water yielded a best fit model featuring a tandem of two water molecules, 2.9 Å and 5.8 Å from sulfur. The model included a S(cys)⋯H-O(w1)H hydrogen-bond of 2.19 Å and of 2.16 Å for H(2)O(w1)⋯H-O(w2)H. One hydrogen-bonding water-sulfur interaction alone was insufficient to fully describe the continuum XAS spectrum. However, density functional theoretical results are convincing that the water-sulfur interaction is weak and should be only transient in water solution. The durable water-sulfur hydrogen bond in aqueous L-cysteine reported here therefore represents a break with theoretical studies indicating its absence. Reconciling the apparent disparity between theory and result remains the continuing challenge.

Vibrational spectra of proximal water in a thermo-sensitive polymer undergoing conformational transition across the lower critical solution temperature

The vibrational spectrum of water near a thermo-sensitive polymer poly(N-isopropylacrylamide) (PNIPAM) undergoing conformational transition through the lower critical solution temperature (LCST) is calculated using molecular dynamics simulations. The characteristic structural features observed at the atomic scale for these proximal water molecules in a solvated polymer chain while undergoing the conformational transition are strongly correlated to their vibrational densities of states. Comparison of the vibrational spectrum below LCST for the proximal water with the vibrational spectrum obtained for bulk water reveals a significant fraction of the hydrogen bonding between the proximal water molecules and the polymer side groups. Hydrogen-bonded bridges of water molecules are formed between two adjacent and alternate monomers. This network of hydrogen bonding results in formation of locally ordered water molecules at temperatures below the LCST. Analysis of the simulation trajectories confirms the presence of a quasi-stable solvation structure near the PNIPAM. The calculated vibrational spectra for proximal water above the LCST suggest significantly reduced hydrogen bonding with the polymer and indicate a reduction in the structural stability of proximal water around a collapsed polymer chain. Systematic trends in the observed peak intensities and frequency shifts at the low- and high-frequency ends of the spectrum can be correlated with the structural and dynamical changes of water molecules below and above the LCST transition, respectively, for various polymer chain lengths. The simulations reveal that, compared to bulk water, the libration bands are blue shifted and OH stretch bands red shifted for water in proximity to PNIPAM with 30 monomer units below the LCST. The simulations suggest that vibrational spectra can be used as a predictive tool for quantifying atomic-scale structural transitions in solvation of thermo-sensitive polymers such as PNIPAM.

Dynamical property of water droplets of different sizes adsorbed onto a poly(methyl methacrylate) surface

A molecular dynamics approach has been employed to study the dynamical behavior of a water droplet adsorbed on a poly(methyl methacrylate) (PMMA) surface. Several sizes of water droplets are considered in order to understand the size influence of the droplet on the dynamical properties of water molecules on the PMMA substrate. The local density profile of water molecules in the droplet upon impact with the PMMA surface is calculated, and the result shows an increase in water penetration with a decrease in the size of the droplet. By examining the velocity field, the regular motion of the water droplet is found both during the equilibrium process and after the droplet reaches the equilibrium state. The dynamical behavior of water molecule is studied by the velocity autocorrelation function (VACF) in different regions for different sizes of water droplets. The result shows that VACFs in different regions are significantly influenced for the droplet with 500 water molecules than for that with 2000 water molecules. Calculations in different regions are made for the vibrational spectrum of the oxygen atom, as well as for hydrogen bond dynamics, the lifetime, and the relaxation time of the hydrogen bond. The changes in the hydrogen bond dynamics are consistent with the change in the distribution of the hydrogen bond angle. We conclude that the dynamical properties of the water molecule are significantly affected by the region relative to the surface but only weakly influenced by the size of the droplet.