Search

Precise prediction of phase diagrams in molecular dynamics simulations is challenging due to the simultaneous need for long time and large length scales and accurate interatomic potentials. We show that thermodynamic integration from low-cost force fields to neural network potentials trained using density-functional theory (DFT) enables rapid first-principles prediction of the solid–liquid phase boundary in the model salt NaCl. We use this technique to compare the accuracy of several DFT exchange–correlation functionals for predicting the NaCl phase boundary and find that the inclusion of dispersion interactions is critical to obtain good agreement with experiment. Importantly, our approach introduces a method to predict solid–liquid phase boundaries for any material at an ab initio level of accuracy, with the majority of the computational cost at the level of classical potentials. [ABSTRACT FROM AUTHOR]

Understanding how nanoparticles deform under compression not only is of scientific importance but also has practical significance in various applications such as tribology, nanoparticle-based probes, and the dry grinding of raw materials. In this study, we conducted compression tests on model brittle glassy nanoparticles using molecular dynamics simulations. We found that during the early stages of plastic deformation, shear bands formed in a similar pattern regardless of the nanoparticle size. However, as the deformation continued, dominant cracks emerged in large nanoparticles while being suppressed in smaller ones. This size-dependent brittle-to-ductile transition can be explained by a simple model based on Griffith's theory. We also investigated the effect of the surface stress state on fracture using thermally tempered nanoparticles. We observed that the presence of compressive surface stress strengthened the nanoparticle by suppressing crack formation, even when a pre-notch was present. On the other hand, tensile surface stress had the opposite effect. Interestingly, nanoparticles with both tensile and compressive surface stress promoted shear deformation, which could potentially compromise the mechanical performance of tempered glass despite delayed crack formation. [ABSTRACT FROM AUTHOR]

Figures 5-8 show density and Si/Al coordination change during the compression-decompression cycle with and without a shear stress component using the two sets of parameters. In the original paper, the short-range interaction parameters for the Si-Li pair in the SHIK potential[1] were used by mistake for the Si-Na pair in the sodium aluminosilicate (NAS) glass. Si and Al coordination changes during the compression-decompression cycle in NAS without [(a) and (b)] and with shear [(c) and (d)] using the two sets of short-range interaction parameters. [Extracted from the article]

Experimental studies have shown that glass systems with high boron content exhibit superior crack resistance under sharp contact loading. However, the underlying mechanism is still not fully understood. In this context, we carried out classical molecular dynamics simulations on sodium aluminosilicate and sodium aluminoborate systems to investigate the effect of boron on the response of glass to nanoindentation. A rigid V-shaped indenter is used to indent the glass sample with a fixed loading rate, during which the indenter interacts with the glass via a repulsive force field. The indenter angle and tip radius are varied to study the effect of indenter sharpness, as what has been done in experiments. These simulated nanoindentation tests reveal how the stress/strain field and the glass structure evolve with deformation underneath the indenter. It was found that a large number of boron atoms in the plastic zone change from three- to fourfold coordination during the loading process, and most of them revert back to the threefold coordination state during the unloading process. Our study shows that this "reversible" boron coordination change plays a critical role in increasing the damage resistance of glass. [ABSTRACT FROM AUTHOR]

Metallic glasses with secondary crystalline phase have been reported to either enhance or degrade the global plasticity. Here, we employed molecular dynamics simulation to investigate a generic glass-forming system, in which crystal-glass composites can be obtained via in situ crystallization. Ex situ composite samples were also constructed by embedding perfect crystals into monolithic glasses in order to delineate the precise influences of the composite microstructure. Our simulations demonstrate both toughening and embrittlement effects due to nanocrystallization. The double-edged-sword role of nanocrystals is mainly controlled by elasticity mismatch between the crystalline phase and the glass matrix. The effects of volume fraction and the size of nanocrystals can be understood in the perspective of cumulative probability of shear bands formation and propagation being interfered by the presence of nanocrystals. [ABSTRACT FROM AUTHOR]

