Effect of twin boundary spacing on deformation behavior of nanotwinned magnesium

The effect of twin spacing and temperature on the deformation behavior of nanotwinned magnesium is investigated using molecular dynamics simulation. The results indicate that there is a pronounced shift in the mechanical behavior of nanotwinned magnesium when twin spacing is smaller than 2.9 nm, and that the yield strength decreases with increasing temperature. The results show that at relatively high temperatures, a strength softening can be observed when twin spacing is larger than 7.8 nm. This study demonstrates that the yield strength is associated with the dislocation storage ability of nanotwinned magnesium and the repulsive force between twin boundaries and dislocations.     

Effect of amorphous lamella on the crack propagation behavior of crystalline Mg/amorphous Mg–Al nanocomposites

The effects of amorphous lamella on the crack propagation behavior in crystalline/amorphous(C/A) Mg/Mg–Al nanocomposites under tensile loading are investigated using the molecular dynamics simulation method. The sample with an initial crack of orientation(1210) [0001] is considered here. For the nano-monocrystal Mg, the crack growth exhibits brittle cleavage. However, for the C/A Mg/Mg–Al nanocomposites, the ‘double hump' behavior can be observed in all the stress–strain curves regardless of the amorphous lamella thickness. The results indicate that the amorphous lamella plays a critical role in the crack deformation, and it can effectively resist the crack propagation. The above mentioned crack deformation behaviors are also disclosed and analyzed in the present work. The results here provide a strategy for designing the high-performance hexagonal-close-packed metal and alloy materials.     

界面结构对Cu/Ni多层膜纳米压痕特性影响的分子动力学模拟

金属多层膜调制周期下降到纳米级时,其力学性质会发生显著改变. Cu-Ni晶格失配度约为2.7%,可以形成共格界面和半共格界面,实验中实现沿[111]方向生长的调制周期为几纳米且具有异孪晶界面结构的Cu/Ni多层膜,其力学性质发生显著改变.本文采用分子动力学方法对共格界面、共格孪晶界面、半共格界面、半共格孪晶界面等四种不同界面结构的Cu/Ni多层膜进行纳米压痕模拟,研究压痕过程中不同界面结构类型的形变演化规律以及位错与界面的相互作用,获取Cu/Ni多层膜不同界面结构对其力学性能的影响特征.计算结果表明,不同界面结构的样品在不同压痕深度时表现出的强化或软化作用机理不同,软化机制主要是由于形成了平行于界面的分位错以及孪晶界面的迁移,强化机制主要是由于界面对位错的限定作用以及失配位错网状结构与孪晶界面迁移时所形成的弓形位错之间的相互作用.     

液态锂在铜的微通道中的流动行为

本文采用分子动力学方法模拟了液态锂在铜的微通道内的流动行为. 通过构建铜(111), (100)和(110)晶面的微通道内壁, 研究了液态锂在流固界面上的微观结构以及在铜微通道中的流动速度分布情况, 并探讨了微通道尺寸对液态锂流动行为的影响. 研究结果表明铜微通道内的液态锂在靠近铜固体壁附近区域呈有序的层状结构分布, 并受铜内壁晶面微观结构的影响. 铜(111)和(100)面内壁附近的液态锂有序层分布结构更明显. 外驱力作用下的液态锂在微通道内的流动速度呈抛物线分布, 流固界面和流动方向对液态锂的流动速度都会产生影响. 液态锂在铜(111)面内壁上流动的速度最大, 且出现了速度滑移; 在铜(110)面内壁上流动速度最小. 通过对不同尺寸的微通道内液态锂流动行为的研究, 发现流动速度的大小随着微通道尺寸的增加而增大, 且最大速度与微通道尺寸呈二次函数关系, 与有关理论计算结果符合得很好.     

Effect of sp-hybridized defects on the oscillatory behavior of carbon nanotube oscillators

Using molecular dynamics simulations, we investigate the oscillatory behaviors of carbon nanotube oscillators containing sp³-hybridized defects formed by hydrogen chemisorption. It is found that the presence of these defects significantly affects the kinetic and potential energies of the nanotube systems, which in turn affects their oscillation periods and frequencies. We have also studied the oscillatory characteristics of the oscillators containing sp³-hybridized Stone-Wales defects. Our results show that it is possible to control the motion of the inner nanotube by introducing sp³-hybridized defects on the outer nanotube, which provides a potential way to tune the oscillatory behavior of nanotube oscillators.     

Molecular dynamic simulation of amorphous carbon and graphite interface

