Evaporation- and surface-induced morphology of symmetric diblock copolymer thin films: a multibody dissipative particle dynamics study

Many-body dissipative particle dynamics simulations of micellization of sodium alkyl sulfates

We present a study of micelle formation in alkyl sulfate surfactants using the simulation method of many-body dissipative particle dynamics (MDPD). We parametrise our model by tuning the intermolecular interactions in order to reproduce experimental values for the chemical potential and density at room temperature. Using this approach, we find that our model shows good agreement with experimental values for the critical micelle concentration (CMC). Furthermore, we show that our model can accurately predict CMC trends, which result from varying properties such as surfactant tail length and the salt concentration. We apply our model to investigate the effect of aggregation number on various micellar properties, such as the shape of individual micelles and the fraction of bound counterions. We show that micelles become aspherical at large aggregation numbers, in line with experimental predictions, and that longer tail surfactants are generally more spherical at all aggregation numbers compared to those which are shorter. We find excellent agreement between our simulations and experimental values for the degree of counterion binding, a factor that is crucial to accurately studying micellar shape, but one that is typically overlooked in the existing literature.

A many-body dissipative particle dynamics parametrisation scheme to study behaviour at air-water interfaces

In this article, we present a general parametrisation scheme for many-body dissipative particle dynamics (MDPD). The scheme is based on matching model components to experimental surface tensions and chemical potentials. This allows us to obtain the correct surface and mixing behaviours of complex, multicomponent systems. The methodology is tested by modelling the behaviour of nonionic polyoxyethylene alkyl ether surfactants at an air/water interface. In particular, the influence of the number of ethylene oxide units in the surfactant head group is investigated. We find good agreement with many experimentally obtained parameters, such as minimum surface area per molecule; and a decrease in the surface tension with increasing surfactant surface density. Moreover, we observe an orientational transition, from surfactants lying directly on the water surface at low surface coverage, to surfactants lying parallel or tilted with respect to the surface normal at high surface coverage. The parametrisation scheme is also extended to cover the zwitterionic surfactant lauryldimethylamine oxide (LDAO), where we provide good predictions for the surface tension at maximum surface coverage. Here, if we exceed this coverage, we are able to demonstrate the spontaneous production of micelles from the surface surfactant layer.

Molecular View of the Distortion and Pinning Force of a Receding Contact Line: Impact of the Nanocavity Geometry

We present a molecular view using many-body dissipative particle dynamics simulations to unravel the pinning phenomenon of a liquid film receding over a solid substrate with a nanocavity. We find that the pinning force and distortion of the pinned contact line vary across different nanocavity shapes. We show that the mechanism of a caterpillar motion, which has previously been proposed for advancing precursor films, persists in a partially pinned receding contact line. Our results also demonstrate a localized clamping effect, which is originated from the variation of the dynamic contact angle along the pinned contact line. The simulation results suggest that the clamping effect can be controlled by the geometry of the nanocavity and hydrophilicity of the underlying substrate.

Highly Oil-Repellent Thermoplastic Boundaries via Surface Delivery of CF3 Groups by Molecular Bottlebrush Additives

We fabricated thermoplastic surfaces possessing extremely limited water and oil wettability without employment of long-chain perfluoroalkyl (LCPFA) substances. Namely, by taking advantage of the structure and behavior of original oleophobic perfluoropolyether (PFPE) methacrylate (PFM) molecular bottlebrush (MBB) additive we obtained polymeric surfaces with oil contact angles well above 80° and surface energy on the level of 10 mN/m. Those angles and surface energies are the highest and the lowest respective values reported to date for any bulk solid flat organic surface not containing LCPFA. We show experimentally and computationally that this remarkable oil repellency is attributed to migration of small quantities of the oleophobic MBB additives to the surface of the thermoplastics. Severe mismatch in the affinity between the densely grafted long side chains of MBB and a host matrix promotes stretching and densification of mobile side chains delivering the lowest surface energy functionalities (CF₃) to the materials' boundary. Our studies demonstrate that PFM can be utilized as an effective low surface energy additive to conventional thermoplastic polymers, such as poly(methyl methacrylate) and Nylon-6. We show that films containing PFM achieve the level of oil repellency significantly higher than that of polytetrafluoroethylene (PTFE), a fully perfluorinated thermoplastic. The surface energy of the films is also significantly lower than that of PTFE, even at low concentrations of PFM additives.

Hydrodynamic Interactions and Entanglements of Polymer Solutions in Many-Body Dissipative Particle Dynamics

Using many-body dissipative particle dynamics (MDPD), polymer solutions with concentrations spanning dilute and semidilute regimes are modeled. The parameterization of MDPD interactions for systems with liquid⁻vapor coexistence is established by mapping to the mean-field Flory⁻Huggins theory. The characterization of static and dynamic properties of polymer chains is focused on the effects of hydrodynamic interactions and entanglements. The coil⁻globule transition of polymer chains in dilute solutions is probed by varying solvent quality and measuring the radius of gyration and end-to-end distance. Both static and dynamic scaling relations for polymer chains in poor, theta, and good solvents are in good agreement with the Zimm theory with hydrodynamic interactions considered. Semidilute solutions with polymer volume fractions up to 0.7 exhibit the screening of excluded volume interactions and subsequent shrinking of polymer coils. Furthermore, entanglements become dominant in the semidilute solutions, which inhibit diffusion and relaxation of chains. Quantitative analysis of topology violation confirms that entanglements are correctly captured in the MDPD simulations.