Molecular dynamics simulations studies of nanoparticles in an isotropic liquid crystal matrix: Single particle behavior and pairwise interactions

We report results of molecular dynamics simulation studies of the behavior of spherical nanoparticles (NPs) in a dense isotropic nematogen matrix comprised of soft spherocylinders (SSCs). The SSCs exhibit a tendency for frustrated planar anchoring at the NP surface that results in a long-range (compared to the size of the NPs and SSCs) reduction in local orientational ordering and increased fluctuations in local orientational ordering compared to the pure isotropic phase of the SSCs. The potential of mean force between two nanoparticles exhibits a novel long-range repulsive tail separated from short-range molecular packing peaks by a shallow local minimum in free energy. The long-range repulsion is caused by NP-induced ordering fluctuations while the shallow minimum results from increased local ordering within the confinement region in between two NPs. The influence of the NPs on local orientational order in the nematogen matrix and the nematogen-induced interaction between NPs are found to depend strongly on the size of the NPs.     

Effect of interacting nanoparticles on the ordered morphology of block copolymer/nanoparticle mixtures

We investigated the effect of hard additives, that is, magnetic nanoparticles (NPs) and metal NPs, on the ordered morphology of block copolymers by varying the NP concentration. To characterize the structural changes of a block copolymer associated with different NP loadings, small-angle X-ray scattering and transmission electron microscopy were performed. Monodisperse maghemite (γ-Fe₂O₃) NPs (7 nm in diameter) and silver (Ag) NPs (6 nm in diameter) with surfaces modified with oleic acids were synthesized, and a cylinder-forming poly(styrene-block-isoprene) diblock copolymer was used as a structure-directing matrix for the NPs. As the NP concentration increased, domains of NP aggregates were observed for both magnetic and metal NPs. In the case of mixtures of cylinder-forming poly(styrene-block-isoprene) and Ag NPs with weak particle–particle interactions, random aggregates of Ag NPs were observed, and the ordered morphology of the block copolymer lost its long-range order with an increase in the NP concentration. However, regular, latticelike aggregates obtained with γ-Fe₂O₃ NPs, because of the strong interparticle interactions, induced an intriguing morphological transformation from hexagonal cylinders to body-centered-cubic spheres via undulated cylinders, whereas the neat block copolymer did not show such a morphological transition over a wide range of temperatures. The interplay between magnetic NPs and the block copolymer was also tested with magnetic NPs of different sizes. © 2006 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 3571–3579, 2006     

Role of steric factors in the mechanism of separation of structural isomers with nematic sorbents

The model of a nematic binary mixture of rigid biaxial molecules is used to calculate the coefficients of selectivity enhancement S _m/p (the ratio between the activity coefficients ofmeta- andpara-substituted benzenes as sorbates at infinite dilution) in the nematogenic matrix of a nematic liquid-crystal sorbent and the parameters of orientation order of the components. It was found that in systems of particles in which interaction between the particles consists in steric repulsion, the Sm/p coefficient is less than 1 and is practically constant in the isotropic phase, whereas in the nematic phase, S _m/p is greater than 1 and increases in parallel with the order parameter.     

Mesophase stability in binary mixtures of monotropic nematogens

Binary phase diagrams were constructed from laterally substituted methyl azo/ester derivatives, namely 4-(4″-substituted phenylazo)-3-methyl phenyl-4″-alkoxy benzoates (In_a–d). In this group of compound the unsubstituted and chloro-substituted derivatives possess the nematic phase monotropically, while the nitro and methyl analogues are enantiotropically nematogenic. The binary phase diagrams constructed were made once from the monotropic nematogens with each other, and another with the enantiotropic nematogens. In both the cases enantiotropic nematic phase was observed that covers a wide range of composition. The mesophase behaviour of all binary mixtures was investigated by differential scanning calorimetry (DSC) and polarised light microscopy (PLM). The nematic phase was exhibited in all binary mixtures. Independent of the alkoxy chain length, the entropy change, ΔS_N–I of the N–I transition of pure components was found to vary irregularly with the anisotropy of polarisability (?α_X) of the polar substituent, X.     

Role of Electrostatic Repulsion on the Viscosity of Bidisperse Silica Suspensions

