Pith. sign in

REVIEW 1 major objections 1 minor 39 references

Energetic coupling between field dislocation mechanics and phase field crystal ignores the incompatible distortion encoding dislocation topology.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · grok-4.3

2026-07-03 19:53 UTC pith:7UCEBIRK

load-bearing objection L2 coupling between FDM and PFC only affects compatible distortions and stays blind to dislocation topology even in general energetic cases. the 1 major comments →

arxiv 2607.01284 v1 pith:7UCEBIRK submitted 2026-07-01 cond-mat.mtrl-sci

On the limits of the energetic coupling between field dislocation mechanics and phase field crystal

Aymane Graini , Jorge Vi\~nals , Manas V. Upadhyay This is my paper
classification cond-mat.mtrl-sci
keywords field dislocation mechanicsphase field crystalenergetic couplingdislocation topologyincompatible distortioncore spreadingvariational analysiscontinuum modeling
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

The paper analyzes the energetic coupling proposed to connect field dislocation mechanics, which handles continuum bodies with dislocations, to the phase field crystal model that captures crystallographic structure. It finds that the coupling, based on penalizing the L2 distance between elastic and configurational distortions, only affects the compatible parts of these fields via divergence-driven terms. This leaves the incompatible components, which carry dislocation topology information, untouched by the coupling. Mechanical boundary conditions thus transmit diffusively into the phase field model rather than elastically, and simulations show the coupling cannot stop unphysical spreading of dislocation cores. The limitation holds even for general energetic couplings.

Core claim

Variational analysis shows that the coupling term in the phase-field evolution equation is divergence-driven and matches only the compatible (curl-free) parts of the distortion fields from FDM and PFC. Its contributions are therefore insensitive to the incompatible (divergence-free) elastic distortion that contains all information on dislocation topology. The configurational nature of the PFC distortion further implies that mechanical boundary conditions are transmitted diffusively from FDM to PFC. Numerical simulations confirm that this coupling cannot prevent the unnatural core spreading in FDM, and the same drawbacks apply to any energetic coupling.

What carries the argument

L2 penalization of the difference between FDM elastic distortion and PFC configurational distortion, whose variational derivative produces a divergence-driven forcing term.

Load-bearing premise

The coupling energy is constructed as an L2 penalization between the FDM elastic distortion and the PFC configurational distortion, leading to a purely divergence-driven variational derivative.

What would settle it

A calculation or simulation where only the incompatible part of the FDM distortion is changed while the compatible part is held fixed, to check if the PFC phase field evolution remains unaffected.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • The coupling term is insensitive to incompatible distortion carrying dislocation topology.
  • Mechanical boundary conditions transmit diffusively rather than elastically.
  • The coupling cannot prevent unnatural core spreading in FDM.
  • These limitations persist in the most general energetic coupling.

Where Pith is reading between the lines

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

  • Direct coupling to incompatibility measures or curl components may be required to link dislocation topology.
  • Non-energetic or hybrid coupling strategies might overcome the diffusive transmission issue.
  • Similar limitations could affect other multi-scale models that rely on energy penalties between continuum and atomistic fields.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

1 major / 1 minor

Summary. The manuscript analyzes the energetic coupling between Field Dislocation Mechanics (FDM) and the Phase Field Crystal (PFC) model, originally proposed via an L² penalization between the FDM elastic distortion and PFC configurational distortion. Variational calculus shows that this term produces only divergence-driven forcing in the PFC evolution equation, coupling exclusively to the compatible (curl-free) distortion components while remaining insensitive to the incompatible (divergence-free) part that encodes dislocation topology. The analysis further concludes that mechanical boundary conditions are transmitted diffusively rather than elastically, that numerical simulations confirm the coupling cannot suppress unnatural core spreading in FDM, and that these limitations persist for arbitrary energetic couplings.

Significance. If the variational results and the extension to the general energetic case are rigorously established, the work provides a clear diagnostic of fundamental limitations in L²-style couplings for reconciling continuum dislocation mechanics with crystallographic models. This could usefully constrain future multiscale modeling efforts in materials science by highlighting the need for couplings that directly engage incompatible distortion fields.

