Screw Dislocation: Theory & Applications

• Screw dislocation is a crystalline defect defined by a parallel Burgers vector and a helical lattice shift, crucial for understanding material behavior.
• Continuum, atomistic, and gauge models show that core regularization and boundary effects determine energy, stress fields, and mobility.
• Screw symmetry enables unique electronic, phononic, and optical modes, offering insights for advanced materials and device applications.

A screw dislocation is a one-dimensional crystalline defect for which the Burgers vector is parallel to the line direction. In anti-plane elasticity, for a straight defect along $e_3$, the only nonzero displacement component is the out-of-plane field $u_z=(b/2\pi)\theta$, with $u_r=u_\theta=0$, so that one circuit around the core produces a net shift $\Delta z=b$. This Volterra description remains the canonical far-field model, but contemporary treatments extend it to lattice variational theories, gauge and gradient regularizations, nonlinear Riemann–Cartan geometry, and symmetry-based descriptions in electronic, phononic, and optical media (Morandotti, 2017, Kobayashi et al., 2024, Sakaguchi et al., 2024).

1. Topological and kinematic characterization

In antiplane shear, the elastic field of a system of screw dislocations is represented by a planar strain field $h=\nabla u\in\mathbb R^2$ satisfying

$$\mathrm{div}\,h=0,\qquad \mathrm{curl}\,h=\sum_{i=1}^n b_i\,\delta_{z_i},$$

where each scalar $b_i\in\{\pm1\}$ is the Burgers modulus and $b_i e_3$ is the Burgers vector of the $i$th line defect. Equivalently, for any loop $\gamma_i$ enclosing exactly one core,

$u_z=(b/2\pi)\theta$0

This circulation law is the two-dimensional expression of the screw character: the defect is encoded not by a smooth displacement potential but by nontrivial winding around the core (Morandotti, 2017).

The continuum Volterra field can be written in cylindrical coordinates as

$u_z=(b/2\pi)\theta$1

The same defect may be represented geometrically by the metric

$u_z=(b/2\pi)\theta$2

which encodes the cut-and-glue operation of a screw dislocation. In this formulation, the single nonzero elastic strain component is $u_z=(b/2\pi)\theta$3 (Mashhadi et al., 2010).

Algebraic formulations replace the continuum branch cut by a covering-space construction. The exact sequence

$u_z=(b/2\pi)\theta$4

is used to realize the defect as a pullback over the punctured plane; the resulting lifted section is an embedded helicoid whose vertical increment after one loop is exactly the Burgers vector. For simple cubic and body-centered cubic lattices, this viewpoint is implemented with free abelian groups and group rings, and the screw displacement is the logarithm of a lattice phase winding (Hamada et al., 2016).

These equivalent descriptions emphasize that a screw dislocation is simultaneously a mechanical singularity, a topological winding defect, and a symmetry defect. That shared structure underlies the otherwise disparate continuum, atomistic, and wave-transport formulations discussed below.

2. Continuum field theories and core regularization

Classical elasticity gives the correct long-range $u_z=(b/2\pi)\theta$5 behavior but leaves the core singular. One major direction therefore replaces the Volterra defect by field theories with intrinsic lengths. In the translation gauge theory of dislocations, the elastic distortion and torsion are

$u_z=(b/2\pi)\theta$6

with Bianchi identity $u_z=(b/2\pi)\theta$7. For a screw dislocation in a functionally graded material whose moduli vary exponentially in $u_z=(b/2\pi)\theta$8 and $u_z=(b/2\pi)\theta$9, consistency requires

$u_r=u_\theta=0$0

which defines the intrinsic gauge length $u_r=u_\theta=0$1. The anti-plane distortion then satisfies the perturbed Helmholtz equation

$u_r=u_\theta=0$2

and the torsion obeys the same operator with source $u_r=u_\theta=0$3. The explicit dislocation density becomes

$u_r=u_\theta=0$4

so the core is nonsingular, the full cylindrical symmetry is broken by gradation, and the fields decay faster in the direction of increasing moduli (Lazar, 2011).

