Polymer composite
materials require softening to reduce their glass
transition temperature and improve processability. To this end, plasticizers
(PLs), which are small organic molecules, are added to the polymer
matrix. The miscibility of these PLs has a large impact on their effectiveness
and, therefore, their interactions with the polymer matrix must be
carefully considered. Many PL characteristics, including their size,
topology, and flexibility, can impact their miscibility and, because
of the exponentially large number of PLs, the current trial-and-error
approach is very ineffective. In this work, we show that using coarse-grained
molecular simulations of a small dataset of 48 PLs, it is possible
to identify topological and thermodynamic descriptors that are proxy
for their miscibility. Using ad-hoc molecular dynamics
simulation setups that are relatively computationally inexpensive,
we establish correlations between the PLs’ topology, internal
flexibility, thermodynamics of aggregation, and degree of miscibility,
and use these descriptors to classify the molecules as miscible or
immiscible. With all available data, we also construct a decision
tree model, which achieves a F1 score of 0.86 ± 0.01 with repeated,
stratified 5-fold cross-validation, indicating that this machine learning
method can be a promising route to fully automate the screening. By
evaluating the individual performance of the descriptors, we show
this procedure enables a 10-fold reduction of the test space and provides
the basis for the development of workflows that can efficiently screen
PLs with a variety of topological features. The approach is used here
to screen for apolar PLs in polyisoprene melts, but similar proxies
would be valid for other polyolefins, while, in cases where polar
interactions drive the miscibility, other descriptors are likely to
be needed.