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.

QM/MM and Classical Molecular Dynamics Simulation of Histidine-Tagged Peptide Immobilization on Nickel Surface

The hybrid quantum mechanics (QM) and molecular mechanics (MM) method is employed to simulate the His-tagged peptide adsorption to ionized region of nickel surface. Based on the previous experiments, the peptide interaction with one Ni ion is considered. In the QM/MM calculation, the imidazoles on the side chain of the peptide and the metal ion with several neighboring water molecules are treated as QM part calculated by “GAMESS”, and the rest atoms are treated as MM part calculated by “TINKER”. The integrated molecular orbital/molecular mechanics (IMOMM) method is used to deal with theQMpart with the transitional metal. By using the QM/MM method, we optimize the structure of the synthetic peptide chelating with a Ni ion. Different chelate structures are considered. The geometry parameters of the QM subsystem we obtained by QM/MM calculation are consistent with the available experimental results. We also perform a classical molecular dynamics (MD) simulation with the experimental parameters for the synthetic peptide adsorption on a neutral Ni(1 0 0) surface. We find that half of the His-tags are almost parallel with the substrate, which enhance the binding strength. Peeling of the peptide from the Ni substrate is simulated in the aqueous solvent and in vacuum, respectively. The critical peeling forces in the two environments are obtained. The results show that the imidazole rings are attached to the substrate more tightly than other bases in this peptide.

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.

QM/MM and Classical Molecular Dynamics Simulation of Histidine-Tagged Peptide Immobilization on Nickel Surface

The hybrid quantum mechanics (QM) and molecular mechanics (MM) method is employed to simulate the His-tagged peptide adsorption to ionized region of nickel surface. Based on the previous experiments, the peptide interaction with one Ni ion is considered. In the QM/MM calculation, the imidazoles on the side chain of the peptide and the metal ion with several neighboring water molecules are treated as QM part calculated by "GAMESS", and the rest atoms are treated as MM part calculated by "TINKER". The integrated molecular orbital/molecular mechanics (IMOMM) method is used to deal with the QM part with the transitional metal. By using the QM/MM method, we optimize the structure of the synthetic peptide chelating with a Ni ion. Different chelate structures are considered. The geometry parameters of the QM subsystem we obtained by QM/MM calculation are consistent with the available experimental results. We also perform a classical molecular dynamics (MD) simulation with the experimental parameters for the synthetic peptide adsorption on a neutral Ni(100) surface. We find that half of the His-tags are almost parallel with the substrate, which enhance the binding strength. Peeling of the peptide from the Ni substrate is simulated in the aqueous solvent and in vacuum, respectively. The critical peeling forces in the two environments are obtained. The results show that the in-tidazole rings are attached to the substrate more tightly than other bases in this peptide.

Molecular Dynamics Simulation of the Bio-Adhesion in Molecular Motors

Molecular dynamics (MD) simulation is employed to study the bio-adhesion in F1 ATP molecular motor. Histidine-peptide is widely used as linkage in micro systems because of its strong binding strength to metals. This paper focuses on the adhesion between a synthetic peptide containing 6xHis-tag (Gly-Gly-Lys-Gly-Gly-Lys-Gly-Gly-His-His-His-His-His-His) and metal substrate, which is used to define the position of the F1 ATP molecular motor on the metal substrate. It is shown that the binding strength between histidine and nickel substrate is the strongest, while that of copper is smaller and that of gold substrate is the smallest. From the result of simulation, we find that the stability of adhesion between histidine and the metal substate result of the ringed structure in histidine.

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.

Characterize dynamic conformational space of human CCR5 extracellular domain by molecular modeling and molecular dynamics simulation

The chemokine receptor CCR5 is the receptor for several chemokines and major coreceptor for R5 human immunodeficiency virus type-1 strains entry into cell. Three-dimensional models of CCR5 were built by using homology modeling approach and 1 ns molecular dynamics (MD) simulation, because studies of site-directed mutagenesis and chimeric receptors have indicated that the N-terminus (Nt) and extracellular loops (ECLs) of CCR5 are important for ligands binding and viral fusion and entry, special attention was focused on disulfide bond function, conformational flexibility, hydrogen bonding, electrostatic interactions, and solvent-accessible surface area of Nt and ECLs of this protein part. We found that the extracellular segments of CCR5 formed a well-packet globular domain with complex interactions occurred between them in a majority of time of MID simulation, but Nt region could protrude from this domain sometimes. The disulfide bond Cys20-Cys269 is essential in controlling specific orientation of Nt region and maintaining conformational integrity of extracellular domain. RMS comparison analysis between conformers revealed the ECL1 of CCR5 stays relative rigid, whereas the ECL2 and Nt are rather flexible. Solvent-accessible surface area calculations indicated that the charged residues within Nt and ECL2 are often exposed to solvent. Integrating these results with available experimental data, a two-step gp120-CCR5 binding mechanism was proposed. The dynamic interaction of CCR5 extracellular domain with gp120 was emphasized. (C) 2004 Elsevier B.V. All rights reserved.

Numerical Simulation of the Flow Field and Concentration Distribution in the Bulk Growth of Silicon Carbide Crystals

The physical vapor transport (PVT) method is being widely used to grow large-size single SiC crystals. The growth process is associated with heat and mass transport in the growth chamber, chemical reactions among multiple species as well as phase change at the crystal/gas interface. The current paper aims at studying and verifying the transport mechanism and growth kinetics model by demonstrating the flow field and species concentration distribution in the growth system. We have developed a coupled model, which takes into account the mass transport and growth kinetics. Numerical simulation is carried out by employing an in-house developed software based on finite volume method. The results calculated are in good agreement with the experimental observation.