The flow behavior of bidisperse aqueous silica suspensions has been studied at different electrolyte concentrations as a function of shear rate, total volume fraction of the particles, and volume ratio of small to large particles. It is shown that the range of the electrostatic repulsion plays an important role in determining the viscosity of the suspension. Binary mixtures of particles of longer range repulsive forces showed higher viscosities than the suspensions of shorter range electrostatic interactions. Bimodal suspensions of long-range interactions showed non-Newtonian behavior over wider ranges of shear due to the deformation of the ionic cloud around the particles, which is larger in these systems. The viscosity of bimodal suspensions used in this study was scaled with respect to the viscosity of the related monosized systems and the viscosity of one bimodal suspension at a fixed total volume fraction of the particles, employing our earlier scaling method. The model normalizes the effect of colloidal forces by introducing a scaling factor that collapses the data into a single curve for bimodal suspensions of a particular size ratio, and it is shown that the model is valid for systems with both short-range and long-range repulsive forces. Copyright 1999 Academic Press.     

Measurement of Colloidal Stability in Solutions of Simple, Nonadsorbing Polyelectrolytes

The stability of a solution of charged polystyrene particles in the presence of nonadsorbing polyelectrolyte macromolecules is measured using optical light scattering. The particles were negatively charged polystyrene latex spheres (0.5–1 μm diameter) while the macromolecules were simulated using negatively charged colloidal silica spheres (5–7 nm diameter). Because of the electrostatic repulsion between the particles, the solution is found to be stable against primary flocculation (irreversible flocculation into a primary energy minima). However, because of long-range attractive depletion forces, reversible secondary flocculation of the particles occurs into a local potential energy minimum. As observed with uncharged macromolecules, the polyelectrolyte first induces flocculation at a critical flocculation concentration (v*), but later restabilizes the system at a critical restabilization concentration (v**). These critical concentrations are found to decrease with decreasing macromolecule size and increasing particle size. The restabilized solutions are found to remain suspended for periods greater than 20 days. Comparison of the measured flocculation and restabilization results to predictions made using a recently developed force-balance model show qualitative agreement.     

Interactions between spherical colloids mediated by a liquid crystal: a molecular simulation and mesoscale study

Monte Carlo simulations and dynamic field theory (DyFT) are used to study the interactions between dilute spherical particles, dispersed in nematic and isotropic phases of a liquid crystal. A recently developed simulation method (expanded ensemble density of states) was used to determine the potential of mean force (PMF) between the two spheres as a function of their separation and size. The PMF was also calculated by a dynamic field theory that describes the evolution of the local tensor order parameter. Both methods reveal an overall attraction between the colloids in the nematic phase; in the isotropic phase, the overall attraction between the colloids is much weaker, whereas the repulsion at short range is stronger. In addition, both methods predict a new topology of the disclination lines, which arises when the particles approach each other. The theory is found to describe the results of simulations remarkably well, down to length scales comparable to the size of the molecules. At separations corresponding to the width of individual molecular layers on the particles' surface, the two methods yield different defect structures. We attribute this difference to the neglect of density inhomogeneities in the DyFT. We also investigate the effects of the size of spherical colloids on their interactions.     

Properties of the discotic-nematic to isotropic transition: influence of short range orientational order

The present investigation concerns the analysis of the influence of short range orientational correlation on the thermodynamic properties of discotic-nematic liquid crystals. Two-site cluster approximation is applied to the orientational molecular coordinates to include the short range orientational correlation. The role of short range orientational order, dispersion interaction, molecular length-to-width ratio and pressure on the thermodynamic and orientational behaviour of discotic nematogens close to the discotic-nematic to isotropic transition are analysed. It is observed that the short range orientational order has a strong influence on the thermodynamic properties and that the transition properties of both the calamitic and discotic mesogens exhibit quite similar behaviour.     

Properties of the discotic-nematic to isotropic transition: influence of short range orientational order

The present investigation concerns the analysis of the influence of short range orientational correlation on the thermodynamic properties of discotic-nematic liquid crystals. Two-site cluster approximation is applied to the orientational molecular coordinates to include the short range orientational correlation. The role of short range orientational order, dispersion interaction, molecular length-to-width ratio and pressure on the thermodynamic and orientational behaviour of discotic nematogens close to the discotic-nematic to isotropic transition are analysed. It is observed that the short range orientational order has a strong influence on the thermodynamic properties and that the transition properties of both the calamitic and discotic mesogens exhibit quite similar behaviour.     

Model for reversible nanoparticle assembly in a polymer matrix

