Sedimentation of concentrated monodisperse colloidal suspensions: role of collective particle interaction forces

The sedimentation velocities and concentration profiles of low-charge, monodisperse hydroxylate latex particle suspensions were investigated experimentally as a function of the particle concentration to study the effects of the collective particle interactions on suspension stability. We used the Kossel diffraction technique to measure the particle concentration profile and sedimentation rate. We conducted the sedimentation experiments using three different particle sizes. Collective hydrodynamic interactions dominate the particle-particle interactions at particle concentrations up to 6.5 vol%. However, at higher particle concentrations, additional collective particle-particle interactions resulting from the self-depletion attraction cause particle aggregation inside the suspension. The collective particle-particle interaction forces play a much more important role when relatively small particles (500 nm in diameter or less) are used. We developed a theoretical model based on the statistical particle dynamics simulation method to examine the role of the collective particle interactions in concentrated suspensions in the colloidal microstructure formation and sedimentation rates. The theoretical results agree with the experimentally-measured values of the settling velocities and concentration profiles.     

Interactions of polarizable media in water: a molecular dynamics approach

We investigate the interactions of polarizable solutes in water as a function of the solute permittivity. A generic and computationally efficient simulation methodology for the investigation of systems involving dielectric discontinuities is introduced. We report results for interactions between two polarizable cylindrical solutes of nanometer dimensions, which demonstrate that the interactions between the solutes strongly depend on the solute permittivity epsilon. For low permittivity, epsilon approximately 1-2, the interactions are dominated by surface tension forces whose origin lies in the formation of a vapor cavity between the two hydrophobic solutes. This effect leads to a drying transition, where the intersolute force changes discontinuously at a specific solute-solute separation. We find that a moderate permittivity, epsilon approximately 20, enhances the solvation of the polarizable objects inhibiting this drying transition. In the limit of moderately high permittivity, the interactions are dominated by solvation forces. These forces are much larger than those calculated using macroscopic models of dielectrics, which consider water as a continuum dielectric medium. Our results emphasize the importance of including the solvent explicitly to investigate dielectric discontinuities and interactions between polarizable media in water.     

Comparison of implicit solvent models for the simulation of protein-surface interactions

Empirical force field-based molecular simulations can provide valuable atomistic-level insights into protein-surface interactions in aqueous solution. While the implicit treatment of solvation effects is desired as a means of improving simulation efficiency, existing implicit solvent models were primarily developed for the simulation of peptide or protein behavior in solution alone, and thus may not be appropriate for protein interactions with synthetic material surfaces. The objective of this research was to calculate the change in free energy as a function of surface-separation distance for peptide-surface interactions using different empirical force field-based implicit solvation models (ACE, ASP, EEF1, and RDIE with the CHARMM 19 force field), and to compare these results with the same calculations conducted using density functional theory (DFT) combined with the self-consistent reaction field (SCRF) implicit solvation model. These comparisons show that distinctly different types of behavior are predicted with each implicit solvation method, with ACE providing the best overall agreement with DFT/SCRF calculations. These results also identify areas where ACE is in need of improvement for this application and provide a basis for subsequent parameter refinement.     

Solvation of coumarin 153 in supercritical fluoroform

We present a study of local density augmentation around an attractive solute (i.e., giving rise to more attractive interaction with the solvent than solvent-solvent interactions) in supercritical fluoroform. This work is based on molecular dynamics simulations of coumarin 153 in supercritical fluoroform at densities both above and below the critical density, ranging from dilute gas-like to liquid-like, at a reduced temperature (T/T(c)) of 1.03. We focused on studying the structure of the solvation shell and the variation of the solute electronic absorption and emission shifts with density. Quantum calculations at the density functional theory (DFT) level were run on the solute in the ground state, and time-dependent DFT calculations were performed in the solute excited state in order to determine the solute-solvent potential parameters. The results obtained for the Stokes shift are in agreement with the experimental measurements. To evaluate local density augmentation from simulations, we used two different definitions, one based on the solvation number and the other derived from solvatochromic shifts. In the former case, the agreement with experimental results is good, while, in the latter case, better agreement is achieved by perturbatively including the induced-dipole contribution to the solvation energy.     

Stability and Flow-Induced Flocculation of Fumed Silica Suspensions in Mixture of Water-Glycerol