Amorphous carbon/graphite interface is modeled by molecular dynamic simulation using a Tersoff-type potential function with the Brenner parameters for in-plane interaction combined with the pair potential function for the interplanar bonding. The interface is created by compressing the amorphous carbon produced in a separate simulation with perfect crystalline graphite terminated to expose (1120) planes. The planar structure and weak interplanar bonding allow the graphitic planes to deform in order to accommodate the bonds formed at the interface, which is consistent with the HRTEM study of the interface. The simulation indicates that the generated interface mostly consists of nearly sp₂ hybridized bonding connecting the two sides. The bonds across the interface when formed are likely to maintain their equilibrium configurations. Due to the large interplanar spacing, many atoms both on the graphite and a-C sides are left unbonded leaving the interface energetically unfavorable with respect to the bulk. These unbonded radicals probably weaken the structural rigidity of the interface providing a fracture path under stress.     

Molecular dynamics study of plastic deformation mechanism in Cu/Ag multilayers

The plastic deformation mechanism of Cu/Ag multilayers is investigated by molecular dynamics(MD) simulation in a nanoindentation process. The result shows that due to the interface barrier, the dislocations pile-up at the interface and then the plastic deformation of the Ag matrix occurs due to the nucleation and emission of dislocations from the interface and the dislocation propagation through the interface. In addition, it is found that the incipient plastic deformation of Cu/Ag multilayers is postponed, compared with that of bulk single-crystal Cu. The plastic deformation of Cu/Ag multilayers is affected by the lattice mismatch more than by the difference in stacking fault energy(SFE) between Cu and Ag. The dislocation pile-up at the interface is determined by the obstruction of the mismatch dislocation network and the attraction of the image force. Furthermore, this work provides a basis for further understanding and tailoring metal multilayers with good mechanical properties, which may facilitate the design and development of multilayer materials with low cost production strategies.     

Effect of the amorphous thin layer on the surface growth of amorphous/crystalline binary multilayer films

The effect of the amorphous thin layer on the surface growth of amorphous/crystalline binary multilayer films has been studied by using a continuum model. It is shown that both the surface roughness and the growth exponent of amorphous/crystalline binary multilayer films decrease with increasing thickness ratio between amorphous and crystalline layers. Our simulations have also revealed, in contrast to the monotonous rise in surface roughness observed in single-layer films grown on flat substrates, the surface growth of a multilayer film consists of two processes: interface smoothing and roughening, namely the film roughness decreases during the growth of amorphous thin layers but increases monotonously during the growth of crystalline thin layers. The observed interface smoothing and roughening can be obviously influenced by the change in the thickness ratio between amorphous and crystalline layers. The rise in thickness ratio between amorphous and crystalline layers enhances the interface smoothing effect but lowers the interface roughening effect and consequently shows a marked smoothing effect on the surface roughness.     

梯度界面层对复合材料有效光学性质的影响

 研究了球形颗粒被无规地浸没在基质中,具有梯度分界面和光滑分界面的球形颗粒复合物界面层对介电响应和光学性质的影响;利用第一原理近似,计算了在稀释极限下具有梯度分布界面层复合物的有效介电常数,讨论了其实部和吸收光谱随外加场频率的变化关系。对结果进行比较,发现梯度分布的界面层使复合材料在等离子体共振频率处的有效介电常数的涨落明显减小,吸收谱线明显宽化并伴有“蓝移”现象,梯度界面层能很好地控制等离子体共振和减小有效介电响应的涨落。     

Amphiphilic peptide binding on crystalline vs. amorphous silica from molecular dynamics simulations

The leucine-lysine amphiphilic peptide LKα14 has been used to study fundamental driving forces in processes such as peptide-surface binding and biomineralization. Here, we employ molecular dynamics (MD) simulations in tandem with replica exchange metadynamics to probe the binding mechanism and thermodynamics of LKα14 on silica. We also investigate the effect that the nature of the silica surface – crystalline vs. amorphous, has on the binding properties and peptide-surface conformations. We find that water adsorbs differently on both surfaces; it forms a denser interfacial layer on the crystalline surface, compared to the amorphous surface. This causes the peptide to bind more strongly on the amorphous surface than the crystalline surface. Cluster analysis shows that the peptide adopts a helical conformation at both surfaces, with a greater distribution of states on the crystalline surface. Peptide binding is primarily through lysine interactions, in line with prior experimental results.     

Molecular dynamics study of effects of nickel coating on torsional behavior of single-walled carbon nanotube

