Design of New Dispersants Using Machine Learning and Visual Analytics.

Artificial intelligence (AI) is an emerging technology that is revolutionizing the discovery of new materials. One key application of AI is virtual screening of chemical libraries, which enables the accelerated discovery of materials with desired properties. In this study, we developed computational models to predict the dispersancy efficiency of oil and lubricant additives, a critical property in their design that can be estimated through a quantity named blotter spot. We propose a comprehensive approach that combines machine learning techniques with visual analytics strategies in an interactive tool that supports domain experts' decision-making. We evaluated the proposed models quantitatively and illustrated their benefits through a case study. Specifically, we analyzed a series of virtual polyisobutylene succinimide (PIBSI) molecules derived from a known reference substrate. Our best-performing probabilistic model was Bayesian Additive Regression Trees (BART), which achieved a mean absolute error of 5.50±0.34 and a root mean square error of 7.56±0.47, as estimated through 5-fold cross-validation. To facilitate future research, we have made the dataset, including the potential dispersants used for modeling, publicly available. Our approach can help accelerate the discovery of new oil and lubricant additives, and our interactive tool can aid domain experts in making informed decisions based on blotter spot and other key properties.

Molecular Simulations of Surfactant Adsorption on Iron Oxide from Hydrocarbon Solvents.

The performance of organic friction modifiers (OFMs) depends on their ability to adsorb onto surfaces and form protective monolayers. Understanding the relationship between OFM concentration in the base oil and the resulting surface coverage is important for improving lubricant formulations. Here, we use molecular dynamics (MD) simulations to study the adsorption of three OFMs─stearic acid (SA), glycerol monoostearate (GMS), and glycerol monooleate (GMO)─onto a hematite surface from two hydrocarbon solvents─n-hexadecane and poly(α-olefin) (PAO). We calculate the potential of mean force of the adsorption process using the adaptive biasing force algorithm, and the adsorption strength increases in the order SA < GMS < GMO. We estimate the minimum area occupied by OFM molecules on the surface using annealing MD simulations and obtained a similar hard-disk area for GMS and GMO but a lower value for SA. Using the MD results, we determine the adsorption isotherms using the molecular thermodynamic theory (MTT), which agree well with one previous experimental data set for SA on hematite. For two other experimental data sets for SA, lateral interactions between surfactant molecules need to be accounted for within the MTT framework. SA forms monolayers with lower surface coverage than GMO and GMS at low concentrations but also has the highest plateau coverage. We validate the adsorption energies from the MD simulations using high-frequency reciprocating rig friction experiments with different concentrations of the OFMs in PAO. For OFMs with saturated tailgroups (SA and GMS), we obtain good agreement between the simulations and the experiments. The results deviate for OFMs containing Z-unsaturated tailgroups (GMO) due to the additional steric hindrance, which is not accounted for in the current simulation framework. This study demonstrates that MD simulations, alongside MTT, are an accurate and efficient tool to predict adsorption isotherms at solid-liquid interfaces.

Toward Automated Tools for Characterization of Molecular Porosity.

The emerging advanced porous materials, e.g. extended framework materials and porous molecular materials, offer an unprecedented level of control of their structure and function. The enormous possibilities for tuning these materials by changing their building blocks mean that, in principle, optimally performing materials for a variety of applications can be systematically designed. However, the process of finding a set of optimal structures for a given application requires computational high-throughput tools to analyze and sieve through many candidate materials. In particular, in the case of porous molecular materials, the analysis and selection of a molecule is one of the key aspects as the structure of the molecule determines the structure of the resulting material, and very often the porosity of the molecule significantly contributes to the porous properties of the resulting material. In this work, we introduce definitions and algorithms to characterize porosity at the molecular level, along with a software implementation of these algorithms. We demonstrate applications of the software tool in the discovery and characterization of porous molecules among ca. 94 million molecules currently enlisted in the PubChem database.

Cluster Glasses of Semiflexible Ring Polymers.

We present computer simulations of concentrated solutions of unknotted nonconcatenated semiflexible ring polymers. Unlike in their flexible counterparts, shrinking involves a strong energetic penalty, favoring interpenetration and clustering of the rings. We investigate the slow dynamics of the centers-of-mass of the rings in the amorphous cluster phase, consisting of disordered columns of oblate rings penetrated by bundles of prolate ones. Scattering functions reveal a striking decoupling of self- and collective motions. Correlations between centers-of-mass exhibit slow relaxation, as expected for an incipient glass transition, indicating the dynamic arrest of the cluster positions. However, self-correlations decay at much shorter time scales. This feature is a manifestation of the fast, continuous exchange and diffusion of the individual rings over the matrix of clusters. Our results reveal a novel scenario of glass formation in a simple monodisperse system, characterized by self-collective decoupling, soft caging, and mild dynamic heterogeneity.

Cluster glasses of ultrasoft particles.