A related regularization arises in second-gradient elasticity for exponentially graded media. There, the anti-plane displacement is written as $u_r=u_\theta=0$5, and $u_r=u_\theta=0$6 solves

$u_r=u_\theta=0$7

with two intrinsic lengths $u_r=u_\theta=0$8 derived from $u_r=u_\theta=0$9. The resulting closed-form fields contain modified Bessel functions $\Delta z=b$0 and $\Delta z=b$1; their short-range contributions cancel the classical $\Delta z=b$2 divergence, while the limit $\Delta z=b$3 recovers the Volterra solution (Davoudi et al., 2010).

A nonlinear geometric formulation uses a Riemann–Cartan manifold $\Delta z=b$4, a plastic frame $\Delta z=b$5, and the torsion relation

$\Delta z=b$6

For an axisymmetric screw density, $\Delta z=b$7. A Helmholtz decomposition of $\Delta z=b$8 and a variational solution of $\Delta z=b$9, followed by elastic embedding into $h=\nabla u\in\mathbb R^2$0, produce stress fields that are finite at $h=\nabla u\in\mathbb R^2$1 and converge to the Volterra $h=\nabla u\in\mathbb R^2$2 field outside the core. In this setting, the Ricci curvature determines the symmetry of the stress field, and the work identifies a stress–curvature duality for screw dislocations (Kobayashi et al., 2024).

Across these formulations, the recurrent result is that the core singularity is not mandatory. It is a feature of the Volterra idealization, not of screw dislocations per se.

3. Renormalized energetics, Peach–Koehler forces, and boundary effects

In bounded domains, screw-dislocation energetics are commonly defined by removing small cores of radius $h=\nabla u\in\mathbb R^2$3 and passing to a renormalized limit. For dislocations at $h=\nabla u\in\mathbb R^2$4, the finite part of the energy can be expressed through the Green function $h=\nabla u\in\mathbb R^2$5 as

$h=\nabla u\in\mathbb R^2$6

with $h=\nabla u\in\mathbb R^2$7. The Peach–Koehler force is then

$h=\nabla u\in\mathbb R^2$8

This separates one-body self-energy, pair interactions, and boundary-image effects in a mathematically explicit way (Morandotti, 2017).

For free boundaries, the leading asymptotic force is attractive. If one dislocation $h=\nabla u\in\mathbb R^2$9 lies at distance $$\mathrm{div}\,h=0,\qquad \mathrm{curl}\,h=\sum_{i=1}^n b_i\,\delta_{z_i},$$0 from $$\mathrm{div}\,h=0,\qquad \mathrm{curl}\,h=\sum_{i=1}^n b_i\,\delta_{z_i},$$1, with nearest boundary point $$\mathrm{div}\,h=0,\qquad \mathrm{curl}\,h=\sum_{i=1}^n b_i\,\delta_{z_i},$$2, then

$$\mathrm{div}\,h=0,\qquad \mathrm{curl}\,h=\sum_{i=1}^n b_i\,\delta_{z_i},$$3

so the dominant term points along the outward normal and has magnitude $$\mathrm{div}\,h=0,\qquad \mathrm{curl}\,h=\sum_{i=1}^n b_i\,\delta_{z_i},$$4. Under overdamped dynamics $$\mathrm{div}\,h=0,\qquad \mathrm{curl}\,h=\sum_{i=1}^n b_i\,\delta_{z_i},$$5, this yields a collision-time estimate

$$\mathrm{div}\,h=0,\qquad \mathrm{curl}\,h=\sum_{i=1}^n b_i\,\delta_{z_i},$$6

for a dislocation initially at boundary distance $$\mathrm{div}\,h=0,\qquad \mathrm{curl}\,h=\sum_{i=1}^n b_i\,\delta_{z_i},$$7. For two nearby dislocations of opposite Burgers vector, the pair-collision time satisfies

$$\mathrm{div}\,h=0,\qquad \mathrm{curl}\,h=\sum_{i=1}^n b_i\,\delta_{z_i},$$8

and, in the leading-order regime discussed in the companion note, scales as $$\mathrm{div}\,h=0,\qquad \mathrm{curl}\,h=\sum_{i=1}^n b_i\,\delta_{z_i},$$9 (Morandotti, 2017, Morandotti, 2017).

