Molecular dynamics study on the evolution of interfacial dislocation network and mechanical properties of Ni-based single crystal superalloys

The evolution of misfit dislocation network at $\gamma /{\gamma }^{\prime }$ phase interface and tensile mechanical properties of Ni-based single crystal superalloys at various temperatures and strain rates are studied by using molecular dynamics (MD) simulations. From the simulations, it is found that with the increase of loading, the dislocation network effectively inhibits dislocations emitted in the γ matrix cutting into the ${\gamma }^{\prime }$ phase and absorbs the matrix dislocations to strengthen itself which increases the stability of structure. Under the influence of the temperature, the initial mosaic structure of dislocation network gradually becomes irregular, and the initial misfit stress and the elastic modulus slowly decline as temperature increasing. On the other hand, with the increase of the strain rate, it almost has no effect on the elastic modulus and the way of evolution of dislocation network, but contributes to the increases of the yield stress and tensile strength. Moreover, tension–compression asymmetry of Ni-based single crystal superalloys is also presented based on MD simulations.     

Dielectrophoresis of nanocolloids: A molecular dynamics study

Dielectrophoresis (DEP), the motion of polarizable particles in non-uniform electric fields, has become an important tool for the transport, separation, and characterization of microparticles in biomedical and nanoelectronics research. In this article we present, to our knowledge, the first molecular dynamics simulations of DEP of nanometer-sized colloidal particles. We introduce a simplified model for a polarizable nanoparticle, consisting of a large charged macroion and oppositely charged microions, in an explicit solvent. The model is then used to study DEP motion of the particle at different combinations of temperature and electric field strength. In accord with linear response theory, the particle drift velocities are shown to be proportional to the DEP force. Analysis of the colloid DEP mobility shows a clear time dependence, demonstrating the variation of friction under non-equilibrium. The time dependence of the mobility further results in an apparent weak variation of the DEP displacements with temperature.     

Effect of thermal annealing on nanoimprinted Cu-Ni alloys using molecular dynamics simulation

The mechanical behaviors of nanoimprinted Cu-Ni alloys before and after annealing were studied using molecular dynamics simulation with a tight-binding potential. The results showed that when the punch is advancing, the punching force obtained from the simulation with a tight-binding potential is lower than with the Morse potential. During and after withdrawing the punch from the specimen, the adhesive phenomena are observed and the large residual stress in the Cu-Ni alloys is induced. During the annealing process, the internal energy of Cu-Ni alloys decreased with increasing the temperature and the component of Cu. In addition, comparing the maximum residual stress in the Cu-Ni alloys with and without annealing treatment, the stress is significantly released after annealing, especially in the higher component of Ni.     

The effect of temperature,defect and strain rate on the mechanical property of multi-layer graphene: Coarse-grained molecular dynamics study

In this work, we investigate the effect of temperature, defect, and strain rate on the mechanical properties of multi-layer graphene using coarse-grained molecular dynamics (CGMD) simulations. The simulation results reveal that the mechanical properties of multi-layer graphene tend to be less sensitive to temperature as the layer increases, but they are sensitive to the distribution and coverage of Stone-Wales (SW) defects. For the same number of defect, there is less decline in the fracture stress and Young's modulus of graphene when the defects have a regular distribution, in contrast to random distribution. In addition, Young's modulus is less influenced by temperature and defect, compared to fracture stress. Both the fracture stress and Young's modulus have little dependence on strain rate.     

氧化硅团簇切削单晶硅粗糙峰的分子动力学模拟研究

应用分子动力学模拟方法研究了氧化硅团簇在不同的切削 深度下切削单晶硅粗糙峰的过程, 考察了切削过程中粗糙峰和氧化硅团簇形态变化、团簇的受力状况、粗糙峰原子配位数和温度分布等. 模拟结果表明: 切削深度小于0.5 nm时, 被去除的材料以原子或者原子簇形式存在, 并黏附在颗粒表面被带走; 当切削深度增大至1 nm时, 材料的去除率增大, 并形成大的切屑. 在切削过程中, 由于压力和温度的升高, 粗糙峰切削区域的单晶硅转变为类似Si-Ⅱ相和Bct5-Si相的过渡结构, 在切削过程后的卸载阶段, 过渡结构由于压力和温度的下降转变为非晶态结构.     

Fe-Cr-Ni合金中点缺陷形成及相互作用的第一性原理研究

