Molecular dynamics simulation of the aggregation of colloidal particles

The spontaneous time evolution of systems containing N colloidal particles (N = 12, 24, 100) in a spherical cell of volume V at a constant volume fraction φ=0.1 was studied by a molecular dynamics method in the NVT ensemble. The starting velocities of the particles are allocated according to the Maxwell distribution at T=273 K.

Pairwise interaction of the particles was specified by molecular, electrostatic and elastic forces. The changes in the potential energy of the systems were calculated during the establishment of dynamic equilibrium. Coagulation takes place at sufficiently high values of the Hamaker constant. The value of the coefficient of Brownian diffusion, which is calculated from the half-time of coagulation, is found to be close to the known value for aqueous dispersions. The inclusion of electrostatic forces prevents coagulation.

The results obtained are in agreement with those obtained using theories of aggregate formation. Some structural characteristics of aggregates and stable systems are discussed.     

Molecular dynamics simulation of the structural configuration of binary colloidal monolayers

Molecular dynamics simulations of binary colloidal monolayers, i.e., monolayers consisting of mixtures of two different particle sizes, are presented. In the simulations, the colloid particles are located at an oil-water interface and interact via an effective dipole-dipole potential. In particular, the influence of the particle ratio on the configurations of the binary monolayers is investigated for two different relative interaction strengths between the particles, and the pair correlation functions corresponding to the binary monolayers are calculated. The simulations show that the binary monolayers can only form two-dimensional crystals for certain particle ratios, for example, 2:1, 6:1, etc., while, for example, for a particle ratio of 7:1 the monolayers are found to be in a disordered, glassy state. The calculations also reveal that in analogy to the Wigner lattice the configurations are very sensitive to the relative interaction strength between the particles but not to the absolute magnitude of the interaction strength, even when particle size effects are taken into account. Consequently, it is argued that a comparison between the calculated configurations and actual binary particle monolayer systems could provide useful information on the relative interaction strength between large and small particles. Possible mechanisms giving rise to disparities in the interaction strength between large and small particles are described briefly.     

Molecular dynamics simulation of amorphous SiO2 nanoparticles

Molecular dynamics simulation of amorphous SiO2 spherical nanoparticles has been carried out in a model with different sizes, 2, 4, and 6 nm, under non-periodic boundary conditions. We use the pair interatomic potentials which have weak Coulomb interaction and Morse type short-range interaction. Models have been obtained by cooling from the melt via molecular dynamics (MD) simulation. Structural properties of amorphous nanoparticles obtained at 350 K have been studied via partial radial distribution functions (PRDFs), mean interatomic distances, coordination numbers, and bond-angle distributions, which are compared with those observed in the bulk. Calculations of the radial density profile in nanoparticles show the tendency of oxygen to concentrate at the surface as observed previously in other amorphous clusters or thin films. Size effects on structure of nanosized models are significant. The calculations show that if the size is larger than 4 nm, amorphous SiO2 nanoparticles have a distorted tetrahedral network structure with the mean coordination number ZSi-O approximately 4.0 and ZO-Si approximately 2.0 like those observed in the bulk. Surface structure, surface energy, and glass transition temperature of SiO2 nanoparticles have been obtained and presented.     

Molecular dynamics simulation with a continuum electrostatic model of the solvent

The accuracy and simplicity of the Poisson-Boltzmann electrostatics model has led to the suggestion that it might offer an efficient solvent model for use in molecular mechanics calculations on biomolecules. We report a successful merger of the Poisson-Boltzmann and molecular dynamics approaches, with illustrative calculations on the small solutes dichloroethane and alanine dipeptide. The algorithm is implemented within the program UHBD. Computational efficiency is achieved by the use of rather coarse finite difference grids to solve the Poisson-Boltzmann equation. Nonetheless, the conformational distributions generated by the new method agree well with reference distributions obtained as Boltzmann distributions from energies computed with fine finite difference grids. The conformational distributions also agree well with the results of experimental measurements and conformational analyses using more detailed solvent models. We project that when multigrid methods are used to solve the finite difference problem and the algorithm is implemented on a vector supercomputer, the computation of solvent electrostatic forces for a protein of modest size will add only about 0.1 s computer time per simulation step relative to a vacuum calculation. © 1995 by John Wiley & Sons, Inc.     

