pith. sign in

arxiv: 2604.21105 · v1 · submitted 2026-04-22 · ❄️ cond-mat.supr-con · cond-mat.dis-nn· cond-mat.mtrl-sci

Healing of topological defects while crystallizing nanocrystals

Pith reviewed 2026-05-09 22:25 UTC · model grok-4.3

classification ❄️ cond-mat.supr-con cond-mat.dis-nncond-mat.mtrl-sci
keywords vortex nanocrystalstopological defectscrystallization dynamicsconfinement effectsdefect healingLangevin simulationssuperconductors
0
0 comments X

The pith

The radial distribution of topological defects in vortex nanocrystals freezes into a stationary profile below the melting line, tuned by vortex properties and confinement.

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

Simulations of vortex nanocrystals with a few hundred vortices reveal how confinement affects the crystallization process and defect healing. The low-temperature structures show defects annealing at the edges, in agreement with experiments on Bi2Sr2CaCu2O8+δ samples. The central finding is that the radial profile of defects becomes stationary once the temperature falls below the melting line, with its shape set by the elasticity of the vortices and the sample boundaries. This stationary state determines the final defect distribution and can be generalized to other confined nanocrystal systems.

Core claim

Vortex nanocrystals formed in field-cooling conditions exhibit healing of topological defects at the edges. The low-temperature radial distribution of topological defects is a stationary profile that freezes at a temperature below the melting line. This profile is tuned by intrinsic properties of the vortex structure and by the confinement effect.

What carries the argument

The stationary radial distribution of topological defects, which freezes below the melting temperature and depends on vortex elasticity and confinement.

If this is right

  • The healing effect at the edges matches quantitative experimental data for different vortex densities and elasticities.
  • Low-temperature structural properties vary with the physical size of the samples.
  • The findings apply to describing physical properties of confined soft condensed matter nanocrystals in general.
  • Changing vortex density or elasticity alters the final radial defect profile.

Where Pith is reading between the lines

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

  • Confinement could be used to engineer specific defect distributions in small crystalline systems.
  • Similar freezing of defect profiles may occur in colloidal or other soft matter nanocrystals under boundary constraints.
  • The model suggests the profile is robust against further temperature decrease once frozen.

Load-bearing premise

The Langevin dynamics model with a few hundred vortices quantitatively reproduces the defect healing seen in real Bi2Sr2CaCu2O8+δ micron-sized samples.

What would settle it

Direct measurement of whether the radial distribution of topological defects in vortex nanocrystals remains constant when the sample is cooled further below the temperature where it first freezes.

Figures

Figures reproduced from arXiv: 2604.21105 by A. B. Kolton, M. I. Dolz, Y. Fasano.

