Calculating Pairwise Similarity of Polymer Ensembles via Earth Mover's Distance.

Synthetic polymers, in contrast to small molecules and deterministic biomacromolecules, are typically ensembles composed of polymer chains with varying numbers, lengths, sequences, chemistry, and topologies. While numerous approaches exist for measuring pairwise similarity among small molecules and sequence-defined biomacromolecules, accurately determining the pairwise similarity between two polymer ensembles remains challenging. This work proposes the earth mover's distance (EMD) metric to calculate the pairwise similarity score between two polymer ensembles. EMD offers a greater resolution of chemical differences between polymer ensembles than the averaging method and provides a quantitative numeric value representing the pairwise similarity between polymer ensembles in alignment with chemical intuition. The EMD approach for assessing polymer similarity enhances the development of accurate chemical search algorithms within polymer databases and can improve machine learning techniques for polymer design, optimization, and property prediction.

Community Resource for Innovation in Polymer Technology (CRIPT): A Scalable Polymer Material Data Structure.

The Community Resource for Innovation in Polymer Technology (CRIPT) data model is designed to address the high complexity in defining a polymer structure and the intricacies involved with characterizing material properties.

Extending BigSMILES to non-covalent bonds in supramolecular polymer assemblies.

As a machine-recognizable representation of polymer connectivity, BigSMILES line notation extends SMILES from deterministic to stochastic structures. The same framework that allows BigSMILES to accommodate stochastic covalent connectivity can be extended to non-covalent bonds, enhancing its value for polymers, supramolecular materials, and colloidal chemistry. Non-covalent bonds are captured through the inclusion of annotations to pseudo atoms serving as complementary binding pairs, minimal key/value pairs to elaborate other relevant attributes, and indexes to specify the pairing among potential donors and acceptors or bond delocalization. Incorporating these annotations into BigSMILES line notation enables the representation of four common classes of non-covalent bonds in polymer science: electrostatic interactions, hydrogen bonding, metal-ligand complexation, and π-π stacking. The principal advantage of non-covalent BigSMILES is the ability to accommodate a broad variety of non-covalent chemistry with a simple user-orientated, semi-flexible annotation formalism. This goal is achieved by encoding a universal but non-exhaustive representation of non-covalent or stochastic bonding patterns through syntax for (de)protonated and delocalized state of bonding as well as nested bonds for correlated bonding and multi-component mixture. By allowing user-defined descriptors in the annotation expression, further applications in data-driven research can be envisioned to represent chemical structures in many other fields, including polymer nanocomposite and surface chemistry.

From well-entangled to partially-entangled wormlike micelles.

We combine mechanical rheometry, diffusing wave spectroscopy (DWS), and small angle neutron scattering (SANS) with a simulation model, the "pointer algorithm", to obtain characteristic lengths and time constants for wormlike micelle (WLM) solutions over a range of salt concentrations encompassing the transition from unentangled to entangled solutions. The solutions contain sodium lauryl ethylene glycol sulfate (SLE1S), cocamidopropyl betaine (CAPB), and NaCl. The pointer algorithm is extended to include relaxation of unentangled micelles, allowing micelle parameters to be extracted from the rheology of partially entangled solutions. DWS provides the data at high frequency needed to determine micelle persistence length accurately. From pointer algorithm fits to rheology, we observe a salt-induced rapid change in micellar length as the solution enters the well-entangled regime and a weaker growth with surfactant concentration consistent with mean-field theory. At a lower surfactant concentration, micelle length and persistence length from SANS are roughly consistent with values from rheology once the lower surfactant concentration used in SANS is accounted for. This is, to our knowledge, the first time that quantitative comparisons of structural features including micelle length are made between rheology and SANS. Finally, scaling laws for micelle diffusion and recombination times indicate that micelle kinetics are reaction controlled leading to mean-field recombination with surrounding micelles over the entire range of concentration of interest except at very low and very high surfactant concentrations where either short micelles or branched micelle clusters are dominant.

Multiscale Modeling of the Effects of Salt and Perfume Raw Materials on the Rheological Properties of Commercial Threadlike Micellar Solutions.

We link micellar structures to their rheological properties for two surfactant body-wash formulations at various concentrations of salts and perfume raw materials (PRMs) using molecular simulations and micellar-scale modeling, as well as traditional surfactant packing arguments. The two body washes, namely, BW-1EO and BW-3EO, are composed of sodium lauryl ethylene glycol ether sulfate (SLEnS, where n is the average number of ethylene glycol repeat units), cocamidopropyl betaine (CAPB), ACCORD (which is a mixture of six PRMs), and NaCl salt. BW-3EO is an SLE3S-based body wash, whereas BW-1EO is an SLE1S-based body wash. Additional PRMs are also added into the body washes. The effects of temperature, salt, and added PRMs on micellar lengths, breakage times, end-cap free energies, and other properties are obtained from fits of the rheological data to predictions of the "Pointer Algorithm" [ Zou , W. ; Larson , R.G. J. Rheol. 2014 , 58 , 1 - 41 ], which is a simulation method based on the Cates model of micellar dynamics. Changes in these micellar properties are interpreted using the Israelachvili surfactant packing argument. From coarse-grained molecular simulations, we infer how salt modifies the micellar properties by changing the packing between the surfactant head groups, with the micellar radius remaining nearly constant. PRMs do so by partitioning to different locations within the micelles according to their octanol/water partition coefficient POW and chemical structures, adjusting the packing of the head and/or tail groups, and by changing the micelle radius, in the case of a large hydrophobic PRM. We find that relatively hydrophilic PRMs with log POW < 2 partition primarily to the head group region and shrink micellar length, decreasing viscosity substantially, whereas more hydrophobic PRMs, with log POW between 2 and 4, mix with the hydrophobic surfactant tails within the micellar core and slightly enhance the viscosity and micelle length, which is consistent with the packing argument. Large and very hydrophobic PRMs, with log POW > 4, are isolated deep inside the micelle, separating from the tails and swelling the radius of the micelle, leading to shorter micelles and much lower viscosities, leading eventually to swollen-droplet micelles.

A hybrid Brownian dynamics/constitutive model for yielding, aging, and rejuvenation in deforming polymeric glasses.

We present a hybrid model for polymeric glasses under deformation that combines a minimal model of segmental dynamics with a beads-and-springs model of a polymer, solved by Brownian dynamics (BD) simulations, whose relaxation is coupled to the segmental dynamics through the drag coefficient of the beads. This coarse-grained model allows simulations that are much faster than molecular dynamics and successfully capture the entire range of mechanical response including yielding, plastic flow, strain-hardening, and incomplete strain recovery. The beads-and-springs model improves upon the dumbbell model for glassy polymers proposed by Fielding et al. (Phys. Rev. Lett., 2012, 108, 048301) by capturing the small elastic recoil seen experimentally without the use of ad hoc adjustments of parameters required in the model of Fielding et al. With appropriate choice of parameters, predictions of creep, recovery, and segmental relaxation are found to be in good agreement with poly(methylmethacrylate) (PMMA) data of Lee et al. (Science, 2009, 323, 231-234). Our model shows dramatic differences in behavior of the segmental relaxation time between extensional creep and steady extension, and between extension and shear. The non-monotonic response of the segmental relaxation time to extensional creep and the small elastic recovery after removal of stress are shown to arise from sub-chains that are trapped between folds, and that become highly oriented and stretched at strains of order unity, connecting the behavior of glassy polymers under creep to that of dilute polymer solutions under fast extensional flows. We are also able to predict the effects of polymer pre-orientation in the parallel or orthogonal direction on the subsequent response to extensional deformation.