Interaction forces between colloidal particles in a solution of like-charged, adsorbing nanoparticles

We have measured the force between a weakly charged micron-sized colloidal particle and flat substrate in the presence of highly charged nanoparticles of the same sign under solution conditions such that the nanoparticles physically adsorb to the colloidal particle and substrate. The objective was to investigate the net effect on the force profile between the microparticle and flat substrate arising from both nanoparticle adsorption and nanoparticles in solution. The experiments used colloidal probe atomic force microscopy (CP-AFM) to measure the force profile between a relatively large (5 μm) colloidal probe glass particle and a planar glass substrate in aqueous solutions at varying concentrations of spherical nanoparticles. At very low nanoparticle concentrations, the primary effect was an increase in the electrostatic repulsion between the surfaces due to adsorption of the more highly charged nanoparticles. As the nanoparticle concentration is increased, a depletion attraction formed, followed by longer-range structural forces at the highest nanoparticle concentrations studied. These results suggest that, depending on their concentration, such nanoparticles can either stabilize a dispersion of weakly-charged colloidal particles or induce flocculation. This behavior is qualitatively different from that in nonadsorbing systems, where the initial effect is the development of an attractive depletion force.     

Molecular dynamics simulation of optically trapped colloidal particles at an oil-water interface

Using molecular dynamics simulations, we calculate the net force on a colloidal particle trapped by an optical tweezer and confined within a particle monolayer which is in motion relative to the trapped particle. The calculations are compared with recent experimental data on polystyrene particles located at an oil-water interface. Good agreement between theory and experiment is obtained over the investigated range of lattice constants for an interaction mechanism between the polystyrene particles which is dominated by an effective dipole-dipole potential. The assumed interaction mechanism is consistent with the formation of surface charge dipoles at the particle-oil interface due to the dissociaton of the hydrophilic sulfate headgroups at the surface of the polystyrene particles. A possible physical mechanism for the formation of the surface charge dipoles, involving a diffuse cloud of fully hydrated counterions, is described, and the fraction of surface groups contributing to the formation of surface charge dipoles is estimated to be of the order of 10(-1) for the present system.     

Molecular dynamics simulation of brushes formed by star polyelectrolytes under theta solvent conditions

Using molecular dynamics simulations, we have studied polyelectrolyte brushes formed by partially or fully charged star polymers tethered on a planar surface under theta solvent conditions. The diagram of states for salt‐free solutions differs in the location of osmotic regime (OB) compared with the respective diagram reported by Borisov and Zhulina. In contrast, simulation results dictate that the OB regime appears for values of the ratio F /α ^?1/2 l_B ^?1 much larger than unity, which is the threshold of counterion localization, with F denoting the branch functionality, α the charged units fraction and l_B the Bjerrum length. The simulation results support that the brush height scaling laws H ～ α ² l_B F ^1.049S ³s ^?1 and H ～ α^0.302 F ^0.23S for the charged Pincus Brush (PB) and osmotic (OB) regimes, respectively, where S is the spacer length and s is the grafted area per star chain. The respective theoretical scaling laws are H ～ α ² l_B F ^1.88S ³s ^?1 and H ～ α ^1/2 F ^0.44S . © 2017 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2017 , 55 , 1110–1117     

A coarse-grained explicit solvent simulation of rheology of colloidal suspensions

We use a simple extension of the dissipative particle dynamics (DPD) model to address the dynamical properties of macrosolutes immersed in complex fluid solvents. In this approach, the solvent particles are still represented as DPD particles, thereby retaining the time and length scale advantages offered by the DPD approach. In contrast, the solute particles are represented as hard particles of the appropriate size. We examine the applicability of this simulation approach to reproduce the correct hydrodynamical characteristics of the mixture. Our results focus on the equilibrium dynamics and the steady-state shear rheological behaviors for a range of volume fractions of the suspension, and demonstrate excellent agreement with many published experimental and theoretical results. Moreover, we are also able to track the glass transition of our suspension and the associated dynamical signatures in both the diffusivities and the rheological properties of our suspension. Our results suggest that the simulation approach can be used as a one-parameter model to examine quantitatively the rheological properties of colloidal suspensions in complex fluid solvents such as polymeric melts and solutions, as well as allied dynamical phenomena such as phase ordering in mixtures of block copolymers and particles.     