Under sharp contact loading, glass deforms elastically and then plastically via densification as well as shear flow. Stress build-up during loading, residual stress build-up after unloading, and possible ultimate cracking all depend sensitively on the competition and interplay between densification and shear flow. The crack resistance was shown in experimental studies to generally improve with the increasing contribution of densification in glasses where shear deformation plays a dominant role under indentation. In this work, the role of densification in deformation behaviors of model metallic glasses under sharp contact loading was studied by 3-D nanoindentation tests using indenters with different sharpness in classical molecular dynamics simulations. Starting from a model metallic glass that favors shear deformation, a Lennard-Jones potential was modified to describe model metallic glasses with different abilities of instantaneous densification under compression and permanent densification after decompression. Our studies show that model metallic glass with a higher densification ability under indentation has less stress build-up and less localized shear deformation during loading, as well as smaller residual stress build-up after unloading. However, our study indicates that both instantaneous and permanent densification need to be tuned for designing damage resistant glasses. [ABSTRACT FROM AUTHOR]

Oxide glass, one of the most transformative materials in the modern world, breaks easily under load due to its brittleness. Using classical molecular dynamics simulations, we prepared amorphous alumina by consolidating glass nanoparticles at room temperature. We showed that consolidated amorphous alumina exhibits work hardening ability, hence deforms homogeneously and fractures via necking under tension, while amorphous alumina obtained from the traditional melt‐quench process fractures catastrophically due to severe shear banding. This finding suggests that if processed properly, amorphous oxides could deform and fracture like ductile metals, which will significantly expand the applications of oxide glasses into new areas where load bearing or mechanical reliability is necessary. [ABSTRACT FROM AUTHOR]

Temperature‐dependent viscosity is critical to decipher two profound questions in condensed matter physics, namely the glass transition and the relaxation of amorphous solids. However, direct measurement of viscosity over a large temperature range is extremely difficult. Here, using classical molecular dynamics (MD) simulations, we report a novel method to calculate the equilibrium viscosity of supercooled liquid both above and below the glass transition temperature (Tg) and to estimate the nonequilibrium viscosity of glass down to room temperature. Based on the shoving model, we derived an analytical formula showing that the shear viscosity in logarithmic scale changes linearly with the shear‐induced variation in shear modulus or potential energy of the glass‐forming system. The shear viscosity as a function of steady‐state potential energy of liquid under different shear strain rates can be directly calculated in MD simulations; together with its equilibrium potential energy, one can extrapolate the zero‐strain‐rate equilibrium viscosity. We verified the proposed model by reliably calculating equilibrium viscosity near Tg of four glass‐forming systems (Kob–Andersen system, silica, Cu45.5Zr45.5Al9, and silicon) with different fragilities. Furthermore, our model can estimate the nonequilibrium viscosity of glass below Tg; the upper‐bound nonequilibrium viscosity of amorphous silica and silicon at room temperature are calculated to be ~1032 and 1025 Pa·s, respectively. [ABSTRACT FROM AUTHOR]

Molecular dynamics simulations are used to investigate the mechanical properties of finite‐width nanocrystalline gold nanoribbons with a thickness of one grain size under tensile strain. We find strong variations in the Young's modulus and flow stress especially for the sample with narrow width. With statistical analysis, Young's modulus has a trend of increase with the decrease of width. Flow stress does not follow this trend. Ductile fracture is observed at large strain in these nanocrystalline nanoribbons, but the fracture strain and mechanism are different from those of thick films as observed in existing experiments. For small widths below 80 nm, fracture is a necking effect, while for larger widths, shear bands and formation of nanovoids due to stress concentration are the main mechanism. Voids are mainly distributed in special GBs that are parallel to the direction of the maximum principal stress of the plastic deformation tensor. We find there is a width regime (about 80–130 nm) for the ultrathin nanoribbon gold, which have good ductility and large fracture strain. [ABSTRACT FROM AUTHOR]