奥氏体不锈钢由于具有良好的力学性能和耐腐蚀性,被广泛应用于核电工业。奥氏体不锈钢受辐照会产生多种点缺陷,进而影响材料力学性能。为了研究304NG奥氏体不锈钢中点缺陷的产生机理,以面心立方Fe-Cr-Ni合金为模型,利用基于密度泛函理论的第一性原理方法研究了多种点缺陷结构的晶格常数,形成能,结合能,扩散迁移及相互作用,并计算了γ-Fe中相应的的性质作为对照.研究结果表明：Fe空位和Fe-Fe <100>哑铃是Fe-Cr-Ni合金中容易形成的点缺陷.基于对这两种点缺陷的迁移能垒分析,发现Fe-Fe <100>哑铃要比空位更易迁移.在Fe-Cr-Ni合金中相邻Fe空位存在吸引作用,会促进孔洞的形成;Fe空位-Fe-Fe <100>哑铃及Fe空位-Fe-Cr <100>哑铃均存在较弱的排斥作用,此外,在第一近邻位置和第二近邻位置的八面体间隙Fe原子都会和空位结合,使体系恢复为完整晶胞.     

Rotational viscosity of a liquid crystal mixture: a fully atomistic molecular dynamics study

Fully atomistic molecular dynamics (MD) simulations at 293, 303 and 313~K have been performed for the four-component liquid crystal mixture, E7, using the software package Material Studio. Order parameters and orientational time correlation functions (TCFs) were calculated from MD trajectories. The rotational viscosity coefficients (RVCs) of the mixture were calculated using the Nemtsov--Zakharov and Fialkowski methods based on statistical-mechanical approaches. Temperature dependences of RVC and density were discussed in detail. Reasonable agreement between the simulated and experimental values was found.     

Effect of Ni and vacancy concentration on initial formation of Cu precipitate in Fe–Cu–Ni ternary alloys by molecular dynamics simulation

In the present work, the effects of Ni atoms and vacancy concentrations(0.1%, 0.5%, 1.0%) on the formation process of Cu solute clusters are investigated for Fe–1.24%Cu–0.62%Ni alloys by molecular dynamics(MD) simulations. The presence of Ni is beneficial to the nucleation of Cu precipitates and has little effect on coarsening rate in the later stage of aging. This result is caused by reducing the diffusion coefficient of Cu clusters and the dynamic migration of Ni atoms. Additionally, there are little effects of Ni on Cu precipitates as the vacancy concentration reaches up to 1.0%,thereby explaining the embrittlement for reactor pressure vessel(RPV) steel. As a result, the findings can hopefully provide the important information about the essential mechanism of Cu cluster formation and a better understanding of ageing phenomenon of RPV steel. Furthermore, these original results are analyzed with a simple model of Cu diffusion, which suggests that the same behavior could be observed in Cu-containing alloys.     

A comprehensive molecular dynamics study of a single polystyrene chain in a good solvent

In this study, molecular characteristics of polystyrene (PS) was calculated measuring its dilute-solution properties in toluene at 288.15 K via molecular dynamics (MD) simulations. The solution models consisted of PS chains with different number of repeating units all of which were in a dilute regime. In order to investigate the compatibility between the polymer and the solvent molecules, interaction energy and Flory-Huggins (FH) interaction parameter were estimated. The simulation results indicate that increasing the chain repeating units enhanced the interaction between the solute and the solvent. Additionally, the chain dimensions were evaluated calculating the radius of gyration (R_g) and end-to-end distance, r0. To determine the dynamic behavior of the chains in the solutions, mean square displacement (MSD) and diffusivity coefficient were calculated. The simulation results indicated that the chain rigidity at low molecular weight and chain flexibility with increasing the molecular weight influenced chains dynamic behavior and diffusivity. Moreover, radial distribution function (RDF) illustrated the effect of steric hindrance of the chains in dilute solution on capturing the solvent molecules. In addition, solution viscosity was calculated by performing non-equilibrium molecular dynamics simulation (NEMD). The obtained results of chain characteristics and viscosity showed a good agreement with experimental results published previously. This agreement confirms the accuracy of the applied simulation method to characterize the dilute solutions and the chains characteristics.     

A molecular dynamics study of the mechanical properties of h-BCN monolayer using a modified Tersoff interatomic potential