A systematic investigation has been performed to relate the effect of glycerol composition to the rheological properties of aqueous suspensions of hydrophilic fumed silica at pH far from the isoelectric point. Steady state/dynamic rheology and electrophoresis measurements are compared to correlate the stability of the suspension with particle-particle and particle-solvent interactions. Although the extent of electrostatic stability is reduced by addition of glycerol, the rheological properties show a transition from a highly flocculated gel to stable dispersions containing no microstructures. This is attributed to a high degree of hydrogen-bonding between glycerol and the Aerosil surface silanol groups. Small dissociation of NaCl and particles reduce the effect of ion exchange and particle bridging mechanisms when the suspensions destabilise in the presence of glycerol. The high viscosity of glycerol is important with respect to the formation of a thick solvation layer around the particles. These parameters give rise to short-range, non-DLVO repulsive solvation forces, which stabilise the dispersion. At intermediate concentrations of glycerol (30–60 wt%) the apparent viscosity increase abruptly and irreversibly as both the extent and time of shearing are increased. The shear rate for the onset of the shear thickening is found to be retarded by decreasing the particle and salt concentration as well as by increasing the glycerol concentration. It is postulated that at intermediate glycerol concentration, where the height of the energy barrier is small, mechanical forces can activate the particles to overcome the energy barrier to enter the region where attractive forces dominate. Here, domination of the hydrodynamic forces to the colloidal forces under the shear results in formation of irreversible gels which does not relax to its initial condition.     

Potential of mean force in confined colloids: integral equations with fundamental measure bridge functions

The potential of mean force for uncharged macroparticles suspended in a fluid confined by a wall or a narrow pore is computed for solvent-wall and solvent-macroparticle interactions with attractive forces. Bridge functions taken from Rosenfeld's density-functional theory are used in the reference hypernetted chain closure of the Ornstein-Zernike integral equations. The quality of this closure is assessed by comparison with simulation. As an illustration, the role of solvation forces is investigated. When the "residual" attractive tails are given a range appropriate to "hard sphere-like" colloids, the unexpected role of solvation forces previously observed in bulk colloids is confirmed in the confinement situation.     

Insight into the C-F bond mechanism of molecular analogs for antibacterial drug design

The activities of biological molecules usually rely on both of intra-molecular and intermolecular interactions between their function groups. These interactions include interonic attraction theory, Van der Waal’s forces and the function of geometry on the individual molecules, whether they are naturally or synthetic. The purpose of this study was to evaluate the antibacterial activity of C-F bond compound using combination of experiments verification and theoretical calculation. We target on the insect natural products from the maggots of Chrysomyis megacephala Fabricius. Based on density functional theory(DFT) and B3LYP method, a theoretical study of the C-F bond on fluoride was designed to explore compounds 2 and 4 antibacterial structure–activity relationship. With the progress in DFT, first-principle calculation based on DFT has gradually become a routine method for drug design, quantum chemistry and other science fields.     

Hybrid density functional theory for homopolymer mixtures confined in a selective nanoslit

By integrating polymer density function theory (DFT) and single-chain molecular simulation, a hybrid DFT is developed for homopolymer mixtures confined in a selective nanoslit. Two weighting functions are adopted separately in the polymer DFT for repulsive and attractive contributions to the excess free energy functional. The theoretical results agree well with simulation data for the density profiles, configurations (tail, loop and train), adsorption amounts, layer thicknesses, and partition coefficients. The polymer-slit interaction is found to have a large effect on the density profiles and partition coefficients but is found to have a small effect on the average sizes and percentages of the configurations. Nearly half of the polymer segments form tails, and the other half form trains. In addition, bridges are observed to form for sufficiently long polymer chains. As the length difference between two polymers increases, the effect of chain connectivity becomes increasingly important.     

Structure, surface excess and effective interactions in polymer nanocomposite melts and concentrated solutions

The Polymer Reference Interaction Site Model (PRISM) theory is employed to investigate structure, effective forces, and thermodynamics in dense polymer-particle mixtures in the one and two particle limit. The influence of particle size, degree of polymerization, and polymer reduced density is established. In the athermal limit, the surface excess is negative implying an entropic dewetting interface. Polymer induced depletion interactions are quantified via the particle-particle pair correlation function and potential of mean force. A transition from (nearly) monotonic decaying, attractive depletion interactions to much stronger repulsive-attractive oscillatory depletion forces occurs at roughly the semidilute-concentrated solution boundary. Under melt conditions, the depletion force is extremely large and attractive at contact, but is proceeded by a high repulsive barrier. For particle diameters larger than roughly five monomer diameters, division of the force by the particle radius results in a nearly universal collapse of the depletion force for all interparticle separations. Molecular dynamics simulations have been employed to determine the depletion force for nanoparticles of a diameter five times the monomer size over a wide range of polymer densities spanning the semidilute, concentrated, and melt regimes. PRISM calculations based on the spatially nonlocal hypernetted chain closure for particle-particle direct correlations capture all the rich features found in the simulations, with quantitative errors for the amplitude of the depletion forces at the level of a factor of 2 or less. The consequences of monomer-particle attractions are briefly explored. Modification of the polymer-particle pair correlations is relatively small, but much larger effects are found for the surface excess including an energetic driven transition to a wetting polymer-particle interface. The particle-particle potential of mean force exhibits multiple qualitatively different behaviors (contact aggregation, steric stabilization, local bridging attraction) depending on the strength and spatial range of the polymer-particle attraction.     