major comments (1)
  1. [Abstract] Abstract (general-case claim): the assertion that 'even in the most general case, an energetic coupling suffers from the same drawbacks' is load-bearing for the central conclusion, yet the provided variational analysis is constructed around the specific L² inner-product penalization; an explicit derivation for arbitrary energetic functionals (without tacit retention of the L² structure or the assumption that the PFC distortion remains curl-free) is required to substantiate the extrapolation.
minor comments (1)
  1. The description of the numerical simulations omits quantitative error measures, convergence data, or explicit statements of boundary conditions, which would allow direct assessment of the core-spreading claim.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful review and constructive feedback. We address the single major comment below.

read point-by-point responses

  1. Referee: [Abstract] Abstract (general-case claim): the assertion that 'even in the most general case, an energetic coupling suffers from the same drawbacks' is load-bearing for the central conclusion, yet the provided variational analysis is constructed around the specific L² inner-product penalization; an explicit derivation for arbitrary energetic functionals (without tacit retention of the L² structure or the assumption that the PFC distortion remains curl-free) is required to substantiate the extrapolation.

    Authors: We acknowledge that the explicit variational calculations are performed for the L² penalization. The general-case statement rests on the structural observation that any local energetic coupling of the form ∫ F(β^FDM − β^PFC) dV produces, upon variation with respect to β^PFC, a forcing term whose leading contribution is the divergence of a tensor obtained from the functional derivative of F. Consequently the PFC evolution equation remains insensitive to the curl-free projection of the distortion difference, independent of the specific form of F. The curl-free character of β^PFC follows directly from its construction within the PFC model (as the symmetrized gradient of the phase field). To make this argument fully rigorous and remove any tacit reliance on the L² structure, we will insert an explicit derivation for a general differentiable functional F in a new subsection of the revised manuscript. revision: yes

Circularity Check

0 steps flagged

No significant circularity; direct variational analysis of stated energy functional

full rationale

The paper applies standard variational calculus to the L2 penalization energy between FDM elastic distortion and PFC configurational distortion (as defined in the 2020 reference being critiqued). The insensitivity to incompatible distortion follows from the divergence structure of the resulting Euler-Lagrange terms, and the general-case extension uses the same structure without introducing fitted parameters, self-referential uniqueness theorems, or ansatzes smuggled via citation. The 2020 citation supplies the object of analysis rather than load-bearing justification for the result. No step reduces the claimed prediction to an input by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The paper introduces no new free parameters or invented entities. It relies on standard variational calculus applied to an energy functional whose coupling term is taken from prior literature.

axioms (2)
  • standard math Variational calculus yields the correct evolution equations from the total energy functional that includes the L2 coupling term
    Invoked to obtain the divergence-driven forcing term in the phase-field equation.
  • domain assumption The configurational distortion field in the PFC model is of a form whose curl-free part can be directly compared to the FDM elastic distortion
    Required for the L2 penalization to be well-defined and for the boundary-condition transmission argument.

pith-pipeline@v0.9.1-grok · 5771 in / 1490 out tokens · 41084 ms · 2026-07-03T19:53:19.109884+00:00 · methodology

0 comments
read the original abstract

This paper investigates the energetic coupling between Field Dislocation Mechanics (FDM) and the Phase Field Crystal (PFC) model proposed in Phys. Rev. B 102, 064109, 2020. While FDM correctly solves the initial boundary value problem of a continuum body with dislocation fields, PFC captures the underlying crystallographic structure. The coupling, which penalizes the $L^2$ distance between elastic distortion from FDM and configurational distortion from PFC in the $L^2$ sense, had been proposed to reconcile dislocation mechanics with crystallography in a single continuum framework. Variational analysis reveals that the coupling term acts as a divergence-driven forcing in the phase-field evolution that matches only the compatible (curl-free) parts of the distortion fields. Consequently, its contributions are insensitive to the incompatible (divergence-free) elastic distortion carrying all the information on dislocation topology. Furthermore, the nature of the configurational distortion causes mechanical boundary conditions to be transmitted diffusively from FDM to PFC rather than elastically. Numerical simulations demonstrate that this coupling cannot prevent the unnatural core spreading in FDM. Finally, it is shown that even in the most general case, an energetic coupling suffers from the same drawbacks, which limits its ability to integrate dislocation mechanics with crystallography.

