Nanoindentation simulations are a helpful complement to experiments. There is a dearth of nanoindentation simulations for bcc metals, partly due to the lack of computationally efficient and reliable ...interatomic potentials at large strains. We carry out indentation simulations for bcc tantalum using three different interatomic potentials and present the defect mechanisms responsible for the creation and expansion of the plastic deformation zone: twins are initially formed, giving rise to shear loop expansion and the formation of sequential prismatic loops. The calculated elastic constants as function of pressure as well as stacking fault energy surfaces explain the significant differences found in the defect structures generated for the three potentials investigated in this study. The simulations enable the quantification of total dislocation length and twinning fraction. The indenter velocity is varied and, as expected, the penetration depth for the first pop-in (defect emission) event shows a strain rate sensitivity m in the range of 0.037–0.055. The effect of indenter diameter on the first pop-in is discussed. A new intrinsic length-scale model is presented based on the profile of the residual indentation and geometrically necessary dislocation theory.
High-power, short-duration, laser-driven, shock compression and recovery experiments on 001 silicon unveiled remarkable structural changes above a pressure threshold. Two distinct amorphous regions ...were identified: (a) a bulk amorphous layer close to the surface and (b) amorphous bands initially aligned with {111} slip planes. Further increase of the laser energy leads to the re-crystallization of amorphous silicon into nanocrystals with high concentration of nano-twins. This amorphization is produced by the combined effect of high magnitude hydrostatic and shear stresses under dynamic shock compression. Shock-induced defects play a very important role in the onset of amorphization. Calculations of the free energy changes with pressure and shear, using the Patel-Cohen methodology, are in agreement with the experimental results. Molecular dynamics simulation corroborates the amorphization, showing that it is initiated by the nucleation and propagation of partial dislocations. The nucleation of amorphization is analyzed qualitatively by classical nucleation theory.
Display omitted
The key to perfect radiation endurance is perfect recovery. Since surfaces are perfect sinks for defects, a porous material with a high surface to volume ratio has the potential to be extremely ...radiation tolerant, provided it is morphologically stable in a radiation environment. Experiments and computer simulations on nanoscale gold foams reported here show the existence of a window in the parameter space where foams are radiation tolerant. We analyze these results in terms of a model for the irradiation response that quantitatively locates such window that appears to be the consequence of the combined effect of two length scales dependent on the irradiation conditions: (i) foams with ligament diameters below a minimum value display ligament melting and breaking, together with compaction increasing with dose (this value is typically ∼5 nm for primary knock on atoms (PKA) of ∼15 keV in Au), while (ii) foams with ligament diameters above a maximum value show bulk behavior, that is, damage accumulation (few hundred nanometers for the PKA's energy and dose rate used in this study). In between these dimensions, (i.e., ∼100 nm in Au), defect migration to the ligament surface happens faster than the time between cascades, ensuring radiation resistance for a given dose-rate. We conclude that foams can be tailored to become radiation tolerant.
The mechanisms of deformation under a nanoindentation in tantalum, chosen as a model body-centered cubic (bcc) metal, are identified and quantified. Molecular dynamics (MD) simulations and ...indentation experiments are conducted for 100, 110 and 111 normals to surface orientations. The simulated plastic deformation proceeds by the formation of nanotwins, which rapidly evolve into shear dislocation loops. It is shown through a dislocation analysis that an elementary twin (three layers) is energetically favorable for a diameter below ∼7nm, at which point a shear loop comprising a perfect dislocation is formed. MD simulations show that shear loops expand into the material by the advancement of their edge components. Simultaneously with this advancement, screw components of the loop cross-slip and generate a cylindrical surface. When opposite segments approach, they eventually cancel by virtue of the attraction between them, forming a quasi-circular prismatic loop composed of edge dislocation segments. This “lasso”-like mechanism by which a shear loop transitions to a prismatic loop is identified for both 001 and 111 indentations. The prismatic loops advance into the material along 〈111〉 directions, transporting material away from the nucleation site. Analytical calculations supplement MD and experimental observations, and provide a framework for the improved understanding of the evolution of plastic deformation under a nanoindenter. Dislocation densities under the indenter are estimated experimentally (∼1.2×1015m−2), by MD (∼7×1015m−2) and through an analytical calculation (2.6–19×1015m−2). Considering the assumptions and simplifications, this agreement is considered satisfactory. MD simulations also show expected changes in pile-up symmetry after unloading, compatible with crystal plasticity.
High Entropy Alloys (HEA) attract attention as possible radiation resistant materials, a feature observed in some experiments that has been attributed to several unique properties of HEA, in ...particular to the disorder-induced reduced thermal conductivity and to the peculiar defect properties originating from the chemical complexity. To explore the origin of such behavior we study the early stages (less than 0.1 ns), of radiation damage response of a HEA using molecular dynamics simulations of collision cascades induced by primary knock-on atoms (PKA) with 10, 20 and 40 keV, at room temperature, on an idealized model equiatomic quinary fcc FeNiCrCoCu alloy, the corresponding “Average Atom” (AA) material, and on pure Ni. We include accurate corrections to describe short-range atomic interactions during the cascade. In all cases the average number of defects in the HEA is lower than for pure Ni, which has been previously used to help claiming that HEA is radiation resistant. However, simulated defect evolution during primary damage, including the number of surviving Frenkel Pairs, and the defect cluster size distributions are nearly the same in all cases, within our statistical uncertainty. The number of surviving FP in the alloy is predicted fairly well by analytical models of defect production in pure materials. All of this indicates that the origin of radiation resistance in HEAs as observed in experiments may not be related to a reduction in primary damage due to chemical disorder, but is probably caused by longer-time defect evolution.
Display omitted
Despite its fundamental importance for a broad range of applications, little is understood about the behaviour of metals during the initial phase of shock compression. Here, we present molecular ...dynamics (MD) simulations of shock-wave propagation through a metal allowing a detailed analysis of the dynamics of high strain-rate plasticity. Previous MD simulations have not seen the evolution of the strain from one- to three-dimensional compression that is observed in diffraction experiments. Our large-scale MD simulations of up to 352 million atoms resolve this important discrepancy through a detailed understanding of dislocation flow at high strain rates. The stress relaxes to an approximately hydrostatic state and the dislocation velocity drops to nearly zero. The dislocation velocity drop leads to a steady state with no further relaxation of the lattice, as revealed by simulated X-ray diffraction.
We use molecular dynamics (MD) simulations to deform single crystal spherical carbon nanoparticles (NP), 4–45 nm diameter, with a hard, flat indenter, compressing along the 001 direction. There is no ...clear amorphization nor phase change in the NP, but there is significant deformation, with bent crystalline planes, and many atoms that retain sp3 coordination, but are no longer recognized as having diamond structure by different structure-identification methods. Machine-learning is used to improve diamond-structure identification. The NP deforms laterally, and volumetric strain is ~0.1 when the uniaxial strain is ~0.5. Poisson's ratio increases with strain, and the elastic limit is reached at 0.2–0.3 strain, at a contact pressure of ~150 GPa. For NPs above 5 nm, dislocations appear and are mostly (1/2) {111} full dislocations, with a few partial dislocations for larger nanoparticles, without twinning. These results agree with the recent observation of plastic deformation in diamond nanopillars. Small NP display elastic modulus, yield stress and hardness increasing with NP size, but NPs with diameter larger than 25 nm display an approximately constant dislocation and dislocation junction density, which leads to a plateau in the hardness versus NP size, at ~150 GPa, close to bulk diamond. Diamond nanoparticles could provide high strength thin coatings, lighter than full-density nanotwinned diamond but with nearly the same strength.
Display omitted
•Indented diamond nanoparticles (NP) show large elastic limits at strains of 0.2‐0.3.•Structure-identification methods including machine-learning discard phase changes.•Dislocations are mostly (1/2)<1 10>{111} full dislocations, screw character.•Elastic modulus and yield strength increase with NP size, then reach a plateau.•NP hardness is ~150 GPa, comparable to bulk diamond.
Hydrogen absorption in Pd thin-films Ramos de Debiaggi, S.; Crespo, E.A.; Braschi, F.U. ...
International journal of hydrogen energy,
05/2014, Letnik:
39, Številka:
16
Journal Article, Conference Proceeding
Recenzirano
Hydrogen absorption isotherms for Pd thin films were modeled at atomistic scale by Monte Carlo (MC) simulation in the TPμN ensamble and by Molecular Dynamics (MD) simulations at 300 K. The ...interaction among atoms was modeled by embedded atom method (EAM) potentials. Simulated samples consisted of monocrystalline nanofilms with different thickness (2–8 nm) and two crystallographic surface orientations, (001) and (111). The isotherms were compared to bulk Pd and a few available experimental results. Instead of the plateau corresponding to the α-β PdH equilibrium in the bulk, the isotherms at nano-films show a two-plateaux behavior: a small one corresponding to a surface–subsurface hydride formation, and a larger one for the subsequent bulk hydride formation. This is strongly correlated with the atomic stress distribution induced within the thin film. The equilibrium pressures at the isotherms depend on the thin-film thickness, with pressure being larger for thicker films. The isotherms of the (001) films display lower equilibrium pressures than those for (111) films.
•Simulation of H absorption isotherms in Pd unsupported thin films.•Absorption isotherms are related to the thickness of the film.•Also to the surface orientation of the film ((111) and (001)).
Dislocations are the primary agents of permanent deformation in crystalline solids. Since the theoretical prediction of supersonic dislocations over half a century ago, there is a dearth of ...experimental evidence supporting their existence. Here we use non-equilibrium molecular dynamics simulations of shocked silicon to reveal transient supersonic partial dislocation motion at approximately 15 km/s, faster than any previous in-silico observation. Homogeneous dislocation nucleation occurs near the shock front and supersonic dislocation motion lasts just fractions of picoseconds before the dislocations catch the shock front and decelerate back to the elastic wave speed. Applying a modified analytical equation for dislocation evolution we successfully predict a dislocation density of 1.5 × 10(12) cm(-2) within the shocked volume, in agreement with the present simulations and realistic in regards to prior and on-going recovery experiments in silicon.