All-atom molecular dynamics simulations of weak polyionic brushes: influence of charge density on the properties of polyelectrolyte chains, brush-supported counterions, and water molecules

Molecular dynamics simulation of salt diffusion in constituting phosphazene-based polymer electrolyte

A growing demand to visualize polymer models in liquid poses a computational challenge in molecular dynamics (MD) simulation, as this requires emerging models under suitable force fields (FFs) to capture the underlying molecular behaviour accurately. In our present study, we have employed TIP3P potential on water and all atomistic optimized potentials for liquid simulations FFs to study the liquid electrolyte behavior of phosphazene-based polymer by considering its potential use in lithium-ion polymer batteries. We have explored the polymer's local structure, chain packing, wettability, and hydrophobic tendencies against the silicon surface using a combination of a pseudocontinuum model in MD simulation, and surface-sensitive sum frequency generation (SFG) vibrational spectroscopy. The finding yields invaluable insights into the molecular architecture of phosphazene. This approach identifies the importance of hydrophobic interactions with air and hydrophilic units with water molecules in understanding the behavior and properties of phosphazene-based polymers at interfaces, contributing to its advancements in materials science. The MD study uniquely captures traces of the polymer-ion linkage, which is observed to become more pronounced with the increase in polymer weight fraction. The theoretical observation of this linkage's influence on lithium-ion diffusion motion offers valuable insights into the fundamental physics governing the behavior of atoms and molecules within phosphazene-based polymer electrolytes in aqueous environments. Further these predictions are corroborated in the molecular-level depiction at the air-aqueous interface, as evidenced from the OH-oscillator strength variation measured by the SFG spectroscopy.The fundamental findings from this study open new avenues for utilizing MD simulation as a versatile methodology to gain profound insights into intermolecular interactions of polymer. It could be useful in the application of biomedical and energy-related research, such as polymer lithium-ion batteries, fuel cells, and organic solar cells.

Strength, number, and kinetics of hydrogen bonds for water confined inside boron nitride nanotubes

Water has shown a myriad of highly interesting properties and behaviors, such as very low friction, phase transition under unexpected conditions, massive property alterations, etc. inside strong nanoconfinements of few-nanometer to sub-nanometer diameters. Water-water hydrogen bonding is one of the most important factors dictating such water behavior and properties inside such strong nanoconfinements. In this paper, we employ Reactive Force Field (ReaxFF) molecular dynamics (MD) simulations for studying multiple facets of such water-water hydrogen bonds (HBs) inside boron-nitride nanotubes (BNNTs) having diameters ranging from a few nanometers to sub-nanometers. First, the strength of the water-water HB interactions, as a function of the HB configuration, is quantified by studying the corresponding PMF (potential of mean force). For water present in extreme confinements (BNNTs with sub-nanometric diameters), we see completely isolated HB basins. On the other hand, for bulk water the HB basin is connected (via a saddle point) to a nearby second PMF well. Therefore, our analysis successfully distinguishes the HB characteristics between the cases of water in extreme confinement and bulk water. Second, we study the kinetics of such water-water HBs: HBs formed by a given pair of water molecules in extreme confinements show a much larger probability of remaining intact once formed or re-forming after they have been broken. Both these results, which shed new light on water-water hydrogen bonding inside strong nanoconfinements, can be explained by the single-file structure formed by the water molecules in extreme BNNT confinements.

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.

Impact of Confinement and Zwitterionic Ligand Chemistry on Ion-Ion Selectivity of Functionalized Nanopores