The recent realization of ferroelectricity in scandium‐ and boron‐substituted AlN thin films has spurred tremendous research interests. Here we established a molecular dynamics simulation framework to model the ferroelectricity of AlN thin films. Through reparameterization of Vashishta potential for AlN, the coercive field strength and the AlN polarization were found to be close to experimental values. Furthermore, we examined the effects of film thickness, temperature, in‐plane strain on polarization‐electric field hysteresis loop, and the thickness‐dependent Curie temperature. Lastly, we incorporated electrodes towards atomic‐level modeling of ferroelectric device, by considering the induced charge at the interface between electrodes and ferroelectric film. We found that low dielectric contrast significantly lowers the coercive field for switching AlN. [ABSTRACT FROM AUTHOR]

Molecular dynamics simulations were used to study the atomic mechanisms of deformation of nanocrystalline gold with 2.65–18 nm in grain size to explore the inverse Hall–Petch effect. Based on the mechanical responses, particularly the flow stress and the elastic-to-plastic transition, one can delineate three regimes: mixed (10–18 nm, dislocation activities and grain boundary sliding), inverse Hall-Petch (5–10 nm, grain boundary sliding), and super-soft (below 5 nm). As the grain size decreases, more grain boundaries present in the nanocrystalline solids, which block dislocation activities and facilitate grain boundary sliding. The transition from dislocation activities to grain boundary sliding leads to strengthening-then-softening due to grain size reduction, shown by the flow stress. It was further found that, samples with large grain exhibit pronounced yield, with the stress overshoot decrease as the grain size decreases. Samples with grain sizes smaller than 5 nm exhibit elastic-perfect plastic deformation without any stress overshoot, leading to the super-soft regime. Our simulations show that, during deformation, smaller grains rotate more and grow in size, while larger grains rotate less and shrink in size. [ABSTRACT FROM AUTHOR]

We designed catalytic nanorotors and investigated the rotational motion and energy conversion efficiency using reactive molecular dynamics in two dimensions. First, a two-arm nanorotor was constructed by decorating a slender beam with catalysts asymmetrically on its two long edges, while fixing the beam center as the rotational axis. Autonomous rotation was observed for the two-arm nanorotor immersing in a fuel environment. Here fuel molecules undergo exothermic combination reaction facilitated by the catalysts. It was found that the angular velocity increases with the catalyst coverage parabolically, while the rotary nanomotor efficiency stays roughly constant. These observations are consistent with a single-collision-momentum-transfer-based propulsion model. Furthermore, multi-arm nanorotors (up to eight arms) were constructed by carving radially distributed arms followed by decorating catalysts. For multi-arm nanorotors, both the angular velocity and the efficiency decrease as the number of arms increases. These behaviors contradict the aforementioned model, which are likely due to the deceleration from secondary collisions between products and the nanorotor arms. Our simulation results show that the optimal design for a nanorotor that maximizes its angular velocity and the motor efficiency is a two-arm nanorotor with nearly full coverage of catalysts. [ABSTRACT FROM AUTHOR]

Using molecular dynamics simulations, model Zr50Cu40Al10 metallic glass (MG) nanopillars were subjected to simulated irradiation processes followed by uniaxial compression tests. As the intensity or dosage of irradiation increases, the plastic deformation mode of the MG nanopillars transits from localized shear banding to homogeneous shear flow. The suppression of shear banding in MG nanopillars is due to irradiation-induced structural disordering. Furthermore, a correlation is found between the average potential energy of MG nanopillars and their deformation modes, common to both irradiation processing and thermal processing. Our results imply that the homogeneous shear flow observed in experimental MG nanopillars carved by focused ion beam may be due to irradiation damage instead of size effect. [ABSTRACT FROM AUTHOR]

We report compression-compression fatigue simulations on model metallic-glass nanowires up to a 4% maximum compressive strain and 100 cycles using the molecular dynamics method. The distribution and irreversibility of the fatigue-induced deformations are characterized by the local shear strain and rise of the system temperature, respectively. Furthermore, the system potential energy, structural signatures in terms of icosahedral centers, and radial distribution functions are obtained to identify structural transformations during cyclic loading. No noticeable changes on the aforementioned isotropic structural signatures have been found during fatigue tests. Based on a binary view of the constituting local structures (stable clusters and floppy clusters) in metallic glasses, we attribute the irreversible deformation to inter-conversion of those floppy clusters, which leads to a constant composition of stable and floppy clusters. Finally, uniaxial compression tests were carried out on the cyclically loaded samples. The mechanical behaviors of fatigue-tested samples are rather similar to the original sample. Upon cyclic loading, unlike macroscopic samples, a metallic-glass nanowire is free of permanent structural damage, and therefore does not exhibit deteriorated mechanical behaviors. The excellent fatigue-resistance behavior of a metallic-glass nanowire might be related to its defect-free structure (no stress concentrators). It appears that defects may play a central role in fatigue damages in experimental metallic-glass systems. [ABSTRACT FROM AUTHOR]

