pith. sign in

arxiv: 2606.12591 · v1 · pith:66SBXBS4new · submitted 2026-06-10 · ❄️ cond-mat.mtrl-sci

Computationally efficient method for determining limiting velocities of edge dislocations in anisotropic crystals

Daniel N. Blaschke This is my paper

Pith reviewed 2026-06-27 08:51 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci
keywords edge dislocationslimiting velocitiesanisotropic crystalsreflection symmetrysecular equationelastic constantsdislocation glide
0
0 comments X

The pith

Reflection symmetry permits an algebraic reduction of the secular equation that gives exact limiting velocities for edge dislocations without numerical slowdown.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

Continuum dislocation theory predicts elastic-energy divergences at geometry-dependent limiting velocities that mark the boundaries between subsonic, transsonic, and supersonic glide. For edge dislocations that possess reflection symmetry yet carry nonzero c'16 or c'26 elastic constants, earlier numerical schemes for locating these velocities became impractically slow. The paper derives a method that exploits the reflection symmetry to reduce the secular equation algebraically. The reduced equation retains the precise roots of the original secular equation while removing the sources of numerical cost and instability. The resulting procedure therefore supplies the velocities required by high-strain-rate strength models in a computationally tractable way.

Core claim

For edge dislocations with reflection symmetry, the secular equation can be algebraically reduced to a form that permits efficient and exact determination of the limiting velocities v_L without additional approximations or numerical instabilities.

What carries the argument

Algebraic reduction of the secular equation that exploits reflection symmetry of the edge-dislocation geometry.

If this is right

  • The subsonic, transsonic, and supersonic glide regimes become directly computable for crystals previously treated as intractable.
  • High-strain-rate material-strength models can incorporate the exact limiting velocities without prohibitive runtime cost.
  • The reduction applies whenever the edge-dislocation geometry satisfies reflection symmetry, independent of the specific values of the elastic constants.
  • No new approximations are introduced, so prior analytic results for simpler symmetries remain unchanged.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same symmetry reduction may apply to other dislocation types or to wave-propagation problems that share analogous secular equations.
  • Implementation inside molecular-dynamics or finite-element codes would allow real-time evaluation of glide regimes during high-rate deformation simulations.
  • A systematic survey of common crystal classes could identify which ones now become tractable for limiting-velocity studies.

Load-bearing premise

Reflection symmetry of the dislocation geometry allows an exact algebraic reduction of the secular equation.

What would settle it

Apply both the new algebraic method and a standard numerical root finder to the same anisotropic crystal with reflection symmetry and nonzero c'16; the two must return identical velocity values to machine precision while the algebraic route uses far fewer operations.

read the original abstract

The continuum-limit theory of dislocations in crystals predicts divergences in the elastic energy at crystal-geometry dependent limiting velocities $v_L$, which separate subsonic, transsonic, and supersonic dislocation glide regimes and are therefore import for material strength models at high strain rates. Although it is known how to calculate those limiting velocities, there is one special case - edge dislocations with reflection symmetry, but non-vanishing elastic constants $c'_{16}$ or $c'_{26}$ - where previous methods have been notoriously slow. In this letter, we address this deficiency by deriving a computationally efficient method for determining the limiting velocities of edge dislocations with reflection symmetry.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript derives a computationally efficient method for calculating the limiting velocities v_L of edge dislocations possessing reflection symmetry in anisotropic crystals. It focuses on the previously slow case where c'_16 or c'_26 are nonzero, by performing an algebraic reduction of the secular equation that is asserted to recover the exact roots of the unreduced problem without approximations or numerical instabilities.

Significance. If the reduction is shown to be exact, the result would remove a known computational bottleneck in high-strain-rate dislocation modeling and material-strength calculations. The explicit claim of parameter-free exactness (no fitted parameters or hidden approximations) is a strength that, if substantiated by direct comparison to the full secular determinant, would make the method immediately usable in existing anisotropic elasticity codes.