Using molecular dynamics (MD) simulations, we investigate the mechanical properties of hexagonal BCN monolayer, a newly synthesized two-dimensional material with an atom ratio of B/C/N = 1:1:1. The Tersoff potential is modified to get good agreement between predicted and measured fracture strengths of graphene. With this modified Tersoff potential, we perform extensive MD simulations to study the effect of temperature, strain rate and vacancy defect on the mechanical properties of h-BCN. It is found that h-BCN is a strong material with fracture strength of 81.4–93.5 GPa, albeit ∼35% lower than that of graphene. Similar to graphene, temperature has strong effect on the mechanical properties of h-BCN. As the temperature increases from 10 K to 1300 K, the fracture strength and strain of h-BCN drops by 55% and 62%, respectively. The strain rate is found to have a moderate effect. When the strain rate increases from 0.00002 to 0.0125 ps⁻¹, the fracture strength and strain of h-BCN increases 6.1% and 12%, respectively. As for the atomic defect, a very small concentration (0.028%) of vacancy in h-BCN is able to cause a 28% reduction in fracture strength and a 35.5% reduction in fracture strain. These findings have significance for its future applications in nanodevices.     

Au-Pd共晶纳米粒子熔化行为的分子动力学研究

采用分子动力学方法结合嵌入原子势, 对Au-Pd共晶纳米粒子的热稳定性进行了模拟研究. 计算结果表明: Au-Pd纳米粒子的熔点明显高于Au单质纳米粒子而低于Pd纳米粒子. 通过计算Lindemann指数发现Au-Pd共晶纳米粒子中的Au原子首先熔化, 然后带动Pd原子的熔化; 熔化所经历的温度区间明显要宽于单质纳米粒子.     

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.     

Evaluating the effect of different test parameters on the tensile mechanical properties of single crystal silver nanowires using molecular dynamics simulation

Nanowires show amazing mechanical properties with respect to their bulk counterpart owing to their very high specific surface and/or interface area and, thus, are widely studied among several researchers. But it is difficult to study the mechanical properties of nanowires at atomistic level, and computational tools provide the required solution. Molecular dynamics simulation studies were carried out in this work to evaluate the mechanical properties of single crystal silver nanowire subjected to tensile deformation under varying wire diameter (4–14 nm), test temperature (100–500 K), and strain velocity (1–6 Å/ps). The simulation were carried out in analogous to real experiment, and the engineering stress and strain were calculated from the simulation result of load and displacement data. The mechanical properties like yield strength and Young’s modulus were calculated from the engineering stress-strain curve. The effect of different test parameters like wire diameter, equilibration temperature, and strain velocity on the mechanical properties were also thoroughly investigated. The result shows that single crystal silver nanowire shows excellent mechanical properties and, thus, can be used as a reinforcing agent to develop ultra-high strength advanced materials for defense and aerospace applications.     

Effects of different parameters on thermal and mechanical properties of aminated graphene/epoxy nanocomposites connected by covalent: A molecular dynamics study

This paper is devoted to studying the thermal and mechanical properties of aminated graphene (AG)/epoxy nanocomposites connected by covalent bond using molecular dynamics (MD) simulation. The effects of crosslinking degree, mass fraction and functionalized graphene (FG) type on AG/epoxy nanocomposites are considered. The elasticity modulus (E), the glass transition temperature (T_g), the coefficient of thermal expansion (CTE) and the interfacial energy (E_int) are also investigated. The MD simulation results indicate that, when the mass fraction of AG is between 1.2% and 3.1% and crosslinking degree reaches about 70%, the E, T_g, E_int and CTE of AG/epoxy nanocomposites are significantly improved compared with those of pure epoxy and graphene/epoxy nanocomposites. The reason is that AG not only possesses some excellent thermodynamic properties of graphene, but also has the function of curing agent to crosslink with epoxy monomer to form the carbon-nitrogen (C–N) covalent bond. A better interfacial interaction between nanoparticles and epoxy is essential in enhancing the thermal and mechanical properties of nanocomposite materials, which will provide a microscopic theoretical basis for the study of epoxy nanocomposites.     

非晶态合金Cu_xPd_(1-x)形成过程合金效应及微观结构转变特性的模拟研究(英文)