We report a dynamic self-assembly system of ‘living’ polymeric chains sustained by chemistry using reactive molecular dynamics simulations. The linear polymeric chains consist of self-assembled nanoparticles connected by metastable linker molecules. As such, the polymeric chains, once assembled, undergo spontaneous dissociation driven by thermodynamics. However, with a continuous supply of linker molecules and the stored chemical energy therein, the polymeric chains can survive and maintain a steady state averaged chain length. These dynamically self-assembled polymeric chains are analogous to biological systems that both are thermodynamically metastable, yet dynamically stable upon continuous influx of matter and energy. [ABSTRACT FROM AUTHOR]

The diamond cutting of transition metal alloys, particularly steel, is severely hindered by accelerated chemical wear of the tool. Recent experimental findings show that the presence of graphene platelets mitigates this problem. However, the specific mechanisms responsible for this wear mitigation are currently unknown. In this paper, molecular dynamics techniques are successfully used to identify these diamond tool wear mitigation mechanisms. A modified embedded atom method force field is first evaluated for its ability to accurately simulate the catalyzed graphitization of diamond in the presence of steel. This force field is then used to simulate nanometric diamond cutting of steel in the presence of a 1 or 3 layered graphene platelet. Coordination analysis of these simulations shows that the presence of the graphene platelet results in 34%–96% reduction in tool wear, when compared to the graphene-free cutting condition. This is attributed to the graphene platelet serving as a physical barrier protecting the tool cutting edge, and also as a sacrificial source for carbon transfer into the workpiece. Other mechanisms, such as platelet cleaving and interlayer sliding, are also observed. The reduction in tool wear reported by the simulations is comparable to the trends observed in prior experiments. [ABSTRACT FROM AUTHOR]

Most existing particle-based polymeric models are non-reactive, and usually constructed to be monodisperse. This deficiency limits the utilization of these models to understand the dependency of polymer properties on polydispersity. Here we report a coarse-grained model with pairwise interactions, yet is reactive and capable of in silico synthesis of polydisperse polymers in a step-wise manner. The polymerization of linear, branched, cross-linked and network polymers can be described by this reactive model. The chain length distribution of the resulting linear polymer agrees well with the Flory-Schulz distribution. The importance of polydispersity in polymer behaviors is highlighted in terms of kinetic (self-diffusivities of linear polymer melts) and mechanical properties (stress-strain responses of linear polymer glasses under uniaxial tension). [ABSTRACT FROM AUTHOR]

In molecular dynamics simulations, the local stress state in the shear band is examined in six different model metallic glasses and one amorphous Si system (also has been perceived as a metallic glass analog) under different loading conditions. For all but the FeP and the amorphous Si systems, the running shear band (RSB) exhibits a liquid-like hydrostatic plus shear stress state. Our results suggest that the liquid feature of a RSB is not due to temperature rise or plastic confinement but due to the disorder driven by flow, which can be offset by strong directionality in bonding, phase segregation, or aging. The knowledge of the liquid-like stress state can be conveniently utilized in experiments to infer the local stress state of the RSB from the global tensile stress for metallic glasses. [ABSTRACT FROM PUBLISHER]