Molecular dynamics simulation of liquid water between two walls

A molecular dynamics simulation, lasting ≈25 ps, has been performed with 150 ST2 water molecules between two quasi-hard repulsive walls, at a temperature of 302 K. A number of static and dynamic properties have been computed as a function of the distance from the walls, showing that water near the walls is in general more “ordered” than in the bulk, and that this bulk water behaves like ordinarv liquid ST2 water.     

Molecular simulation of interaction between passivated gold nanoparticles in supercritical CO2

Molecular dynamics simulations have been performed to study the potential of mean force (PMF) between passivated gold nanoparticles (NPs) in supercritical CO(2) (scCO(2)). The nanoparticle model consists of a 140 atom gold nanocore and a surface self-assembled monolayer, in which two kinds of fluorinated alkanethiols were considered. The molecular origin of the thermodynamics interaction and the solvation effect has been comprehensively studied. The simulation results demonstrate that increasing the solvent density and ligand length can enhance the repulsive feature of the free energy between the passivated Au nanoparticles in scCO(2), which is in good agreement with previous experimental results. The interaction forces between the two passivated NPs have been decomposed to reveal various contributions to the free energy. It was revealed that the interaction between capping ligands and the interaction between the capping ligands and scCO(2) solvent molecules cooperatively determine the total PMF. A thermodynamic entropy-energy analysis for each PMF contribution was used to explain the density dependence of PMF in scCO(2) fluid. Our simulation study is expected to provide a novel microscopic understanding of the effect of scCO(2) solvent on the interaction between passivated Au nanoparticles, which is helpful to the dispersion and preparation of functional metal nanoparticles in supercritical fluids.     

Probing colloidal forces between a Si3N4 AFM tip and single nanoparticles of silica and alumina

The atomic force microscope (AFM) has been used to measure surface forces between silicon nitride AFM tips and individual nanoparticles deposited on substrates in 10(-4) and 10(-2) M KCl solutions. Silica nanoparticles (10 nm diameter) were deposited on an alumina substrate and alumina particles (5 to 80 nm diameter) were deposited on a mica substrate using aqueous suspensions. Ionic concentrations and pH were used to manage attractive substrate-particle electrostatic forces. The AFM tip was located on deposited nanoparticles using an operator controlled offset to achieve stepwise tip movements. Nanoparticles were found to have a negligible effect on long-range tip-substrate interactions, however, the forces between the tip and nanoparticle were detectable at small separations. Exponentially increasing short-range repulsive forces, attributed to the hydration forces, were observed for silica nanoparticles. The effective range of hydration forces was found to be 2-3 nm with the decay length of 0.8-1.3 nm. These parameters are in a good agreement with the results reported for macroscopic surfaces of silica obtained using the surface force apparatus suggesting that hydration forces for the silica nanoparticles are similar to those for flat silica surfaces. Hydration forces were not observed for either alumina substrates or alumina nanoparticles in both 10(-4) M KCl solution at pH 6.5 and 10(-2) M KCl at pH 10.2. Instead, strong attractive forces between the silicon nitride tip and the alumina (nanoparticles and substrate) were observed.     

Adhesion forces between hybrid colloidal particles and concanavalin A

