🔬 materials science

Crystal forming ability of amorphous refractory metals under nanoindentation: a molecular dynamics study

This molecular dynamics study reveals that amorphous refractory metals (V, Nb, Mo, Ta, W) undergo shear-driven bulk nucleation and growth to transform into bcc crystals during nanoindentation, with their crystal forming ability decreasing from V to W and scaling inversely with indentation velocity, indicating that the transformation rate is governed by cohesive bond strength rather than thermodynamic driving force.

Original authors: Prashant Dwivedi, Alberto Fraile, Tomas Polcar

Published 2026-07-30

📖 7 min read🧠 Deep dive

CC BY 4.0

Original authors: Prashant Dwivedi, Alberto Fraile, Tomas Polcar

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Secret Life of Metals: When Glass Gets Shaken

Imagine a world where metal doesn't just bend or break, but can actually turn into glass. You might think of glass as something fragile, like a window pane, but in the world of advanced materials, "metallic glass" is a super-strong, super-tough substance. Unlike normal metals, which are made of tiny, orderly crystals (like soldiers standing in perfect rows), metallic glass is a chaotic mess of atoms, frozen in place just like the atoms in a liquid. This disorder gives it incredible strength, but it also makes it a bit unstable, like a house of cards waiting for a breeze.

Scientists have long known that if you heat this glassy metal, it will eventually melt back into a liquid and then try to re-form into those orderly crystals. But what happens if you don't use heat? What happens if you just push on it? This is the question of "mechanical devitrification." It's like asking if you can turn a chaotic pile of sand into a perfect sandcastle just by poking it with a stick. For years, researchers have wondered: if you poke different types of glassy metals, which ones will snap back into their crystal form the easiest? And does the speed of your poke matter? This paper dives deep into that mystery, using powerful computer simulations to watch these atomic transformations happen in real-time, revealing that the answer depends on how tightly the atoms hold hands with each other.

The Great Atomic Shake-Up

In this study, the researchers acted like cosmic chefs, but instead of cooking food, they were cooking metals. They took five specific "refractory" metals—Vanadium (V), Niobium (Nb), Molybdenum (Mo), Tantalum (Ta), and Tungsten (W)—which are known for being incredibly tough and heat-resistant. First, they melted these metals down in a computer simulation and then froze them so fast (at rates of about $10^{13}$ Kelvin per second!) that the atoms didn't have time to line up. Instead, they got stuck in a messy, glassy state.

Once they had their five different glassy metals, they put them to the test. Using a virtual "nano-indenter"—think of it as a super-tiny, perfectly round marble the size of a virus—they pressed down on the surface of each glass. They did this at different speeds, ranging from a slow $0.075$ Å ps $^{-1}$ to a faster $0.30$ Å ps $^{-1}$ . The goal was to see if the pressure of the poke would force the chaotic atoms to snap back into their original, orderly body-centered cubic (bcc) crystal structure.

The Results: Who Crystallizes First?

The simulations showed that yes, the pressure worked! In every single case, the glassy metal started to turn back into a crystal right under the tip of the indenter. It wasn't a slow, uniform change; instead, tiny crystal "seeds" popped up in random spots, grew bigger, and then crashed into each other to form a solid crystal patch.

But here is the big discovery: not all metals were equally eager to change. The researchers created a score called "Crystal Forming Ability" (CFA) to measure how fast the transformation happened per inch of depth the indenter pushed. The results were a clear ranking:

Vanadium (V) was the fastest, turning into crystal most easily.
Molybdenum (Mo) came in second.
Niobium (Nb) was third.
Tantalum (Ta) was fourth.
Tungsten (W) was the slowest, the most stubborn, and the hardest to force into a crystal.

The order was V > Mo > Nb > Ta > W.

The Speed Factor: A Surprising Constant

One might guess that if you poke faster, the metal has less time to change, so it would transform less. And indeed, if you poke faster, the transformation happens deeper down. However, the researchers found something fascinating about the rate of change. When they looked at the math, they discovered that the product of the Crystal Forming Ability and the speed of the poke ( $CFA \cdot v$ ) was almost the same for all speeds.