We review a minimalist's reactive force field, reactive state summation (RSS) potential. The essence of RSS potential scheme is to model each reactive state by individual non-reactive force fields, then modulate each term by a reaction-coordinate-dependent weight function, finally sum together to obtain the reactive potential. Compared with existing reactive potentials, RSS potential is easier to formulate and parameterise and is computationally efficient, at the expense of lesser accuracy. Thus, RSS potential can be regarded as a ‘reactive Lennard-Jones’ potential. Three exemplary RSS potentials are described in the context of their respective chemical systems: RSS-nitrogen for modelling detonation, RSS-carbon for modelling pyrolysis of activated carbon and RSS-fuel-catalyst for modelling catalytic chemical reaction. [ABSTRACT FROM AUTHOR]

Tensile fracture of metallic glasses due to shear-band-to-crack transition was studied by molecular dynamics. Such transition is difficult to characterize experimentally as it occurs within nanoseconds in the buried nanometer-thin shear band. We show that the sample fractures via shear band cavitation under high initial tensile stress, or with long sample length. The critical thermomechanical conditions for shear band cavitation, strongly influenced by shear, were identified. An analytical model was developed, leading to a length scale and a time scale for shear-band-to-crack transition. [ABSTRACT FROM AUTHOR]

Abstract: We report the dynamic self-assembly of chemically propelled nanoscale building blocks (NBBs) under confinement using reactive molecular dynamics simulations. The NBBs are hollow rectangles with one side open (tail) to expose the catalysts inside to the fuel environment. Adjacent NBBs can self assembly into tail–tail or tail–head (the side opposite to the tail) pairs. It was found that chemical propulsion forces unlock many tail–tail pairs but not tail–head pairs. Although the distributions of both tail–tail and tail–head pairs are orientation-dependent with respect to the confining wall, the effect of chemistry seems to be orientation-independent. [Copyright &y& Elsevier]

Abstract: Simulations of a monatomic model amorphous matrix embedded with approximately 37% of a body-centered cubic phase demonstrate mechanisms by which nanocrystallites can alter the mechanical response of metallic glass. Three effects affect the resulting ductility: (i) the presence of weak amorphous–crystalline interfaces, (ii) the fraction of nanocrystallites oriented to prevent twinning relative to the loading stress, and (iii) the shear-induced growth and dissolution of the nanocrystallites when they are impinged by shear bands. While the first effect dominates in these simulations due to system size limitations, the third effect appears to be crucial for understanding the ductility of experimental samples. These simulations indicate that shear-induced growth of existing nanocrystallites, rather than nucleation of new crystalline regions, may account for the observed enhancement in ductility. [Copyright &y& Elsevier]

Abstract: Simulated nanoindentation tests on a three-dimensional model of a binary metallic glass-forming alloy reveal how the stress field and material structure interact to control deformation beneath the indenter. Initially, the stress field follows the Hertzian solution with fluctuations due to heterogeneities. Homogeneous or localized plastic deformation arises, depending on the processing history of the material. In the case of localized deformation, the first shear band initiates sub-surface, relaxing the local shear stress. The subsequent shear band morphology is observed to be rate-dependent. Deformation induced changes in material structure are characterized in terms of short-range ordering and free volume, with the former providing significantly greater signal-to-noise. [Copyright &y& Elsevier]

Abstract: Nanoindentation was simulated in a two dimensional model metallic glass thin film using molecular dynamics. Strain localization was observed in simulations where the system was sufficiently structurally relaxed prior to deformation. Indentation simulations utilized an atomic indenter that adhered to the surface or a frictionless indenter. Boundary conditions were varied to constrain the film or allow the film to relax in-plane to examine the effect on shear band formation. The most constrained system, i.e. that with the atomic indenter and the constrained boundaries exhibited the highest hardness. [Copyright &y& Elsevier]

The strength of nanocrystalline (NC) metal has been found to be sensitive to strain rate. Here, by molecular dynamics simulation, we explore the strain rate effects on apparent Young's modulus, flow stress and grain growth of NC gold with small size. The simulation results indicate that the apparent Young's modulus of NC gold decreases with the decrease of strain rate, especially for strain rates above 1 ns−1. The rearrangement of atoms near grain boundaries is a response to the decrease of apparent Young's modulus. Indeed, the flow stress is also sensitive to the strain rate and decreases following the strain rate's decrease. This can be found from the change of strain rate sensitivity and activation volume with the strain rate. Temperature has little effect on the activation volume of NC gold with small grain size, but has an obvious effect on that of relatively large grain size (such as 18 nm) under low strain rate (0.01 ns−1). Finally, grain growth in the deformation process is found to be sensitive to strain rate and the critical size for grain growth increases following the decrease of strain rate. [ABSTRACT FROM AUTHOR]