The same framework also yields confinement results under prescribed boundary data. When one imposes tangential boundary conditions $b_i\in\{\pm1\}$0 with total circulation $b_i\in\{\pm1\}$1, the limiting finite energy

$b_i\in\{\pm1\}$2

diverges as any $b_i\in\{\pm1\}$3, where $b_i\in\{\pm1\}$4 is the minimum distance to the boundary or to another core. Consequently, minimizers stay strictly in the interior and remain well separated (Morandotti, 2017).

Zero-temperature lattice models sharpen this picture at the atomistic level. A variational anti-plane model on the triangular lattice establishes existence of a global minimizer with net Burgers vector $b_i\in\{\pm1\}$5, despite non-coercivity of the energy in the displacement variable. A related analysis proves locally stable multi-dislocation states provided dislocations are sufficiently far from one another and from the boundary, and shows that linear elasticity remains accurate beyond a few lattice spacings from the cores (Hudson et al., 2013, Hudson et al., 2014).

4. Core structure, cross-slip, and mobility in crystalline metals

In HCP Mg, the perfect screw $b_i\in\{\pm1\}$6 dislocation dissociates on the basal plane into two Shockley partials bounding an $b_i\in\{\pm1\}$7-type stacking fault of width $b_i\in\{\pm1\}$8 Å in DFT or $b_i\in\{\pm1\}$9 Å in EAM. Density-functional calculations give two nearly degenerate basal-dissociated cores and a finite-size-scaled cross-slip barrier

$b_i e_3$0

A line-tension model based on the DFT Peierls potential gives a jog-pair formation enthalpy $b_i e_3$1–$b_i e_3$2 eV at zero applied $b_i e_3$3, or $b_i e_3$4–$b_i e_3$5 eV at $b_i e_3$6 MPa. Molecular dynamics further shows jerky prismatic glide at low temperature, the appearance of plateau segments above $b_i e_3$7 K, and a $b_i e_3$8–$b_i e_3$9 reduction of cross-slip CRSS when a $i$0-Å hard sphere is placed on the basal plane at $i$1 K and $i$2 K (Itakura et al., 2015).

In fcc Al, real-space orbital-free DFT indicates that the energetic core size of a perfect screw is $i$3, substantially larger than the $i$4–$i$5 estimate obtained from displacement fields alone. Upon relaxation, the perfect screw dissociates into two Shockley partials with separations $i$6 Å and $i$7 Å. The relaxed core energy per unit length is nearly linear in applied strain,

$i$8

with $i$9 eV/Å and $\gamma_i$0 eV/Å for equi-triaxial volumetric strain. In discrete dislocation networks, the resulting core forces include $\gamma_i$1, $\gamma_i$2, and line-tension-type terms, and the type-I contribution can remain at least $\gamma_i$3 of the Peach–Koehler force out to $\gamma_i$4–$\gamma_i$5 nm (Das et al., 2017).

For the $\gamma_i$6 screw dislocation in $\gamma_i$7-Fe, ab initio calculations resolve an additional core field beyond the Volterra screw field. Its symmetry is that of a biaxial dilatation with

$\gamma_i$8

Including this field yields converged core energies $\gamma_i$9 meV/Å for the easy core and $u_z=(b/2\pi)\theta$00 meV/Å for the hard core when $u_z=(b/2\pi)\theta$01 Å. The same core field strengthens short-range interactions, increasing the dipole passing stress by up to $u_z=(b/2\pi)\theta$02 for $u_z=(b/2\pi)\theta$03 Å and increasing carbon binding by $u_z=(b/2\pi)\theta$04 eV for sites within $u_z=(b/2\pi)\theta$05–$u_z=(b/2\pi)\theta$06 Å of the dislocation (Clouet et al., 2011).