Think of it like a conveyor belt. If the belt moves twice as fast, the items on it move twice as far in the same time, but the rate at which they get processed stays the same. In these simulations, the transformation rate per second was almost independent of how fast the indenter was moving. This suggests that the process is driven by the mechanical work being done, not just by how much time the atoms have to think about it.

Why Do Some Metals Resist?

So, why is Tungsten (W) so stubborn while Vanadium (V) is so eager? The team tested several theories. They looked at the "thermodynamic driving force"—basically, how much the metal wants to be a crystal because it's a lower energy state. Surprisingly, this didn't match the results. In fact, Tungsten had the strongest desire to be a crystal, yet it was the hardest to force into that shape.

Instead, the answer lay in cohesive bond strength. This is a measure of how tightly the atoms hold onto each other. The metals that were hardest to transform (like Tungsten) had the strongest bonds. It's like trying to rearrange a group of people holding hands: if they are holding hands loosely (Vanadium), it's easy to shuffle them into a line. If they are gripping each other with super-strong glue (Tungsten), it takes a massive amount of effort to break that grip and reorganize them. The study showed that the resistance to turning back into a crystal tracks directly with how strong the atomic bonds are, not with how much energy is released by the change.

The Role of Shear: The "Twist" That Does It

The researchers also looked at how the atoms moved just before they turned into crystals. They found that the atoms didn't just get squished by pressure (hydrostatic pressure). Instead, the key was shear strain—a kind of twisting or sliding motion. The atoms that eventually turned into crystals were the ones that had already been "slipping" or moving in a chaotic, non-uniform way (called non-affine displacement) before the crystal formed.

This confirms that the transformation is a "shear-associated" event. It's not just about being squeezed; it's about being twisted and sheared. The atoms that were already sliding past their neighbors were the ones most likely to snap into the new crystal order.

The Final Picture: A Population of Grains

Finally, the study tracked what happened to the tiny crystal grains as they grew. They found that the process involves three stages:

Nucleation: Tiny crystal seeds appear.
Growth: These seeds get bigger.
Coalescence: The seeds run into each other and merge.

The researchers found that for most of these metals, growth was the main driver, accounting for between 64% and 93% of the volume increase before the crystals peaked. However, the final size of the grains depended on how many seeds survived and how much they merged. For example, Vanadium and Molybdenum had lots of seeds that merged heavily, while Tantalum and Niobium had fewer seeds that grew into larger, distinct grains.

The Bottom Line

This paper didn't just watch metals turn into crystals; it mapped out the rules of the game. It showed that for these five refractory metals, the ease of turning glass back into crystal is a race against the strength of their own atomic bonds. The stronger the bond, the harder it is to force the change, regardless of how much energy the change would release. It also proved that this transformation is a mechanical dance driven by shear and sliding, not just by heat or simple pressure. While these results come from computer simulations (which are incredibly detailed but still models), they provide a clear, unified picture of how stress can rewrite the atomic structure of some of the world's toughest metals.

Technical Summary: Crystal Forming Ability of Amorphous Refractory Metals Under Nanoindentation

Problem Statement
Amorphous refractory metal coatings offer high hardness and chemical inertness but possess a metastable nature that renders them susceptible to mechanically induced crystallisation (devitrification) under contact loading. While the glass-forming ability (GFA) of multicomponent alloys is well-characterized by thermal descriptors, no comparable normalized metric exists for the inverse propensity: the ease with which a deformed, single-element amorphous metal reverts to its crystalline ground state under mechanical work. Furthermore, the variation in this propensity across the body-centered cubic (bcc) refractory series (V, Nb, Mo, Ta, W) and the specific structural and mechanical variables governing the transformation pathway remain unestablished.

Methodology
The study employs large-scale molecular dynamics (MD) simulations to investigate the nanoindentation response of amorphous V, Nb, Mo, Ta, and W.