Effective interactions in lysozyme aqueous solutions: a small-angle neutron scattering and computer simulation study

We report protein-protein structure factors of aqueous lysozyme solutions at different pH and ionic strengths, as determined by small-angle neutron scattering experiments. The observed upturn of the structure factor at small wavevectors, as the pH increases, marks a crossover between two different regimes, one dominated by repulsive forces, and another one where attractive interactions become prominent, with the ensuing development of enhanced density fluctuations. In order to rationalize such experimental outcome from a microscopic viewpoint, we have carried out extensive simulations of different coarse-grained models. We have first studied a model in which macromolecules are described as soft spheres interacting through an attractive r(-6) potential, plus embedded pH-dependent discrete charges; we show that the uprise undergone by the structure factor is qualitatively predicted. We have then studied a Derjaguin-Landau-Verwey-Overbeek (DLVO) model, in which only central interactions are advocated; we demonstrate that this model leads to a protein-rich/protein-poor coexistence curve that agrees quite well with the experimental counterpart; experimental correlations are instead reproduced only at low pH and ionic strengths. We have finally investigated a third, "mixed" model in which the central attractive term of the DLVO potential is imported within the distributed-charge approach; it turns out that the different balance of interactions, with a much shorter-range attractive contribution, leads in this latter case to an improved agreement with the experimental crossover. We discuss the relationship between experimental correlations, phase coexistence, and features of effective interactions, as well as possible paths toward a quantitative prediction of structural properties of real lysozyme solutions.     

Solvent density effects on the solvation behavior and configurational structure of bare and passivated 38-atom gold nanoparticle in supercritical ethane

In exploring the effects of solvent density on the mode and the degree of solvation of the bare and passivated 38-atom gold particle in supercritical ethane, we have extended the molecular dynamics simulations of the system, reported previously,(34) to cover a range of isotherms in the T > T(c) regime, where T(c) is the critical temperature of the solvent. Consonant with our previous observations, the modes of solvation of the bare and the passivated particle, deduced from the radial distribution of the solvent about the metal core center of mass, are found to be vastly different from each other at all solvent densities: while the molecules solvating the bare particle form a well-defined, two-region layer around it, those solvating the passivated particle are loosely dispersed in the passivating layer. For the bare particle, the degree of solvation (vartheta) as a function of solvent density passes through a maximum occurring in the close vicinity of the critical point, consistent with our previous results and in agreement with Debenedetti's theoretical analysis,(22,23) which predicts a solvation enhancement effect in the critical region for systems where the unlike solvent/solute interaction is much stronger than the solvent/solvent interaction. Taking the degree of solvation (vartheta) as a measure of solvent quality, we have investigated how the solvent quality would vary along the solvent-density isotherms. In the solvent-density regime rho > rho(c), the solvent quality is found to be a decreasing function of the density as a result of progressive dominance of the excluded volume effect over the attractive particle/solvent interactions. The particle/solvent affinity is greatly reduced in the presence of the passivating layer, resulting in considerable shrinkage of the good-solvent-quality domain in the supercritical regime. The solvent environment and the presence of the passivating chains produce significant disorder in the equilibrium structure assumed by the nanoparticle core.     

Experimental investigation of electrostatic particle-particle interactions in optoelectronic tweezers

This paper reports experimental and theoretical investigation of electrostatic attraction and repulsion of microparticles in an optoelectronic tweezers (OET). When we manipulate dielectric particles suspended in a fluid using OET, the electrostatic interactions of the polarized particles occur, limiting the effective manipulation of microparticles using a light-induced dielectrophoresis. In this study, we first demonstrate the electrostatic particle-particle interactions in the OET device using a liquid crystal display. At the same time, the experimental investigation of the dipole interactions between two spherical particles has been performed using the OET device. On the basis of the point-dipole model, simulation studies on the dipole forces acting on the particles and their trajectories by the forces are also performed. The experimental results show good agreement with the previously reported numerical studies as well as the results of our simulation studies.     

Density functional theory of homopolymer mixtures confined in a slit

A density functional theory (DFT) is developed for polymer mixtures with shorted-ranged attractive interparticle interactions confined in a slit. Different weighting functions are used separately for the repulsive part and the attractive part of the excess free energy functional by applying the weighted density approximation. The predicted results by DFT are in good agreement with the corresponding simulation data indicating the reliability of the theory. Furthermore, the center-of-mass profiles and the end-to-end distance distributions are obtained by the single chain simulation; the predictions also agree well with simulation data. The results reveal that both the attraction of the slit wall and the temperature has stronger effect on longer chains than on shorter ones because the intrasegment correlation of chains increases with increasing chain length.     