The effect of thermal-activation on dislocation processes at an atomistic crack-tip

The effects of thermal activation on the dislocation emission from an atomistic crack tip are discussed, Molecular dynamics simulations at different constant temperatures are carried out to investigate the thermal effects. The simulated results show that the processes of the partial dislocation generation and emission are temperature dependent. As the temperature increases, the incipient duration of the partial dislocation nucleation becomes longer, the critical stress intensity factor for partial dislocation emission is reduced and, at the same loading level, more dislocations are emitted. The dislocation velocity moving away from the crack tip and the separations of partial dislocations are apparently not temperature dependent. The simulated results also show that, as the temperature increases, the stress distribution along the crack increases slightly. Therefore stress softening at the crack tip induced by thermal activation does not exist in the present simulation. A simple model is proposed to evaluate the relation of the critical stress intensity factor versus temperature. The obtained relation is in good agreement with our molecular dynamics results.

Molecular dynamics studies of brittle fracture of SiO2

By applying for molecular dynamics (MD) simulation and Griffith fracture criterion, the brittle behavior of crack extension of mode I type is investigated. The critical stress intensity factor (SIF)K-Ic(MD) of crack extension is calculated, and the evolution of atoms near crack tip is observed. It is found that K-Ic(MD) is in good agreement with the Griffith ftacture criterion K-Ic(Griffith).

Effect Of Water On Brittle Fracture Of Sio2 By Molecular Dynamics Study

The influence of water on the brittle behavior of beta-cristobalite is studied by means of molecular dynamics (MD) simulation With the TTAM potential. Crack extension of mode 1 type is observed as the crack opening is filled LIP With water. The critical stress intensity factor K-lc(MD) is used to characterize the crack extension of MD simulation. The surface energy of SiO2 covered with layers of water is calculated at temperature of 300 K. Based oil the Griffith fracture criterion, the critical stress intensity factor K-lc(Griffith) is calculated, and it is in good agreement with that of MD simulation. (C) 2008 Elsevier B.V. All rights reserved.

Temporal characteristic simulation for stimulated emission depletion microscopy

The imaging technology of stimulated emission depletion (STED) utilizes the nonlinearity relationship between the fluorescence saturation and the excited state stimulated depletion. It implements three-dimensional (3D) imaging and breaks the diffraction barrier of far-field light microscopy by restricting fluorescent molecules at a sub-diffraction spot. In order to improve the resolution which attained by this technology, the computer simulation on temporal behavior of population probabilities of the sample was made in this paper, and the optimized parameters such as intensity, duration and delay time of the STED pulse were given.

Motion simulation of annular polishing

Based on the Coulomb friction model, the frictional motion model of workpiece relating to the polishing pad was presented in annular polishing. By the dynamic analysis software, the model was simulated and analysed. The conclusions from the results were that the workpiece did not rotate steadily. When the angular velocity of ring and the direction were the same as that of the polishing pad, the angular velocity of workpiece hoicked at the beginning and at the later stage were the same as that of the polishing pad before contacting with the ring. The angular velocity of workpiece vibrated at the moment of contacting with the ring. After that the angular velocity of workpiece increased gradually and fluctuated at a given value, while the angular velocity of ring decreased gradually and also fluctuated at a given value. Since the contact between the workpiece and the ring was linear, their linear velocities and directions should be the same. But the angular velocity of workpiece was larger than that of the polishing pad on the condition that the radius of the workpiece was less than that of the ring. This did not agree with the pure translation principle and the workpiece surface could not be flat, either. Consequently, it needed to be controlled with the angular velocity of ring and the radii of the ring and the workpiece, besides friction to make the angular velocity of workpiece equal to that of the polishing pad for obtaining fine surface flatness of the workpiece. Copyright © 2007 Inderscience Enterprises Ltd.}

"神光-Ⅱ"装置第九路靶场终端光学组件的研制

针对“神光-Ⅱ”装置第九路系统主激光瞄准精度小于等于30μm和大焦斑辐照均匀性优于10％的要求，提出了靶场终端光学组件的设计结构。应用有限元法对组件关键机械元件和ICF靶室整体进行动静态分析，优化了设计参数。同时与聚焦透镜配合进行数值分析列阵透镜，确定了单元数、曲率和厚度以及单元长和宽等参数。经过实验测试，主激光瞄准精度达到28．9μm，大焦斑辐照的形状为1000μm×500μm，均匀性为12．0％。

Numerical simulation of nc-Si : H/c-Si heterojunction solar cells

AMPS simulator, which was developed by Pennsylvania State University, has been used to simulate photovoltaic performances of nc-Si:H/c-Si solar cells. It is shown that interface states are essential factors prominently influencing open circuit voltages (V-OC) and fill factors (FF) of these structured solar cells. Short circuit current density (J(SC)) or spectral response seems more sensitive to the thickness of intrinsic a-Si:H buffer layers inserted into n(+)-nc-Si:H layer and p-c-Si substrates. Impacts of bandgap offset on solar cell performances have also been analyzed. As DeltaE(C) increases, degradation of VOC and FF owing to interface states are dramatically recovered. This implies that the interface state cannot merely be regarded as carrier recombination centres, and impacts of interfacial layer on devices need further investigation. Theoretical maximum efficiency of up to 31.17% (AM1.5,100mW/cm(2), 0.40-1.1mum) has been obtained with BSF structure, idealized light-trapping effect(R-F=0, R-B=1) and no interface states.