Sample Preparation: Amorphous configurations were generated via melt-quenching from bcc crystals at effective cooling rates of $1.4 \text{--} 4.6 \times 10^{13} \text{ K s}^{-1}$ to 300 K. Three independent glass replicas were created for each element using distinct random seeds.
Potential Validation: Selected interatomic potentials (EAM, Finnis-Sinclair, and angular-dependent potentials) were validated against ab initio liquid radial distribution functions (RDFs) to ensure the simulated amorphous structures inherited correct short- and medium-range topology from the liquid precursors.
Nanoindentation Protocol: Spherical nanoindentation was performed at 300 K with a tip radius of 100 Å at velocities of 0.075, 0.15, and 0.30 Å ps $^{-1}$ . Ten seeded trajectories were analyzed for each glass replica and velocity combination.
Analysis: Local bcc environments were identified using polyhedral template matching. The transformation kinetics were quantified via a "Crystal Forming Ability" (CFA), defined as the maximum slope of the sigmoidal bcc fraction versus indentation depth curve. A case-control analysis compared transforming atoms against depth-matched amorphous controls to identify precursor deformation states.

Key Contributions and Results

Operational Crystal Forming Ability (CFA):
The study introduces CFA as a depth-normalized metric for transformation resistance. At the reference velocity, the CFA ranking is V > Mo > Nb > Ta > W. This indicates that Vanadium transforms most readily, while Tungsten requires the greatest indentation advance per unit of transformed fraction.
- Velocity Dependence: CFA decreases with increasing indenter velocity following a power law $CFA \propto v^{-m}$ , with a mean exponent $m \approx 1.08$ . Consequently, the product $CFA \cdot v$ is nearly constant, suggesting the time-domain transformation rate is almost independent of indenter velocity across the tested range.
Transformation Mechanism:
The amorphous-to-bcc transition proceeds via bulk nucleation, growth, and coalescence.
- Shear Bias: Atoms that subsequently transform carry significantly higher non-affine displacement ( $D^2_{min}$ ) and local shear strain in their precursor state compared to atoms that remain amorphous. Hydrostatic pressure shows only a minor negative contrast, indicating the pathway is shear-associated and displacive rather than pressure-driven.
- Kinetics: The transformation involves the formation of compact bcc grains that grow and impinge. Grain growth contributes 64% to 93% of the pre-peak volume increase, while coalescence determines the final microstructure.
Drivers of Resistance:
- Cohesive Bond Strength vs. Thermodynamics: Resistance to devitrification correlates perfectly with cohesive energy ( $E_{coh}$ ), with W (highest $E_{coh}$ ) being the most resistant. This contradicts a simple thermodynamic driving force model, as the estimated bulk driving force is actually largest for W. The results suggest that resistance tracks the ease of local atomic rearrangement in more weakly bound metals rather than the thermodynamic driving force.
- Early Order Predictors: The initial bcc-like fraction at a depth of 0.1 nm ( $\phi_{0.1}$ ) strongly predicts both the sharpness of the depth response ( $n_h$ ) and the half-transformation depth ( $h_{50}$ ).
- Activation Work: A dimensionless activation work ( $\Pi_W$ ) based on bulk modulus and indentation depth orders the series as Mo < V < Ta < Nb < W. The peak persistent nucleus density ( $\Pi_n$ ) scales with activation work as $\Pi_n \propto \Pi_W^{-0.479}$ , statistically consistent with a $-1/2$ exponent.
Microstructural Evolution:
The terminal microstructure is governed by a population balance between nucleation density and coalescence. Elements like Mo and V exhibit dense persistent populations with strong coalescence, whereas Ta and Nb form fewer grains that grow to larger terminal sizes. The final grain size does not strictly follow the CFA ranking.

Significance and Claims
The paper establishes a common, normalized basis for comparing stress-driven devitrification across bcc refractory metals. It demonstrates that:

Mechanical loading can drive multistage ordering in single-element glasses, distinct from thermal evolution.
The transformation is a shear-biased process where resistance is governed by cohesive bond strength rather than thermodynamic driving force.
The operational CFA metric, combined with early structural order and dimensionless activation work, provides a consistent framework for linking early structure, mechanical loading, and microstructure evolution.

The authors explicitly note that these findings are specific to the five-element series and the interatomic potentials used. The scaling relations are presented as internally consistent empirical regularities for this specific set of potentials and conditions, rather than universal laws derived from first principles. The study acknowledges limitations regarding the ultrafast quench rates and high indentation velocities inherent to MD simulations, which differ significantly from experimental timescales.