Figures

Figures reproduced from arXiv: 2607.01284 by Aymane Graini, Jorge Vi\~nals, Manas V. Upadhyay.

Figure 1
Figure 1. Figure 1: (a) Order parameter field in the presence of a single edge dislocation at the center of a 2D hexagonal lattice. (b) Corresponding dislocation density αxz. (d) Contour plot of the penalty-induced forcing for this configuration. In both, the gray circle shows the extent of the core. (c) A horizontal slice along y = 0 of the penalty forcing whereas (e) show the behavior of the xx of Q and U e . The vertical s… view at source ↗
Figure 2
Figure 2. Figure 2: Immobile edge dislocation dipole in Climb configuration. (a): Order parameter ψ at t = 0 and the dislocation density tensor α. (b) relaxation of ∥U e − Q∥ with increasing cpen. (c), corresponding evolution of ∥∇ · (U e − Q)∥. At t = 0, in the absence of coupling (cpen = 0), the elastic distortion U e and the configurational lattice distortion Q exhibit different spatial distributions( [PITH_FULL_IMAGE:fig… view at source ↗
Figure 3
Figure 3. Figure 3: a [PITH_FULL_IMAGE:figures/full_fig_p012_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Fpen at t = 0 and t = ∞ for different cpen (a) t = 0 (b) t = tf (c) t = 0 (d) t = tf [PITH_FULL_IMAGE:figures/full_fig_p014_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Vertical profiles at x = L/2 of the distortion fields for cpen = 2, compared with the uncoupled elastic solution for the same defect configuration: (a) xx component at t = 0; and at t = tf in (b). (c) yy component at t = 0; and at t = tf in (d). 5.2 Evolving dislocation core in a periodic domain without macroscopic loading Before investigating dislocation motion under an applied loading, we first examine t… view at source ↗
Figure 7
Figure 7. Figure 7: d [PITH_FULL_IMAGE:figures/full_fig_p014_7.png] view at source ↗
Figure 6
Figure 6. Figure 6: Horizontal profiles at y = 3H/8 of the distortion fields for cpen = 2, compared with the uncoupled elastic solution for the same defect configuration: (a) xy component at t = 0; and at t = tf in (b). (c) yx component at t = 0; and at t = tf in (d). As we activate the coupling with cpen = 2, the behavior of the system remains qualitatively similar (Figures 7.b and 7.e). The spreading of the core is still no… view at source ↗
Figure 8
Figure 8. Figure 8: a [PITH_FULL_IMAGE:figures/full_fig_p015_8.png] view at source ↗
Figure 7
Figure 7. Figure 7: Evolution of the dislocation density fields [PITH_FULL_IMAGE:figures/full_fig_p016_7.png] view at source ↗
Figure 9
Figure 9. Figure 9: a [PITH_FULL_IMAGE:figures/full_fig_p016_9.png] view at source ↗
Figure 8
Figure 8. Figure 8: Initial configuration of the dislocation annihilation simulation. [PITH_FULL_IMAGE:figures/full_fig_p017_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: In (a) and (b),Time evolution of the dislocation core positions for two values of B.Dashed lines track the PFC density αe and solid lines the FDM density α. (c) shows the snapshot of FDM cores for different values of B and cpen, (d) shows the PFC cores for the same parameters. spatially uniform correction U corr(x, t) = U (x, t) + dt Seff : Σ − U (t) + U p (t) + cpen Seff : Qsym(t) [PITH_FULL_IMAGE:figure… view at source ↗
Figure 10
Figure 10. Figure 10: a [PITH_FULL_IMAGE:figures/full_fig_p017_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: i [PITH_FULL_IMAGE:figures/full_fig_p018_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: a [PITH_FULL_IMAGE:figures/full_fig_p018_12.png] view at source ↗
Figure 11
Figure 11. Figure 11: Evolution of a dislocation dipole during PFC-driven glide and annihilation: [PITH_FULL_IMAGE:figures/full_fig_p019_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: (a): Evolution of the coordinate of the cores for different values of cpen. (b) to (d) same evolution for fixed values of cpen and different applied macroscopic stresses τ . Reference PFC uncoupled motion shown in solid line. 19 [PITH_FULL_IMAGE:figures/full_fig_p019_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Plastic and PFC dissipation in the uncoupled case (a) and coupled case (b) for both PFC and FDM driven cases. Shaded part and dashed black line marks the annihilation start (cores are less than a0 distance). 21 [PITH_FULL_IMAGE:figures/full_fig_p021_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Decay of Fpen in penalty only case for for an immobile dipole in climb. For numerical validation. C.2 Dislocation Transport: Edge Dislocation in 2D In the framework of Field Dislocation Mechanics (FDM), the evolution of the dislocation density tensor α is governed by the transport equation: α˙ + ∇ × α × v d  = 0. (76) We consider a 2D setting in the (x, y) plane, assuming invariance along the z-axis and … view at source ↗
Figure 15
Figure 15. Figure 15: Time evolution of αe and the Peach-Koehler (pk) fit for B in eq. 48 References [1] A. Acharya. A model of crystal plasticity based on the theory of continuously distributed dislo￾cations. Journal of the Mechanics and Physics of Solids, 49(4):761–784, Apr. 2001. [2] A. Acharya. Microcanonical Entropy and Mesoscale Dislocation Mechanics and Plasticity. Journal of Elasticity, 104(1-2):23–44, Aug. 2011. [3] A… view at source ↗

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 · 39 canonical work pages

  1. [1]

    A. Acharya. A model of crystal plasticity based on the theory of continuously distributed dislo- cations. Journal of the Mechanics and Physics of Solids, 49(4):761–784, Apr. 2001

    work page 2001

  2. [2]

    A. Acharya. Microcanonical Entropy and Mesoscale Dislocation Mechanics and Plasticity.Journal of Elasticity, 104(1-2):23–44, Aug. 2011

    work page 2011

  3. [3]

    Acharya, L

    A. Acharya, L. Angheluta, and J. Viñals. Elasticity versus phase field driven motion in the phase field crystal model.Modelling and Simulationin Materials Science and Engineering, 30(6):064005, Sept. 2022

    work page 2022

  4. [4]

    Acharya and A

    A. Acharya and A. Roy. Size effects and idealized dislocation microstructure at small scales: Predictions of a Phenomenological model of Mesoscopic Field Dislocation Mechanics: Part I. Journal of the Mechanics and Physics of Solids, 54(8):1687–1710, Aug. 2006

    work page 2006

  5. [5]

    Acharya and J

    A. Acharya and J. Viñals. Field dislocation mechanics and phase field crystal models.Physical Review B, 102(6):064109, Aug. 2020

    work page 2020

  6. [6]

    P. M. Anderson, J. P. Hirth, and J. Lothe.Theory of Dislocations. Cambridge University Press, Jan. 2017

    work page 2017

  7. [7]

    Arora, R

    A. Arora, R. Arora, and A. Acharya. Mechanics of micropillar confined thin film plasticity.Acta Materialia, 238:118192, Oct. 2022

    work page 2022

  8. [8]

    R. Arora. Computational Approximation of Mesoscale Field Dislocation Mechanics at Finite Deformation. Thesis, Carnegie Mellon University, Mar. 2019

    work page 2019

  9. [9]

    Arora, X

    R. Arora, X. Zhang, and A. Acharya. Finite element approximation of finite deformation dislo- cation mechanics. Computer Methods in Applied Mechanics and Engineering, 367:113076, Aug. 2020

    work page 2020

  10. [10]

    Berry, M

    J. Berry, M. Grant, and K. R. Elder. Diffusive atomistic dynamics of edge dislocations in two dimensions. Physical Review E, 73(3):031609, Mar. 2006. 31

    work page 2006

  11. [11]

    Bradbury, R

    J. Bradbury, R. Frostig, P. Hawkins, M. J. Johnson, C. Leary, D. Maclaurin, G. Necula, A. Paszke, J. VanderPlas, S. Wanderman-Milne, and Q. Zhang. JAX: Composable transformations of Python+NumPy programs, 2018

    work page 2018

  12. [12]

    P. Y. Chan, G. Tsekenis, J. Dantzig, K. A. Dahmen, and N. Goldenfeld. Plasticity and Dislocation Dynamics in a Phase Field Crystal Model.Physical Review Letters, 105(1):015502, June 2010

    work page 2010

  13. [13]

    S. M. Cox and P. C. Matthews. Exponential Time Differencing for Stiff Systems.Journal of Computational Physics, 176(2):430–455, Mar. 2002

    work page 2002

  14. [14]

    M. C. Cross and P. C. Hohenberg. Pattern formation outside of equilibrium.Reviews of modern physics, 65(3):851, 1993

    work page 1993

  15. [15]

    Devincre, V

    B. Devincre, V. Pontikis, Y. Brechet, G. Canova, M. Condat, and L. Kubin. Three-Dimensional Simulations of Plastic Flow in Crystals. In M. Mareschal and B. L. Holian, editors,Microscopic Simulations of Complex Hydrodynamic Phenomena, pages 413–423. Springer US, Boston, MA, 1992

    work page 1992

  16. [16]

    K. R. Elder and M. Grant. Modeling elastic and plastic deformations in nonequilibrium processing using phase field crystals.Physical Review E, 70(5):051605, Nov. 2004

    work page 2004

  17. [17]

    K. R. Elder, M. Katakowski, M. Haataja, and M. Grant. Modeling Elasticity in Crystal Growth. Physical Review Letters, 88(24):245701, June 2002

    work page 2002

  18. [18]

    K. R. Elder, N. Provatas, J. Berry, P. Stefanovic, and M. Grant. Phase-field crystal modeling and classical density functional theory of freezing.Physical Review B, 75(6):064107, Feb. 2007

    work page 2007

  19. [19]

    Eymard, T

    R. Eymard, T. Gallouët, and R. Herbin. Finite volume methods. In Solution of Equation in R (Part 3), Techniques of Scientific Computing (Part 3), volume 7 ofHandbook of Numerical Analysis, pages 713–1018. Elsevier, Jan. 2000

    work page 2000

  20. [20]

    StrongStabilityPreservingRunge-KuttaandMultistep Time Discretizations

    S.Gottlieb, D.Ketcheson, andC.-W.Shu. StrongStabilityPreservingRunge-KuttaandMultistep Time Discretizations. WORLD SCIENTIFIC, 2011

    work page 2011

  21. [21]

    Halperin

    B. Halperin. Physics of defects. In Physics of Defects, number XXXV in Proceedings of the Les Houches Summer School, pages 814–857. North Holland Publishing Company, r. balian, m. kleman, and j.-p. poirier edition, 1981

    work page 1981

  22. [22]

    Heinonen, C

    V. Heinonen, C. V. Achim, J. M. Kosterlitz, S.-C. Ying, J. Lowengrub, and T. Ala-Nissila. Con- sistent Hydrodynamics for Phase Field Crystals. Physical Review Letters, 116(2):024303, Jan. 2016

    work page 2016

  23. [23]

    DYNAMICALTHEORYOFDISLOCATIONS

    A.M.Kosevich. DYNAMICALTHEORYOFDISLOCATIONS. SovietPhysicsUspekhi, 7(6):837, June 1965

    work page 1965

  24. [24]

    G. F. Mazenko. Vortex Velocities in the $\mathit{O}(\mathit{n})$ Symmetric Time-Dependent Ginzburg-Landau Model. Physical Review Letters, 78(3):401–404, Jan. 1997

    work page 1997

  25. [25]

    T. Mura. Continuous distribution of moving dislocations.Philosophical Magazine, 8(89):843–857, May 1963

    work page 1963

  26. [26]

    A. Roy. Finite element approximation of field dislocation mechanics.Journal of the Mechanics and Physics of Solids, 53(1):143–170, Jan. 2005

    work page 2005

  27. [27]

    Skaugen, L

    A. Skaugen, L. Angheluta, and J. Viñals. Dislocation dynamics and crystal plasticity in the phase-field crystal model.Physical Review B, 97(5):054113, Feb. 2018. 32

    work page 2018

  28. [28]

    Skaugen, L

    A. Skaugen, L. Angheluta, and J. Viñals. Separation of Elastic and Plastic Timescales in a Phase Field Crystal Model.Physical Review Letters, 121(25):255501, Dec. 2018

    work page 2018

  29. [29]

    Skogvoll, L

    V. Skogvoll, L. Angheluta, A. Skaugen, M. Salvalaglio, and J. Viñals. A phase field crystal theory of the kinematics of dislocation lines.Journal of the Mechanics and Physics of Solids, 166:104932, Sept. 2022

    work page 2022

  30. [30]

    Skogvoll, J

    V. Skogvoll, J. Rønning, M. Salvalaglio, and L. Angheluta. A unified field theory of topological defects and non-linear local excitations.npj Computational Materials, 9(1):1–13, July 2023

    work page 2023

  31. [31]

    Skogvoll, M

    V. Skogvoll, M. Salvalaglio, and L. Angheluta. Hydrodynamic phase field crystal approach to interfaces, dislocations, and multi-grain networks, Sept. 2022

    work page 2022

  32. [32]

    Stefanovic, M

    P. Stefanovic, M. Haataja, and N. Provatas. Phase-Field Crystals with Elastic Interactions. Physical Review Letters, 96(22):225504, June 2006

    work page 2006

  33. [33]

    M. V. Upadhyay and J. Viñals. Coupling Phase Field Crystal and Field Dislocation Mechanics for a consistent description of dislocation structure and elasticity.European Journal of Mechanics - A/Solids, 108:105419, Nov. 2024

    work page 2024

  34. [34]

    van Leer

    B. van Leer. Towards the ultimate conservative difference scheme. V. A second-order sequel to Godunov’s method.Journal of Computational Physics, 32(1):101–136, July 1979

    work page 1979

  35. [35]

    S. N. Varadhan, A. J. Beaudoin, A. Acharya, and C. Fressengeas. Dislocation transport using an explicit Galerkin/least-squares formulation. Modelling and Simulation in Materials Science and Engineering, 14(7):1245–1270, Oct. 2006

    work page 2006

  36. [36]

    Sur l’equilibre des cirps elastiques multiplement connexes.Ann

    Voltera. Sur l’equilibre des cirps elastiques multiplement connexes.Ann. Ecole Norm. Super, 24:401–517, 1907

    work page 1907

  37. [37]

    Wijnen, R

    J. Wijnen, R. H. J. Peerlings, J. P. M. Hoefnagels, and M. G. D. Geers. A discrete slip plane model for simulating heterogeneous plastic deformation in single crystals.International Journal of Solids and Structures, 228:111094, Oct. 2021

    work page 2021

  38. [38]

    Zhang, A

    X. Zhang, A. Acharya, N. J. Walkington, and J. Bielak. A single theory for some quasi-static, supersonic, atomic, and tectonic scale applications of dislocations.Journal of the Mechanics and Physics of Solids, 84:145–195, Nov. 2015

    work page 2015

  39. [39]

    S. J. Zhou, D. L. Preston, P. S. Lomdahl, and D. M. Beazley. Large-scale molecular dynamics simulations of dislocation intersection in copper.Science (New York,N.Y.), 279(5356):1525–1527, Mar. 1998. 33

    work page 1998