本文根据量子Sutton-Chen多体势,采用分子动力学方法对含有50000个原子大系统液态二元合金Cu_xPd_(1-x)(CuPd的原子半径比为1.14)在快速凝固过程中的微观结构转变特性进行模拟研究.运用Honeycutt-Andersen(HA)键型指数法和原子成团类型指数法(CTIM)分析了液态和固态的微观结构特性.研究结果发现:在7.73×10~(13)K/S冷却速率下,Cu_xPd_(1-x)合金形成以1551、1541和1431三种键型为主的非晶态结构;系统以1551键型和由1551键型构成的(12 0 12 0)二十面体团族在所有的键型和团簇中占主导地位,并且在液态合金Cu_xPd_(1-x)微观结构转变中起着关键性作用.通过分析键型、团簇和平均原子体积,我们发现液态合金Cu_xPd_(1-x)的玻璃转变温度是573K.同时还发现,原子的平均配位数的变化与1551,1441,1661键型的变化趋势相当接近,这反映出体系对称性结构的变化规律与配位数的变化有关.     

碳纳米锥力学特性的分子动力学研究

结合原子间短程作用势(Brenner势)和长程作用势(Lennard-Jones势),利用分子动力学方法对各种锥角的碳纳米锥进行拉伸和压缩实验,获得其载荷-应变关系曲线、受拉/压载荷极限、应变极限和构形演变等力学特性,并与等量原子组成的碳纳米管进行比较研究.研究结果表明,等量碳原子组成的碳纳米锥的受拉/压载荷极限随着锥角的增大先是增大后减小,受拉/压应变极限则随着锥角的增大而增大.与碳纳米锥相比,等量碳原子组成的碳纳米管的受拉/压载荷极限和应变极限显得既不突出也不逊色.在受压构形演化方面,与碳纳米管丰富的径向屈曲/扭转/侧向屈曲组合形变不同,112.88°和83.62°锥角的碳纳米锥受压沿轴向完美内陷,而60.0°和38.94°锥角的碳纳米锥受压发生侧向屈曲.     

Theoretical study on liquid crystal cyanobiphenyls: Phase stability and phase behavior

This article presents a theoretical study on liquid crystalline materials in homologous series of 4'-n-alkyl-4-cyanobiphenyl (nCB) with propyl (3CB), pentyl (5CB), and heptyl (7CB) groups. The atomic net charge and dipole moment components at each atomic center have been evaluated using the complete neglect differential overlap (CNDO/2) method. The modified Rayleigh–Schrodinger perturbation theory along with the multicentered-multipole expansion method has been employed to evaluate the long-range intermolecular interactions, while a ‘6-exp’ potential function has been assumed for short-range interactions. Further, these interaction energy values have been used as input to calculate the translational entropy, and free energy of nCB (n=3, 5, and 7) molecules during the stacking, and in-plane interactions. The observed results have been correlated with the mesogenic behavior and phase stability based on the thermodynamic parameters introduced in this article. Further, an attempt has been made to elucidate the flexibility of a configuration at a particular temperature, which has a direct relation with phase transition property of the molecules.     

Influence of line defects on relaxation properties of graphene: A molecular dynamics study

The relaxation properties of single layer graphene sheets containing line defects were investigated using molecular dynamics simulation with AIROBE bond-order interatomic potential. The dynamic evolution of graphene sheets during relaxation condition was analyzed. The simulation results show that the single layer graphene sheets are not perfectly flat in an ideal state, and the graphene sheet shows a significant corrugations at the verge of sheet. The graphene sheet is bent with the line defects at the end of the sheet, and the extent of this bend also increases with the increase of the defect number. Furthemore, the graphene sheet transforms into a paraboloid with the line defects at the middle of the sheet.     

KDP晶体固-液界面吸附行为的分子模拟研究

通过构建晶体表面-KDP分子界面吸附结构模型, 采用分子动力学和密度泛函计算方法研究KDP分子在(001)和(010)面吸附的物理化学过程, 考察了温度对物理吸附行为的影响. 研究表明: KDP晶体表面的吸附过程和生长习性主要由化学吸附主导, 化学吸附能的计算表明[K-O₈]基元在(001)界面的结合能是(010)界面结合能的2.86倍; 在饱和温度附近, [H₂PO₄]^-阴离子在KDP界面的物理结合能随温度的变化呈现振荡特征, 溶液中有较多的离子团簇形成, 溶液变得很不稳定; 当温度从323 K降低至308 K时, 水分子在界面的结合能总体呈下降趋势, 而KDP分子在界面的吸附能总体呈上升趋势, 脱水过程是水分子和[H₂PO₄]^-阴离子在固液界面边界层竞争吸附的结果. 研究结果对确足晶体生长界面动力学过程发展和完善晶体生长理论有重要意义. 关键词： 分子动力学 双层结构模型 结合能