Silicon shows potential as an anode material in lithium ion batteries due to its high specific capacity, yet its considerable volume expansion during lithiation leads to fracture and pulverization. Unfortunately, neither the atomic-level structural evolution, nor the mechanical behavior of the anode during lithiation and delithiation cycles is well understood. Interestingly, the lithiation process of a-Si provides an interesting continuum from open-structured network glass to densely-packed atomic glass, which could be used to obtain useful insights regarding commonalities in glasses. Here atomic level simulation has been used to investigate one cycle of lithiation and delithiation of amorphous silicon electrode, using Grand Canonical Monte Carlo (GCMC) and molecular dynamics (MD) simulations. The atomic level structural transformation and damage accumulation of the anode during cycling has been systematically analyzed, as well as their mechanical responses in compact tension tests. There appears to be a ductile-brittle-ductile transition for the amorphous silicon anode during both the lithiation and delithiation cycle. In other words, amorphous silicon is particularly vulnerable at intermediate lithiation. The fracture behavior of lithiated silicon was found to correlate to the Poisson's ratio, due to variations in bond covalency and structural disorder. Image 1 [ABSTRACT FROM AUTHOR]

Understanding the dynamics of shear band propagation in metallic glasses remains elusive due to the limited temporal and spatial scales accessible in experiments. In micron-scale molecular dynamics simulations on two model metallic glasses, we studied the propagation of a dominant shear band under uniaxial tension with a macroscopic strain of 3-5%. For both materials, the shear band can be intersonic with a propagation speed exceeding their respective shear wave speeds. The propagation exhibits intrinsic instability that manifests itself as microbranching and considerable fluctuations in velocity. The shear strain singularity ahead of propagating shear band tip scales as 1/ r (r is the distance away from the tip), independent of the macroscopic tensile strain. In addition, we studied the intersection of two shear bands under uniaxial tension, during which path deflection, speed slowing-down, and temperature rise at the junction region were observed. The dynamics of propagating shear band shown here indicate that shear band in metallic glasses can be viewed as shear crack under the framework of weakly nonlinear fracture mechanics theory. [Display omitted] [ABSTRACT FROM AUTHOR]

We investigated the wear formation during a single-asperity sliding against a rigid substrate using molecular dynamics simulations. Two distinct wear mechanisms were observed: low-load atomic wear (isolated debris atoms or clusters) and high-load plastic wear (collective debris formation from plastic flow). The atomic wear rate depends on the normal stress exponentially, in agreement with the mechanically assisted bond-rupture theory. However, the plastic wear rate exhibits linear dependency on the normal stress. It was found that the asperity-substrate adhesion reduces the critical normal stress for the atomic-to-plastic wear transition, suppressing the atomic wear mode. Finally, a wear mechanism map of atomic/plastic wear was constructed in the domain of normal stress and adhesion, which is consistent with existing simulation and experimental results. This wear map could resolve the recent controversy on whether Archard׳s linear law is applicable for low-load tip wear. [ABSTRACT FROM AUTHOR]

[Display omitted] Al 2 O 3 –SiO 2 binary glasses with structural and chemical heterogeneities were prepared by consolidating glassy nanoparticles with different compositions using classical molecular dynamics simulations. Consolidated glasses show both excellent ductility and enhanced flow strength. It was found that the structural heterogeneities such as over-coordinated network formers and neighboring oxygen induced during consolidation serve as plasticity carriers to increase the ductility. On the other hand, the enhanced yield strength in consolidated glasses is due to the chemical heterogeneity inherited from the starting glassy nanoparticles, which does not compromise the ductility. Furthermore, apparent work hardening behavior appears upon cold work in the consolidated glasses, with an increase in yield strength from ~3.3 to ~6.4 GPa after 40% cold work. Consolidation with glassy nanoparticles could be a viable way to synthesize strong, damage resistant and transparent oxide glasses that cannot be obtained through the traditional melt-quench process. [ABSTRACT FROM AUTHOR]

