Brownian dynamics simulation of the crystallization dynamics of charged colloidal particles

Crystal formation process of charged colloidal particles is investigated using Brownian dynamics (BD) simulations. The particles are assumed to interact with the pair-additive repulsive Yukawa potential. The time evolution of crystallization process and the crystal structure during the simulation are characterized by means of the radial distribution functions (RDF) and mean square displacement (MSD). The simulations show that when the interaction is featured with long-range, particles can spontaneously assemble into body-centered-cubic (BCC) arrays at relatively low particle number density. When the interaction is short-ranged, with increasing the number density particles become trapped into a stagnant disordered configuration before the crystallization could be actualized. The simulations further show that as long as the trapped configurations are bypassed, the face-centered-cubic (FCC) structures can be achieved and are actually more stable than BCC structures. (C) 2010 Elsevier Inc. All rights reserved.

Molecular Dynamics Simulation Of Adhesion Forces In A Dipalmitoylphosphatidylcholine Membrane

Adhesion forces of Dipalmitoylphosphatidylcholine ( DPPC) membrane in the gel phase are investigated by molecular dynamics ( MD) simulation. In the simulations, individual DPPC molecules are pulled out of DPPC membranes with different rates and we get the maximum adhesion forces of DPPC membrane. We find that the maximum adhesion forces increase with pull rate, from about 400 to 700 pN when pull rates are from 0.001 to 0.03 nm/ps. We analyze the relationship between pull rate and adhesion forces of different origins using Brownian dynamics and notice that viscosity of solvent plays an important role in adhesion forces. Then we simulate the motion of a single DPPC molecule in solvent and it elucidates that the maximum drag force is almost linear with respect to the pull rate. We use Stokes' relation to describe the motion of a single DPPC molecule and deduce the effective length of a DPPC molecule. Conformational analyses indicate that the free energy variation of a DPPC molecule inside and outside of the DPPC membrane is an essential part of adhesion energy.

Computer simulation on the collision-sticking dynamics of two colloidal particles in an optical trap

Collisions of a particle pair induced by optical tweezers have been employed to study colloidal stability. In order to deepen insights regarding the collision-sticking dynamics of a particle pair in the optical trap that were observed in experimental approaches at the particle level, the authors carry out a Brownian dynamics simulation. In the simulation, various contributing factors, including the Derjaguin-Landau-Verwey-Overbeek interaction of particles, hydrodynamic interactions, optical trapping forces on the two particles, and the Brownian motion, were all taken into account. The simulation reproduces the tendencies of the accumulated sticking probability during the trapping duration for the trapped particle pair described in our previous study and provides an explanation for why the two entangled particles in the trap experience two different statuses. (c) 2007 American Institute of Physics.

Molecular dynamics simulation of adhesion forces in a dipalmitoylphosphatidylcholine membrane

Adhesion forces of Dipalmitoylphosphatidylcholine ( DPPC) membrane in the gel phase are investigated by molecular dynamics ( MD) simulation. In the simulations, individual DPPC molecules are pulled out of DPPC membranes with different rates and we get the maximum adhesion forces of DPPC membrane. We find that the maximum adhesion forces increase with pull rate, from about 400 to 700 pN when pull rates are from 0.001 to 0.03 nm/ps. We analyze the relationship between pull rate and adhesion forces of different origins using Brownian dynamics and notice that viscosity of solvent plays an important role in adhesion forces. Then we simulate the motion of a single DPPC molecule in solvent and it elucidates that the maximum drag force is almost linear with respect to the pull rate. We use Stokes' relation to describe the motion of a single DPPC molecule and deduce the effective length of a DPPC molecule. Conformational analyses indicate that the free energy variation of a DPPC molecule inside and outside of the DPPC membrane is an essential part of adhesion energy.

Towards an Understanding of the Influence of Sedimentation on Colloidal Aggregation by Peclet Number

The Peclet number is a useful index to estimate the importance of sedimentation as compared to the Brownian motion. However, how to choose the characteristic length scale for the Peclet number evaluation is rather critical because the diffusion length increases as the square root of the time whereas the drifting length is linearly related to time. Our Brownian dynamics simulation shows that the degree of sedimentation influence on the coagulation decreases when the dispersion volume fraction increases. Therefore using a fixed length, such as the diameter of particle, as the characteristic length scale for Peclet number evaluation is not a good choice when dealing with the influence of sedimentation on coagulation. The simulations demonstrated that environmental factors in the coagulation process, such as dispersion volume fraction and size distribution, should be taken into account for more reasonable evaluation of the sedimentation influence.

重力对溶液悬浮粒子聚集过程影响的布朗动力学模拟

采用布朗动力学的计算机模拟方法,研究重力因素对于稀溶液中悬浮粒子聚集过程的影响.通过在计算机模拟程序中加入和排除重力因素的影响,对不同重力条件下的粒子团数量随时间的变化曲线进行对比研究,得到了重力对溶液中的粒子团总数和不同大小的粒子团数量随时间变化的影响规律,可以总结为:在聚集阶段初期,重力不影响粒子的聚集,而在聚集阶段后期,重力加快了粒子的聚集.同时,从动力学分析的角度出发,对重力如何影响悬浮溶液中粒子聚集过程的机制也进行了更加深入的探讨. 

流体动力学作用对粒子聚集影响的布朗动力学模拟

