Search

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 propose a method to calculate the equilibrium contact angle of heterogeneous 3-phase solid/fluid/fluid systems using molecular dynamics simulations. The proposed method, which combines the phantom-wall method [F. Leroy and F. Müller-Plathe, J. Chem. Phys. 133, 044110 (2010)] and Bennett's acceptance ratio approach [C. H. Bennett, J. Comput. Phys. 22, 245 (1976)], is able to calculate the solid/fluid surface tension relative to the solid surface energy. The calculated relative surface tensions can then be used in Young's equation to estimate the equilibrium contact angle. A fluid droplet is not needed for the proposed method, in contrast to the situation for direct simulations of contact angles. In addition, while prior free-energy based methods for contact angles mainly focused on the wetting of fluids in coexistence with their vapor on solid surfaces, the proposed approach was designed to study the contact angles of fluid mixtures on solid surfaces above the fluid saturation pressures. Using the proposed approach, the contact angles of binary Lennard-Jones fluid mixtures on a non-polar solid substrate were calculated at various interaction parameters and the contact angle of water in equilibrium with CO2 on a hydrophilic polar silica surface was obtained. For both nonpolar and polar systems, the calculated contact angles from the proposed method were in agreement with those obtained from the geometry of a cylindrical droplet. The computational cost of the proposed method was found to be comparable to that of simulations that use fluid droplets, but the new method provides a way to calculate the contact angle directly from Young's equation without ambiguity. [ABSTRACT FROM AUTHOR]

The phase equilibria of water with two higher hydrocarbons, n-hexane and n-butane, have been calculated for a range of temperatures. The relatively simple pairwise potentials become increasingly inaccurate for mixture predictions as the asymmetry in molecular size of unlike molecules grows.

Research was conducted to examine a new model for the n-alkane homologous series which has been created and parametrized to the vapor-liquid coexistence characteristics for a range of chain lengths. The Buckingham exponential-6 potential was used to describe the nonbonded interaction energy. The model was compared to the TraPPE and NERD models which showed that all three models reproduce adequately experimental saturated liquid densities. Results for the mixture behavior of four mixtures indicate that simulations using united-atom models fit to pure component data can be utilized to determine the phase behavior of complex non-polar mixtures with reasonable accuracy.

The phase equilibria of binary mixtures of water with methane or ethane has been determined. The Monte Carlo method was used in computing for Henry's Law constant of the two hydrocarbons in water over a wide range of temperatures. Results show that as the pressure increases the mutual solubility also increases. The mixture is then best described by the compositions of the two coexisting phases. The computed results show good agreement with the experimental data. However, the exp-6 model is more accurate than the MSPC/E model near the pure water critical point.

Triblock Janus colloids, which are colloidal spheres decorated with attractive patches at each pole, have recently generated significant interest as potential building blocks for functional materials. Their inherent anisotropy is known to induce self-assembly into open structures at moderate temperatures and pressures, where they are stabilized over close-packed crystals by entropic effects. We present a numerical investigation of the equilibrium phases of triblock Janus particles with many different patch geometries in three dimensions, using Monte Carlo simulations combined with free energy calculations. In all cases, we find that the free energy difference between crystal polymorphs is less than 0.2 kBT per particle. By varying the patch fraction and interaction range, we show that large patches stabilize the formation of structures with four bonds per patch over those with three. This transition occurs abruptly above a patch fraction of 0.30 and has a strong dependence on the interaction range. Furthermore, we find that a short interaction range favors four bonds per patch, with longer range increasingly stabilizing structures with only three bonds per patch. By quantifying the effect of patch geometry on the stability of the equilibrium crystal structures, we provide insights into the fundamental design rules for constructing complex colloidal crystals. [ABSTRACT FROM AUTHOR]

A fixed-point charge potential model for water was developed. The model characterizes the vapor-liquid coexistence properties namely saturated vapor and liquid densities, vapor pressures and liquid structure of water. The Buckingham exponential-6 potential model was utilized to describe the nonpolar interactions. Critical temperature was recorded at within 0.2%, critical density at 8% and saturated liquid densities to an average of 2.5%. The model may be useful in understanding environmental remediation, biological systems and separation phenomena.