Solids exhibit transverse shrinkage when they are stretched, except auxetics that abnormally demonstrate lateral expansion instead. Graphene possesses the unique normal-auxeticity (NA) transition when it is stretched along the armchair direction but not along the zigzag direction. Here we report on the anisotropic temperature-dependent NA transitions in strained graphene using molecular dynamics simulations. The critical strain where the NA transition occurs increases with respect to an increase in the tilt angle deviating from armchair direction upon uniaxial loading. The magic angle for the NA transition is 10.9°, beyond which the critical strain is close to fracture strain. In addition, the critical strain decreases with an increasing temperature when the tilt angle is smaller than the NA magic angle. Our results shed lights on the unprecedented nonlinear dimensional response of graphene to the large mechanical loading at various temperatures. Image 1 [ABSTRACT FROM AUTHOR]

Understanding the fuzz growth mechanism on tungsten (W) surface under the helium (He) particles irradiation is an important guide for structural design of the first wall materials. In this work, molecular dynamics simulation was used to explore the process of He bubble formation, growth and rupture under irradiation to understand the initial stage of fuzz nanostructure formation. It was found that He clusters had the most movable ability at the surface temperature around 800∼1800 K and thus big bubbles were favorable to be aggregated. We simulated the growth, rupture and merging of He bubbles to reproduce the process of protrusion's formation on the surface which induced the formation of fuzz nanostructure. It was found that the shallow bubble burst led to the pinhole in surface layer. The dislocation slip in the growth process of deep bubbles was the main factor for the configurations of protrusion on surface. It mainly pushes the W above He bubble along the <111> directions, while the different surfaces have different relation with the <111> directions. Thus, the different morphology of protrusion induced by the bubble appeared on different surfaces as observed in experiments. [ABSTRACT FROM AUTHOR]

The paradox of biomimetic research is to perform bio-functionality, usually associated with sophisticated structures optimized by nature, with minimal structural complexity for the ease of fabrication. Here we show that a three-particle trimer can exhibit kinesin-like autonomous walk on a track via reactive molecular dynamics simulations. The autonomous motion is due to imbalanced transitions resulting from exothermic catalytic reactions, and the spatial asymmetry from the track. This molecular design can be realized by reproducing the particle–particle interactions in functionalized nano- or colloidal particles. Our results open up the possibility of fabricating bio-mimetic nano-systems in a minimalist approach. [ABSTRACT FROM AUTHOR]

Low-cycle fatigue fracture of metallic glass nanowires was investigated using molecular dynamics simulations. The nanowires exhibit work hardening or softening, depending on the applied load. The structural origin of the hardening/softening response was identified as the decrease/increase of the tetrahedral clusters, as a result of the non-hardsphere nature of the glass model. The fatigue fracture is caused by shear banding initiated from the surface. The plastic-strain-controlled fatigue tests show that the fatigue life follows the Coffin–Manson relation. Such power-law form originates from plastic-strain-dependent microscopic damage accumulation. Lastly, the effect of a notch on low-cycle fatigue of nanowires in terms of failure mode and fatigue life was also discussed. [ABSTRACT FROM AUTHOR]

Atomic wear is generally understood as a mechanically-assisted and thermally-activated bond-breaking chemical process. In other words, debris formation is proportional to time just like a chemical reaction (instead of sliding distance), the contact area, and strongly temperature-dependent. Here using molecular dynamics, we show in a model glassy tribo-system at relatively high sliding speeds that the atomic wear depends weakly on the temperature and not on the sliding speed and increases sub-linearly with the contact area. We propose a different wear formulation for this system by considering that the debris-formation process at relatively high sliding speed is shear-assisted and athermally-activated. • Atomic wear in a model glass depends weakly on temperature & sliding speed. • Atomic wear in the tribo-system can be described by the effective temperature. • Effective temperature is a concept used in non-equilibrium thermodynamics. • An 'athermal' wear equation is derived to improve quantification of wear. [ABSTRACT FROM AUTHOR]