The clustering of nanoparticles (NPs) in solutions and polymer melts depends sensitively on the strength and directionality of the NP interactions involved, as well as the molecular geometry and interactions of the dispersing fluids. Since clustering can strongly influence the properties of polymer-NP materials, we aim to better elucidate the mechanism of reversible self-assembly of highly symmetric NPs into clusters under equilibrium conditions. Our results are based on molecular dynamics simulations of icosahedral NP with a long-ranged interaction intended to mimic the polymer-mediated interactions of a polymer-melt matrix. To distinguish effects of polymer-mediated interactions from bare NP interactions, we compare the NP assembly in our coarse-grained model to the case where the NP interactions are purely short ranged. For the "control" case of NPs with short-ranged interactions and no polymer matrix, we find that the particles exhibit ordinary phase separation. By incorporating physically plausible long-ranged interactions, we suppress phase separation and qualitatively reproduce the thermally reversible cluster formation found previously in computations for NPs with short-ranged interactions in an explicit polymer-melt matrix. We further characterize the assembly process by evaluating the cluster properties and the location of the self-assembly transition. Our findings are consistent with a theoretical model for equilibrium clustering when the particle association is subject to a constraint. In particular, the density dependence of the average cluster mass exhibits a linear concentration dependence, in contrast to the square root dependence found in freely associating systems. The coarse-grained model we use to simulate NP in a polymer matrix shares many features of potentials used to model colloidal systems. The model should be practically valuable for exploring factors that control the dispersion of NP in polymer matrices where explicit simulation of the polymer matrix is too time consuming.     

Flexible polymers in solution in the isotropic phase of a nematogenic liquid crystal. Viscosity and Kerr effect

The behaviour of flexible polymers in solution in the isotropic phase of a nematogenic solvent was studied using viscosity and electric birefringence measurements. The results obtained show the existence of an intraction between the polymer and the local orientational order for polymers of low molecular weight, which disappears for those of high molecular weight. This change in behaviour can be explained by comparing the parameters characterizing the local orientational order to those characterizing the polymeric coil.     

Hydrodynamic Characteristics and Adsorption Particularity of Nanobiological Feedstock Along the Bed Height in a Novel Chromatography Column

Expanded bed adsorption (EBA) is a practical method for the separation of nanoparticulates. In order to analysis the local hydrodynamic and adsorption behavior of nanoparticle (NP)-based biological feedstock, a modified Nano Biotechnology Group EBA column with a 26-mm inner diameter was used to withdraw liquid from different axial positions of the column. Fabricated egg albumin (EA) NPs with an average size of 70?nm were employed as a model system and viral size/charge mimic to assess the relationship between hydrodynamic and adsorption performance of NPs at the different column regions. The effects of influential factors, including flow velocity and initial concentration of NPs, on NP hydrodynamic behavior and adsorption kinetics along the bed height were investigated. NP hydrodynamic studies confirmed that non-uniform behavior dominated the system and a decreasing trend of liquid mixing/dispersion with increase of bed height was observed in this column. The results demonstrated an increase in the mixing/dispersion at certain bed heights with the increase in both the velocity and feed initial concentration. Breakthrough curves were measured at various column points to determine the adsorption performance [dynamic binding capacity (DBC) and yield] in different bed positions/zones. Yield and DBC of NPs were improved along the bed height, whereas liquid velocity had the opposite effect. Increasing the initial concentration of NPs enhanced only the DBC. Separation of EA NPs under optimal conditions was 87?%, which is an excellent result for a one-pass frontal chromatography method.     

On modifying properties of polymeric melts by nanoscopic particles

We study geometric and energetic factors that partake in modifying properties of polymeric melts via inserting well‐dispersed nanoscopic particles (NP). Model systems are cis‐1,4‐polybutadiene melts including a single atomic clusters of size varied in the range 10–150 atoms (3–7 Å in radius; 0.1–1.5% v/v). We modify the interactions between the chains and the particle by tuning attractive van der Waals interactions. Using molecular dynamics, we study equilibrium fluctuations and dynamical properties at the interface. The NPs move in the polymer matrix in two different regimes corresponding to trapped and free diffusion, depending on the NP size. Furthermore, degree of crowding around the NP by the polymer chains is quantified. Effect of NP size and interaction strength both on volume and volumetric fluctuations is manifested in mechanical properties, quantified here by bulk modulus, K. Tuning NP size and nonbonded interactions results in ～15% enhancement in K by addition of a maximum of 1.5% v/v NP. © 2012 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys, 2012     

