pith. sign in

arxiv: 2512.03833 · v3 · pith:3TLE23LOnew · submitted 2025-12-03 · ⚛️ physics.atom-ph

Shell formation and two-dimensional nanofriction in three-dimensional ion Coulomb crystals

Pith reviewed 2026-05-21 18:13 UTC · model grok-4.3

classification ⚛️ physics.atom-ph
keywords Coulomb crystalsnanofrictionshell formationPeierls-Nabarro potentialmolecular dynamicsPaul trapsrotational barrier2D friction
0
0 comments X

The pith

Small changes in ion number can alter the outer-shell rotational barrier in three-dimensional Coulomb crystals by up to a factor of 60.

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

This paper examines how three-dimensional ion Coulomb crystals in linear Paul traps form concentric shells that act as curved interfaces for two-dimensional nanofriction studies. Molecular-dynamics simulations map shell formation versus ion number N and trap aspect ratio, then define a Peierls-Nabarro-type potential using the outer-shell rotation angle as a collective coordinate against a static inner core. The central result is that the height of the rotational energy barrier depends strongly on N, with a one-ion change shifting the barrier by up to a factor of seven and changes of only a few ions producing shifts up to a factor of sixty. These large variations arise from the competition among inter-shell interactions, the outer shell's mechanical response, and trap confinement, which in turn produce pinned, stick-slip, and smooth-sliding regimes plus occasional hysteresis. A reader would care because the setup supplies an atomically precise, experimentally tunable model for nanoscale friction that could guide the design of nanorotors or low-friction mechanical elements.

Core claim

In self-organized three-dimensional ion Coulomb crystals the outer shell rotates against a static inner core under a Peierls-Nabarro-type potential whose barrier is governed by the system-dependent interplay between inter-shell interaction, outer-shell response and trap confinement. Changing N by one can alter the effective rotational barrier by up to a factor of approximately seven, while changes by only a few ions produce variations up to a factor of approximately sixty. Dynamical simulations with applied torques show pinned, stick-slip and smooth-sliding regimes whose depinning thresholds depend on ion number, inner-shell geometry and trap aspect ratio; some configurations exhibit torque-

What carries the argument

The rotation angle of the outer shell relative to the inner core, used as the collective coordinate to compute the Peierls-Nabarro-type potential.

Load-bearing premise

The inner core remains perfectly static while the outer shell rotates as a rigid body, allowing the rotation angle to serve as a clean collective coordinate for the Peierls-Nabarro-type potential.

What would settle it

Directly measuring the torque needed to rotate the outer shell for two crystals whose ion numbers differ by one and checking whether the observed barrier heights differ by the predicted factor of up to seven.

Figures

Figures reproduced from arXiv: 2512.03833 by L.-A. R\"uffert, T. E. Mehlst\"aubler.