Membranes incorporating zwitterionic chemistries have recently emerged as promising candidates for facilitating challenging ion-ion separations. Transport of ions in such membranes predominantly occurs in hydrated nanopores lined with zwitterionic monomers. To shed light on the physics of ion-ion selectivity underlying such materials, we conducted molecular dynamics simulations of sodium halide transport in model nanopores grafted with sulfobetaine methacrylate molecules. Our results reveal that in both functionalized and unfunctionalized nanopores smaller ions prefer to reside near the pore center, while the larger ions tend to reside near the pore walls. An enhancement in the selective transport of larger anions is observed within the unfunctionalized nanopores relative to that in salt-in-water solutions. Upon functionalization of the nanopores with zwitterions (ZIs), the disparities in the anionic distribution profiles within the pores coupled with differences in the anion-ZI interactions result in a slowdown of larger anions relative to smaller anions. Increasing the ZI grafting density exacerbates these effects, further promoting the selective transport of smaller anions. Our results suggest that selectivity toward large anions can be realized by using nanoporous membranes with ZI content that is high enough to facilitate ion/water partitioning into the pores while preserving the characteristic tendency of the unfunctionalized pores to facilitate faster transport of the larger anions. On the other hand, selectivity toward smaller anions can be achieved by targeting ZI content within the pores that is high enough to significantly slow down the transport of large anions but not high enough to hinder the partitioning of ions/water molecules into the pore due to steric effects.

Design, preparation, and characterization of lubricating polymer brushes for biomedical applications

The lubrication modification of biomedical devices significantly enhances the functionality of implanted interventional medical devices, thereby providing additional benefits for patients. Polymer brush coating provides a convenient and efficient method for surface modification while ensuring the preservation of the substrate's original properties. The current research has focused on a "trial and error" method to finding polymer brushes with superior lubricity qualities, which is time-consuming and expensive, as obtaining effective and long-lasting lubricity properties for polymer brushes is difficult. This review summarizes recent research advances in the biomedical field in the design, material selection, preparation, and characterization of lubricating and antifouling polymer brushes, which follow the polymer brush development process. This review begins by examining various approaches to polymer brush design, including molecular dynamics simulation and machine learning, from the fundamentals of polymer brush lubrication. Recent advancements in polymer brush design are then synthesized and potential avenues for future research are explored. Emphasis is placed on the burgeoning field of zwitterionic polymer brushes, and highlighting the broad prospects of supramolecular polymer brushes based on host-guest interactions in the field of self-repairing polymer brush applications. The review culminates by providing a summary of methodologies for characterizing the structural and functional attributes of polymer brushes. It is believed that a development approach for polymer brushes based on "design-material selection-preparation-characterization" can be created, easing the challenge of creating polymer brushes with high-performance lubricating qualities and enabling the on-demand creation of coatings. STATEMENT OF SIGNIFICANCE: Biomedical devices have severe lubrication modification needs, and surface lubrication modification by polymer brush coating is currently the most promising means. However, the design and preparation of polymer brushes often involves "iterative testing" to find polymer brushes with excellent lubrication properties, which is both time-consuming and expensive. This review proposes a polymer brush development process based on the "design-material selection-preparation-characterization" strategy and summarizes recent research advances and trends in the design, material selection, preparation, and characterization of polymer brushes. This review will help polymer brush researchers by alleviating the challenges of creating polymer brushes with high-performance lubricity and promises to enable the on-demand construction of polymer brush lubrication coatings.

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.

Specific Ion and Electric Field Controlled Diverse Ion Distribution and Electroosmotic Transport in a Polyelectrolyte Brush Grafted Nanochannel