Figure 1
Figure 1. Figure 1: FIG. 1. Typical temperature-history protocol used in our sim [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Total density of defects in the vortex nanocrystal at [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. (a) Example of a snapshot of a nanocrystal with [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Examples of Delaunay triangulations of the snapshots of vortex nanocrystals obtained at the lowest temperature, [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Radial density of defects in the vortex nanocrystal [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Dependence of the phenomenological healing length [PITH_FULL_IMAGE:figures/full_fig_p006_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Central density of topological defects as a function [PITH_FULL_IMAGE:figures/full_fig_p007_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Radial density of topological defects as a function of [PITH_FULL_IMAGE:figures/full_fig_p008_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Variation with temperature of the total density of [PITH_FULL_IMAGE:figures/full_fig_p008_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Sequence of snapshots while field-cooling a vortex nanocrystal nucleated at 16 G in a sample with 50 [PITH_FULL_IMAGE:figures/full_fig_p009_10.png] view at source ↗
read the original abstract

Understanding the role of confinement while crystallizing nanocrystals is very relevant for predicting their structure and physical properties. With this aim we perform Langevin dynamics simulations of nanocrystals of the model system of few hundred vortices nucleated in micron-sized superconductors. We study the crystallization dynamics and the low-temperature structural properties of vortex nanocrystals nucleated in field-cooling conditions when changing vortex density or elasticity of the system and physical size of the samples. The low-temperature snapshots obtained in simulations present a healing effect at the edges that is in quantitative agreement with experimental data in Bi2Sr2CaCu2O8+{\delta} micron-sized samples. We show that the low-temperature radial distribution of topological defects is a stationary profile frozen at a temperature below the melting line tuned by intrinsic properties of the vortex structure and on the confinement effect. These findings on the dynamics and spatial profile of topological defects can be applied to describe the physical properties of confined soft condensed matter nanocrystals in general.

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 reports Langevin dynamics simulations of a few hundred point vortices in confined 2D geometries modeling vortex nanocrystals in superconductors under field-cooling. By varying vortex density, elasticity, and sample size, the authors study crystallization dynamics and low-temperature structural properties, focusing on edge healing of topological defects. They claim quantitative agreement between simulated low-T snapshots and experimental observations in Bi2Sr2CaCu2O8+δ micron-sized samples, concluding that the radial defect distribution forms a stationary profile frozen below the melting line and tuned by intrinsic vortex properties plus confinement. The results are proposed to generalize to confined soft condensed matter nanocrystals.

Significance. If the reported quantitative agreement holds under detailed scrutiny and the idealized model captures the dominant physics, the work would usefully demonstrate how confinement induces stationary defect profiles during crystallization in vortex systems. The controlled variation of density, elasticity, and size isolates their contributions, which is a methodological strength. The link to Bi2212 data and the suggestion of broader applicability to soft-matter nanocrystals could inform structure prediction in confined systems, provided the simulation-to-experiment mapping is robustly documented.

major comments (2)
  1. [Abstract] Abstract: The central claim of 'quantitative agreement' between simulated low-temperature edge healing and Bi2Sr2CaCu2O8+δ data is load-bearing for interpreting the defect profile as an intrinsic stationary state; however, no error bars, overlap metrics, radial-distribution-function comparisons, or fitting details are supplied, preventing assessment of whether the agreement is statistically meaningful or merely qualitative.
  2. [Model and Results] Model and Results sections: The overdamped 2D Langevin dynamics for point vortices omits pinning centers and 3D layered effects known to exist in real BSCCO nanocrystals, yet the manuscript asserts that the low-T profile is tuned solely by intrinsic vortex properties and confinement; without a sensitivity test or justification for these omissions, the claim that the profile is model-independent does not follow.
minor comments (2)
  1. [Abstract] Abstract: The phrase 'few hundred vortices' should specify the exact range (e.g., 100–500) and how results scale with N to allow reproducibility.
  2. Notation: The symbol for elasticity (likely a spring constant or shear modulus) is introduced without an explicit equation reference, making it unclear how it enters the Langevin equation.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive assessment of our work and the constructive comments on the abstract and model assumptions. We address the major comments point by point below.

read point-by-point responses
  1. Referee: [Abstract] Abstract: The central claim of 'quantitative agreement' between simulated low-temperature edge healing and Bi2Sr2CaCu2O8+δ data is load-bearing for interpreting the defect profile as an intrinsic stationary state; however, no error bars, overlap metrics, radial-distribution-function comparisons, or fitting details are supplied, preventing assessment of whether the agreement is statistically meaningful or merely qualitative.

    Authors: We agree that the abstract would benefit from more explicit reference to the quantitative aspects. In the manuscript, the agreement is shown via direct overlay of the radial defect density profiles from simulations and experiments, with the simulated profiles falling within the experimental scatter. We will revise the abstract to state that the agreement is within the variability observed in both datasets and add a note on the comparison method. Additionally, we will include error bars in the relevant figures and a brief description of the overlap in the revised manuscript. revision: yes

  2. Referee: [Model and Results] Model and Results sections: The overdamped 2D Langevin dynamics for point vortices omits pinning centers and 3D layered effects known to exist in real BSCCO nanocrystals, yet the manuscript asserts that the low-T profile is tuned solely by intrinsic vortex properties and confinement; without a sensitivity test or justification for these omissions, the claim that the profile is model-independent does not follow.

    Authors: The 2D point-vortex model is standard for studying vortex crystallization in thin films and captures the essential physics of defect healing under confinement for the parameters relevant to BSCCO nanocrystals. Pinning is neglected because the samples are high-quality with minimal pinning effects at the fields used, and 3D effects are secondary in the thin limit. We do not assert model-independence but rather that the stationary profile emerges from the competition between vortex-vortex interactions, elasticity, and boundary conditions, as demonstrated by varying these parameters in our simulations. We will add a dedicated subsection discussing the model approximations and supporting references to justify the omissions, showing that they do not alter the qualitative and quantitative features reported. revision: partial

Circularity Check

0 steps flagged

No circularity: simulation-derived stationary profile is independent of inputs

full rationale

The paper runs overdamped Langevin dynamics on a few hundred point vortices in confined 2D geometry, varying only density, elasticity, and sample size. It reports that low-T edge healing and the radial defect distribution emerge as stationary profiles below the melting line. This is presented as a direct numerical observation, not a fit to data or a self-citation chain. The quantitative agreement with Bi2212 experiments is stated as an outcome, not an input that forces the profile by construction. No equation or step equates a claimed prediction back to a fitted parameter or prior self-result; the derivation remains self-contained within the simulation protocol.

Axiom & Free-Parameter Ledger

3 free parameters · 1 axioms · 0 invented entities

The central claim rests on the assumption that standard Langevin dynamics with tunable vortex density, elasticity, and sample size captures the essential physics; no new entities are introduced.

free parameters (3)
  • vortex density
    Varied to study crystallization dynamics
  • elasticity of the system
    Varied to study crystallization dynamics
  • physical size of samples
    Varied to study confinement effects
axioms (1)
  • domain assumption Langevin dynamics accurately models vortex motion and interactions in the superconducting system
    Invoked for all reported simulation results

pith-pipeline@v0.9.0 · 5475 in / 1262 out tokens · 23969 ms · 2026-05-09T22:25:04.722286+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

49 extracted references

  1. [1]

    J. P. Sethna, M. K. Bierbaum, K. A. Dahmen, C. P. Goodrich, J. R. Greer, L. X. Hayden, J. P. Kent-Dobias, E. D. Lee, D. B. Liarte, X. Ni, K. N. Quinn, A. Raju, D. Z. Rocklin, A. Shekhawat, and S. Zapperi, Annual Review Materials Research47, 217 (2017)

  2. [2]

    V. J. Anderson and H. N. W. Lekkerkerker, Nature416, 811 (2002)

  3. [3]

    Taloni, M

    A. Taloni, M. Vodret, G. Costantini, and S. Zapperi, Na- ture Reviews Materials3, 211–224 (2018)

  4. [4]

    D. R. Nelson,Defects and geometry in condensed mat- ter Physics, 1st ed. (Cambridge University Press, Cam- bridge, United Kingdom, 2002)

  5. [5]

    Hillert, Acta Metallographica13, 227–238 (1965)

    M. Hillert, Acta Metallographica13, 227–238 (1965)

  6. [6]

    Moretti, M.-C

    P. Moretti, M.-C. Miguel, and S. Zapperi, Phys. Rev. B 72, 014505 (2005)

  7. [7]

    Salvalaglio, R

    M. Salvalaglio, R. Backofen, K. R. Elder, and A. Voigt, Phys. Rev. Mater.2, 053804 (2018)

  8. [8]

    Salvalaglio, A

    M. Salvalaglio, A. Voigt, and K. R. Elder, npj Computa- tional Materials5, 48 (2019)

  9. [9]

    C. Qiu, M. Punke, Y. Tian, Y. Han, S. Wang, Y. Su, M. Salvalaglio, X. Pan, D. J. Srolovitz, and J. Han, Sci- ence385, 980 (2024)

  10. [10]

    N. J. Petch, J. Iron Steel Res. Int.19, 25 (1953)

  11. [11]

    E. O. Hall, Proc. Phys. Soc., Sec. B64, 747 (1954)

  12. [12]

    Benoit–Mar´ echal, I

    L. Benoit–Mar´ echal, I. Nitschke, A. Voigt, and M. Sal- valaglio, Mechanics of Materials198, 105114 (2024)

  13. [13]

    Rudolph, AIP Conference Proceedings916(2007)

    P. Rudolph, AIP Conference Proceedings916(2007)

  14. [14]

    Reichhardt and C

    C. Reichhardt and C. J. Olson Reichhardt, Rep. Prog. Phys.80, 026501 (2017)

  15. [15]

    P. M. Chaikin and T. C. Lubensky,Principles of Con- densed Matter Physics(Cambridge University Press, 1995)

  16. [16]

    Menghini, Y

    M. Menghini, Y. Fasano, and F. de la Cruz, Phys. Rev. B65, 064510 (2002)

  17. [17]

    Peeters, B

    F. Peeters, B. Baelus, and M. Miloˇ sevi´ c, Physica E: Low-dimensional Systems and Nanostructures18, 312 (2003), 23rd International Conference on Low Temper- ature Physics (LT23)

  18. [18]

    Baus and C

    M. Baus and C. F. Tejero, eds., Topological defects, inEquilibrium Statistical Physics:Phases of Matter and Phase Transitions(Springer Berlin Heidelberg, Berlin, Heidelberg, 2008) pp. 307–313

  19. [19]

    van der Beek, S

    C. van der Beek, S. Demirdis, M. Konczykowski, Y. Fasano, N. Cejas Bolecek, H. Pastoriza, D. Colson, and F. Rullier-Albenque, Physica B: Condensed Matter 407, 1746 (2012), proceedings of the International Work- shop on Electronic Crystals (ECRYS-2011)

  20. [20]

    J. B. Llorens, L. Embon, A. Correa, J. D. Gonz´ alez, E. Herrera, I. Guillam´ on, R. F. Luccas, J. Azpeitia, F. J. Mompe´ an, M. Garc´ ıa-Hern´ andez, C. Munuera, J. A. S´ anchez, Y. Fasano, M. V. Miloˇ sevi´ c, H. Suderow, and Y. Anahory, Phys. Rev. Res.2, 013329 (2020)

  21. [21]

    J. A. S´ anchez, R. C. Maldonado, M. L. Amig´ o, G. Nieva, A. Kolton, and Y. Fasano, Phys. Rev. B107, 094508 (2023)

  22. [22]

    J. Puig, J. A. S´ anchez, E. Herrera, I. Guillam´ on, Z. Pribulov´ a, J. Kacmarcik, H. Suderow, A. B. Kolton, and Y. Fasano, Phys. Rev. B110, 024108 (2024). 11

  23. [23]

    G. R. Berdiyorov, B. J. Baelus, M. V. Miloˇ sevi´ c, and F. M. Peeters, Phys. Rev. B68, 174521 (2003)

  24. [24]

    I. V. Grigorieva, W. Escoffier, J. Richardson, L. Y. Vin- nikov, S. Dubonos, and V. Oboznov, Phys. Rev. Lett.96, 077005 (2006)

  25. [25]

    H. J. Zhao, V. R. Misko, F. M. Peeters, V. Oboznov, S. V. Dubonos, and I. V. Grigorieva, Phys. Rev. B78, 104517 (2008)

  26. [26]

    Turlo and T

    V. Turlo and T. Rupert, Acta Materialia151(2017)

  27. [27]

    J. D. Gonz´ alez, O. J. Gaona, and J. Barba-Ortega, Jour- nal of Low Temperature Physics198, 123 (2020)

  28. [28]

    A. O. Haridy, A. A. Sheekhoon, and V. V. Kresin, Phys. Rev. B113, L041404 (2026)

  29. [29]

    Bubeck, C

    R. Bubeck, C. Bechinger, S. Neser, and P. Leiderer, Phys. Rev. Lett.82, 3364 (1999)

  30. [30]

    S. M. Reimann and M. Manninen, Rev. Mod. Phys.74, 1283 (2002)

  31. [31]

    Moshchalkov, L

    V. Moshchalkov, L. Gielen, C. Strunk, R. Jonckheere, X. Qiu, C. Haesendonck, and Y. Bruynseraede, Nature 373, 319 (1995)

  32. [32]

    A. Geim, I. Grigorieva, S. Dubonos, J. Lok, J. Maan, A. Filippov, and F. M. Peeters, Nature390, 259 (1997)

  33. [33]

    J. J. Palacios, Phys. Rev. B58, R5948 (1998)

  34. [34]

    V. A. Schweigert, F. M. Peeters, and P. S. Deo, Phys. Rev. Lett.81, 2783 (1998)

  35. [35]

    Y. M. Wang, M. S. Fuhrer, A. Zettl, S. Ooi, and T. Tamegai, Phys. Rev. Lett.86, 3626 (2001)

  36. [36]

    M. I. Dolz, Y. Fasano, N. R. Cejas Bolecek, H. Pastoriza, V. Mosser, M. Li, and M. Konczykowski, Phys. Rev. Lett. 115, 137003 (2015)

  37. [37]

    Blatter, M

    G. Blatter, M. V. Feigel’man, V. B. Geshkenbein, A. I. Larkin, and V. M. Vinokur, Rev. Mod. Phys.66, 1125 (1994)

  38. [38]

    J. Puig, F. El´ ıas, J. A. S´ anchez, R. C. Maldonado, G. Rumi, G. Nieva, P. Pedrazzini, A. B. Kolton, and Y. Fasano, Communications Materials3, 32 (2022)

  39. [39]

    Y. M. Wang, A. Zettl, S. Ooi, and T. Tamegai, Phys. Rev. B65, 184506 (2002)

  40. [40]

    N. R. Cejas Bolecek, M. I. Dolz, A. Kolton, H. Pastor- iza, C. J. v. d. Beek, M. Konczykowski, M. Menghini, G. Nieva, and Y. Fasano, Journal of Low Temperature Physics179, 35 (2015)

  41. [41]

    N. R. C. Bolecek, M. I. Dolz, H. Pastoriza, M. Kon- czykowski, C. J. van der Beek, A. B. Kolton, and Y. Fasano, Phys. Rev. B96, 024507 (2017)

  42. [42]

    M. Dolz, Y. Fasano, H. Pastoriza, V. Mosser, M. Li, and M. Konczykowski, Physical Review B90, 144507 (2014)

  43. [43]

    Reichhardt, J

    C. Reichhardt, J. Groth, C. J. Olson, S. B. Field, and F. Nori, Phys. Rev. B54, 16108 (1996)

  44. [44]

    M. P. Allen and D. J. Tildesley,Computer Simulation of Liquids, 2nd ed. (Oxford University Press, Oxford, 2017)

  45. [45]

    Fasano, M

    Y. Fasano, M. D. Seta, M. Menghini, H. Pastoriza, and F. de la Cruz, Proceedings of the National Academy of Sciences of the United States of America102, 3898 (2005)

  46. [46]

    Fasano, M

    Y. Fasano, M. De Seta, M. Menghini, H. Pastoriza, and F. de la Cruz, Solid State Communications128, 51 (2003)

  47. [47]

    Hull and D

    D. Hull and D. J. Bacon,Introduction to dislocations, 5th ed. (Butterworth-Heinemann, Elsevier, Oxford, 2011)

  48. [48]

    Cruz-Garc´ ıa, J

    A. Cruz-Garc´ ıa, J. Puig, R. M. Besana, A. B. Kolton, and Y. Fasano, Phys. Rev. B113, 024112 (2026)

  49. [49]

    Fasano, J

    Y. Fasano, J. Herbsommer, and F. de la Cruz, Physica Status Solidi (b)215, 563 (1999)