Search

Colloidal crystals are often prepared by evaporation from solution, and there is considerable interest to link the processing conditions to the crystal morphology and quality. Here, we study the evaporation-induced assembly of colloidal crystals using massive-scale nonequilibrium molecular dynamics simulations. We apply a recently developed machine-learning technique to characterize the assembling crystal structures with unprecedented microscopic detail. In agreement with previous experiments and simulations, faster evaporation rates lead to earlier onset of crystallization and more disordered surface structures. Surprisingly, we find that collective rearrangements of the bulk crystal during later stages of drying reduce the influence of the initial surface structure, and the final morphology is essentially independent of the evaporation rate. Our structural analysis reveals that the crystallization process is well-described by two time scales, the film drying time and the crystal growth time, with the latter having an unexpected dependence on the evaporation rate due to equilibrium thermodynamic effects at high colloid concentrations. These two time scales may be leveraged to control the relative influence of equilibrium and nonequilibrium growth mechanisms, suggesting a route to rapidly process colloidal crystals while also removing defects. Our analysis additionally reveals that solvent-mediated interactions play a critical role in the crystallization kinetics and that commonly used implicit-solvent models do not faithfully resolve nonequilibrium processes such as drying. [ABSTRACT FROM AUTHOR]

Molecular dynamic simulations are used to investigate the structural effects of treating a glassy polymer thin film with solvents of varying quality and subsequently evaporating the solvent. Both a monodisperse film and a polydisperse film are studied for poor to good solvent conditions, including the limit in which the polymer film is fully dissolved. In agreement with previous studies, the dissolved polymer-solvent mixtures form a polymer-rich skin on top of the forming film during evaporation. In the case of the polydisperse films, a segregation of the lower molecular weight polymer to the film interface is observed. We provide a detailed, systematic analysis of the interface structure and properties during and after evaporation. We find that for non-dissolved films, the surface width of the film after solvent evaporation is enhanced compared to the case without solvent. Our results show that due to the kinetic arrest of the surface structure, the increased surface width is preserved after solvent evaporation for both mono- and polydisperse films. We conclude that it is important to take poor solvent effects into account for the surface morphology of already formed thin glassy films, an effect which is often neglected. [ABSTRACT FROM AUTHOR]

We perform hybrid molecular dynamics simulations to study the flow behavior of rigid colloids dispersed in a dilute polymer solution. The underlying Newtonian solvent and the ensuing hydrodynamic interactions are incorporated through multiparticle collision dynamics, while the constituent polymers are modeled as bead-spring chains, maintaining a description consistent with the colloidal nature of our system. We study the cross-stream migration of the solute particles in slit-like channels for various polymer lengths and colloid sizes and find a distinct focusing onto the channel center under specific solvent and flow conditions. To better understand this phenomenon, we systematically measure the effective forces exerted on the colloids. We find that the migration originates from a competition between viscoelastic forces from the polymer solution and hydrodynamically induced inertial forces. Our simulations reveal a significantly stronger fluctuation of the lateral colloid position than expected from thermal motion alone, which originates from the complex interplay between the colloid and polymer chains. [ABSTRACT FROM AUTHOR]

We study the structure and dynamics of semidilute solutions of semiflexible polymers at rest and under shear using hybrid molecular dynamics simulations that take hydrodynamic interactions into account. We show that the polymer center-of-mass diffusion coefficient significantly decreases with increasing chain stiffness at fixed monomer density. The zero-shear viscosity shows a corresponding increase due to the intermolecular interactions of stiffer chains. We apply steady shear flow to the polymer solutions and show that at high shear rates the flow properties become almost independent of polymer stiffness. We characterize the polymer conformations under shear and find that in this regime polymers are elongated and aligned along the flow direction, but semiflexible polymers surprisingly form an increased number of hairpin-like conformations due to increased tumbling. [ABSTRACT FROM AUTHOR]

Stratification in binary colloidal mixtures was investigated using implicit-solvent molecular dynamics simulations. For large particle size ratios and film Péclet numbers greater than unity, smaller colloids migrated to the top of the film, while big colloids were pushed to the bottom, creating an "inverted" stratification. This peculiar behavior was observed in recent simulations and experiments conducted by Fortini et al. [ Phys. Rev. Lett. 2016 , 116 , 118301 ]. To rationalize this behavior, particle size ratios and drying rates spanning qualitatively different Péclet number regimes were systematically studied, and the dynamics of the inverted stratification were quantified in detail. The stratified layer of small colloids was found to grow faster and to larger thicknesses for larger size ratios. Interestingly, inverted stratification was observed even at moderate drying rates where the film Péclet numbers were comparable to unity, but the thickness of the stratified layer decreased. A model based on dynamical density functional theory is proposed to explain the observed phenomena. [ABSTRACT FROM AUTHOR]

Simulation of the Rouse model in flow underlies a great variety of numerical investigations of polymer dynamics, in both entangled melts and solutions and in dilute solution. Typically a simple explicit stochastic Euler method is used to evolve the Rouse model. Here we compare this approach to an operator splitting method, which splits the evolution operator into stochastic linear and deterministic nonlinear parts and takes advantage of an analytical solution for the linear Rouse model in terms of the noise history. We show that this splitting method has second-order weak convergence, whereas the Euler method has only first-order weak convergence. Furthermore, the splitting method is unconditionally stable, in contrast to the limited stability range of the Euler method. Similar splitting methods are applicable to a broad class of problems in stochastic dynamics in which noise competes with ordering and flow to determine steady-state order parameter structures. [ABSTRACT FROM AUTHOR]