Controlling ion distribution inside a charged nanochannel is central to using such channels in diverse applications. Here, we show the possibility of using a charged polyelectrolyte (PE) brush-grafted nanochannel for triggering diverse nanoscopic ion distribution and nanofluidic electroosmotic transport by controlling the valence and size of the counterions (that screen the charges of the PE brushes) and the strength of an externally applied axial electric field. We atomistically simulate separate cases of fully charged polyacrylic acid (PAA) brush functionalized nanochannels with Na⁺, Cs⁺, Ca²⁺, Ba²⁺, and Y³⁺ counterions screening the PE charges. Four key findings emerge from our simulations. First, we find that the counterions with a greater valence and a smaller size prefer to remain localized inside the brush layer. Second, for the case where there is an added chloride salt with the same cation (as the screening counterions), there are more coions (Cl^- ions) in the brush-free bulk than counterions (for counterions Na⁺, Ca²⁺, Ba²⁺, Y³⁺): this is a manifestation of the overscreening (OS) of the PE brush layer. Contrastingly, the number of Cs⁺ ions remain higher than the Cl^- ions inside the brush-free bulk, ensuring that there is no OS effect for this case. Third, large applied electric field enables a few Na⁺, Cs⁺, and Ba²⁺ counterions to leave the brush layer and to go to the bulk: this makes the OS of the PE brush layer disappear for the cases of PE brushes being screened by the Na⁺ and Ba²⁺ ions. On the other hand, no such electric-field-mediated disappearance of OS is observed for the cases of Ca²⁺ and Y³⁺ screening counterions; we attribute this to the firm attachment of these counterions to the negatively charged monomers. Free energy associated with a counterion binding to a PE chain corroborates this diversity in the counterion-specific response to the applied electric field. Finally, we demonstrate that such diverse ion distributions, along with specific electric-field-strength-dependent ion properties, lead to (1) electroosmotic (EOS) transport in nanochannels grafted with PAA brushes screened with Cs⁺ ions to be always counterion dominated, (2) EOS transport in nanochannels grafted with PAA brushes screened with Ca²⁺ and Y³⁺ ions to be always coion-dominated, and (3) EOS transport in nanochannels grafted with PAA brushes screened with Na⁺ and Ba²⁺ ions to be coion dominated for smaller electric fields and counterion dominated for larger electric fields.

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

The configuration of densely grafted charged polyelectrolyte (PE) brushes is strongly dictated by the properties and behavior of the counterions that screen the PE brush charges and the solvent molecules (typically water) that solvate the brush molecules and these screening counterions. Only recently, efforts have been made to study the PE brushes atomistically, thereby shedding light on the properties of brush-supported ions and water molecules. However, even for such efforts, there are limitations associated with using a generic definition to estimate certain properties of water and ions inside the brush layer. For example, water-water hydrogen bonds (HBs) will behave differently for locations outside and inside the brush layer, given the fact that the densely closely grafted PE brush molecules create a soft nanoconfinement where the water connectivity becomes highly disrupted: therefore, using the same definition to quantify the HBs inside and outside the brush layer will be unwise. In this paper, we address this limitation by employing an unsupervised machine learning (ML) approach to predict the water-water hydrogen bonding inside a cationic PE brush layer modeled using all-atom molecular dynamics (MD) simulations. The ML method, which relies on a clustering approach and uses the equilibrium coordinates of the water molecules (obtained from the all-atom MD simulations) as the input, is capable of identifying the structural modification of water-water HBs (revealed through appropriate clustering of the data) inside the PE brush layer induced soft nanoconfinement. Such capabilities would not have been possible by using a generic definition of the HBs. Our calculations lead to four key findings: (1) the clusters formed inside and outside the brush layer are structurally similar; (2) the margin of the cluster is shorter inside the PE brush layer confirming the possible disruption of the HBs inside the PE brush layer; (3) the average "hydrogen-acceptor-oxygen-donor-oxygen" angle that defines the HB is reduced for the HBs formed inside the brush layer; (4) the use of the generic definition (definition usable for characterizing the HBs in brush-free bulk) leads to an overprediction of the number of HBs formed inside the PE brush layer.

Effect of water nanoconfinement on the dynamic properties of paramagnetic colloidal complexes