Figure 1
Figure 1. Figure 1: FIG. 1: Experimental image of a three-dimensional (3D), Doppler-cooled Coulomb crystal of approximately 200 [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Model system of two-dimensional (2D) nanofric [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: (a) Number of shells as different colored regions with respect to the number of ions [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Rotation of the outer shell (blue wireframe) around the [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: Angle dependent modulations of the Peierls-Nabarro-type potential of the outer shell for selected crystal [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7: Normalized effective energy barrier [PITH_FULL_IMAGE:figures/full_fig_p009_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: shows Ω(t)/Ωmax over time, of the outer shell of the α = 1.52, N = 44 configuration. Different driving torques are applied and the system is propagated un￾til a steady state has been reached, before the data is recorded. The oscillations in Ω(t)/Ωmax are indicative of stick-slip motion, while the periodicities in the graphs reflect full revolutions of the outer shell. For increasing driving torques, the im… view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9: Sliding efficiency [PITH_FULL_IMAGE:figures/full_fig_p011_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10: Normalized average angular velocity of each [PITH_FULL_IMAGE:figures/full_fig_p012_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11: Schematic illustration of multidimensional fric [PITH_FULL_IMAGE:figures/full_fig_p012_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12: Deltahedron shapes based on the ground states of ion crystals with trapping parameter [PITH_FULL_IMAGE:figures/full_fig_p017_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: FIG. 13: (a) Histogram of a 100-ion crystal with the iteration variable [PITH_FULL_IMAGE:figures/full_fig_p017_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: FIG. 14: Configurations of the inner shell ions for selected crystals. Orange lines mark the projection of the ion [PITH_FULL_IMAGE:figures/full_fig_p019_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: FIG. 15: Radial distance [PITH_FULL_IMAGE:figures/full_fig_p020_15.png] view at source ↗
read the original abstract

Self-organized three-dimensional (3D) ion Coulomb crystals in linear Paul traps naturally form concentric shells that provide a curved, atomically resolved interface for studying two-dimensional (2D) nanofriction. Building on earlier studies of one-dimensional nanofriction and orientational melting in 2D ion crystals, we extend friction studies from linear chains and planar rings to 3D shell structures. Using molecular-dynamics simulations, we map shell formation as a function of ion number N and trap aspect ratio and obtain a simple scaling relation that can aid ion-number estimation in experiments. We compute a Peierls-Nabarro-type potential for rotating the outer shell against a static inner core, using the rotation angle as a collective coordinate. Changing N by one can alter the effective rotational barrier by up to a factor of ~7, while changes by only a few ions can lead to variations up to a factor of ~60. Combining geometric commensurability analysis with energy decomposition, we show that the barrier is governed by a system-dependent interplay between inter-shell interaction, outer-shell response and trap confinement. Dynamical simulations with applied torques reveal pinned, stick-slip and smooth-sliding regimes with depinning thresholds that depend on ion number, inner-shell geometry and trap aspect ratio. Some configurations show hysteresis due to torque-induced metastable states. We further find that spatially varying coupling to the inner-core corrugation can create coexisting fast and slow domains within the rotating outer shell, realizing multidimensional friction in which intra-shell shear and inter-shell nanofriction act simultaneously. Our results establish self-organized ion Coulomb crystals as model systems for 2D nanofriction and suggest routes toward ion-based nanorotors, torque sensors and ultra-low-friction nanomechanical systems.

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 / 3 minor

Summary. The manuscript uses molecular-dynamics simulations to map shell formation in three-dimensional ion Coulomb crystals confined in linear Paul traps as a function of ion number N and trap aspect ratio, deriving a scaling relation for experimental ion-number estimation. It computes a Peierls-Nabarro-type potential for rotating the outer shell against a static inner core using the rotation angle as collective coordinate, reporting that single-ion changes in N can alter the effective rotational barrier by factors up to ~7 and changes by a few ions by up to ~60. Dynamical torque simulations reveal pinned, stick-slip, and smooth-sliding regimes whose depinning thresholds depend on N, inner-shell geometry, and aspect ratio, together with hysteresis from metastable states and coexisting fast/slow domains arising from spatially varying inter-shell coupling.

Significance. If the numerical results are robust, the work provides a concrete platform for studying two-dimensional nanofriction on curved, atomically resolved interfaces and extends prior studies of one-dimensional chains and planar rings. The MD-derived scaling relation and regime diagrams constitute a useful contribution, and the reported extreme sensitivity of barriers to small changes in N offers a route to tunable nanofriction with possible applications to nanorotors and torque sensors. The authors receive credit for obtaining these relations directly from trajectories and energy decompositions without re-use of prior fitted parameters.

major comments (2)
  1. [Abstract and Methods] Abstract and Methods section: no information is supplied on the precise form of the inter-ion potential (Coulomb plus trap), integration timestep, equilibration protocol, or statistical uncertainties on the reported barrier variations (factors of ~7 and ~60). These omissions are load-bearing because the central claims of N-sensitivity and the existence of distinct dynamical regimes rest entirely on the simulation outputs.
  2. [Energy decomposition and dynamical torque simulations] Energy-decomposition and torque-simulation paragraphs: the Peierls-Nabarro potential is constructed by treating the rotation angle as the sole collective coordinate while holding the inner core perfectly static and the outer shell rigid. If torque-driven or thermal displacements occur in the inner shell, or if the outer shell experiences intra-shell shear or radial breathing, the extracted barrier heights and the claimed interplay between inter-shell interaction, outer-shell response, and trap confinement would incorporate additional degrees of freedom not accounted for in the reported values.
minor comments (3)
  1. [Abstract] The explicit functional form of the 'simple scaling relation' for shell formation is not stated in the abstract or early text; providing the relation (or at least its leading dependence on N and aspect ratio) would improve immediate accessibility.
  2. [Figures] Figure captions for the regime diagrams should list the specific values of N, aspect ratio, and temperature used, together with any averaging procedure over trajectories.
  3. [Dynamical simulations] A few sentences describing the hysteresis mechanism could be clarified by indicating whether the metastable states are identified from energy minima or from long-time trajectory statistics.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and for the constructive comments, which help improve the clarity and robustness of our presentation. We address each major comment below and will revise the manuscript accordingly.

read point-by-point responses
  1. Referee: [Abstract and Methods] Abstract and Methods section: no information is supplied on the precise form of the inter-ion potential (Coulomb plus trap), integration timestep, equilibration protocol, or statistical uncertainties on the reported barrier variations (factors of ~7 and ~60). These omissions are load-bearing because the central claims of N-sensitivity and the existence of distinct dynamical regimes rest entirely on the simulation outputs.

    Authors: We agree that these details are necessary for full reproducibility and to underpin the central claims. In the revised manuscript we will expand the Methods section to specify the inter-ion potential (Coulomb repulsion plus the linear Paul trap harmonic confinement), the integration timestep, the equilibration protocol (including annealing schedule and relaxation criteria), and quantitative estimates of statistical uncertainties on the barrier heights obtained from multiple independent trajectories and energy decompositions. revision: yes

  2. Referee: [Energy decomposition and dynamical torque simulations] Energy-decomposition and torque-simulation paragraphs: the Peierls-Nabarro potential is constructed by treating the rotation angle as the sole collective coordinate while holding the inner core perfectly static and the outer shell rigid. If torque-driven or thermal displacements occur in the inner shell, or if the outer shell experiences intra-shell shear or radial breathing, the extracted barrier heights and the claimed interplay between inter-shell interaction, outer-shell response, and trap confinement would incorporate additional degrees of freedom not accounted for in the reported values.

    Authors: The rigid-inner-core and rigid-outer-shell approximation is adopted deliberately to isolate the rotational Peierls-Nabarro barrier as a function of the collective rotation angle. The manuscript already reports an energy decomposition that explicitly includes outer-shell response and trap-confinement contributions, and the subsequent dynamical torque simulations allow all degrees of freedom to evolve freely, reproducing the pinned, stick-slip and smooth-sliding regimes. We will add a dedicated paragraph in the revised text that quantifies the energetic cost of inner-shell displacements and intra-shell shear, thereby clarifying the regime of validity of the reported barriers while preserving the central physical conclusions. revision: partial

Circularity Check

0 steps flagged

No significant circularity; results from direct MD and energy decomposition

full rationale

The paper derives shell formation, scaling relations, and rotational barriers directly from molecular-dynamics trajectories and explicit energy decompositions that vary the rotation angle as a collective coordinate under the stated static-inner-core assumption. No load-bearing step reduces a claimed prediction or uniqueness result to a parameter fitted from the same data and re-used, nor to a self-citation chain whose content is itself unverified. The reported factor-of-7 and factor-of-60 barrier variations are simulation outputs, not quantities defined by construction from prior fits or ansatzes imported via the authors' own earlier work. The derivation chain therefore remains self-contained against external simulation benchmarks.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The central claims rest on classical molecular-dynamics assumptions for Coulomb interactions in a harmonic trap plus the modeling choice that shells can be treated as rotatable rigid bodies for the purpose of computing a collective-coordinate potential. No new particles or forces are postulated.

free parameters (1)
  • trap aspect ratio
    Used to map shell formation; its specific values are chosen to produce the reported shell structures and are not derived from first principles.
axioms (2)
  • domain assumption Ions interact via pure Coulomb repulsion in a static harmonic trap potential
    Standard for Paul-trap ion-crystal simulations; invoked implicitly throughout the shell-formation and rotation analysis.
  • domain assumption Classical molecular dynamics sufficiently captures the relevant dynamics at the simulated temperatures
    Required for all reported energy barriers and sliding regimes.

pith-pipeline@v0.9.0 · 5858 in / 1616 out tokens · 79765 ms · 2026-05-21T18:13:27.687681+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

70 extracted references · 70 canonical work pages

  1. [1]

    The complete algorithm is explained in more detail in Appendix A

    Multiple SA runs are executed until the same con- figuration has been found repeatedly, indicating that the structural ground state has been found. The complete algorithm is explained in more detail in Appendix A. While SA is a powerful algorithm to explore the energy landscape and reliably converges on low-energy config- urations, it cannot guarantee tha...

  2. [2]

    Holmberg and A

    K. Holmberg and A. Erdemir, Influence of tribology on global energy consumption, costs and emissions, Friction 5, 263 (2017)

  3. [3]

    Binnig, C

    G. Binnig, C. F. Quate, and Ch. Gerber, Atomic Force Microscope, Phys. Rev. Lett.56, 930 (1986)

  4. [4]

    Martin, C

    Y. Martin, C. C. Williams, and H. K. Wickramasinghe, Atomic force microscope–force mapping and profiling on a sub 100-˚A scale, J. Appl. Phys.61, 4723 (1987)

  5. [5]

    C. M. Mate, G. M. McClelland, R. Erlandsson, and S. Chiang, Atomic-scale friction of a tungsten tip on a graphite surface, Phys. Rev. Lett.59, 1942 (1987)

  6. [6]

    Dienwiebel, G

    M. Dienwiebel, G. S. Verhoeven, N. Pradeep, J. W. M. Frenken, J. A. Heimberg, and H. W. Zandbergen, Super- lubricity of Graphite, Phys. Rev. Lett.92, 126101 (2004)

  7. [7]

    Xu, T.-B

    L. Xu, T.-B. Ma, Y.-Z. Hu, and H. Wang, Vanishing stick–slip friction in few-layer graphenes: The thickness effect, Nanotechnology22, 285708 (2011)

  8. [8]

    Andersson and A

    D. Andersson and A. S. De Wijn, Understanding the fric- tion of atomically thin layered materials, Nat. Commun. 11, 420 (2020)

  9. [9]

    Vanossi, N

    A. Vanossi, N. Manini, and E. Tosatti, Static and dy- namic friction in sliding colloidal monolayers, Proc. Natl. 15 Acad. Sci. U.S.A.109, 16429 (2012)

  10. [10]

    Mandelli, A

    D. Mandelli, A. Vanossi, and E. Tosatti, Stick-slip nanofriction in trapped cold ion chains, Phys. Rev. B 87, 195418 (2013)

  11. [11]

    Bylinskii, D

    A. Bylinskii, D. Gangloff, and V. Vuleti´ c, Tuning friction atom-by-atom in an ion-crystal simulator, Science348, 1115 (2015)

  12. [12]

    Kiethe, R

    J. Kiethe, R. Nigmatullin, T. Schmirander, D. Kalincev, and T. E. Mehlst¨ aubler, Nanofriction and motion of topo- logical defects in self-organized ion Coulomb crystals, New J. Phys.20, 123017 (2018)

  13. [13]

    Brazda, A

    T. Brazda, A. Silva, N. Manini, A. Vanossi, R. Guerra, E. Tosatti, and C. Bechinger, Experimental Observation of the Aubry Transition in Two-Dimensional Colloidal Monolayers, Phys. Rev. X8, 011050 (2018)

  14. [14]

    J. I. Cirac and P. Zoller, Quantum Computations with Cold Trapped Ions, Phys. Rev. Lett.74, 4091 (1995)

  15. [15]

    Blatt and D

    R. Blatt and D. Wineland, Entangled states of trapped atomic ions, Nature453, 1008 (2008)

  16. [16]

    Monroe and J

    C. Monroe and J. Kim, Scaling the Ion Trap Quantum Processor, Science339, 1164 (2013)

  17. [17]

    Arnold, E

    K. Arnold, E. Hajiyev, E. Paez, C. H. Lee, M. D. Barrett, and J. Bollinger, Prospects for atomic clocks based on large ion crystals, Phys. Rev. A92, 032108 (2015)

  18. [18]

    Keller, T

    J. Keller, T. Burgermeister, D. Kalincev, J. Kiethe, and T. E. Mehlst¨ aubler, Evaluation of trap-induced system- atic frequency shifts for a multi-ion optical clock at the 10-19 level, J. Phys. Conf. Ser.723, 012027 (2016)

  19. [19]

    Keller, H

    J. Keller, H. N. Hausser, I. M. Richter, T. Nord- mann, N. M. Bhatt, J. Kiethe, H. Liu, E. Benkler, B. Lipphardt, S. D¨ orscher, K. Stahl, J. Klose, C. Lis- dat, M. Filzinger, N. Huntemann, E. Peik, and T. E. Mehlst¨ aubler, High-accuracy multi-ion spectroscopy with mixed-species Coulomb crystals, J. Phys. Conf. Ser. 2889, 012050 (2024)

  20. [20]

    P. O. Schmidt, T. Rosenband, C. Langer, W. M. Itano, J. C. Bergquist, and D. J. Wineland, Spectroscopy Using Quantum Logic, Science309, 749 (2005)

  21. [21]

    Champenois, M

    C. Champenois, M. Marciante, J. Pedregosa-Gutierrez, M. Houssin, M. Knoop, and M. Kajita, Ion ring in a lin- ear multipole trap for optical frequency metrology, Phys. Rev. A81, 043410 (2010)

  22. [22]

    Islam, E

    R. Islam, E. Edwards, K. Kim, S. Korenblit, C. Noh, H. Carmichael, G.-D. Lin, L.-M. Duan, C.-C. Joseph Wang, J. Freericks, and C. Monroe, Onset of a quantum phase transition with a trapped ion quantum simulator, Nat. Commun.2, 377 (2011)

  23. [23]

    J. W. Britton, B. C. Sawyer, A. C. Keith, C.-C. J. Wang, J. K. Freericks, H. Uys, M. J. Biercuk, and J. J. Bollinger, Engineered two-dimensional Ising interactions in a trapped-ion quantum simulator with hundreds of spins, Nature484, 489 (2012)

  24. [24]

    J. G. Bohnet, B. C. Sawyer, J. W. Britton, M. L. Wall, A. M. Rey, M. Foss-Feig, and J. J. Bollinger, Quantum spin dynamics and entanglement generation with hun- dreds of trapped ions, Science352, 1297 (2016)

  25. [25]

    R¨ uffert, E

    L.-A. R¨ uffert, E. A. Dijck, L. Timm, J. R. C. L´ opez- Urrutia, and T. E. Mehlst¨ aubler, Domain formation and structural stabilities in mixed-species Coulomb crystals induced by sympathetically cooled highly charged ions, Phys. Rev. A110, 063110 (2024)

  26. [26]

    D. H. E. Dubin and T. M. O’Neil, Trapped nonneutral plasmas, liquids, and crystals (the thermal equilibrium states), Rev. Mod. Phys.71, 87 (1999)

  27. [27]

    Poindron, J

    A. Poindron, J. Pedregosa-Gutierrez, and C. Champ- enois, Thermal bistability in laser-cooled trapped ions, Phys. Rev. A104, 043116 (2021)

  28. [28]

    Baldovin, G

    M. Baldovin, G. Vallet, G. Hagel, E. Trizac, and C. Champenois, Self-diffusion in a strongly coupled non- neutral plasma, Phys. Rev. A109, 043116 (2024)

  29. [29]

    Zaris, W

    J. Zaris, W. Johnson, A. Shankar, J. J. Bollinger, and S. E. Parker, Numerical simulations of three-dimensional ion crystal dynamics in a penning trap using the fast multipole method, J. Plasma Phys.91, E53 (2025)

  30. [30]

    S. L. Gilbert, J. J. Bollinger, and D. J. Wineland, Shell- Structure Phase of Magnetically Confined Strongly Cou- pled Plasmas, Phys. Rev. Lett.60, 2022 (1988)

  31. [31]

    R. W. Hasse and J. P. Schiffer, The structure of the cylin- drically confined Coulomb lattice, Ann. Phys.203, 419 (1990)

  32. [32]

    R. W. Hasse and V. V. Avilov, Structure and Madelung energy of spherical Coulomb crystals, Phys. Rev. A44, 4506 (1991)

  33. [33]

    Drewsen, C

    M. Drewsen, C. Brodersen, L. Hornekær, J. S. Hangst, and J. P. Schiffer, Large Ion Crystals in a Linear Paul Trap, Phys. Rev. Lett.81, 2878 (1998)

  34. [34]

    Radzvilaviˇ cius and E

    A. Radzvilaviˇ cius and E. Anisimovas, Topological defect motifs in two-dimensional Coulomb clusters, J. Phys.: Condens. Matter23, 385301 (2011)

  35. [35]

    Drewsen, Ion Coulomb crystals, Physica B460, 105 (2015)

    M. Drewsen, Ion Coulomb crystals, Physica B460, 105 (2015)

  36. [36]

    L. Duca, N. Mizukami, E. Perego, M. Inguscio, and C. Sias, Orientational Melting in a Mesoscopic System of Charged Particles, Phys. Rev. Lett.131, 083602 (2023)

  37. [37]

    Bonitz, V

    M. Bonitz, V. Golubnychiy, A. Filinov, and Y. Lozovik, Single-electron control of Wigner crystallization, Micro- electron. Eng.63, 141 (2002)

  38. [38]

    Golubnychiy, P

    V. Golubnychiy, P. Ludwig, A. Filinov, and M. Bonitz, Controlling intershell rotations in mesoscopic electron clusters, Superlattices Microstruct.34, 219 (2003)

  39. [39]

    Kiesenhofer, H

    D. Kiesenhofer, H. Hainzer, A. Zhdanov, P. C. Holz, M. Bock, T. Ollikainen, and C. F. Roos, Controlling two-dimensional Coulomb crystals of more than 100 ions in a monolithic radio-frequency trap, PRX Quantum4, 020317 (2023)

  40. [40]

    J. Ahn, Z. Xu, J. Bang, P. Ju, X. Gao, and T. Li, Ul- trasensitive torque detection with an optically levitated nanorotor, Nat. Nanotechnol.15, 89 (2020)

  41. [41]

    Ohira, T

    R. Ohira, T. Mukaiyama, and K. Toyoda, Breaking ro- tational symmetry in a trapped-ion quantum tunneling rotor, Phys. Rev. A101, 022106 (2020)

  42. [42]

    J. Shao, W. Zhu, X. Zhang, and Y. Zheng, Molecu- lar rotors with designed polar rotating groups possess mechanics-controllable wide-range rotational speed, npj Comput. Mater.6, 185 (2020)

  43. [43]

    Singhania, S

    A. Singhania, S. Kalita, P. Chettri, and S. Ghosh, Ac- counts of applied molecular rotors and rotary motors: Recent advances, Nanoscale Adv.5, 3177 (2023)

  44. [44]

    Paul, Electromagnetic traps for charged and neutral particles, Rev

    W. Paul, Electromagnetic traps for charged and neutral particles, Rev. Mod. Phys.62, 531 (1990)

  45. [45]

    Okada, M

    K. Okada, M. Wada, T. Takayanagi, S. Ohtani, and H. A. Schuessler, Characterization of ion Coulomb crystals in a linear Paul trap, Phys. Rev. A81, 013420 (2010)

  46. [46]

    Rahman and J

    A. Rahman and J. P. Schiffer, Structure of a One- Component Plasma in an External Field: A Molecular- Dynamics Study of Particle Arrangement in a Heavy-Ion Storage Ring, Phys. Rev. Lett.57, 1133 (1986)

  47. [47]

    H. Totsuji, Static and Dynamic Properties of Strongly- 16 Coupled Classical One-Component Plasmas: Numerical Experiments on Supercooled Liquid State and Simulation of Ion Plasma in the Penning Trap, inStrongly Coupled Plasma Physics, edited by F. J. Rogers and H. E. Dewitt (Springer US, Boston, MA, 1987) pp. 19–33

  48. [48]

    Dubin and T

    D. Dubin and T. O’Neil, Computer simulation of ion clouds in a Penning trap, Phys. Rev. Lett.60, 511 (1988)

  49. [49]

    Birkl, S

    G. Birkl, S. Kassner, and H. Walther, Multiple-shell structures of laser-cooled 24Mg+ ions in a quadrupole storage ring, Nature357, 310 (1992)

  50. [50]

    Prandtl, Ein Gedankenmodell zur kinetischen Theorie der festen K¨ orper, Z

    L. Prandtl, Ein Gedankenmodell zur kinetischen Theorie der festen K¨ orper, Z. Angew. Math. Mech.8, 85 (1928)

  51. [51]

    G. A. Tomlinson, A molecular theory of friction, Philos. Mag. Ser. 77, 905 (1929)

  52. [52]

    Frenkel and T

    Y. Frenkel and T. Kontorova, On the theory of plastic deformation and twinning, Zh. Eksp. Teor. Fiz.8, 1340 (1938), in Russian

  53. [53]

    Aubry, The New Concept of Transitions by Break- ing of Analyticity in a Crystallographic Model, inSoli- tons and Condensed Matter Physics, Vol

    S. Aubry, The New Concept of Transitions by Break- ing of Analyticity in a Crystallographic Model, inSoli- tons and Condensed Matter Physics, Vol. 8, edited by M. Cardona, P. Fulde, H.-J. Queisser, A. R. Bishop, and T. Schneider (Springer Berlin Heidelberg, Berlin, Heidel- berg, 1978) pp. 264–277

  54. [54]

    Peyrard and S

    M. Peyrard and S. Aubry, Critical behaviour at the tran- sition by breaking of analyticity in the discrete Frenkel- Kontorova model, J. Phys. C: Solid State Phys.16, 1593 (1983)

  55. [55]

    Kiethe, R

    J. Kiethe, R. Nigmatullin, D. Kalincev, T. Schmiran- der, and T. E. Mehlst¨ aubler, Probing nanofriction and Aubry-type signatures in a finite self-organized system, Nat. Commun.8, 15364 (2017)

  56. [56]

    D. A. Gangloff, A. Bylinskii, and V. Vuleti´ c, Kinks and nanofriction: Structural phases in few-atom chains, Phys. Rev. Res.2, 013380 (2020)

  57. [57]

    Voisin, F

    C. Voisin, F. Renard, and J.-R. Grasso, Long term fric- tion: From stick-slip to stable sliding, Geophys. Res. Lett.34(2007)

  58. [58]

    Gourdon and J

    D. Gourdon and J. N. Israelachvili, Transitions between smooth and complex stick-slip sliding of surfaces, Phys. Rev. E68, 021602 (2003)

  59. [59]

    Drummond and J

    C. Drummond and J. Israelachvili, Dynamic phase transi- tions in confined lubricant fluids under shear, Phys. Rev. E63, 041506 (2001)

  60. [60]

    Benassi, A

    A. Benassi, A. Vanossi, and E. Tosatti, Nanofriction in cold ion traps, Nat. Commun.2, 236 (2011)

  61. [61]

    Braiman, J

    Y. Braiman, J. Baumgarten, J. Jortner, and J. Klafter, Symmetry-breaking transition in finite Frenkel- Kontorova chains, Phys. Rev. Lett.65, 2398 (1990)

  62. [62]

    D. H. E. Dubin, Theory of structural phase transitions in a trapped Coulomb crystal, Phys. Rev. Lett.71, 2753 (1993)

  63. [63]

    J.-F. Han, B. Liu, and W.-S. Duan, The superlubricity of the special material with hexagonal symmetry in a two- dimensional Frenkel–Kontorova model, Indian J. Phys. 94, 521 (2020)

  64. [64]

    Cang-Long, D

    W. Cang-Long, D. Wen-Shan, Y. Yang, and C. Jian-Min, Application of Two-Dimensional Frenkel–Kontorova Model to Nanotribology, Commun. Theor. Phys.54, 112 (2010)

  65. [65]

    K. Pyka, J. Keller, H. L. Partner, R. Nigmatullin, T. Burgermeister, D. M. Meier, K. Kuhlmann, A. Ret- zker, M. B. Plenio, W. H. Zurek, A. Del Campo, and T. E. Mehlst¨ aubler, Topological defect formation and spontaneous symmetry breaking in ion Coulomb crystals, Nat. Commun.4, 2291 (2013)

  66. [66]

    Kirkpatrick, C

    S. Kirkpatrick, C. D. Gelatt, and M. P. Vecchi, Optimiza- tion by Simulated Annealing, Science220, 671 (1983)

  67. [67]

    Caracciolo, A

    S. Caracciolo, A. K. Hartmann, S. Kirkpatrick, and M. Weigel, Simulated annealing, optimization, searching for ground states (2023)

  68. [68]

    D. J. Berkeland, J. D. Miller, J. C. Bergquist, W. M. Itano, and D. J. Wineland, Minimization of ion micro- motion in a Paul trap, J. Appl. Phys.83, 5025 (1998)

  69. [69]

    Kalincev, L

    D. Kalincev, L. S. Dreissen, A. P. Kulosa, C.-H. Yeh, H. A. F¨ urst, and T. E. Mehlst¨ aubler, Motional heating of spatially extended ion crystals, Quantum Sci. Technol. 6, 034003 (2021)

  70. [70]

    Frausto-Sol´ ıs, E

    J. Frausto-Sol´ ıs, E. Li˜ n´ an-Garc´ ıa, and G. Santamar´ ıa- Bonfil, Tuned simulated annealing based on Boltzmann and Bose–Einstein distribution applied to MAXSAT problem, J. Asian Sci. Res.4, 14 (2014). Appendix A: Simulated Annealing To explore potential ground state configurations, we employ a simulated annealing (SA) algorithm [65, 66, 69] that ite...