We present molecular dynamics (MD) simulations results for dense fluids of ultrasoft, fully penetrable particles. These are a binary mixture and a polydisperse system of particles interacting via the generalized exponential model, which is known to yield cluster crystal phases for the corresponding monodisperse systems. Because of the dispersity in the particle size, the systems investigated in this work do not crystallize and form disordered cluster phases. The clustering transition appears as a smooth crossover to a regime in which particles are mostly located in clusters, isolated particles being infrequent. The analysis of the internal cluster structure reveals microsegregation of the big and small particles, with a strong homo-coordination in the binary mixture. Upon further lowering the temperature below the clustering transition, the motion of the clusters' centers-of-mass slows down dramatically, giving way to a cluster glass transition. In the cluster glass, the diffusivities remain finite and display an activated temperature dependence, indicating that relaxation in the cluster glass occurs via particle hopping in a nearly arrested matrix of clusters. Finally we discuss the influence of the microscopic dynamics on the transport properties by comparing the MD results with Monte Carlo simulations.

Static and dynamic contributions to anomalous chain dynamics in polymer blends.

By means of computer simulations, we investigate the relaxation of the Rouse modes in a simple bead-spring model for non-entangled polymer blends. Two different models are used for the fast component, namely fully flexible and semiflexible chains. The latter, which incorporate intramolecular barriers with bending and torsion terms, are semiflexible in the sense that static intrachain correlations are strongly non-gaussian at all length scales. The dynamic asymmetry in the blend is strongly enhanced with decreasing temperature, inducing confinement effects on the fast component. The dynamics of the Rouse modes show very different trends for the two models of the fast component. For the fully flexible case, the relaxation times exhibit a progressive deviation from Rouse scaling on increasing the dynamic asymmetry. This anomalous effect has a dynamic origin. It is not related to particular static features of the Rouse modes, which indeed are identical to those of the fully flexible homopolymer, and are not modified by the dynamic asymmetry in the blend. On the contrary, in the semiflexible case the relaxation times approximately exhibit the same scaling behaviour as the amplitudes of the modes. This suggests that the origin of the anomalous dynamic scaling for semiflexible chains confined in the blend is essentially of static nature. We discuss the implications of these observations for the applicability of theoretical approaches to chain dynamics in polymer blends.

From caging to Rouse dynamics in polymer melts with intramolecular barriers: a critical test of the mode coupling theory.

By means of computer simulations and solution of the equations of the mode coupling theory (MCT), we investigate the role of the intramolecular barriers on several dynamic aspects of nonentangled polymers. The investigated dynamic range extends from the caging regime characteristic of glass-formers to the relaxation of the chain Rouse modes. We review our recent work on this question, provide new results, and critically discuss the limitations of the theory. Solutions of the MCT for the structural relaxation reproduce qualitative trends of simulations for weak and moderate barriers. However, a progressive discrepancy is revealed as the limit of stiff chains is approached. This disagreement does not seem related with dynamic heterogeneities, which indeed are not enhanced by increasing barrier strength. It is not connected either with the breakdown of the convolution approximation for three-point static correlations, which retains its validity for stiff chains. These findings suggest the need of an improvement of the MCT equations for polymer melts. Concerning the relaxation of the chain degrees of freedom, MCT provides a microscopic basis for time scales from chain reorientation down to the caging regime. It rationalizes, from first principles, the observed deviations from the Rouse model on increasing the barrier strength. These include anomalous scaling of relaxation times, long-time plateaux, and nonmonotonous wavelength dependence of the mode correlators.

The role of intramolecular barriers on the glass transition of polymers: Computer simulations versus mode coupling theory.

We present computer simulations of a simple bead-spring model for polymer melts with intramolecular barriers. By systematically tuning the strength of the barriers, we investigate their role on the glass transition. Dynamic observables are analyzed within the framework of the mode coupling theory (MCT). Critical nonergodicity parameters, critical temperatures, and dynamic exponents are obtained from consistent fits of simulation data to MCT asymptotic laws. The so-obtained MCT lambda-exponent increases from standard values for fully flexible chains to values close to the upper limit for stiff chains. In analogy with systems exhibiting higher-order MCT transitions, we suggest that the observed large lambda-values arise form the interplay between two distinct mechanisms for dynamic arrest: general packing effects and polymer-specific intramolecular barriers. We compare simulation results with numerical solutions of the MCT equations for polymer systems, within the polymer reference interaction site model (PRISM) for static correlations. We verify that the approximations introduced by the PRISM are fulfilled by simulations, with the same quality for all the range of investigated barrier strength. The numerical solutions reproduce the qualitative trends of simulations for the dependence of the nonergodicity parameters and critical temperatures on the barrier strength. In particular, the increase in the barrier strength at fixed density increases the localization length and the critical temperature. However the qualitative agreement between theory and simulation breaks in the limit of stiff chains. We discuss the possible origin of this feature.

Dynamic arrest in polymer melts: competition between packing and intramolecular barriers.

We present molecular dynamics simulations of a simple model for polymer melts with intramolecular barriers. We investigate structural relaxation as a function of the barrier strength. Dynamic correlators can be consistently analyzed within the framework of the mode coupling theory of the glass transition. Control parameters are tuned in order to induce a competition between general packing effects and polymer-specific intramolecular barriers as mechanisms for dynamic arrest. This competition yields unusually large values of the so-called mode coupling theory exponent parameter and rationalizes qualitatively different observations for simple bead-spring and realistic polymers. The systematic study of the effect of intramolecular barriers presented here also establishes a fundamental difference between the nature of the glass transition in polymers and in simple glass formers.