The effects of nickel coating on the torsional behaviors of single-walled carbon nanotubes (SWCNTs) subject to torsion motion are investigated using the molecular dynamics (MD) simulation method. The simulation results show that regardless of chirality, defect or radius, nickel coating can considerably enhance the critical torque of SWCNTs. However, by comparing the critical torsion angle between nickel-coated SWCNTs and corresponding pristine SWCNTs, it is found that nickel coating in small-radius nanotubes does induce a reduction in the critical torsion angle. The results also show that the structural failure of nickel coated imperfect (9,0) SWCNT occurs at an obviously higher critical torque in comparison with uncoated (9,0) SWCNT with a vacancy defect. Furthermore, we also find that the critical torque of a short nickel coated SWCNT is bigger than that of a long one, while the critical torsion angle for a short tube is smaller than that for a long one.     

Numerical simulation of deformation behavior of Al particles impacting on Al substrate and effect of surface oxide films on interfacial bonding in cold spraying

In this study, a comprehensive examination of the deformation behavior of Al particles impacting on Al substrate was conducted by using the Arbitrary Lagrangian Eulerian (ALE) method to clarify the deposition characteristics of Al powder and the effect of surface oxide films in cold spraying. It was found that the deformation behavior of Al particles is different from that of Cu particles under the same impact conditions owing to its lower density and thus less kinetic energy upon impact. The results indicated that a higher velocity was required for Al particles to reach the same compression ratio as that of Cu particles. On the other hand, the numerical results showed that the oxide films at particle surfaces influenced the deformation and bonding condition between the particle and substrate. The inclusions of the crushed oxide films at the interfaces between the depostied particles inhibit the deformation.     

Ab initio study of the structure and dynamical properties of crystalline ice

We investigated the structural and dynamical properties of a tetrahedrally coordinated crystalline ice from first principles based on density functional theory within the generalized gradient approximation with the projected augmented wave method. First, we report the structural behaviour of ice at finite temperatures based on the analysis of radial distribution functions obtained by molecular dynamics simulations. The results show how the ordering of the hydrogen bonding breaks down in the tetrahedral network of ice with entropy increase, in agreement with the neutron diffraction data. We also calculated the phonon spectra of ice in a 3× 1× 1 supercell using the direct method. So far, due to the direct method used in this calculation, the phonon spectra are obtained without taking into account the effect of polarization arising from dipole–dipole interactions of water molecules, which is expected to yield the splitting of longitudinal and transverse optic modes at the Γ point. The calculated longitudinal acoustic velocities from the initial slopes of the acoustic mode are in reasonable agreement with the neutron scattering data. Analysis of the vibrational density of states shows the existence of a boson peak at low energy of the translational region, a characteristic common to amorphous systems.     

First-principles molecular dynamics simulation of biased electrode/solution interface

First-principles molecular dynamics simulations have been carried out for water in contact with Pt(1 1 1) surface. To apply negative bias potential to the water/Pt interface, excess electrons were added to our slab model using the recently developed computational scheme called “effective screening medium (ESM)”. Water molecules located away from the surface reoriented themselves to screen the electric field, but they responded differently near the surface. Water molecules nearest to the surface, forming a distinct layered structure with the hydrogen atom directed to the surface, increased the density with increasing field. On these bases, we discuss microscopic aspects of the electric double layer.     

金属Ni原子对石墨烯生长结构的影响

基于Material Studio软件平台,利用分子动力学方法,对Ni原子与石墨烯层状结构相互作用和晶体结构变化过程进行模拟分析,得到如下结论：低浓度Ni原子会吸附在石墨烯表面层沿边缘生长,活性从中心向边缘逐渐降低,高浓度的Ni原子会溶解到内层石墨烯中.当石墨烯层数增加,附着在表层石墨烯的Ni原子生长排列范围扩大,且在石墨烯表面形成的点阵排列被破坏,附着在内层石墨烯的Ni原子比表层石墨烯Ni原子排列更散乱,同时石墨烯生长结构逐渐出现弯曲;随着层数增多和Ni原子浓度增加,石墨烯的拉伸强度也随之增加,石墨烯生长缺陷的偏转角度也随之增大.通过计算以上结构的径向分布函数(RDF),验证了石墨烯长程有序到短程有序的结构变化过程.     

Effect of deformation on the electrochemical behavior of hot-dip galvanized steel sheets