Alternative methods for determining the critical micelle concentration (cmc) are investigated using canonical and grand canonical Monte Carlo simulations of a lattice surfactant model. A common measure of the cmc is the "free" (unassociated) surfactant concentration in the presence of micellar aggregates. Many prior simulations of micellizing systems have observed a decrease in the free surfactant concentration with overall surfactant loading for both ionic and nonionic surfactants, contrary to theoretical expectations from mass-action models of aggregation. In the present study, we investigate a simple lattice nonionic surfactant model in implicit solvent, for which highly reproducible simulations are possible in both the canonical (NVT) and grand canonical (μVT) ensembles. We confirm the previously observed decrease of free surfactant concentration at higher overall loadings and propose an algorithm for the precise calculation of the excluded volume and effective concentration of unassociated surfactant molecules in the accessible volume of the solution. We find that the cmc can be obtained by correcting the free surfactant concentration for volume exclusion effects resulting from the presence of micellar aggregates. We also develop an improved method for determination of the cmc based on the maximum in curvature for the osmotic pressure curve determined from μVT simulations. Excellent agreement in cmc and other micellar properties between NVT and μVT simulations of different system sizes is observed. The methodological developments in this work are broadly applicable to simulations of aggregating systems using any type of surfactant model (atomistic/coarse grained) or solvent description (explicit/implicit). [ABSTRACT FROM AUTHOR]

