🔬 physics

Competition between dislocation emission and cleavage at a blunt notch from the exact elastic image force on a mixed dislocation

This paper derives an exact closed-form expression for the elastic image force on a mixed dislocation near a blunt elliptic notch to construct a material-specific phase diagram that determines whether a crystal fails by cleavage or dislocation emission based on the interplay between its intrinsic properties and the flaw's tip radius.

Original authors: Anirudh Sai Ram Bommidi, Anirudh Udupa

Published 2026-10-05

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Anirudh Sai Ram Bommidi, Anirudh Udupa

Original paper licensed under CC BY 4.0 (https://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

In the microscopic world of solid materials, there is a fundamental choice that determines whether a piece of metal or crystal will bend and stretch or shatter into pieces. This decision happens at the very tip of a tiny flaw, a microscopic crack or notch hidden within the material. When stress is applied, the material faces a fork in the road. One path leads to cleavage, a brittle fracture where atomic bonds snap cleanly, causing the material to break apart suddenly. The other path leads to dislocation emission, a process where the material yields by sliding layers of atoms past one another, allowing it to deform plastically without breaking. For decades, scientists have understood this competition by imagining the flaw as a perfectly sharp point, an idealized mathematical edge where the stress becomes infinite. However, real-world flaws are never perfectly sharp; they always have a tiny, rounded tip, much like the rounded end of a needle rather than a razor's edge. The shape of this tip matters, yet the rules governing whether a material bends or breaks have traditionally been written for the sharp ideal, leaving a gap in our understanding of how real, rounded flaws behave.

Researchers at the Indian Institute of Technology Madras have now closed that gap by calculating exactly how the roundness of a flaw tip changes the outcome. They replaced the theoretical sharp crack with a finite, rounded void and traced the forces acting on a single defect within the crystal lattice. Their work reveals that the competition between bending and breaking is not a fixed property of the material alone but is a joint property of the material and the sharpness of the flaw it contains. By deriving a precise mathematical description of the forces at play, they mapped out exactly when a material will switch from being brittle to being ductile as the tip of a flaw becomes rounder.

The study focused on six different crystals, including common metals like copper and iron, as well as silicon and salt. The researchers calculated the balance of forces at the tip of a rounded flaw. On one side, there is the force trying to pull the material apart, which grows stronger as the tip becomes sharper. On the other side, there is a force that encourages the material to slide and deform, which depends on the specific arrangement of atoms in the crystal. The team found that for some materials, like iron, the balance shifts dramatically as the tip rounds out. When the tip of an iron flaw is very sharp, the material snaps cleanly. But if the tip is rounded just a little bit more—specifically, if the radius of the curve is about twice the width of a single atom's step—the material suddenly prefers to bend and slide instead of breaking. This critical point, where the behavior flips, was calculated to be at a radius of 2.03 times the size of the atomic step for iron.

For other materials, the story is different. Copper and sodium were found to be ductile, or bendable, even when the flaws are quite sharp, provided the tip is larger than the size of a single atom. Their calculations showed that these materials would only break if the flaw were sharper than the atomic scale itself, a condition that rarely exists in the real world. In contrast, materials like silicon, lithium fluoride, and beryllium remained brittle across the entire range of tip sizes the researchers examined. No matter how much the tip was rounded, these crystals preferred to snap rather than bend. The researchers also noted that for beryllium, the geometry of the atomic layers is so specific that rounding the tip actually makes it harder for the material to slide, keeping it firmly in the brittle category.

This work provides a complete, analytical map of how fracture behavior changes with the sharpness of a flaw. It moves beyond the old assumption that a material is simply "brittle" or "ductile" and shows that these labels depend on the geometry of the damage. The researchers did not rely on computer simulations that approximate the answer; instead, they derived an exact solution that works for any degree of rounding. They confirmed their results by checking them against known solutions for flat surfaces and by ensuring the math held up under extreme conditions. The findings suggest that the classic view of material failure is only a special case of a broader reality. A material's fate is not written in stone; it is written in the shape of the crack. By understanding the precise radius at which a material switches from breaking to bending, scientists can better predict how structures will fail and perhaps design materials that are more resistant to catastrophic fracture. The study concludes that the sharp-crack classification is just one edge of a much larger landscape, where the interplay between the crystal's internal structure and the flaw's external shape dictates whether a material will hold together or fall apart.

Technical Summary: Competition between Dislocation Emission and Cleavage at a Blunt Notch

Problem Statement
The ductile–brittle transition in crystalline materials is classically determined by a competition between two load thresholds: the stress required for Griffith cleavage and the stress required for dislocation emission from a crack tip. This classification, established by Rice and Thomson, relies on the idealization of an atomically sharp crack. However, real flaws possess finite tip radii ( $\rho$ ), which render the stress concentration finite rather than singular. Consequently, the force balance governing dislocation emission becomes a function of notch acuity. While prior numerical studies have addressed finite tip radii, a closed-form analytical resolution of the competition between emission and cleavage at a blunt notch—specifically accounting for the exact elastic image force on a mixed dislocation—has been lacking. This work addresses the fundamental question of whether the sharp-crack ductile–brittle classification is robust to crack-tip radius and, if not, identifies the critical tip radius at which a crystal's fracture character changes.

Methodology
The study replaces the atomically sharp crack with a traction-free elliptic void of semi-major axis $a_e$ and semi-minor axis $b_e$ , defining a tip radius $\rho = b_e^2/a_e$ . The analysis employs the complex-potential and conformal-mapping framework of Muskhelishvili.

Exact Elastic Image Force: The core technical advancement is the derivation of the exact elastic image force on a straight mixed dislocation near the elliptic void. The perturbation construction generates three infinite series for the edge component of the image force. Unlike previous approaches that truncate these series, the authors sum them in exact closed form, yielding finite rational expressions (Eqs. 29–31). This eliminates truncation error entirely and reduces the evaluation cost from $O(N)$ to $O(1)$ for any exterior dislocation position and elliptic aspect ratio. The screw component is derived independently in closed form.
Force Balance and Criteria: The total glide force is the sum of the exact image force and the resolved remote-tension field (evaluated exactly, not via $K$ $K$ -field asymptotes).
- Emission Criterion: A discrete core-cutoff criterion is retained. Emission is deemed spontaneous (athermal) when the unstable equilibrium distance $d_{eq}$ falls within the material-specific core cutoff radius $r_c$ (where $r_c$ is a multiple of the Burgers vector $b$ ). The emission threshold $K_{Id}(\rho)$ is defined as the load where $d_{eq} = r_c$ .
- Cleavage Criterion: Cleavage is governed by the Creager–Paris leading-order opening stress at the notch tip, $\sigma_{max} = 2K_I/\sqrt{\pi\rho}$ . The cleavage threshold $K_{Ic}(\rho)$ is reached when this stress attains the theoretical strength $\sigma_{th} = \sqrt{2E\gamma_s/b}$ .
Material Parameters: The framework is applied to six crystals (Cu, Na, LiF, Si, Be, Fe) using crystallographic data (cleavage and slip systems) and elastic constants from Rice and Thomson. The analysis accounts for the specific edge–screw decomposition of the Burgers vector for each system.

Key Contributions

Closed-Form Solution: The paper provides the first closed-form analytical resolution of the emission–cleavage competition at a finite tip radius. By summing the perturbation series for the edge image force exactly, it avoids the numerical truncation errors that plague series-based methods near the void surface.
Notch-Acuity Phase Diagram: The work constructs a phase diagram in the $(\text{material}, \rho/b)$ space. It demonstrates that fracture mode is a joint property of the crystal and the acuity of the flaw, rather than a property of the crystal alone.
Analytical Efficiency: The closed-form expressions allow for the computation of thresholds across a wide range of tip radii ( $\rho/b \in [0.05, 100]$ ) for multiple materials without solving singular integral equations at each load step, a requirement of prior Peierls-based numerical treatments.

Results
The study computes the critical blunting radius $\rho^\star$ where the emission threshold $K_{Id}(\rho)$ and cleavage threshold $K_{Ic}(\rho)$ cross ( $K_{Id} = K_{Ic}$ ).

Iron (Fe): Exhibits a critical crossing at $\rho^\star = 2.03b$ . Iron is predicted to cleave below this radius and emit dislocations (ductile) above it.
Copper (Cu) and Sodium (Na): Show formal crossings at $\rho^\star = 0.68b$ and $0.98b$ , respectively. These values lie below the nominal continuum-resolution scale ( $\rho \sim b$ ). Consequently, the model predicts both materials to be ductile for all physically resolvable blunt notches ( $\rho \gtrsim b$ ).
LiF, Silicon (Si), and Beryllium (Be): No crossing occurs within the swept range ( $\rho/b \in [0.05, 100]$ ). These materials remain brittle across the entire range, although blunting increases the load required for failure.
Non-monotonicity: The emission threshold $K_{Id}(\rho)$ is non-monotonic for most crystals, passing through a minimum, whereas the cleavage threshold $K_{Ic}(\rho)$ rises monotonically as $\rho^{1/2}$ . This behavior creates the conditions for the crossing. Beryllium is an exception due to its orthogonal slip geometry, where $K_{Id}$ increases monotonically.

Significance and Claims
The authors claim that this work resolves the competition between dislocation emission and cleavage at a blunt notch in closed form, extending the classical Rice–Thomson sharp-crack classification. The paper asserts that the sharp-crack classification is recovered as a limiting edge (atomically sharp) of the derived phase diagram, but the $\rho \to 0$ limit of the current finite-radius criteria does not reproduce the Rice–Thomson result directly due to the different scaling of the cleavage threshold.

The significance lies in establishing that the ductile–brittle character of a material is not intrinsic but depends on the geometry of the flaw. Specifically, for materials like iron, blunting alone can switch the failure mode from brittle cleavage to ductile emission. The authors note that while the framework complements Peierls-based nucleation physics (which model the slip zone continuously), it retains the discrete core-cutoff criterion to enable exact finite-radius solutions. The results are presented as a "no-ledge" model, acknowledging that ledge forces could alter specific values (particularly for sodium) but do not change the fundamental geometric dependence of the competition. The work does not propose new experimental protocols but provides an analytical benchmark against which numerical and atomistic simulations can be compared.

Technical Summary: Competition between Dislocation Emission and Cleavage at a Blunt Notch

More like this