The main purpose of employing pre-coated steel sheets is to minimize corrosion of steel. However, coatings can be severely affected by forming processes. In the forming processes, due to the different modes of deformation, the strain levels are different and so can affect the properties of the coatings to a varying degree. Special attention has to be paid to the influence of deformation conditions on the performance of the coating, as regards protection against corrosion. The adhesion of the coating must remain good, and the surface should not be damaged during forming. Therefore, it is necessary to study the behavior of its corrosion resistance against the deformation.In this work, effect of strain path on the corrosion behavior of hot-dip galvanized steel sheets has been studied. Corrosion behavior of hot-dip galvanized steel sheets at various strain levels has been evaluated under four different strain paths namely, biaxial, plane strain, uniaxial (drawing) and tensile modes. The sheets were deformed by a limiting dome height test (LDH) set-up. A correlation between the degree of deformation and the loss in extent of corrosion protection offered by the coating has been established by carrying out electrochemical studies such as potentiodynamic polarization and impedance spectroscopy (EIS) in 3.5% NaCl solution. The present study shows that increase in deformation increases the extent of delamination of the coating for all the modes of deformation. The severity of deformation on delamination, however, has been found to vary in the order of, tensile < uniaxial < plane < biaxial.     

非晶Ag晶化过程中不同类型晶核结构的识别与跟踪

采用分子动力学模拟研究了非晶Ag的等温晶化过程,通过原子轨迹逆向追踪法分析了不同类型晶体团簇的结构遗传与组态演化.在团簇类型指数法的基础上,根据基本团簇种类与联结方式不同,提出了一种可区分fcc单晶、多晶与混晶团簇的分析方法.在非晶Ag等温晶化过程中,基于团簇结构的连续遗传性特征,发展了一种可区分fcc单晶、多晶与混晶晶胚与晶核的结构分析技术.结果发现:不论临界尺寸还是几何构型,不同类型的晶核结构都存在差异,其中fcc单晶临界尺寸最小,多晶次之,混晶最大; fcc单晶与多晶壳层原子中有少量hcp和bcc原子,而混晶壳层则全部为非晶类原子,并且fcc单晶、多晶与混晶的临界晶核都不是球型结构.     

Dynamic and spatial behavior of a corrugated interface in the driven lattice gas model

The spatiotemporal behavior of an initially corrugated interface in the two-dimensional driven lattice gas (DLG) model with attractive nearest-neighbors interactions is investigated via Monte Carlo simulations. By setting the system in the ordered phase, with periodic boundary conditions along the external field axis. i.e. horizontal, and open along the vertical directions respectively, an initial interface was imposed, that consists in a series of sinusoidal profiles with amplitude A₀ and wavelength λ set parallel to the applied driving field axis. We studied the dynamic behavior of its statistical width or roughness W(t), defined as the root mean square of the interface position. We found that W(t) decays exponentially for all λ and lattice longitudinal sizes L_x, i.e., the lattice side that runs along the axis of the external field. We determined its relaxation time τ, and found that depends on λ as a power law τ∝λ^p, where p depends on the temperature and L_x. At low T’s (T?T_c(E)) and large L_x, p approaches to p=3/2. At intermediate T’s (T<T_c(E)), p decreases up to p≈1, and is free of finite effects. This indicates that the interface stabilizes faster than in the equilibrium model, i. e. the Ising lattice gas (E=0) where p=3. At higher T’s p increases for T?T_c(E), and the finite size dependence is recovered. Also, if T is fixed, p increases with L_x until it saturates at large values of it, while this regime is vanishing at T?T_c(E). In this way, the dynamic relaxation process of a sinusoidal interface is improved by the external driving field with respect to its equilibrium counterpart, if the system is set in an intermediate temperature stage far from T_c(E) and in a lattice with a sufficiently large longitudinal side. The behavior of τ was also investigated as a function of E and in the intermediate stage T<T_c(E). It was found that τ decreases exponentially with E in the interval 0<E?1, while for higher fields it remains constant. The exponential decay depends on the wavelength of the initial profile.In order to study the spatial evolution of the profiles, we evaluated the structure factor of the interface, and the Fourier coefficients corresponding to the same wave vector of the initial profile. The obtained results allowed us to conclude that the spatial evolution of the profile maintains its initial wavelength, does not travel along the external field axis, and its shape is preserved over all the relaxation process.     

模型二元有序合金固液界面结构的分子动力学研究

采用分子动力学方法,研究了模型二元有序合金体系的平衡界面结构和界面处原子的扩散行为.计算结果表明,该二元有序合金的固液界面属于光滑界面.由于固体中同时存在结构和化学有序,从而导致界面处的原子结构与单质以及异质固液界面的结构明显不同.在界面法向方向上,粒子数密度呈复杂的波动行为,并延伸到液体中约30 (A).对界面层的二维结构分析表明,固液转变层部分原子形成了二维固体团簇.从固体到液体,扩散系数从零逐渐增加到一个饱和值.在界面处附近,平行于界面方向的扩散系数明显比垂直于界面方向的大.