Effect of stiffness on the phase behavior of cubic lattice chains

Gran canonica Monte Carlo (GCMC) simulazioni assistite da tecniche di riponderazione istogramma sono stati utilizzati per studiare l'effetto della flessibilità catena sul comportamento di soluzione fase di cubi catene reticolari corti con 4-32 segmenti. Ciò è stato fatto variando un parametro di rigidità gradualmente fino alla media calcolata end-to-end distanza avvicinato la lunghezza totale. Per entrambe le catene flessibili e rigide si è riscontrato che la temperatura critica, ottenuta tramite mista analisi dei campi di dimensioni finite, aumentata lunghezza della catena e la densità critica trasferisce a valori più bassi, in accordo con le osservazioni sperimentali. Il estrapolato lungh…

Directed Assembly of Soft Colloids through Rapid Solvent Exchange

We studied the directed assembly of soft nanoparticles through rapid micromixing of polymers in solution with a nonsolvent. Both experiments and computer simulations were performed to elucidate the underlying physics and to investigate the role of various process parameters. In particular, we discovered that no external stabilizing agents or charged end groups are required to keep the colloids separated from each other when water is used as the nonsolvent. Furthermore, the size of the nanoparticles can be reliably tuned through the mixing rate and the ratio between polymer solution and nonsolvent. Our results demonstrate that this mechanism is highly promising for the mass fabrication of un…

Phase Equilibria of Lattice Polymers from Histogram Reweighting Monte Carlo Simulations