Hybrid particles of poly(methyl methacrylate) and carboxymethylcellulose, PMMA/CMC, were attached to atomic force microscopy cantilevers and probed against concanavalin A (ConA) films formed either on Si wafers or on CMC substrate. Regardless of the substrate, the approach curves showed different inclinations, indicating that the probe first touches a soft surface and then a hard substrate. The distance corresponding to the soft layer was estimated as 20 +/- 10 nm and was attributed to the CMC layers attached to the hybrid particles surfaces. Probing PMMA/CMC particles against ConA adsorbed onto Si wafers yielded retract curves with a sawlike pattern. The average range of adhesion forces (maximum pull-off distance) and mean adhesion force were estimated as 100 +/- 40 nm and -11 +/- 7 nN, respectively, evidencing multiple adhesions between CMC sugar residues and ConA. However, upon probing against ConA adsorbed onto CMC substrates, the mean pull-off distance and mean adhesion force were reduced to 37 +/- 18 nm and -3 +/- 1 nN, respectively, indicating that the ConA molecules immobilized onto CMC films are less available to interact with the hybrid particle than the ConA molecules adsorbed onto Si wafers. Another set of experiments, where PMMA/CMC particle probed against ConA-covered Si wafers in the presence of mannose, showed that the addition of mannose led to a considerable decrease in the mean adhesion force from -11 +/- 7 to -3 +/- 1 nN. Two hypotheses have been considered to explain the effect caused by mannose addition. The first suggested the desorption of ConA from the substrate so that the hybrid particle would probe bare Si wafer (weak adhesion). The second proposed the adsorption of mannose onto the ConA layer so that mannose layer would probe against another mannose layer, leading to low adhesion forces. In situ ellipsometry and capillary electrophoresis have been applied to check the hypotheses.     

Molecular dynamics simulation on the interaction between polymer inhibitors and anhydrite surface

Investigation on the microscopic interaction between polymer inhibitors and calcium sulfate will be helpful for understanding its scale inhibition mechanism and can provide a theoretical guidance to developing new scale inhibitors. In this work, molecular dynamics simulations with COMPASS force field have been performed to simulate the interaction between hydrolyzed polymaleic anhydride (HPMA), polyaspartic acid (PASP), polyepoxysuccinic acid (PESA), polyacrylic acid (PAA) and the (001) and (020) surfaces of anhydrite (AD) crystal with and without water. The results show that the sequence of binding energies between four polymer inhibitors and AD (001) and (020) with water is PESA > PASP > HPMA > PAA. The binding energy of the same polymer inhibitor on AD (001) is smaller than that on AD (020). Water molecules weaken the deformations of HPMA and PAA but aggravate those of PASP and PESA. Natural bond orbital (NBO) charges of the repeat units of polymer inhibitors were calculated by B3LYP/6‐31G* method. The Coulomb interaction is formed between the O atoms of polymer inhibitors and the Ca atoms of AD crystal. The system of polymer–AD is mainly contributed from the non‐bonding interaction. Polymer inhibitors do not interact directly with AD crystal, but indirectly through the interactions between inhibitor–H₂O and H₂O–AD, i.e. water molecules participate in scale inhibition of polymer inhibitors to AD crystal. Water molecules cannot be ignored when the interaction models are constructed, i.e. solvent effect cannot be ignored. Copyright © 2015 John Wiley & Sons, Ltd.     

金属晶体间微观静摩擦特性的分子动力学模拟研究

以光滑干摩擦接触平面为对象,利用金属晶体间的强体积效应特征,建立了简化计算静摩擦力的界面势能模型.根据第一性原理的方法模拟得出界面分子势能的变化,通过界面分子势能计算出静摩擦力大小,并将数据结果通过通用黏附能量函数计算出的静摩擦力大小进行验证,也将计算结果与超高真空原子力显微镜试验结果进行对比.最后拟合出最大静摩擦力与法向载荷的线性函数关系,得出摩擦力的数值为真实接触面积的函数,并与法向载荷成正比的结论.从微观上对同种金属材料间库伦摩擦定律进行验证与研究.     

Molecular dynamics simulation study of interaction between model rough hydrophobic surfaces

We study some aspects of hydrophobic interaction between molecular rough and flexible model surfaces. The model we use in this work is based on a model we used previously (Eun, C.; Berkowitz, M. L. J. Phys. Chem. B 2009, 113, 13222-13228), when we studied the interaction between model patches of lipid membranes. Our original model consisted of two graphene plates with attached polar headgroups; the plates were immersed in a water bath. The interaction between such plates can be considered as an example of a hydrophilic interaction. In the present work, we modify our previous model by removing the charge from the zwitterionic headgroups. As a result of this procedure, the plate character changes: it becomes hydrophobic. By separating the total interaction (or potential of mean force, PMF) between plates into the direct and the water-mediated interactions, we observe that the latter changes from repulsive to attractive, clearly emphasizing the important role of water as a medium. We also investigate the effect of roughness and flexibility of the headgroups on the interaction between plates and observe that roughness enhances the character of the hydrophobic interaction. The presence of a dewetting transition in a confined space between charge-removed plates confirms that the interaction between plates is strongly hydrophobic. In addition, we notice that there is a shallow local minimum in the PMF in the case of the charge-removed plates. We find that this minimum is associated with the configurational changes that flexible headgroups undergo as the two plates are brought together.     