Effect of attractions on the structure of polymer solutions confined between surfaces: a density functional approach

A density functional theory is presented to study the effect of attractions on the structure of polymer solutions confined between surfaces. The polymer molecules have been modeled as a pearl necklace of freely jointed hard spheres and the solvent as hard spheres, both having Yukawa-type attractions and the mixture being confined between attractive Yukawa-type surfaces. The present theory treats the ideal gas free energy functional exactly and uses weighted density approximation for the hard chain and hard sphere contributions to the excess free energy functional. The attractive interactions are calculated using the direct correlation function obtained from the polymer reference interaction site model theory along with the mean spherical approximation closure. The theoretical predictions on the density profiles of the polymer and the solvent molecules are found to agree quite well with the Monte Carlo simulation results for varying densities, chain lengths, wall separations, and different sets of interaction potentials.     

Measurement of dynamical forces between deformable drops using the atomic force microscope. I. Theory

Recent experimental developments have enabled the measurement of dynamical forces between two moving liquid drops in solution using an atomic force microscope (AFM). The drop sizes, interfacial tension, and approach velocities used in the experiments are in a regime where surface forces, hydrodynamics, and drop deformation are all significant. A detailed theoretical model of the experimental setup which accounts for surface forces, hydrodynamic interactions, droplet deformation, and AFM cantilever deflection has been developed. In agreement with experimental observations, the calculated force curves show pseudo-constant compliance regions due to drop flattening, as well as attractive pull-off forces due mainly to hydrodynamic lubrication forces.     

Anomalous potentials from inverse analyses of interfacial polydisperse attractive colloidal fluids

This paper investigates effects of using monodisperse inverse analyses to extract particle-particle and particle-surface potentials from simulated interfacial colloidal fluids of polydisperse attractive particles. Effects of polydispersity are investigated as functions of particle concentration and attractive well depth and range for van der Waals and depletion potentials. Forward Monte Carlo simulations are used to generate particle distribution functions for polydisperse interfacial colloidal fluids from which inverted potentials are obtained using an inverse Ornstein-Zernike analysis and an inverse Monte Carlo simulation method. Attractive potentials are successfully recovered for monodisperse colloidal fluids, but polydispersity that is unaccounted for in inverse analyses produces (1) apparent softening of strong forces, (2) anomalous repulsive and attractive interactions, and (3) aphysical particle overlaps. This investigation provides insights into the role of polydispersity in altering the equilibrium structure and corresponding inverted potentials of attractive colloidal fluids near surfaces. These findings should assist the design and interpretation of optical microscopy experiments involving interfacial colloidal fluids similar to the simulated experiments reported here.     

1,3,5-trinitro-1,3,5-triazine decomposition and chemisorption on Al(111) surface: first-principles molecular dynamics study

We have investigated the decomposition and chemisorption of a 1,3,5-trinitro-1,3,5-triazine (RDX) molecule on Al(111) surface using molecular dynamics simulations, in which interatomic forces are computed quantum mechanically in the framework of the density functional theory (DFT). The real-space DFT calculations are based on higher-order finite difference and norm-conserving pseudopotential methods. Strong attractive forces between oxygen and aluminum atoms break N-O and N-N bonds in the RDX and, subsequently, the dissociated oxygen atoms and NO molecules oxidize the Al surface. In addition to these Al surface-assisted decompositions, ring cleavage of the RDX molecule is also observed. These reactions occur spontaneously without potential barriers and result in the attachment of the rest of the RDX molecule to the surface. This opens up the possibility of coating Al nanoparticles with RDX molecules to avoid the detrimental effect of oxidation in high energy density material applications.     

Application of a two-length-scale field theory to the solvation of neutral and charged molecules

We develop a continuous self-consistent theory of solute-water interactions that allows determination of the hydrophobic layer around a solute molecule of any geometry, with an explicit account of solvent structure described by its correlation function. We compute the mean solvent density profile n(r) surrounding the solute molecule as well as its solvation free energy deltaG. We compare the two-length-scale field theory to the numerical data of Monte Carlo simulations found in the literature for spherical molecules and discuss the possibility of self-consistent adjustment of the free parameters of the theory. In the framework of this approach, we compute the solvation free energies of alkane molecules and the free energy of interaction of two spheres of radius R separated by the distance D. We describe the general setting of the self-consistent account of electrostatic interactions in the framework of our model where the water is considered not as a continuous medium but as a gas of dipoles. We analyze the limiting cases where the proposed theory coincides with the electrostatics of a continuous medium.