major comments (2)
  1. [derivation of the reduced secular equation] The central claim that the algebraic reduction preserves the exact set of real positive roots when c'_16 or c'_26 are nonzero must be demonstrated explicitly. The manuscript should show, either analytically or by direct numerical comparison for several values of these constants, that the reduced secular equation yields identical roots to the unreduced 6×6 (or 3×3) determinant; any missed coupling term would either drop legitimate limiting velocities or introduce extraneous ones.
  2. [reflection symmetry reduction] The reflection-symmetry assumption is used to decouple displacement components, yet the handling of the c'_16 and c'_26 terms that couple those components is not shown to leave the characteristic polynomial unchanged. A concrete step-by-step substitution into the Stroh matrix (or equivalent) followed by determinant evaluation is required to confirm that no roots are lost or added.
minor comments (2)
  1. The abstract states the method is 'computationally efficient' but provides no timing benchmarks or operation-count comparison against the standard method; a short table or sentence quantifying the speedup would strengthen the letter.
  2. Notation for the rotated elastic constants (c'_ij) should be defined at first use and kept consistent with standard conventions in anisotropic elasticity literature.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading and for highlighting the importance of explicit verification of the reduction. We address the two major comments below and have revised the manuscript accordingly to include the requested demonstrations.

read point-by-point responses

  1. Referee: [derivation of the reduced secular equation] The central claim that the algebraic reduction preserves the exact set of real positive roots when c'_16 or c'_26 are nonzero must be demonstrated explicitly. The manuscript should show, either analytically or by direct numerical comparison for several values of these constants, that the reduced secular equation yields identical roots to the unreduced 6×6 (or 3×3) determinant; any missed coupling term would either drop legitimate limiting velocities or introduce extraneous ones.

    Authors: We agree that the original manuscript asserted exactness via algebraic reduction but did not include explicit verification. The revised manuscript now contains both an analytical proof (via direct substitution showing the reduced polynomial factors identically from the full secular determinant) and numerical comparisons for several nonzero values of c'_16 and c'_26. These comparisons confirm that all real positive roots match to machine precision, with no dropped or extraneous roots introduced. revision: yes

  2. Referee: [reflection symmetry reduction] The reflection-symmetry assumption is used to decouple displacement components, yet the handling of the c'_16 and c'_26 terms that couple those components is not shown to leave the characteristic polynomial unchanged. A concrete step-by-step substitution into the Stroh matrix (or equivalent) followed by determinant evaluation is required to confirm that no roots are lost or added.

    Authors: The revised manuscript now includes a detailed step-by-step substitution of the reflection-symmetric Stroh matrix elements, explicitly retaining the c'_16 and c'_26 coupling terms through the determinant evaluation. This shows that the characteristic polynomial for the relevant roots is unchanged, confirming that no legitimate limiting velocities are lost and no extraneous roots are added. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained from standard anisotropic elasticity

full rationale

The paper derives an algebraic reduction of the secular equation for limiting velocities of edge dislocations exploiting reflection symmetry. This is presented as a direct consequence of continuum anisotropic elasticity theory without any fitted parameters, self-citations as load-bearing premises, or renaming of prior results. The central claim is that the reduction preserves exact roots without approximations, and no evidence in the provided abstract or description shows the output reducing to the input by construction. The method is independent of any self-referential steps.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract supplies no information on free parameters, background axioms, or new entities; all such quantities remain unknown.

pith-pipeline@v0.9.1-grok · 5631 in / 986 out tokens · 20357 ms · 2026-06-27T08:51:15.093962+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

39 extracted references · 4 linked inside Pith

  1. [1]

    A dislocation- based multi-rate single crystal plasticity model

    B. L. Hansen, I. J. Beyerlein, C. A. Bronkhorst, E. K. Cerreta, and D. Dennis-Koller, “A dislocation- based multi-rate single crystal plasticity model”,Int. J . Plast.44(2013) 129–146

    2013

  2. [2]

    Simulation of shock wave propagation in single crystal and polycrystalline aluminum

    J. T. Lloyd, J. D. Clayton, R. Becker, and D. L. McDowell, “Simulation of shock wave propagation in single crystal and polycrystalline aluminum”,Int. J . Plast.60(2014) 118–144

    2014

  3. [3]

    Coupling continuum dislocation transport with crystal plasticity for application to shock loading conditions

    D. J. Luscher, J. R. Mayeur, H. M. Mourad, A. Hunter, and M. A. Kenamond, “Coupling continuum dislocation transport with crystal plasticity for application to shock loading conditions”,Int. J . Plast. 76(2016) 111–129

    2016

  4. [4]

    Elastic precursor wave decay in shock-compressed aluminum over a wide range of tem- perature

    R. A. Austin, “Elastic precursor wave decay in shock-compressed aluminum over a wide range of tem- perature”,J . Appl. Phys.123(2018) 035103

    2018

  5. [5]

    Dynamic flow stress of pure polycrys- talline aluminum: Pressure-shear plate impact experiments and extension of dislocation-based model- ing to large strains

    B. Zuanetti, D. J. Luscher, K. Ramos, C. Bolme, and V. Prakash, “Dynamic flow stress of pure polycrys- talline aluminum: Pressure-shear plate impact experiments and extension of dislocation-based model- ing to large strains”,J . Mech. Phys. Solids146(2021) 104185

    2021

  6. [6]

    The mechanics and physics of high- speed dislocations: a critical review

    B. Gurrutxaga-Lerma, J. Verschueren, A. P. Sutton, and D. Dini, “The mechanics and physics of high- speed dislocations: a critical review”,Int. Mater . Rev.66(2021) 215–255

    2021

  7. [7]

    Relativistic effect inducing drag on fast-moving dislocation in discrete system

    S. Kim, H. Kim, K. Kang, and S. Y. Kim, “Relativistic effect inducing drag on fast-moving dislocation in discrete system”,Int. J . Plast.126(2020) 102629. 5

    2020

  8. [8]

    Direct comparison between experiments and dislocation dynamics simulations of high rate deformation of single crystal copper

    S. Akhondzadeh, M. Kang, R. B. Sills, K. Ramesh, and W. Cai, “Direct comparison between experiments and dislocation dynamics simulations of high rate deformation of single crystal copper”,Acta Materialia 250(2023) 118851

    2023

  9. [9]

    A discrete dislocation dynamics framework for modeling polycrystal plasticity with hardening

    T. N. Tak, A. Prakash, I. Samajdar, A. A. Benzerga, and P. Guruprasad, “A discrete dislocation dynamics framework for modeling polycrystal plasticity with hardening”,Int. J . Solids Struct.281(2023) 112442

    2023

  10. [10]

    Unified crystal plasticity model for fcc metals: From quasistatic to shock loading

    C. Ye, G. Liu, K. Chen, J. Liu, J. Hu, Y. Yu, Y. Mao, and Y. Shen, “Unified crystal plasticity model for fcc metals: From quasistatic to shock loading”,Phys. Rev.B107(2023) 024105

    2023

  11. [11]

    Learning dislocation dynamics mobility laws from large-scale MD simulations

    N. Bertin, V. V. Bulatov, and F. Zhou, “Learning dislocation dynamics mobility laws from large-scale MD simulations”,npj Comp. Mater .10(2024) 192

    2024

  12. [12]

    Dislocation activity in aluminum at ultra-high strain rates: Atomistic investigation and continuum modeling

    K. D. Manukhina, V. S. Krasnikov, D. S. Voronin, and A. E. Mayer, “Dislocation activity in aluminum at ultra-high strain rates: Atomistic investigation and continuum modeling”,Comput. Mater . Sci.244 (2024) 113269

    2024

  13. [13]

    How to determine limiting velocities of dislocations in anisotropic crystals

    D. N. Blaschke, “How to determine limiting velocities of dislocations in anisotropic crystals”,J . Phys.: Cond. Mat.33(2021) 503005,arXiv:2107.01220 [cond-mat.mtrl-sci]

    arXiv 2021

  14. [14]

    The transient motion of a dislocation with a ramp-like core

    X. Markenscoff and L. Ni, “The transient motion of a dislocation with a ramp-like core”,J . Mech. Phys. Solids49(2001) 1603–1619

    2001

  15. [15]

    Uniformly-moving non-singular dislocations with ellipsoidal core shape in anisotropic media

    Y.-P. Pellegrini, “Uniformly-moving non-singular dislocations with ellipsoidal core shape in anisotropic media”,J . Micromech. Molec. Phys.3(2018) 1840004,arXiv:1808.10272 [physics.class-ph]

    Pith/arXiv arXiv 2018

  16. [16]

    Quantum statistical theory of dislocation mobility in discrete lattices

    B. Gurrutxaga-Lerma, “Quantum statistical theory of dislocation mobility in discrete lattices”,Phys. Rev. Mater .9(2025) 123605

    2025

  17. [17]

    Dislocations faster than the speed of sound

    P. Gumbsch and H. Gao, “Dislocations faster than the speed of sound”,Science283(1999) 965–968

    1999

  18. [18]

    Atomistic simulations of dis- location mobility in Al, Ni and Al/Mg alloys

    D. L. Olmsted, L. G. Hector Jr., W. A. Curtin, and R. J. Clifton, “Atomistic simulations of dis- location mobility in Al, Ni and Al/Mg alloys”,Mod. Simul. Mater . Sci. Eng.13(2005) 371, arXiv:cond-mat/0412324

    Pith/arXiv arXiv 2005

  19. [19]

    Dislocation kinematics: a molecular dynamics study in Cu

    E. Oren, E. Yahel, and G. Makov, “Dislocation kinematics: a molecular dynamics study in Cu”,Mod. Simul. Mater . Sci. Eng.25(2017) 025002

    2017

  20. [20]

    Limiting velocities and transonic dislocations in Mg

    K. Dang, D. N. Blaschke, S. Fensin, and D. J. Luscher, “Limiting velocities and transonic dislocations in Mg”,Comput. Mater . Sci.215(2022) 111786,arXiv:2205.11687 [cond-mat.mtrl-sci]

    arXiv 2022

  21. [21]

    Resonance with surface waves induces forbidden velocity bands in dislocation glide

    T. Duong and M. J. Demkowicz, “Resonance with surface waves induces forbidden velocity bands in dislocation glide”,J . Mech. Phys. Solids180(2023) 105422

    2023

  22. [22]

    Exploring the relation between transonic dislocation glide and stacking fault width in FCC metals

    K. R. Jones, K. Dang, D. N. Blaschke, S. J. Fensin, and A. Hunter, “Exploring the relation between transonic dislocation glide and stacking fault width in FCC metals”,Mod. Simul. Mater . Sci. Eng.33 (2025) 025020,arXiv:2409.10705 [cond-mat.mtrl-sci]

    arXiv 2025

  23. [23]

    Transonic dislocation propagation in diamond

    K. Katagiri, T. Pikuz, L. Fang,et al., “Transonic dislocation propagation in diamond”,Science382(2023) 69–72,arXiv:2303.04370 [cond-mat.mtrl-sci]

    arXiv 2023

  24. [24]

    Observing formation and evolution of dislocation cells during plastic deformation

    A. Zelenika, A. A. Cretton, F. Frankus,et al., “Observing formation and evolution of dislocation cells during plastic deformation”,Sci. Rep.15(2025) 8655

    2025

  25. [25]

    Individual dislocation identification in dark-field X-ray mi- croscopy

    S. Borgi, G. Winther, and H. F. Poulsen, “Individual dislocation identification in dark-field X-ray mi- croscopy”,J . Appl. Crystallogr .58(2025) 813–821

    2025

  26. [26]

    Dislocation drag and its influence on elastic precursor decay

    D. N. Blaschke and D. J. Luscher, “Dislocation drag and its influence on elastic precursor decay”,Int. J . Plast.144(2021) 103030,arXiv:2101.10497 [cond-mat.mtrl-sci]

    arXiv 2021

  27. [27]

    Anisotropic continuum theory of lattice defects

    D. J. Bacon, D. M. Barnett, and R. O. Scattergood, “Anisotropic continuum theory of lattice defects”, Prog. Mater . Sci.23(1980) 51–262

    1980

  28. [28]

    J. P. Hirth and J. Lothe,Theory of Dislocations, second ed., (New York: Wiley, 1982)

    1982

  29. [29]

    Elastic surface waves in anisotropic crystals: a simplified method for calculating Rayleigh velocities using dislocation theory

    D. M. Barnett, J. Lothe, K. Nishioka, and R. J. Asaro, “Elastic surface waves in anisotropic crystals: a simplified method for calculating Rayleigh velocities using dislocation theory”,J . Phys. F: Met. Phys.3 (1973) 1083–1096. 6

    1973

  30. [30]

    Uniformly moving dislocations in anisotropic media

    R. Bullough and B. A. Bilby, “Uniformly moving dislocations in anisotropic media”,Proc. Phys. Soc.B67 (1954) 615–624

    1954

  31. [31]

    Dynamical behavior of dislocations in anisotropic media

    L. J. Teutonico, “Dynamical behavior of dislocations in anisotropic media”,Phys. Rev.124(1961) 1039– 1045

    1961

  32. [32]

    The transient motion of a screw dislocation in an anisotropic medium

    X. Markenscoff and L. Ni, “The transient motion of a screw dislocation in an anisotropic medium”,J . Elast.14(1984) 93–95

    1984

  33. [33]

    Nonuniform motion of an edge dislocation in an anisotropic solid. I

    X. Markenscoff and L. Q. Ni, “Nonuniform motion of an edge dislocation in an anisotropic solid. I”, Quart. Appl. Math.41(1984) 475–494

    1984

  34. [34]

    Clarifying the definition of ‘transonic’ screw dislocations

    D. N. Blaschke, J. Chen, S. Fensin, and B. Szajewski, “Clarifying the definition of ‘transonic’ screw dislocations”,Phil. Mag.101(2021) 997–1018,arXiv:2008.13760 [cond-mat.mtrl-sci]

    arXiv 2021

  35. [35]

    D. N. Blaschke,PyDislocDyn, 2018–2026, #C18073,URL:https://github.com/dblaschke-LANL/ PyDislocDyn

    2018

  36. [36]

    PyDislocDyn: A Python code for calculating dislocation drag and other crystal proper- ties

    D. N. Blaschke, “PyDislocDyn: A Python code for calculating dislocation drag and other crystal proper- ties”,JOSS10(2025) 9309,arXiv:2509.02900 [cond-mat.mtrl-sci]

    arXiv 2025

  37. [37]

    Line tension of a dislocation moving through an anisotropic crys- tal

    D. N. Blaschke and B. A. Szajewski, “Line tension of a dislocation moving through an anisotropic crys- tal”,Phil. Mag.98(2018) 2397–2424,arXiv:1711.10555 [cond-mat.mtrl-sci]

    Pith/arXiv arXiv 2018

  38. [38]

    The phonon-dislocation interaction and its role in dislocation dragging and thermal re- sistivity

    V. I. Alshits, “The phonon-dislocation interaction and its role in dislocation dragging and thermal re- sistivity”, inElastic Strain Fields and Dislocation Mobility, V. L. Indenbom and J. Lothe, eds., vol. 31 ofModern Problems in Condensed Matter Sciences, pp. 625–697, (Elsevier, 1992)

    1992

  39. [39]

    Velocity dependent dislocation drag from phonon wind and crystal geometry

    D. N. Blaschke, “Velocity dependent dislocation drag from phonon wind and crystal geometry”,J . Phys. Chem. Solids124(2019) 24–35,arXiv:1804.01586 [cond-mat.mtrl-sci]. 7

    Pith/arXiv arXiv 2019