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.

Molecular Dynamics Simulation of the Uniaxial Tensile Deformation of Nanocrystalline Copper

纳米材料是由尺度在1～100 nm的微小颗粒组成的体系,由于它具有独特的性能而备受关注.本文简要地回顾了分子动力学在纳米材料研究中的应用,并运用它模拟了平均晶粒尺寸从1.79～5.38nm的纳米晶体的力学性质.模拟结果显示:随着晶粒尺寸的减小,系统与晶粒内部的原子平均能量升高,而晶界上则有所下降;纳米晶体的弹性模量要小于普通多晶体,并随着晶粒尺寸的减小而减小;纳米晶铜的强度随着晶粒的减小而减小,显示了反常的Hall-Petch效应;纳米晶体的塑性变形主要是通过晶界滑移与运动,以及晶粒的转动来实现的;位错运动起着次要的、有限的作用;在较大的应变下(约大于5%),位错运动开始起作用;这种作用随着晶粒尺寸的增加而愈加明显.

The Influence of Grain Size and Temperature on the mechanical Deformation of Nanocrystalline Materials: Molecular Dynamics Simulation

Nanocrystalline (nc) materials are characterized by a typical grain size of 1-100nm. The uniaxial tensile deformation of computer-generated nc samples, with several average grain sizes ranging from 5.38 to 1.79nm, is simulated by using molecular dynamics with the Finnis-Sinclair potential. The influence of grain size and temperature on the mechanical deformation is studied in this paper. The simulated nc samples show a reverse Hall-Petch effect. Grain boundary sliding and motion, as well as grain rotation are mainly responsible for the plastic deformation. At low temperatures, partial dislocation activities play a minor role during the deformation. This role begins to occur at the strain of 5%, and is progressively remarkable with increasing average grain size. However, at elevated temperatures no dislocation activity is detected, and the diffusion of grain boundaries may come into play.

Molecular dynamics simulation on dislocation emission process

A correlative reference model for a computer simulation of molecular dynamics is proposed in this paper. Based on this model, a flexible displacement boundary scheme is naturally introduced and the dislocations emitted from a crack tip are presumed to continuously pass through the border of an inner discrete atomic region to pile up at an outer continuum region. The simulations for a Mo crystal show that the interaction between a crack and emitted dislocations results in the decrease in local stress intensity factor gradually.

Molecular Dynamics Simulation of Microstructure of Nanocrystalline Copper

The microstructure of computer generated nanocrystalline coppers is simulated by using molecular dynamics with the Finnis-Sinclair potential, analysed by means of radial distribution functions, coordination number, atomic energy and local crystalline order. The influence of the grain size on the nanocrystalline structure is studied. The results reveal that as the grain size is reduced, the grain boundary shows no significant structural difference, but the grain interior becomes more disordered, and their structural difference diminishes gradually; however, the density and the atomic average energy of the grain boundary present different tendencies from those of the grain interior.

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.

Forced Dissociation of Selectin-ligand Complexes Using Steered Molecular Dynamics Simulation

Selectin-ligand interactions are crucial to such biological processes as inflammatory cascade or tumor metastasis. How transient formation and dissociation of selectin-ligand bonds in blood flow are coupled to molecular conformation at atomic level, however, has not been well understood. In this study, steered molecular dynamics (SMD) simulations were used to elucidate the intramolecular and intermolecular conformational evolutions involved in forced dissociation of three selectin-ligand systems: the construct consisting of P-selectin lectin (Lec) and epidermal growth factor (EGF)-like domains (P-LE) interacting with synthesized sulfoglycopeptide or SGP-3, P-LE with sialyl Lewis X (sLeX), and E-LE with sLeX. SMD simulations were based on newly built-up force field parameters including carbohydrate units and sulfated tyrosine(s) using an analogy approach. The simulations demonstrated that the complex dissociation was coupled to the molecular extension. While the intramolecular unraveling in P-LESGP-3 system mainly resulted from the destroy of the two anti-parallel sheets of EGF domain and the breakage of hydrogen-bond cluster at the Lec-EGF interface, the intermolecular dissociation was mainly determined by separation of fucose (FUC) from Ca2+ ion in all three systems. Conformational changes during forced dissociations depended on pulling velocities and forces, as well as on how the force was applied. This work provides an insight into better understanding of conformational changes and adhesive functionality of selectin-ligand interactions under external forces.

Cross-relaxation dynamics on anomalous saturation processes in low pressure supersonic cw HF chemical laser amplifier

It is assumed that both translational and rotational nonequilibrium cross-relaxations play a role simultaneoulsy in low pressure supersonic cw HF chemical laser amplifier. For two-type models of gas flow medium with laminar and turbulent flow diffusion mixing, the expressions of saturated gain spectrum are derived respectively, and the numerical calculations are performed as well. The numerical results show that turbulent flow diffusion mixing model is in the best agreement with the experimental result.