Electronic-structure effects also control alloying trends in bcc Mo. For a/2$u_z=(b/2\pi)\theta$07 screw cores, Ta and W raise both the complex formation energy and the zero-stress Peierls barrier, whereas Os, Ir, and Pt lower them. At the first-nearest-neighbor site, the reported values are $u_z=(b/2\pi)\theta$08 eV and $u_z=(b/2\pi)\theta$09 eV for Ta, $u_z=(b/2\pi)\theta$10 eV and $u_z=(b/2\pi)\theta$11 eV for W, $u_z=(b/2\pi)\theta$12 eV and $u_z=(b/2\pi)\theta$13 eV for Os, $u_z=(b/2\pi)\theta$14 eV and $u_z=(b/2\pi)\theta$15 eV for Ir, and $u_z=(b/2\pi)\theta$16 eV and $u_z=(b/2\pi)\theta$17 eV for Pt. The fitted stress-dependent barrier for pure Mo is

$u_z=(b/2\pi)\theta$18

with analogous parameter changes under solute addition (Zhou et al., 2024).

At larger scales, recent atomistic simulations in Nb, Mo, dilute Nb–Mo alloys, and equiatomic NbMo at $u_z=(b/2\pi)\theta$19 K over strain rates from $u_z=(b/2\pi)\theta$20 to $u_z=(b/2\pi)\theta$21 and line lengths from $u_z=(b/2\pi)\theta$22 to $u_z=(b/2\pi)\theta$23 nm show that cross-kinks form not only in concentrated alloys but also in pure BCC metals. High-rate depinning proceeds predominantly via vacancy–interstitial cluster formation, whereas low-rate, long-line configurations exhibit lateral cross-kink migration, three-dimensional forward–backward cross-slip, and prismatic loop formation (Chakraborty et al., 18 Feb 2026).

5. Screw symmetry as a generator of electronic, phononic, and optical states

In a three-dimensional strong $u_z=(b/2\pi)\theta$24 topological insulator, a screw dislocation can bind one-dimensional helical electronic modes. For weak indices $u_z=(b/2\pi)\theta$25, the Ran–Zhang–Vishwanath criterion is

$u_z=(b/2\pi)\theta$26

Thus, for Burgers vector $u_z=(b/2\pi)\theta$27, odd $u_z=(b/2\pi)\theta$28 gives gapless dislocation modes and even $u_z=(b/2\pi)\theta$29 does not. In the continuum Dirac description, the defect is encoded by the twisted boundary condition

$u_z=(b/2\pi)\theta$30

and the bound states have dispersion $u_z=(b/2\pi)\theta$31. Transport calculations separate $u_z=(b/2\pi)\theta$32; under periodic boundary conditions in $u_z=(b/2\pi)\theta$33 and $u_z=(b/2\pi)\theta$34, odd $u_z=(b/2\pi)\theta$35 gives $u_z=(b/2\pi)\theta$36 at low $u_z=(b/2\pi)\theta$37, while even $u_z=(b/2\pi)\theta$38 gives $u_z=(b/2\pi)\theta$39. A Zeeman gap on the side surfaces suppresses $u_z=(b/2\pi)\theta$40 without affecting $u_z=(b/2\pi)\theta$41 (Sakaguchi et al., 2024).

A closely related symmetry mechanism appears in phononics, but without requiring a nontrivial bulk topological invariant. In an HCP phononic crystal with a $u_z=(b/2\pi)\theta$42 screw dislocation, the defect preserves a local screw symmetry

$u_z=(b/2\pi)\theta$43

Modes can be labeled by Bloch momentum and screw eigenvalue, with the selection condition $u_z=(b/2\pi)\theta$44. Within the bulk band gap $u_z=(b/2\pi)\theta$45–$u_z=(b/2\pi)\theta$46 kHz, exactly two defect modes appear for each $u_z=(b/2\pi)\theta$47, with angular indices $u_z=(b/2\pi)\theta$48. Their transverse displacement fields have $u_z=(b/2\pi)\theta$49 and $u_z=(b/2\pi)\theta$50 exactly $u_z=(b/2\pi)\theta$51 out of phase, and near $u_z=(b/2\pi)\theta$52 they are described by

$u_z=(b/2\pi)\theta$53

Propagation direction and helicity are therefore locked by the screw symmetry itself (Zhou et al., 2024).