Hydrodynamic Characteristics and Adsorption Particularity of Nanobiological Feedstock Along the Bed Height in a Novel Chromatography Column

Expanded bed adsorption (EBA) is a practical method for the separation of nanoparticulates. In order to analysis the local hydrodynamic and adsorption behavior of nanoparticle (NP)-based biological feedstock, a modified Nano Biotechnology Group EBA column with a 26-mm inner diameter was used to withdraw liquid from different axial positions of the column. Fabricated egg albumin (EA) NPs with an average size of 70 nm were employed as a model system and viral size/charge mimic to assess the relationship between hydrodynamic and adsorption performance of NPs at the different column regions. The effects of influential factors, including flow velocity and initial concentration of NPs, on NP hydrodynamic behavior and adsorption kinetics along the bed height were investigated. NP hydrodynamic studies confirmed that non-uniform behavior dominated the system and a decreasing trend of liquid mixing/dispersion with increase of bed height was observed in this column. The results demonstrated an increase in the mixing/dispersion at certain bed heights with the increase in both the velocity and feed initial concentration. Breakthrough curves were measured at various column points to determine the adsorption performance [dynamic binding capacity (DBC) and yield] in different bed positions/zones. Yield and DBC of NPs were improved along the bed height, whereas liquid velocity had the opposite effect. Increasing the initial concentration of NPs enhanced only the DBC. Separation of EA NPs under optimal conditions was 87 %, which is an excellent result for a one-pass frontal chromatography method.

     

Modeling Heterogeneity in UO2 Nanoparticles Using X-ray Absorption Spectroscopy

EXAFS provides the capability to interrogate nanoparticle (NP) structure in atomistic detail without relying on long-range crystallinity. There is a limitation in that EXAFS provides averaged structural information, making it difficult to separate a small amount of heterogeneous structure from bulk. In this work, models were developed to extract surface-specific information from conventional EXAFS measurements collected on UO₂ NPs of varying size. Specifically, the surface terminating species of UO₂ NPs was determined from comparison of coordination numbers with geometric models while the origin of static disorder was interrogated from user-defined simulations. Results show that the degree of oxygenation on the NP surface does not significantly deviate from bulk surface and that static disorder is highly enhanced in NP surface layers but cannot be attributed to surface relaxation effects alone. The approach described herein has the potential to be adapted to a range of inorganic NP systems to interrogate surface structure.     

Single-chain compaction of long duplex DNA by cationic nanoparticles: modes of interaction and comparison with chromatin

The compaction of long duplex DNA by cationic nanoparticles (NP) used as a primary model of histone core particles has been investigated. We have systematically studied the effect of salt concentration, particle size, and particle charge by means of single-molecule observations-fluorescence microscopy (FM) and transmission electron microscopy (TEM)-and molecular dynamics (MD) simulations. We have found that the large-scale DNA compaction is progressive and proceeds through the formation of beads-on-a-string structures of various morphologies. The DNA adsorbed amount per particle depends weakly on NP concentration but increases significantly with an increase in particle size and is optimal at an intermediate salt concentration. Three different complexation mechanisms have been identified depending on the correlation between DNA and NPs in terms of geometry, chain rigidity, and electrostatic interactions: free DNA adsorption onto NP surface, DNA wrapping around NP, and NP collection on DNA chain.     

Surface molecular imprinting on polypropylene fibers for rhodamine B selective adsorption

The interaction between silver nanoparticles (Ag NPs) of different surface charge and surfactants relevant to the laundry cycle has been investigated to understand changes in speciation, both in and during transport from the washing machine. Ag NPs were synthesized to exhibit either a positive or a negative surface charge in solution conditions relevant for the laundry cycle (pH 10 and pH 7). These particles were characterized in terms of size and surface charge and compared to commercially laser ablated Ag NPs. The surfactants included anionic sodium dodecylbenzenesulfonate (LAS), cationic dodecyltrimethylammoniumchloride (DTAC) and nonionic Berol 266 (Berol). Surfactant-Ag NP interactions were studied by means of dynamic light scattering, Raman spectroscopy, zeta potential, and Quartz Crystal Microbalance. Mixed bilayers of CTAB and LAS were formed through a co-operative adsorption process on positively charged Ag NPs with pre-adsorbed CTAB, resulting in charge reversal from positive to negative zeta potentials. Adsorption of DTAC on negatively charged synthesized Ag NPs and negatively charged commercial Ag NPs resulted in bilayer formation and charge reversal. Weak interactions were observed for nonionic Berol with all Ag NPs via hydrophobic interactions, which resulted in decreased zeta potentials for Berol concentrations above its critical micelle concentration. Differences in particle size were essentially not affected by surfactant adsorption, as the surfactant layer thicknesses did not exceed more than a few nanometers. The surfactant interaction with the Ag NP surface was shown to be reversible, an observation of particular importance for hazard and environmental risk assessments.     