Molecular dynamics simulation in the grand canonical ensemble

An extended system Hamiltonian is proposed to perform molecular dynamics (MD) simulation in the grand canonical ensemble. The Hamiltonian is similar to the one proposed by Lynch and Pettitt (Lynch and Pettitt, J Chem Phys 1997, 107, 8594), which consists of the kinetic and potential energies for real and fractional particles as well as the kinetic and potential energy terms for material and heat reservoirs interacting with the system. We perform a nonlinear scaling of the potential energy parameters of the fractional particle, as well as its mass to vary the number of particles dynamically. On the basis of the equations of motion derived from this Hamiltonian, an algorithm has been proposed for MD simulation at constant chemical potential. The algorithm has been tested for the ideal gas, for the Lennard-Jones fluid over a wide range of temperatures and densities, and for water. The results for the low-density Lennard-Jones fluid are compared with the predictions from a truncated virial equation of state. In the case of the dense Lennard-Jones fluid and water our predicted results are compared with the results reported using other available methods for the calculation of the chemical potential. The method is also applied to the case of vapor-liquid coexistence point predictions.     

Molecular dynamics simulation of cellobiose in water

The conformational behavior of cellobiose was studied by molecular dynamics simulation in a periodic box of waters. Several different initial conformations were used and the results compared with equivalent vacuum simulations. The average positions and rms fluctuations within single torsional conformations of cellobiose were affected only slightly by the solvent. However, water damped local torsional librations and transitions. The conformational energies of the solute and their fluctuations were also sensitive to the presence of solvent. Intramolecular hydrogen bonding was weakened relative to that observed in vacuo due to competition with solvating waters. All cellobiose hydroxyl groups participated in intermolecular hydrogen bonds with water, with approximately eight hydrogen bonds formed per glucose ring. The hydrogen bonding was predominantly between water hydrogens and solute hydroxyl oxygens. Intermolecular hydrogen bonding to ring and bridge oxygens was seldom present. The diffusion coefficients of both water and solute agree closely with experimental values. Water interchanged rapidly between the solvating first shell and the bulk on the picosecond time scale. © 1993 John Wiley & Sons, Inc.     

Dissipative particle dynamics simulation on a ternary system with nanoparticles, double-hydrophilic block copolymers, and solvent

Dissipative particle dynamics (DPD) simulations are performed to study the aggregation of hydrophobic nanoparticles in the presence of double-hydrophilic block copolymer (DHBC). A single compact spherical nanoparticle aggregate is formed in the absence of DHBC. The response of the aggregate to a continuous increase in the concentration of DHBC has been investigated in detail. We observe the evolvement from single spherical aggregate, through single ellipsoidal aggregate, single platelike aggregate, single long and curly rod, dispersed aggregates, then to hexagonally packed cylinders, and ultimately to ordered lamellar structures upon slow addition of DHBC chains. However, when nanoparticles and DHBCs are added into the system simultaneously at the beginning of simulation, we only obtain single spherical aggregate, dispersed aggregates, hexagonally packed cylinders, and ordered lamellar structures at different concentrations of DHBC. Phase diagrams of structures against concentration of DHBC are presented for these two methods, and the stabilities of structures obtained with the two methods are compared.     

Molecular dynamics simulation study of the water-mediated interaction between zwitterionic and charged surfaces

We calculated the potential of mean force (PMF) for the interaction between a model zwitterionic bilayer and a model charged bilayer. To understand the role of water, we separated the PMF into two components: one due to direct interaction and the other due to water-mediated interaction. In our calculations, we observed that water-mediated interaction is attractive at larger distances and repulsive at shorter. The calculation of the entropic and enthalpic contributions to the solvent-mediated components of the PMF showed that attraction is entropically dominant, while repulsion is dominated by the enthalpy.