Histogram-reweighting Monte Carlo simulations were used to obtain polymer / solvent phase diagrams for lattice homopolymers of chain lengths up to r=1000 monomers. The simulation technique was based on performing a series of grand canonical Monte Carlo calculations for a small number of state points and combining the results to obtain the phase behavior of a system over a range of temperatures and densities. Critical parameters were determined from mixed-field finite-size scaling concepts by matching the order parameter distribution near the critical point to the distribution for the three-dimensional Ising universality class. Calculations for the simple cubic lattice (coordination number z…

Rapid Production of Internally Structured Colloids by Flash Nanoprecipitation of Block Copolymer Blends.

Colloids with internally structured geometries have shown great promise in applications ranging from biosensors to optics to drug delivery, where the internal particle structure is paramount to performance. The growing demand for such nanomaterials necessitates the development of a scalable processing platform for their production. Flash nanoprecipitation (FNP), a rapid and inherently scalable colloid precipitation technology, is used to prepare internally structured colloids from blends of block copolymers and homopolymers. As revealed by a combination of experiments and simulations, colloids prepared from different molecular weight diblock copolymers adopt either an ordered lamellar morph…

Effect of Stiffness on the Micellization Behavior of Model H4T4 Surfactant Chains

The micellization behavior of a series of model surfactants, all with four head and tail groups (H4T4) but with different degrees of chain stiffness, was studied using grand canonical Monte Carlo simulations on a cubic lattice. The critical micelle concentration, micellar size, and thermodynamics of micellization were examined. In all cases investigated, the critical micelle concentration was found to increase with increasing temperature as observed for nonionic surfactants in apolar or slightly polar solvents. At a fixed reduced temperature and increasing chain stiffness, in agreement with previous observations, it was found that the critical micelle concentration decreased and the average…

On the Stability of Polymeric Nanoparticles Fabricated through Rapid Solvent Mixing.

We study the stability of polymeric nanoparticles fabricated through the rapid mixing of polymers in a good solvent with a poor solvent that is miscible with the good solvent. In previous experiments where water was used as the poor solvent, a negative surface charge was measured on the precipitated nanoparticles, which led to the long-time stability of the dispersion. It was argued that these charges originate presumably from either water or hydroxide adsorption at the hydrophobic nanoparticle surface or from impurities in the feed streams that preferentially adsorb on the precipitated nanoparticles. To elucidate the origin of this stabilization mechanism, we performed experiments wherein …

Stratification Dynamics in Drying Colloidal Mixtures

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. Th…

Evaporation-induced assembly of colloidal crystals

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 durin…

Micellization in Model Surfactant Systems

Formation of micelles in model lattice surfactant systems was studied by a novel methodology based on grand-canonical Monte Carlo simulations. The methodology involves combining free-energy information from a series of simulations in small systems by histogram reweighting. The solution osmotic pressure as a function of overall volume fraction of surfactant shows a sharp break at the critical micelle concentration (cmc) at sufficiently low temperatures. Studies in larger systems at appropriate values of the surfactant chemical potential are used to investigate the size distribution of micellar aggregates. The methodology allows for a clear distiction between micellization and macroscopic pha…

Axial dispersion of Brownian colloids in microfluidic channels

Controlled production of patchy particles from the combined effects of nanoprecipitation and vitrification

Using molecular dynamics simulations, we study a simple and scalable method for fabricating patchy nanoparticles via the assembly of binary polymer blends under a rapid solvent exchange. Patchiness can be achieved by incorporating a glassy component, which kinetically traps the particle morphology along the path to the equilibrium configuration. Our simulations reveal that the number of surface patches increases for larger nanoparticles and for more asymmetric blend ratios, while the size distribution of the patches remains rather uniform. Other than multi-patch nanoparticles, Janus structures have been obtained for small nanoparticles. Further, ribbon structures with elongated surface doma…

Structured Nanoparticles from the Self-Assembly of Polymer Blends through Rapid Solvent Exchange

Molecular dynamics simulations were performed to study systematically the rapid mixing of a polymer blend in solution with a miscible nonsolvent. In agreement with experiments, we observe that polymers self-assemble into complex nanoparticles, such as Janus and core-shell particles, when the good solvent is displaced by the poor solvent. The emerging structures can be predicted on the basis of the surface tensions between the polymers as well as between the polymers and the surrounding liquid. Furthermore, the size of the nanoparticles can be independently tuned through the mixing rate and the polymer concentration in the feed stream; meanwhile, the composition of the nanoparticles can be c…

Efficient mesoscale hydrodynamics: Multiparticle collision dynamics with massively parallel GPU acceleration

Abstract We present an efficient open-source implementation of the multiparticle collision dynamics (MPCD) algorithm that scales to run on hundreds of graphics processing units (GPUs). We especially focus on optimizations for modern GPU architectures and communication patterns between multiple GPUs. We show that a mixed-precision computing model can improve performance compared to a fully double-precision model while still providing good numerical accuracy. We report weak and strong scaling benchmarks of a reference MPCD solvent and a benchmark of a polymer solution with research-relevant interactions and system size. Our MPCD software enables simulations of mesoscale hydrodynamics at lengt…

Multi-scale simulations of polymeric nanoparticle aggregation during rapid solvent exchange.

Using a multi-scale approach which combines both molecular dynamics (MD) and kinetic Monte Carlo (KMC) simulations, we study a simple and scalable method for fabricating charge-stabilized nanoparticles through a rapid solvent exchange, i.e., Flash NanoPrecipitation (FNP). This multi-scale approach is based on microscopic information from MD simulations and uses a KMC algorithm to access macroscopic length- and time scales, which allows direct comparison with experiments and quantitative predictions. We find good agreement of our simulation results with the experiments. In addition, the model allows us to understand the aggregation mechanism on both microscopic and macroscopic levels and det…

Erratum: “Evaporation-induced assembly of colloidal crystals” [J. Chem. Phys. 149, 094901 (2018)]

Phase behavior of low-functionality, telechelic star block copolymers.

We apply state-of-the-art, Grand Canonical Monte Carlo simulations to determine the self-organization and phase behavior of solutions of block copolymer stars. The latter consist of f AB-block copolymers with N monomers each, which contain a solvophilic block A and solvophobic block B, and which are tethered on a common center on their A-side. We vary the degree of polymerization N and the relative composition of the block copolymer arms and investigate the interplay between macrophase and microphase separation in the system. Preliminary results of the effect of increasing the number of arms, f of the stars are also presented.

Self-assembly scenarios of block copolymer stars

We examine the self-organization scenarios of star-shaped AB-block copolymers, consisting of a solvophilic A-block and a solvophobic B-block, in which f such blocks are chemically anchored on a common centre on their A-parts, leaving the B-blocks exposed on their exterior. We employ a lattice model and we perform Grand Canonical Monte Carlo simulations for the case f = 6, varying thereby the percentage of attractive monomers as well as the concentration of stars. In agreement with previous studies on the low-functionality case f = 3 [F. Lo Verso, A.Z. Panagiotopoulos, and C.N. Likos, Phys. Rev. E 79, 010401(R) (2009)], we find that when the majority of monomers in the star are attractive, m…

Stratification in Drying Polymer–Polymer and Colloid–Polymer Mixtures

Drying polymer-polymer and colloid-polymer mixtures were studied using Langevin dynamics computer simulations. Polymer-polymer mixtures vertically stratified into layers, with the shorter polymers enriched near the drying interface and the longer polymers pushed down toward the substrate. Colloid-polymer mixtures stratified into a polymer-on-top structure when the polymer radius of gyration was comparable to or smaller than the colloid diameter, and a colloid-on-top structure otherwise. We also developed a theoretical model for the drying mixtures based on dynamical density functional theory, which gave excellent quantitative agreement with the simulations for the polymer-polymer mixtures a…

Self-Assembly of Polymer Blends and Nanoparticles through Rapid Solvent Exchange.

Molecular dynamics simulations were performed to study the fabrication of polymeric colloids containing inorganic nanoparticles (NPs) via the flash nanoprecipitation (FNP) technique. During this process, a binary polymer blend, initially in a good solvent for the polymers, is rapidly mixed with NPs and a poor solvent for the polymers that is miscible with the good solvent. The simulations reveal that the polymers formed Janus particles with NPs distributed either on the surface of the aggregates, throughout their interior, or aligned at the interface between the two polymer domains, depending on the NP-polymer and NP-solvent interactions. The loading and surface density of NPs can be contro…