In amorphous optical media, the same defect geometry acts as an effective gauge field for paraxial propagation. With metric $u_z=(b/2\pi)\theta$54, the paraxial envelope obeys

$u_z=(b/2\pi)\theta$55

equivalently with gauge potential $u_z=(b/2\pi)\theta$56. The orbital angular momentum shifts as $u_z=(b/2\pi)\theta$57, allowing vortex annihilation or generation, while the elasto-optic birefringence induces polarization precession and an optical Hall effect (Mashhadi et al., 2010).

In GaN, exact restoration of the screw-dislocation group algebra yields band-connectivity and optical-selection rules. For the $u_z=(b/2\pi)\theta$58 threading screw, the eigenvalues of $u_z=(b/2\pi)\theta$59 satisfy

$u_z=(b/2\pi)\theta$60

and the band-flow rule under $u_z=(b/2\pi)\theta$61 is $u_z=(b/2\pi)\theta$62. Electric-dipole transitions obey

$u_z=(b/2\pi)\theta$63

The computed recombination coefficients are $u_z=(b/2\pi)\theta$64 and $u_z=(b/2\pi)\theta$65 at $u_z=(b/2\pi)\theta$66 eV and $u_z=(b/2\pi)\theta$67 K, with $u_z=(b/2\pi)\theta$68 by four orders of magnitude. The reported mechanism is a spiraling piezoelectric potential that spatially separates electrons and holes and suppresses radiative recombination (Xie et al., 27 Jan 2026).

These results directly contradict the common reduction of screw dislocations to purely mechanical defects. In the cited systems, the defect line is a symmetry-bearing channel that can bind, filter, or suppress excitations.

6. Unifying themes across scales and models

Several themes recur across otherwise different theories. First, the core is environment-dependent rather than universal. Gauge theory introduces $u_z=(b/2\pi)\theta$69 and an effective inverse length $u_z=(b/2\pi)\theta$70 in graded media; second-gradient elasticity introduces $u_z=(b/2\pi)\theta$71 and $u_z=(b/2\pi)\theta$72; nonlinear geometry replaces a Dirac core by a smooth density $u_z=(b/2\pi)\theta$73; atomistic calculations in Al assign an energetic core radius of $u_z=(b/2\pi)\theta$74; and $u_z=(b/2\pi)\theta$75-Fe requires an additional biaxial core field beyond the Volterra solution (Lazar, 2011, Davoudi et al., 2010, Das et al., 2017, Clouet et al., 2011).

Second, the far field is robust while the near field is model-sensitive. The variational lattice analyses show that linear elasticity is essentially correct outside a few lattice spacings, and the nonlinear geometric construction recovers the Volterra $u_z=(b/2\pi)\theta$76 field away from the core. By contrast, cross-slip barriers, strain-dependent core energies, short-range solute binding, and symmetry-filtered waveguiding are all controlled by the detailed core structure (Hudson et al., 2014, Kobayashi et al., 2024, Itakura et al., 2015).

Third, screw symmetry itself is operational. In phononics it protects $u_z=(b/2\pi)\theta$77 helical modes; in topological insulators it determines whether a dislocation binds conducting channels through $u_z=(b/2\pi)\theta$78; in GaN it organizes band connectivity and dipole selection rules; and in optics it enters as an Aharonov–Bohm-type phase shift. A plausible implication is that screw dislocations should be treated not only as defects to be regularized, but also as one-dimensional symmetry objects that can be exploited or mitigated depending on context (Zhou et al., 2024, Sakaguchi et al., 2024, Xie et al., 27 Jan 2026).

Finally, boundaries and external fields are not secondary perturbations. Free surfaces attract screw dislocations; prescribed boundary data can confine them; applied strain modifies core energy in Al; hard-sphere obstacles catalyze cross-slip in Mg at low temperature; and ferromagnetic coating can isolate dislocation transport in a topological insulator by gapping side-surface states (Morandotti, 2017, Das et al., 2017, Itakura et al., 2015, Sakaguchi et al., 2024).

Within this modern view, the screw dislocation is best understood as a topological line defect whose mechanics, energetics, and transport signatures are inseparable from intrinsic length scales, local symmetry, and environmental coupling.