Dynamical simulation of dipolar Janus colloids: Equilibrium structure and thermodynamics

The static microstructures and thermodynamics of a colloidal dispersion of dipolar Janus (DJ) particles-that is, dipolar spheres in which each hemisphere is specified by a different charge interaction-have been investigated through simulation. DJ particles are modeled at a high level of detail with pairwise potentials represented as a sum of a spherically symmetric soft repulsion and an orientation-dependent electrostatic component using continuous potentials. The latter is important because it allows for the use of conventional molecular dynamics simulations, and is in contrast to the patch model and dipolar hard sphere model, which are discontinuous and therefore do not. The electrostatics are represented through a rigorous pointwise (PW) covering of two different hemispheres filled by points of corresponding charge. An isotropic coarse-graining (CG) of the PW models serves as a limit of the structure wherein the orientations of the DJ particles can be pairwise averaged. Over the range of volume fractions and DJ charge densities studied-consistent with reversible structures absent of long-range correlations-the CG model agrees well with the PW model with respect to equilibrium structure (isotropic pair correlation) and ensemble free energy. Time-dependent relaxation simulations of the PW model suggest that chain structures are not expected in liquid phases in contrast to that which has been observed for point dipole models of simple polar fluids.     

In situ coarsening study of inverse micelle-prepared Pt nanoparticles supported on γ-Al2O3: pretreatment and environmental effects

The thermal stability of inverse micelle prepared Pt nanoparticles (NPs) supported on nanocrystalline γ-Al(2)O(3) was monitored in situ under different chemical environments (H(2), O(2), H(2)O) via extended X-ray absorption fine-structure spectroscopy (EXAFS) and ex situ via scanning transmission electron microscopy (STEM). Drastic differences in the stability of identically synthesized NP samples were observed upon exposure to two different pre-treatments. In particular, exposure to O(2) at 400 °C before high temperature annealing in H(2) (800 °C) was found to result in the stabilization of the inverse micelle prepared Pt NPs, reaching a maximum overall size after moderate coarsening of ～1 nm. Interestingly, when an analogous sample was pre-treated in H(2) at ～400 °C, a final size of ～5 nm was reached at 800 °C. The beneficial role of oxygen in the stabilization of small Pt NPs was also observed in situ during annealing treatments in O(2) at 450 °C for several hours. In particular, while NPs of 0.5 ± 0.1 nm initial average size did not display any significant sintering (0.6 ± 0.2 nm final size), an analogous thermal treatment in hydrogen leads to NP coarsening (1.2 ± 0.3 nm). The same sample pre-dosed and annealed in an atmosphere containing water only displayed moderate sintering (0.8 ± 0.3 nm). Our data suggest that PtO(x) species, possibly modifying the NP/support interface, play a role in the stabilization of small Pt NPs. Our study reveals the enhanced thermal stability of inverse micelle prepared Pt NPs and the importance of the sample pre-treatment and annealing environment in the minimization of undesired sintering processes affecting the catalytic performance of nanosized particles.     

Random nematic structures in the absence of inherent frustrations

We study the impact of anisotropic nanoparticles (NPs) on the nematic liquid crystal (LC) order. Within a mesoscopic Flory–Huggins-type approach we have estimated regimes where LC–NP mixtures are essentially homogenous. Using a lattice Lebwohl–Lasher type approach we have also studied the impact of anisotropic NPs on LC ordering. We analysed the cases where the orientations of NPs are either frozen-in or could be varied, to which we refer as random field mixtures and annealed mixtures, respectively. In the latter case we have demonstrated the existence of qualitatively different regimes. In particular, we determined the concentration regime, where LC configurations resembling a transparent nematic phase could be observed. Such domain patterns are stabilised by NPs hindering the annihilation of topological defects in LC order.