Grafted Dipolar Chains: Dipoles and Restricted Freedom Lead to Unexpected Hairpins

Learning glass transition temperatures via dimensionality reduction with data from computer simulations: Polymers as the pilot case

Machine learning methods provide an advanced means for understanding inherent patterns within large and complex datasets. Here, we employ the principal component analysis (PCA) and the diffusion map (DM) techniques to evaluate the glass transition temperature (Tg) from low-dimensional representations of all-atom molecular dynamic simulations of polylactide (PLA) and poly(3-hydroxybutyrate) (PHB). Four molecular descriptors were considered: radial distribution functions (RDFs), mean square displacements (MSDs), relative square displacements (RSDs), and dihedral angles (DAs). By applying Gaussian Mixture Models (GMMs) to analyze the PCA and DM projections and by quantifying their log-likelihoods as a density-based metric, a distinct separation into two populations corresponding to melt and glass states was revealed. This separation enabled the Tg evaluation from a cooling-induced sharp increase in the overlap between log-likelihood distributions at different temperatures. Tg values derived from the RDF and MSD descriptors using DM closely matched the standard computer simulation-based dilatometric and dynamic Tg values for both PLA and PHB models. This was not the case for PCA. The DM-transformed DA and RSD data resulted in Tg values in agreement with experimental ones. Overall, the fusion of atomistic simulations and DMs complemented with the GMMs presents a promising framework for computing Tg and studying the glass transition in a unified way across various molecular descriptors for glass-forming materials.

Effects of Amino Acid Side-Chain Length and Chemical Structure on Anionic Polyglutamic and Polyaspartic Acid Cellulose-Based Polyelectrolyte Brushes

We used atomistic molecular dynamics (MD) simulations to study polyelectrolyte brushes based on anionic α,L-glutamic acid and α,L-aspartic acid grafted on cellulose in the presence of divalent CaCl₂ salt at different concentrations. The motivation is to search for ways to control properties such as sorption capacity and the structural response of the brush to multivalent salts. For this detailed understanding of the role of side-chain length, the chemical structure and their interplay are required. It was found that in the case of glutamic acid oligomers, the longer side chains facilitate attractive interactions with the cellulose surface, which forces the grafted chains to lie down on the surface. The additional methylene group in the side chain enables side-chain rotation, enhancing this effect. On the other hand, the shorter and more restricted side chains of aspartic acid oligomers prevent attractive interactions to a large degree and push the grafted chains away from the surface. The difference in side-chain length also leads to differences in other properties of the brush in divalent salt solutions. At a low grafting density, the longer side chains of glutamic acid allow the adsorbed cations to be spatially distributed inside the brush resulting in a charge inversion. With an increase in grafting density, the difference in the total charge of the aspartic and glutamine brushes disappears, but new structural features appear. The longer sides allow for ion bridging between the grafted chains and the cellulose surface without a significant change in main-chain conformation. This leads to the brush structure being less sensitive to changes in salt concentration.

Conformational transitions and helical structures of a dipolar chain in external electric fields

The conformational behavior of a single dipolar chain in a uniform electric field is investigated by molecular dynamics simulations. The dipolar chain is modeled as a backbone bead-on-spring chain of equally charged beads, each connected by a rigid spring with an oppositely charged side bead that can freely rotate around the backbone bead. In the strong coupling regime, when the dipolar chain is in the globular state due to a strong electrostatic correlational attraction, the application of an electric field causes the chain swelling and elongation along the field direction. In the weak coupling regime, a qualitatively new regime is found when the swollen dipolar chain shrinks along the field direction adopting flattened conformations due to the field-induced anisotropy of the chain rigidity and the head-to-tail attraction of the dipoles orienting along the field lines. A novel helical conformation is detected for low-polar media and strong electric fields. An increasing rigidity of the backbone chain leads to some stabilization of the helical conformation and the formation of double and triple helices as well as flat spread springs. Fine tuning of the interplay between dipolar and volume interactions by external electric fields induces re-orientation of rod-like dipolar chains in dilute solutions. The obtained results can provide new ways to control dipolar polymer conformations and design materials with responsive properties.

Branched versus linear lactide chains for cellulose nanoparticle modification: an atomistic molecular dynamics study

We studied the structure of brushes consisting of branched oligolactide (OLA) chains grafted onto the surface of cellulose nanoparticles (CNPs) in polylactide (PLA) and compared the outcomes to the case of grafting linear OLA chains using atomistic molecular dynamics simulations. The systems were considered in a melt state. The branched model OLA chains comprised one branching point and three branches, while the linear OLA chains examined had a molecular weight similar to the branched chains. It was shown that free branches of the branched OLA chains tend to fold back toward the CNPs due to dipole-dipole interactions within the grafted layer, in contrast to the well-established behavior of the grafted uncharged branched chains. This result, however, is in qualitative agreement with the conformational behavior known for linear OLA chains. At the same time, no significant difference in the effectiveness of covering the filler surface with grafted branched or linear OLA chains was found. In terms of the expelling ability of the grafted chains and the interaction between PLA and CNP or OLA, the linear chains were broadly similar (sparse grafting) or better (intermediate or dense grafting) compared to the branched ones. Thus, the grafted lactide chains with a linear architecture, rather than their branched counterpart, may be preferable for the covalent modification of cellulose nanoparticles.

The Structure of Dipolar Polymer Brushes and Their Interaction in the Melt. Impact of Chain Stiffness

By using the numerical lattice Scheutjens-Fleer self-consistent field (SF-SCF) method we have studied the effect of the restricted flexibility of grafted chains on the structure and mutual interaction of two opposing planar conventional and A-type dipolar brushes. Brushes are immersed in the solvent consisting of chains similar to the grafted ones. The increase of the chain rigidity enhances the segregation of grafted chains in a A-type dipolar brush into two populations: backfolded chains with terminal monomers near the grafting surface and chains with the ends at the brush periphery. The fraction of backfolded chains grows by an increase of the Kuhn segment length. It is shown that two opposite A-type dipolar brushes from semi-rigid chains are attracted to each other at short distances. The attraction becomes more pronounced and begins at larger distances for more rigid chains with the same brush characteristics: polymerization degree, grafting density, and dipole moments of monomer units. This attraction is connected with the dipole-dipole interactions between chains of oncoming brushes with oppositely directed dipoles penetrating deeply into each other upon contact. This effect of the chain rigidity is opposite to that for conventional brushes without dipoles in the chains. For such brushes, an increase in the chain rigidity leads to the enhanced repulsion between them.