The damage evaluation of rigid particle filled polymer

In this paper, the rigid particle filled polymer is studied in the hope to understand the real damage mechanisms. Two damage parameters were introduced and measured. One is the macro-damage of the materials calculated from the modulus measured, another is micro-damage describing the interfacial debonding or the percentage of the particle debonded from the matrix. The damage rate of the macro damage decreases, while the micro damage increases with the applied stress.

Molecular Dynamics Simulation of Thermal and Mechanical Properties of Polyimide-Carbon-Nanotube Composites

An aromatic polyimide and its mixture with randomly distributed carbon nanotubes (NTs) are simulated by using molecular dynamics, repeated energy minimization and cooling processes. The glass transition temperatures are identified through volume-temperature curves. Stress-strain curves, Young's moduli, densities and Poisson ratios are computed at different temperatures. It is demonstrated that the carbon NT reduces the softening effects of temperature on mechanical properties and increases the ability to resist deformation.

Molecular Dynamics Simulation of Barnacle Cement

Barnacle cement is an underwater adhesive that is used for permanent settlement. Its main components are insoluble protein complexes that have not been fully studied. In present article, we chose two proteins of barnacle cement for study, 36-KD protein and Mrcp-100K protein. In order to investigate the characteristic of above two proteins, we introduced the method of molecular modeling. And the simulation package GROMACS was used to simulate the behavior of these proteins. In this article, before the simulations, we introduce some theories to predict the time scale for polymer relaxation. During the simulation, we mainly focus on two properties of these two proteins: structural stability and adhesive force to substrate. First, we simulate the structural stability of two proteins in water, and then the stability of 36-KD protein in seawater environment is investigated.We find that the stability varies in the different environments. Next, to study adhesive ability of two proteins, we simulate the process of peeling the two proteins from the substrate (graphite). Then, we analyze the main reasons of these results. We find that hydrogen bonds in proteins play an important role in the protein stability. In the process of the peeling, we use Lennard–Jones 12-6 potential to calculate the van der Waals interactions between proteins and substrate.

Dynamics of collapsing fluidized beds and its application in the simulation of pulsed fluidized beds

A new method to study collapsing process of fluidized bed was proposed. The method is based on the analysis of the pressure variation during collapse. A model is proposed to describe the pressure variation on any location of the bed during collapse. Three kind of particles were fluidized by air and the pressure variation after shutting down the gas supply was measured by pressure transducers and simulated by the proposed model. The simulated results were in good agreement with the experimental data. The parameters of the fluidized bed, such as bubble fraction alpha(b), particle fraction in dense phase alpha(p), bubble velocity u(b) and slip velocity of particle u(i), can be estimated from the measured pressure data by using the proposed model. The model was also used to calculate the pressure variation in pulsed fluidized beds with low frequency.

Forced extension of P-selectin construct using steered molecular dynamics

P-selectin, a 70-nm-long cellular adhesive molecule, possesses elastic and extensible properties when neutrophils roll over the activated endotheliam of blood vessel in inflammatory reaction. Transient formation and dissociation of P-selectin/ligand bond on applied force of blood flow induces the extension of P-selectin and relevant ligands. Steered molecular dynamics simulations were performed to stretch a single P-selectin construct consisting of a lectin (Lec) domain and an epithelial growth factor (EGF)-like domain, where P-selectin construct was forced to extend in water with pulling velocities of 0.005-0.05 nm/ps and with constant forces of 1000-2500 pN respectively. Resulting force-extension profiles exhibited a dual-peak pattern on various velocities, while both plateaus and shoulders appeared in the extension-time profiles on various forces. The force or extension profiles along stretching pathways were correlated to the conformational changes, suggesting that the structural collapses of P-selectin Lec/EGF domains were mainly attributed to the burst of hydrogen bonds within the major beta sheet of EGF domain and the disruptions of two hydrophobic cores of Lee domain. This work furthers the understanding of forced dissociation of P-selectin/ligand bond.

A Molecular Dynamics Simulation: the Effect of Finite Size on the Thermal Conductivity in a Single Crystal Silicon

Non-equilibrium molecular dynamics (NEMD) simulations are performed to calculate thermal conductivity. The environment-dependent interatomic potential (EDIP) potential on crystal silicon is adopted as a model system. The issues are related to nonlinear response, local thermal equilibrium and statistical averaging. The simulation results by non-equilibrium molecular dynamics show that the calculated thermal conductivity decreases almost linearly as the film thickness reduced at the nanometre scale. The effect of size on the thermal conductivity is also obtained by a theoretic analysis of the kinetic theory and formulas of the heat capacity. The analysis reveals that the contributions of phonon mean free path (MFP) and phonon number in a finite cell to thermal conductivity are very important.