Abstract: Molecular dynamics simulations were carried out to investigate the wear formation of a truncated-cone-shaped tip sliding against a rigid flat substrate. Under a constant normal load, the worn volume of the tip increases sub-linearly with the sliding distance, which holds true for three distinct debris definitions (geometry-based, material-transfer-based and displacement-based). Under a constant normal stress, the tip height reduces linearly with the sliding distance, which can be conveniently used to calculate the wear rate. It was found that the tip height loss per sliding distance (ΔH/ΔL) is a function of the normal stress and the tapering angle of the tip, but independent of the sliding speed and the contact area. Above a critical stress, ΔH/ΔL increases linearly with the stress while the debris is generated from plastic flow. ΔH/ΔL is almost zero below the critical stress while the debris is generated from atom-by-atom attrition. Our results show that Archard's linear wear law is not applicable to single-asperity tip wear in the plastic wear regime. Finally, a new wear rate formula is given for tip wear as a function of the load, the contact area and the tapering angle of the tip. [Copyright &y& Elsevier]

Metallic glasses, like many other glasses and ceramics, generally have zero tensile ductility. Distinct from common wisdom which seeks to toughen brittle glasses with ductile phases, we show that composites made from two brittle glasses can be ductile, strong and even harden under tension as a result of microstructure and stiffness contrast. Using molecular dynamics simulations, we designed and tested composites consisting of alternating wavy nanofilaments of two brittle glasses with different stiffness. The composites exhibit failure strain over 40% and strain hardening modulus over 2 GPa. The tensile ductility is determined by both the flexibility of nanofilaments and the propensity of crack deflection, with which a ductile composite design map can be constructed. The strain hardening is due to filament alignment analogous to linear polymer. The brittle-brittle composite design strategy proposed here is not material-specific and should be applicable to, for instance, ultra-hard materials as constituents in achieving new damage-resistance composites. [Display omitted] [ABSTRACT FROM AUTHOR]

The changes of flow stress in voids-free, and HV and FV films are related to the dislocation density in the strain processes. Following the decrease of relative stacking faults, the decrease of flow stress in FV films is faster than that in HV film. Nanocrystalline ultrathin metallic thin films have various applications for which mechanical properties and durability are important, such as electrodes and electrical interconnects. We investigated the effect of voids on mechanical behavior of thin nanocrystalline gold films using molecular dynamics. Different kinds of voids were considered. We found that the nature and geometry of the void was important for the mechanical properties in this model ductile material. It was found that in the half voids films the flow stress decreased following the decrease of mass density with a linear relation. For the full voids films, the strain–stress curve was different from that of voids-free film and half voids films. The plastic deformation could be activated with smaller amount of dislocations. Following the strain increase over the strain of ultimate tensile strength, the flow stress decreased. This implied that the voids which penetrated through the films would weaken the ductility of films. The deformation main mechanism was grain boundaries sliding with voids coalescence. On the contrary, if the void depth was less than the thickness of film, the ductility of film wasn't influenced apparently, though the ultimate tensile strength decreased with the increase of voids size and decrease of mass density. [ABSTRACT FROM AUTHOR]

Surface diffusion is an interesting and practically important phenomenon in many areas of chemistry and physics. One emerging problem is that of Li wetting on high-Z plasma-facing metals, especially tungsten, for application in tokamak reactors. We report construction of a Li-W binary force field to describe the interface between liquid Li and solid W. We show that this force field can simulate well the atomic processes of liquid Li diffusing across the W surface, as compared with results from first-principles calculations. We find that Li atoms have different wetting behaviors on different surfaces of W. Diffusion is the fastest on the (110) surface and the slowest on the (100) surface. Diffusion of second Li layer on the (110) surface is activated at relatively low temperatures above 262 K. In addition, the diffusion rate of single-layer Li changes with the coverage rate of Li on (110) surface. High diffusion coefficients at coverage rate of θ = 1/6, 1/3 and 1 are found to be related to the 2D phase transition of Li on the (110) surface. [ABSTRACT FROM AUTHOR]