Probing the limits of metal plasticity with molecular dynamics simulations.

Ordinarily, the strength and plasticity properties of a metal are defined by dislocations-line defects in the crystal lattice whose motion results in material slippage along lattice planes. Dislocation dynamics models are usually used as mesoscale proxies for true atomistic dynamics, which are computationally expensive to perform routinely. However, atomistic simulations accurately capture every possible mechanism of material response, resolving every "jiggle and wiggle" of atomic motion, whereas dislocation dynamics models do not. Here we present fully dynamic atomistic simulations of bulk single-crystal plasticity in the body-centred-cubic metal tantalum. Our goal is to quantify the conditions under which the limits of dislocation-mediated plasticity are reached and to understand what happens to the metal beyond any such limit. In our simulations, the metal is compressed at ultrahigh strain rates along its [001] crystal axis under conditions of constant pressure, temperature and strain rate. To address the complexity of crystal plasticity processes on the length scales (85-340 nm) and timescales (1 ns-1µs) that we examine, we use recently developed methods of in situ computational microscopy to recast the enormous amount of transient trajectory data generated in our simulations into a form that can be analysed by a human. Our simulations predict that, on reaching certain limiting conditions of strain, dislocations alone can no longer relieve mechanical loads; instead, another mechanism, known as deformation twinning (the sudden re-orientation of the crystal lattice), takes over as the dominant mode of dynamic response. Below this limit, the metal assumes a strain-path-independent steady state of plastic flow in which the flow stress and the dislocation density remain constant as long as the conditions of straining thereafter remain unchanged. In this distinct state, tantalum flows like a viscous fluid while retaining its crystal lattice and remaining a strong and stiff metal.

Atomistic insights into metal hardening.

For millennia, humans have exploited the natural property of metals to get stronger or harden when mechanically deformed. Ultimately rooted in the motion of dislocations, mechanisms of metal hardening have remained in the cross-hairs of physical metallurgists for over a century. Here, we performed atomistic simulations at the limits of supercomputing that are sufficiently large to be statistically representative of macroscopic crystal plasticity yet fully resolved to examine the origins of metal hardening at its most fundamental level of atomic motion. We demonstrate that the notorious staged (inflection) hardening of metals is a direct consequence of crystal rotation under uniaxial straining. At odds with widely divergent and contradictory views in the literature, we observe that basic mechanisms of dislocation behaviour are the same across all stages of metal hardening.

Mesoscale Particle-Based Model of Electrophoretic Deposition.

We present and evaluate a semiempirical particle-based model of electrophoretic deposition using extensive mesoscale simulations. We analyze particle configurations in order to observe how colloids accumulate at the electrode and arrange into deposits. In agreement with existing continuum models, the thickness of the deposit increases linearly in time during deposition. Resulting colloidal deposits exhibit a transition between highly ordered and bulk disordered regions that can give rise to an appreciable density gradient under certain simulated conditions. The overall volume fraction increases and falls within a narrow range as the driving force due to the electric field increases and repulsive intercolloidal interactions decrease. We postulate ordering and stacking within the initial layer(s) dramatically impacts the microstructure of the deposits. We find a combination of parameters, i.e., electric field and suspension properties, whose interplay enhances colloidal ordering beyond the commonly known approach of only reducing the driving force.

On coarse projective integration for atomic deposition in amorphous systems.

Direct molecular dynamics simulation of atomic deposition under realistic conditions is notoriously challenging because of the wide range of time scales that must be captured. Numerous simulation approaches have been proposed to address the problem, often requiring a compromise between model fidelity, algorithmic complexity, and computational efficiency. Coarse projective integration, an example application of the "equation-free" framework, offers an attractive balance between these constraints. Here, periodically applied, short atomistic simulations are employed to compute time derivatives of slowly evolving coarse variables that are then used to numerically integrate differential equations over relatively large time intervals. A key obstacle to the application of this technique in realistic settings is the "lifting" operation in which a valid atomistic configuration is recreated from knowledge of the coarse variables. Using Ge deposition on amorphous SiO2 substrates as an example application, we present a scheme for lifting realistic atomistic configurations comprised of collections of Ge islands on amorphous SiO2 using only a few measures of the island size distribution. The approach is shown to provide accurate initial configurations to restart molecular dynamics simulations at arbitrary points in time, enabling the application of coarse projective integration for this morphologically complex system.

Extreme Metastability of Diamond and its Transformation to the BC8 Post-Diamond Phase of Carbon.

Diamond possesses exceptional physical properties due to its remarkably strong carbon-carbon bonding, leading to significant resilience to structural transformations at very high pressures and temperatures. Despite several experimental attempts, synthesis and recovery of the theoretically predicted post-diamond BC8 phase remains elusive. Through quantum-accurate multimillion atom molecular dynamics (MD) simulations, we have uncovered the extreme metastability of diamond at very high pressures, significantly exceeding its range of thermodynamic stability. We predict the post-diamond BC8 phase to be experimentally accessible only within a narrow high pressure-temperature region of the carbon phase diagram. The diamond to BC8 transformation proceeds through premelting followed by BC8 nucleation and growth in the metastable carbon liquid. We propose a double-shock compression pathway for BC8 synthesis, which is currently being explored in experiments at the National Ignition Facility.

Atomistics of vapour-liquid-solid nanowire growth.

Vapour-liquid-solid route and its variants are routinely used for scalable synthesis of semiconducting nanowires, yet the fundamental growth processes remain unknown. Here we employ atomic-scale computations based on model potentials to study the stability and growth of gold-catalysed silicon nanowires. Equilibrium studies uncover segregation at the solid-like surface of the catalyst particle, a liquid AuSi droplet, and a silicon-rich droplet-nanowire interface enveloped by heterogeneous truncating facets. Supersaturation of the droplets leads to rapid one-dimensional growth on the truncating facets and much slower nucleation-controlled two-dimensional growth on the main facet. Surface diffusion is suppressed and the excess Si flux occurs through the droplet bulk which, together with the Si-rich interface and contact line, lowers the nucleation barrier on the main facet. The ensuing step flow is modified by Au diffusion away from the step edges. Our study highlights key interfacial characteristics for morphological and compositional control of semiconducting nanowire arrays.

Size effects on the mechanical behavior of nanoporous Au.

Recent nanomechanical tests on submicron metal columns and wires have revealed a dramatic increase in yield strength with decreasing sample size. Here, we demonstrate that nanoporous metal foams can be envisioned as a three-dimensional network of ultrahigh-strength nanowires, thus bringing together two seemingly conflicting properties: high strength and high porosity. Specifically, we characterized the size-dependent mechanical properties of nanoporous gold using a combination of nanoindentation, column microcompression, and molecular dynamics simulations. We find that nanoporous gold can be as strong as bulk Au, despite being a highly porous material, and that the ligaments in nanoporous gold approach the theoretical yield strength of Au.