The anomalous behavior of confined water at the nanoscale has remarkable implications in a number of nanotechnological applications. In this work, we analyze the effect of water self-diffusion on the dynamic properties of a solvated gadolinium-based paramagnetic complex, typically used for contrast enhancement in magnetic resonance imaging. In particular, we examine the effect of silica-based nanostructures on water behavior in the proximity of the paramagnetic complex via atomistic simulations, and interpret the resulting tumbling dynamics in the light of the local solvent modification based on the Lipari-Szabo formalism and of the fractional Stokes-Einstein relation. It is found that the local water confinement induces an increased "stiffness" on the outer sphere of the paramagnetic complex, which eventually reduces its tumbling properties. These model predictions are found to explain well the relaxivity enhancement observed experimentally by confining paramagnetic complexes into porous nanoconstructs, and thus offer mechanistic guidelines to design improved contrast agents for imaging applications.

Study on the Mechanism of the Reversible Color Change of Polyacrylic Acid Modified Gold Nanoparticles Responding to pH

In view of various explanations regarding the pH response of the nanocomposite of gold nanoparticles (AuNPs) modified with polyacrylic acid (PAA) molecules in reported literature, in this work, AuNPs with a size of 20 nm saturatedly loaded with PAA molecules (AuNPs-PAAs) were used to investigate the following aspects of this issue. We investigated the effects of pH on the stability of AuNPs-PAAs in the presence of salt, CTAB, poly (sodium styrenesulfonate) (PSS), ethanol, and free PAA, respectively. Common techniques were undertaken to evaluate the stability, including UV-Vis spectroscopy, Zeta potential analysis, and TEM. The results show that AuNPs-PAAs could respond to pH variations, having a reversible aggregation-to-disaggregation, accompanying their Zeta potential change. The proposed corresponding mechanism was that this reversible change was attributes to the net charge variation of AuNPs-PAAs induced by a reversible protonation-to-deprotonation of PAA rather than the conformational change. It was found that salt, CTAB, PSS, and free PAA could strengthen the dispersity of AuNPs-PAAs, even though their absolute Zeta potential values were decreased to small values or dropped to nearly zero. This abnormal phenomenon was explained by solvation. It was also found that AuNPs-PAAs have an opposite pH response in aqueous and ethanol solutions, justifying the solvation effect. All these results revealed the conformational stability of PAAs immobilized on AuNPs. The methods and the findings of this investigation give some new insights to understand the pH-response of AuNPs-PAAs composites and the design of AuNPs-PAAs-based functional sensors.

Overscreening, Co-Ion-Dominated Electroosmosis, and Electric Field Strength Mediated Flow Reversal in Polyelectrolyte Brush Functionalized Nanochannels

Controlling the direction and strength of nanofluidic electrohydrodyanmic transport in the presence of an externally applied electric field is extremely important in a number of nanotechnological applications. Here, we employ all-atom molecular dynamics simulations to discover the possibility of changing the direction of electroosmotic (EOS) liquid flows by merely changing the electric field strength in a nanochannel functionalized with polyelectrolyte (PE) brushes. In exploring this, we have uncovered three facets of nanoconfined PE brush behavior and resulting EOS transport. First, we identify the onset of an overscreening effect: such overscreening refers to the presence of more counterions (Na⁺) within the brush layer than needed to neutralize the negative brush charges. Accordingly, as a consequence of the overscreening, in the bulk liquid outside the brush layer, there is a greater number of co-ions (Cl^-) than counterions in the presence of an added salt (NaCl). Second, this specific ion distribution ensures that the overall EOS flow is along the direction of motion of the co-ions. Such co-ion-dictated EOS transport directly contradicts the notion that EOS flow is always dictated by the motion of the counterions. Finally, for large-enough electric fields, the brush height reduces significantly, causing some of the excess overscreening-inducing counterions to squeeze out of the PE brush layer into the brush-free bulk. As a result, the overscreening effect disappears and the number of co-ions and counterions outside the PE brush layer become similar. Despite that there is an EOS transport, this EOS transport, unlike the standard EOS transport that occurs due to the imbalance of the co-ions and counterions, occurs since a larger residence time of the water molecules in the first solvation shell of the counterions (Na⁺) ensures a water transport in the direction of motion of the counterions. The net effect is the reversal of the direction of the EOS transport by merely changing the strength of the electric field.