We study the effective interactions and phase behavior of star polymer-colloid mixtures through theory and Monte Carlo simulations. We extend previous theoretical approaches for calculating the effective star-colloid pair potential to take into account attractive contributions, which become significant at worsening solvent conditions. In order to assess the validity of our simulation and theory, we compute the effective interactions via virtual move parallel tempering Monte Carlo simulations using a microscopic bead-spring model for the star polymer and achieve excellent agreement. Finally, we perform grand canonical Monte Carlo simulations of the coarse-grained systems to study the effect of solvent quality on the phase behavior. [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 investigate the thermodynamic behavior of a model polymer-grafted nanoparticle (GNP) system on a fine lattice, using grand canonical Monte Carlo simulations, to compare and contrast the validity of two different models for GNPs: "nanoparticle amphiphiles" versus "patchy particles." In the former model, continuous self-assembly processes are expected to dominate the system, whereas the latter are characterized by first-order phase separation into novel equilibrium phases such as "empty liquids." We find that, in general, considering GNPs as amphiphiles within the framework of a recent mean-field theory [Pryamtisyn et al., J. Chem. Phys. 131, 221102 (2009)] provides a qualitatively accurate description of the thermodynamics of GNP systems, revealing either first-order phase separation into two isotropic phases or continuous self-assembly. Our model GNPs display no signs of empty liquid formation, suggesting that these nanoparticles do not provide a route to such phases. [ABSTRACT FROM AUTHOR]

The mean ionic activity coefficients of aqueous NaCl solutions of varying concentrations at 298.15 K and 1 bar have been obtained from molecular dynamics simulations by gradually turning on the interactions of an ion pair inserted into the solution. Several common non-polarizable water and ion models have been used in the simulations. Gibbs-Duhem equation calculations of the thermodynamic activity of water are used to confirm the thermodynamic consistency of the mean ionic activity coefficients. While the majority of model combinations predict the correct trends in mean ionic activity coefficients, they overestimate their values at high salt concentrations. The solubility predictions also suffer from inaccuracies, with all models underpredicting the experimental values, some by large factors. These results point to the need for further ion and water model development. [ABSTRACT FROM AUTHOR]

Molecular dynamics and Monte Carlo simulations have been carried out to obtain thermodynamic and transport properties of the binary mixture H2O+NaCl at temperatures from T = 298 to 473 K. In particular, vapor pressures, liquid densities, viscosities, and vapor-liquid interfacial tensions have been obtained as functions of pressure and salt concentration. Several previously proposed fixedpoint-charge models that include either Lennard-Jones (LJ) 12-6 or exponential-6 (Exp6) functional forms to describe non-Coulombic interactions were studied. In particular, for water we used the SPC and SPC/E (LJ) models in their rigid forms, a semiflexible version of the SPC/E (LJ) model, and the Errington-Panagiotopoulos Exp6 model; for NaCl, we used the Smith-Dang and Joung-Cheatham (LJ) parameterizations as well as the Tosi-Fumi (Exp6) model. While none of the model combinations are able to reproduce simultaneously all target properties, vapor pressures are well represented using the SPC plus Joung-Cheathem model combination, and all LJ models do well for the liquid density, with the semiflexible SPC/E plus Joung-Cheatham combination being the most accurate. For viscosities, the combination of rigid SPC/E plus Smith-Dang is the best alternative. For interfacial tensions, the combination of the semiflexible SPC/E plus Smith-Dang or Joung-Cheatham gives the best results. Inclusion of water flexibility improves the mixture densities and interfacial tensions, at the cost of larger deviations for the vapor pressures and viscosities. The Exp6 water plus Tosi-Fumi salt model combination was found to perform poorly for most of the properties of interest, in particular being unable to describe the experimental trend for the vapor pressure as a function of salt concentration. [ABSTRACT FROM AUTHOR]

We employ extensive computer simulations to study the flow behavior of spherical, nanoscale colloids in a viscoelastic solvent under Poiseuille flow. The systems are confined in a slit-like microfluidic channel, and viscoelasticity is introduced explicitly through the inclusion of polymer chains on the same length scale as the dispersed solute particles. We systematically study the effects of flow strength and polymer concentration, and identify a regime in which the colloids migrate to the centerline of the microchannel, expelling the polymer chains to the sides. This behavior was recently identified in experiments, but a detailed understanding of the underlying physics was lacking. To this end, we provide a detailed analysis of this phenomenon and discuss ways to maximize its effectiveness. The focusing mechanism can be exploited to separate and capture particles at the sub-micrometer scale using simple microfluidic devices, which is a crucial task for many biomedical applications, such as cell counting and genomic mapping. [ABSTRACT FROM AUTHOR]

In this contribution, we develop a coarse-graining methodology for mapping specific block copolymer systems to bead-spring particle-based models. We map the constituent Kuhn segments to Lennard-Jones particles, and establish a semi-empirical correlation between the experimentally determined Flory-Huggins parameter χ and the interaction of the model potential. For these purposes, we have performed an extensive set of isobaric-isothermal Monte Carlo simulations of binary mixtures of Lennard-Jones particles with the same size but with asymmetric energetic parameters. The phase behavior of these monomeric mixtures is then extended to chains with finite sizes through theoretical considerations. Such a top-down coarse-graining approach is important from a computational point of view, since many characteristic features of block copolymer systems are on time and length scales which are still inaccessible through fully atomistic simulations. We demonstrate the applicability of our method for generating parameters by reproducing the morphology diagram of a specific diblock copolymer, namely, poly(styrene-b-methyl methacrylate), which has been extensively studied in experiments. [ABSTRACT FROM AUTHOR]

It has been recently suggested that hydroxide ions can be formed in the electrolyte of molten carbonate fuel cells when water vapor is present. The hydroxide can replace carbonate in transporting electrons across the electrolyte, thereby reducing the CO2 separation efficiency of the fuel cell although still producing electricity. In this work, we obtain the equilibrium concentration of hydroxide in five molten alkali carbonate salts from molecular simulations. The results reveal that there can be a substantial amount of hydroxide in the electrolyte at low partial pressures of CO2. In addition, we find that the equilibrium concentration of molecular water dissolved in the electrolyte is over two orders of magnitude higher than that of CO2. Increasing the size and polarizability (or in other words reducing the "hardness") of the cations present in the electrolyte can reduce the hydroxide fraction, but at the cost of lowering ionic conductivity. [ABSTRACT FROM AUTHOR]

We use the weighted histogram analysis method [S. Kumar, D. Bouzida, R. H. Swendsen, P. A. Kollman, and J. M. Rosenberg, J. Comput. Chem. 13, 1011 (1992)] to calculate the free energy surface of the ST2 model of water as a function of density and bond-orientational order. We perform our calculations at deeply supercooled conditions (T = 228.6 K, P = 2.2 kbar; T = 235 K, P = 2.2 kbar) and focus our attention on the region of bond-orientational order that is relevant to disordered phases. We find a first-order transition between a low-density liquid (LDL, ρ ≈ 0.9 g/cc) and a high-density liquid (HDL, ρ ≈ 1.15 g/cc), confirming our earlier sampling of the free energy surface of this model as a function of density [Y. Liu, A. Z. Panagiotopoulos, and P. G. Debenedetti, J. Chem. Phys. 131, 104508 (2009)]. We demonstrate the disappearance of the LDL basin at high pressure and of the HDL basin at low pressure, in agreement with independent simulations of the system's equation of state. Consistency between directly computed and reweighted free energies, as well as between free energy surfaces computed using different thermodynamic starting conditions, confirms proper equilibrium sampling. Diffusion and structural relaxation calculations demonstrate that equilibration of the LDL phase, which exhibits slow dynamics, is attained in the course of the simulations. Repeated flipping between the LDL and HDL phases in the course of long molecular dynamics runs provides further evidence of a phase transition. We use the Ewald summation with vacuum boundary conditions to calculate long-ranged Coulombic interactions and show that conducting boundary conditions lead to unphysical behavior at low temperatures. [ABSTRACT FROM AUTHOR]

Thin films of water under nanoscopic confinement are prevalent in natural and manufactured materials. To investigate the equilibrium and dynamic behavior of water in such environments, we perform molecular dynamics simulations of water confined between atomistically detailed hydrophobic plates at T = 298 K for pressures (-0.1) <= P <= 1.0 GPa and plate separations of 0.40 <= d <= 0.80 nm. From these simulations, we construct an expanded P-d phase diagram for confined water, and identify and characterize a previously unreported confined monolayer ice morphology. We also study the decompression-induced sublimation of bilayer ice in a d = 0.6 nm slit, employing principal component analysis to synthesize low-dimensional embeddings of the drying trajectories and develop insight into the sublimation mechanism. Drying is observed to proceed by the nucleation of a bridging vapor cavity at one corner of the crystalline slab, followed by expansion of the cavity along two edges of the plates, and the subsequent recession of the remaining promontory of bilayer crystal into the bulk fluid. Our findings have implications for the understanding of diverse phenomena in materials science, nanofluidics, and protein folding and aggregation. [ABSTRACT FROM AUTHOR]

Coarse-grained models of poly(ethylene oxide) oligomer-grafted nanoparticles are established by matching their structural distribution functions to atomistic simulation data. Coarse-grained force fields for bulk oligomer chains show excellent transferability with respect to chain lengths and temperature, but structure and dynamics of grafted nanoparticle systems exhibit a strong dependence on the core-core interactions. This leads to poor transferability of the core potential to conditions different from the state point at which the potential was optimized. Remarkably, coarse graining of grafted nanoparticles can either accelerate or slowdown the core motions, depending on the length of the grafted chains. This stands in sharp contrast to linear polymer systems, for which coarse graining always accelerates the dynamics. Diffusivity data suggest that the grafting topology is one cause of slower motions of the cores for short-chain oligomer-grafted nanoparticles; an estimation based on transition-state theory shows the coarse-grained core-core potential also has a slowing-down effect on the nanoparticle organic hybrid materials motions; both effects diminish as grafted chains become longer. [ABSTRACT FROM AUTHOR]

The diffusivity and structural relaxation characteristics of oligomer-grafted nanoparticles have been investigated with simulations of a previously proposed coarse-grained model at atmospheric pressure. Solvent-free, polymer-grafted nanoparticles as well as grafted nanoparticles in a melt were compared to a reference system of bare (ungrafted) particles in a melt. Whereas longer chains lead to a larger hydrodynamic radius and lower relative diffusivity for grafted particles in a melt, bulk solvent-free nanoparticles with longer chains have higher relative diffusivities than their short chain counterparts. Solvent-free nanoparticles with short chains undergo a glass transition as indicated by a vanishing diffusivity, diverging structural relaxation time and the formation of body-centered-cubic-like order. Nanoparticles with longer chains exhibit a more gradual increase in the structural relaxation time with decreasing temperature and concomitantly increasing particle volume fraction. The diffusivity of the long chain nanoparticles exhibits a minimum at an intermediate temperature and volume fraction where the polymer brushes of neighboring particles overlap, but must stretch to fill the interparticle space. [ABSTRACT FROM AUTHOR]

Dissipative particle dynamics simulations were used to study the effects of mixing time, solute solubility, solute and diblock copolymer concentrations, and copolymer block length on the rapid coprecipitation of polymer-protected nanoparticles. The simulations were aimed at modeling Flash NanoPrecipitation, a process in which hydrophobic solutes and amphiphilic block copolymers are dissolved in a water-miscible organic solvent and then rapidly mixed with water to produce composite nanoparticles. A previously developed model by Spaeth et al. [J. Chem. Phys. 134, 164902 (2011)] was used. The model was parameterized to reproduce equilibrium and transport properties of the solvent, hydrophobic solute, and diblock copolymer. Anti-solvent mixing was modeled using time-dependent solvent-solute and solvent-copolymer interactions. We find that particle size increases with mixing time, due to the difference in solute and polymer solubilities. Increasing the solubility of the solute leads to larger nanoparticles for unfavorable solute-polymer interactions and to smaller nanoparticles for favorable solute-polymer interactions. A decrease in overall solute and polymer concentration produces smaller nanoparticles, because the difference in the diffusion coefficients of a single polymer and of larger clusters becomes more important to their relative rates of collisions under more dilute conditions. An increase in the solute-polymer ratio produces larger nanoparticles, since a collection of large particles has less surface area than a collection of small particles with the same total volume. An increase in the hydrophilic block length of the polymer leads to smaller nanoparticles, due to an enhanced ability of each polymer to shield the nanoparticle core. For unfavorable solute-polymer interactions, the nanoparticle size increases with hydrophobic block length. However, for favorable solute-polymer interactions, nanoparticle size exhibits a local minimum with respect to the hydrophobic block length. Our results provide insights on ways in which experimentally controllable parameters of the Flash NanoPrecipitation process can be used to influence aggregate size and composition during self-assembly. [ABSTRACT FROM AUTHOR]

We have developed explicit- and implicit-solvent models for the flash nanoprecipitation process, which involves rapid coprecipitation of block copolymers and solutes by changing solvent quality. The explicit-solvent model uses the dissipative particle dynamics (DPD) method and the implicit-solvent model uses the Brownian dynamics (BD) method. Each of the two models was parameterized to match key properties of the diblock copolymer (specifically, critical micelle concentration, diffusion coefficient, polystyrene melt density, and polyethylene glycol radius of gyration) and the hydrophobic solute (aqueous solubility, diffusion coefficient, and solid density). The models were simulated in the limit of instantaneous mixing of solvent with antisolvent. Despite the significant differences in the potentials employed in the implicit- and explicit-solvent models, the polymer-stabilized nanoparticles formed in both sets of simulations are similar in size and structure; however, the dynamic evolution of the two simulations is quite different. Nanoparticles in the BD simulations have diffusion coefficients that follow Rouse behavior (D ∝ M-1), whereas those in the DPD simulations have diffusion coefficients that are close to the values predicted by the Stokes-Einstein relation (D ∝ R-1). As the nanoparticles become larger, the discrepancy between diffusion coefficients grows. As a consequence, BD simulations produce increasingly slower aggregation dynamics with respect to real time and result in an unphysical evolution of the nanoparticle size distribution. Surface area per polymer of the stable explicit-solvent nanoparticles agrees well with experimental values, whereas the implicit-solvent nanoparticles are stable when the surface area per particle is roughly two to four times larger. We conclude that implicit-solvent models may produce questionable results when simulating nonequilibrium processes in which hydrodynamics play a critical role. [ABSTRACT FROM AUTHOR]

Nonlinear dimensionality reduction techniques can be applied to molecular simulation trajectories to systematically extract a small number of variables with which to parametrize the important dynamical motions of the system. For molecular systems exhibiting free energy barriers exceeding a few kBT, inadequate sampling of the barrier regions between stable or metastable basins can lead to a poor global characterization of the free energy landscape. We present an adaptation of a nonlinear dimensionality reduction technique known as the diffusion map that extends its applicability to biased umbrella sampling simulation trajectories in which restraining potentials are employed to drive the system into high free energy regions and improve sampling of phase space. We then propose a bootstrapped approach to iteratively discover good low-dimensional parametrizations by interleaving successive rounds of umbrella sampling and diffusion mapping, and we illustrate the technique through a study of alanine dipeptide in explicit solvent. [ABSTRACT FROM AUTHOR]

We perform histogram-reweighting grand canonical Monte Carlo simulations of the Lennard-Jones fluid confined between two parallel hard walls and determine the vapor-liquid critical and coexistence properties in the range of σ≤H≤6σ and 10σ≤Lx,Ly≤28σ, where H is the wall separation, Lx=Ly is the system size and σ is the characteristic length. By matching the probability distribution of the ordering operator, P(M), to the three-dimensional (3D) and two-dimensional (2D) Ising universality classes according to the mixed-field finite-size scaling approach, we establish a “phase diagram” in the (H,L) plane, showing the boundary between four types of behavior: 3D, quasi-3D, quasi-2D, and 2D. In order to facilitate 2D critical point calculation, we present a four-parameter analytical expression for the 2D Ising universal distribution. We show that the infinite-system-size critical points obtained by extrapolation from the apparent 3D and 2D critical points have only minor differences with each other. In agreement with recent reports in the literature [Jana et al., J. Chem. Phys. 130, 214707 (2009)], we find departure from linearity in the relationship between critical temperature and inverse wall separation, as well as nonmonotonic dependence of the critical density and the liquid density at coexistence upon wall separation. Additional studies of the ST2 model of water show similar behavior, which suggests that these are quite general properties of confined fluids. [ABSTRACT FROM AUTHOR]

We use molecular dynamics simulations over microsecond time scales to study the micellization behavior of recently proposed continuum-space, coarse grained surfactant models. In particular, we focus on the MARTINI model by Marrink et al. [J. Phys. Chem. B 111, 7812 (2007)] and a model by Shinoda et al. [Soft Matter 4, 2454 (2008)]. We obtain the critical micelle concentration (cmc) and equilibrium aggregate size distributions at low surfactant loadings. We present evidence justifying modest extrapolations for determining the cmc at low temperatures, where significant sampling difficulties remain. The replica exchange method provides only modest improvements of sampling efficiency for these systems. We find that the two coarse grained models significantly underpredict experimental cmc near room temperature for zwitterionic surfactants, but are closer to measured values for nonionic ones. The aggregation numbers for both zwitterionic and nonionic surfactants are near those observed experimentally, but the temperature dependence of the cmc is incorrect in both cases, because of the use of an unstructured solvent. Possible refinements to the models to bring them into quantitative agreement with experiment are discussed. [ABSTRACT FROM AUTHOR]

Recent experimental results demonstrated that polymer grafted nanoparticles in solvents display self-assembly behavior similar to the microphase separation of block copolymers and other amphiphiles. We present a mean-field theory and complementary computer simulations to shed light on the parametric underpinnings of the experimental observations. Our theory suggests that such self-assembled structures occur most readily when the nanoparticle size is comparable to the radius of gyration of the polymer brush chains. Much smaller particle sizes are predicted to yield uniform particle dispersions, while larger particles are expected to agglomerate due to phase separation from the solvent. Selected aspects of our theoretical predictions are corroborated by computer simulations. [ABSTRACT FROM AUTHOR]

A lattice model for diblock surfactants that incorporates orientational bonding has been developed for studying self-assembly in dilute solutions. Using grand canonical Monte Carlo simulations with histogram reweighting and mixed field finite size scaling, we examine the effect of amphiphile architecture on phase transitions and distinguish between first order transitions that create a disordered liquid phase and higher order transitions that indicate the formation of finite sized aggregates. As the solution temperature increases, we find that the critical micelle concentration for the orientational bonding model surfactants reaches a minimum value at a temperature that can be controlled by varying the number of bonding orientations between the solvophobic surfactant monomers and the implicit solvent. This trend is qualitatively similar to experimental data for ionic and nonionic surfactants in aqueous solutions. A comparable dependence on temperature is observed in the limit of amphiphile solubility for phase separating systems. None of the model surfactants considered here undergo both a first and a higher order transition over the range of densities and temperatures examined. [ABSTRACT FROM AUTHOR]