采用布朗动力学模拟方法,研究了流体动力学作用对稀溶液中悬浮粒子聚集过程的影响.模拟中忽略了一个粒子同时与多个粒子碰撞聚集的可能,引入了前人有关两粒子间流体动力学作用影响的研究成果.模拟结果证实了流体动力学的作用在比较大的幅度上减缓了粒子的聚集过程,是导致粒子聚集速率的实验值低于Smoluchowski理论值的重要原因之一.另外,在分别加入和排除重力作用,以及考虑和忽略粒子间流体动力学作用在内的各种条件下模拟了粒子的聚集过程,得到了两种因素相互耦合作用时各因素对粒子聚集过程影响的结果,并从动力学的角度对这些因素的影响机制进行了相应的讨论.

Molecular dynamics simulation of the role of dislocations in microcrack healing

The molecular dynamics method is used to simulate microcrack healing during heating or/and under compressive stress. A centre microcrack in Cu crystal would be sealed under compressive stress or by heating. The role of compressive stress and heating in crack healing was additive. During microcrack healing, dislocation generation and motion occurred. When there were pre-existing dislocations around the microcrack, the critical temperature or compressive stress necessary for microcrack healing would decrease, and, the higher the number of dislocations, the lower the critical temperature or compressive stress. The critical temperature necessary for microcrack healing depended upon the orientation of the crack plane. For example, the critical temperature for the crack along the (001) plane was the lowest, i.e. 770K.

Dislocation dynamics of strain relaxation in epitaxial layers

Many experimental observations have clearly shown that dislocation interaction plays a crucial role in the kinetics of strain relaxation in epitaxial thin films. A set of evolution equations are presented in this article. The key feature of the equations

Formation of Fivefold Deformation Twins in Nanocrystalline Face-Centered-Cubic Copper Based on Molecular Dynamics Simulations

Fivefold deformation twins were reported recently to be observed in the experiment of the nanocrystalline face-centered-cubic metals and alloys. However, they were not predicted previously based on the molecular dynamics (MD) simulations and the reason was thought to be a uniaxial tension considered in the simulations. In the present investigation, through introducing pretwins in grain regions, using the MD simulations, the authors predict out the fivefold deformation twins in the grain regions of the nanocrystal grain cell, which undergoes a uniaxial tension. It is shown in their simulation results that series of Shockley partial dislocations emitted from grain boundaries provide sequential twining mechanism, which results in fivefold deformation twins. (c) 2006 American Institute of Physics.

Molecular Dynamics Simulation of Peeling a DNA Molecule on Substrate

Molecular dynamics (MD) simulations are performed to study adhesion and peeling of a short fragment of single strand DNA (ssDNA) molecule from a graphite surface. The critical peel-off force is found to depend on both the peeling angle and the elasticity of ssDNA. For the short ssDNA strand under investigation, we show that the simulation results can be explained by a continuum model of an adhesive elastic band on substrate. The analysis suggests that it is often the peak value, rather than the mean value, of adhesion energy which determines the peeling of a nanoscale material.

Correlative reference model and molecular dynamics simulation of dislocation emission process

A correlative reference model for computer molecular dynamics simulations is proposed. Based on this model, a flexible displacement boundary scheme is introduced and the dislocations emitted from a crack tip can continuously pass through the border of the inner discrete atomic region and pile up at the outer continuum region. The effect of the emitted dislocations within the plastic zone on the inner atomistic region can be clearly demonstrated. The simulations for a molybdinum crystal show that a full dislocation in a bcc crystal is dissociated into three partial dislocations and interaction between the crack and the emitted dislocations results in gradual decrease of the local stress intensity factor.

Molecular dynamics simulation of hydrogen enhancing dislocation emission

The interactive pair potential between Al and H is obtained based on the ab initio calculation and the Chen-Mobius 3D lattice inversion formula. By utilizing the pair potentials calculated, the effects of hydrogen on the dislocation emission from crack tip have been studied. The simulated result shows that hydrogen can reduce the cohesive strength for Al single crystal, and then the critical stress intensity factor for partial dislocation emission decreases from 0.11 MPa root m (C-H = 0) to 0.075 MPa root m (C-H=0.72%) and 0.06 MPa root m (C-H = 1.44%). This indicates thar hydrogen can enhance the dislocation emission. The simulation also shows that atoms of hydrogen can gather and turn into small bubbles, resulting in enhancement of the equilibrium vacancy concentration.

Get Acquaited with Rarefied Gas Dynamics

以通俗易懂的方式介绍了空气动力学当气体间断分子效应显著时发展起来的特殊分支--稀薄气体动力学。讨论了非平衡现象与稀薄气体动力学的关系。通过与8速度气体模型的间断Boltzmann方程的对比,解释了Boltzmann方程碰撞项的物理意义和数学困难,简要综述了其一般解法。讨论了分子在物体表面的反射和问题的边界条件,着重介绍了直接模拟Monte Carlo(DSMC)方法和为克服低速稀薄流动(如MEMS中流动)中模拟困难的信息保存(IP)方法。

Molecular Dynamics Modeling of Diffusion Bonding

Molecular dynamics simulations on diffusion bonding of Cu-Ag showed that the thickness of the interfacial region depended on the stress. The interfacial region became amorphous during diffusion bonding, and it would normally transform from amorphous into crystalline structure when the structure was cooled to the room temperature.