Search NASA⌕ Search

Engineering topics

Reichhardt, C. J. O. (ORCID:0000000234875089)

Publications and source records attributed to Reichhardt, C. J. O. (ORCID:0000000234875089).

Depinning, melting, and sliding of generalized Wigner crystals in Moiré systems

We numerically examine the depinning, sliding, and melting of commensurate and incommensurate Wigner crystals on two-dimensional hexagonal periodic substrates near fillings of 1 / 3 , 1 / 2 , and 2 / 3 to model the dynamics of generalized Wigner crystals in moiré heterostructures. At low temperatures where thermal fluctuations are irrelevant, for commensurate fillings we find a strongly pinned state that depins into a sliding crystal, while at incommensurate fillings the depinning threshold is strongly reduced. For fillings below 1 / 2 , the depinning occurs in a two-step process. Above the first depinning threshold, there is an extended range of drives where the conduction occurs via the sliding of antikinks along the charge stripes, followed by a second threshold where all the charges begin to slide. At finite temperatures, the external driving reduces the effective melting temperature at commensurate fillings and enhances creep at incommensurate fillings. We show that depinning into different sliding states, such as moving fluids or moving crystals, produces nonlinear features in the transport curves. We also show that transport is asymmetric on either side of a commensurate filling due to the different dynamics for interstitials versus holes in the commensurate structure. A variety of sliding states should be accessible for generalized Wigner crystals at finite temperatures even for low drives. If experiments can realize sliding dynamics in moiré systems, it would open another class of commensurate and incommensurate phases for study. Additionally, the sliding dynamics and antikink flow could provide unique functionalities for charge transport in moiré heterostructures. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Stripe and bubble ratchets on asymmetric substrates

We show that various nonmonotonic ratchet effects can arise when mesophase pattern-forming systems, which exhibit anisotropic crystal, stripe, and bubble regimes, are coupled to one-dimensional asymmetric substrates under ac driving. The patterns emerge in the absence of a substrate when the ratio of attraction to repulsion is varied for particles with competing short-range attraction and long-range repulsion potentials. In the presence of the substrate, we find that the ratchet efficiency varies nonmonotonically with increasing attraction, depending upon how well the mesophase morphology matches the substrate spacing and periodicity. For strong repulsion, there is a weak but finite ratchet effect, while at intermediate attraction, there is a robust ratchet effect in which the system forms stripes aligned with the substrate symmetry direction. For strong attraction, large bubbles appear that have weak or no ratchet effects when the bubble width exceeds the substrate lattice spacing, causing the bubble to be only weakly coupled to the substrate. For very strong attraction, small bubbles form and undergo a strong ratchet effect with an efficiency that oscillates as a function of ac drive amplitude. Although the small bubble regime is a strongly correlated regime in which each bubble contains many particles, the system behaves as if it is in a single-particle regime since the bubble width is much smaller than the substrate lattice spacing. We map out the different rectification phases as a function of the pattern morphology, substrate strength, and ac drive amplitude. The pronounced ratchet effects that we observe in some regimes can be exploited for pattern sorting in hard and soft matter systems. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Skyrmion soliton motion on periodic substrates by atomistic and particle-based simulations

Abstract We compare the dynamical behavior of magnetic skyrmions interacting with square and triangular defect arrays just above commensuration using both atomistic and particle-based models. Under applied drives, the initial motion is a kink traveling through the pinned skyrmion lattice. For the square defect array, both models agree well and show a regime in which the soliton motion is locked along 45°. The atomistic model also produces locking of a soliton along 30°, which is absent in the particle-based model. For the triangular defect array, the atomistic model exhibits 30° soliton motion over a wide region of external current values. In contrast, the particle-based model gives 45° soliton motion over a small range of external driving force values. The difference arises because the nondeforming particle model facilitates meandering skyrmion orbits while the deformable atomistic model enables stronger skyrmion-skyrmion interactions that reduce the meandering. Our results indicate that soliton motion through pinned skyrmion lattices on a periodic substrate is a robust effect.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Reversible to irreversible transitions for ac driven skyrmions on periodic substrates

Abstract Using atomistic simulations, we investigate the dynamical behavior of magnetic skyrmions in dimer and trimer molecular crystal arrangements, as well as bipartite lattices at 3/2 and 5/2 fillings, under ac driving over a square array of anisotropy defects. For low ac amplitudes, at all fillings reversible motion appears in which the skyrmions return to their original positions at the end of each ac drive cycle and the diffusion is zero. We also identify two distinct irreversible regimes. The first is a translating regime in which the skyrmions form channels of flow in opposing directions and translate by one substrate lattice constant per ac drive cycle. The translating state appears in the dimer and trimer arrangements, and produces pronounced peaks in the diffusivity in the direction perpendicular to the external drive. For larger ac amplitudes, we find chaotic irreversible motion in which the skyrmions can randomly exchange places with each other over time, producing long-time diffusive behavior both parallel and perpendicular to the ac driving direction.

36 MATERIALS SCIENCE↗

Sliding dynamics for bubble phases on periodic modulated substrates

We analyze a bubble-forming system composed of particles with competing long-range repulsive and short-range attractive interactions driven over a quasi-one-dimensional periodic substrate. We find various pinned and sliding phases as a function of substrate strength and drive amplitude. When the substrate is weak, a pinned bubble phase appears that depins elastically into a sliding bubble lattice. For stronger substrates, we find anisotropic bubbles, disordered bubbles, and stripe phases. Plastic depinning occurs via the hopping of individual particles from one bubble to the next in a pinned bubble lattice, and as the drive increases, there is a transition to a state where all of the bubbles are moving but are continuously shedding and absorbing individual particles. This is followed at high drives by a moving bubble lattice in which the particles can no longer escape their individual bubbles. The transition between the plastic and elastic sliding phases can be detected via signatures in the velocity-force curves, differential conductivity, and noise. When the bubbles shrink due to an increase in the attractive interaction term, they fit better inside the pinning troughs and become more strongly pinned, leading to a reentrant pinning phase. For weaker attractive terms, the size of the bubbles becomes greater than the width of the pinning troughs and the depinning becomes elastic with a reduced depinning threshold. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Activated vortex lattice transition in a superconductor with combined sixfold and twelvefold anisotropic interactions

Abstract Numerical simulations are used to examine the transition dynamics between metastable and ground state vortex lattice phases in a system with combined sixfold and twelvefold contributions to the vortex–vortex interactions. The system is initially annealed using a twelvefold anisotropy, yielding domains of two different orientations and separated by grain boundaries. The vortex–vortex interaction is then suddenly changed to a sixfold anisotropy, rendering the twelvefold state metastable. Applying a drive that mimics an oscillating magnetic field causes the metastable state to decay, indicated by the structure factor that evolves from twelve to six peaks. The results fit the behavior seen in recent small-angle neutron scattering studies of the vortex lattice in MgB 2 . At higher drive amplitudes, the decay exhibits a two step process in which the initial fast decrease is followed by a slower regime where avalanches or bursts are correlated with dislocation annihilation events. The results are compared to other types of metastable systems with quenched disorder that decay under a periodic external drive.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Spontaneous skyrmion conformal lattice and transverse motion during dc and ac compression

Abstract We use atomistic-based simulations to investigate the behavior of ferromagnetic skyrmions being continuously compressed against a rigid wall under dc and ac drives. The compressed skyrmions can be annihilated close to the wall and form a conformal crystal with both a size and a density gradient, making it distinct from conformal crystals observed previously for superconducting vortices and colloidal particles. For both dc and ac driving, the skyrmions can move transverse to the compression direction due to a combination of density and size gradients. Forces in the compression direction are converted by the Magnus force into transverse motion. Under ac driving, the amount of skyrmion annihilation is reduced and we find a skyrmion Magnus ratchet pump. We also observe shear banding in which skyrmions near the wall move up to twice as fast as skyrmions further from the wall. When we vary the magnitude of the applied drive, we find a critical current above which the skyrmions are completely annihilated during a time scale that depends on the magnitude of the drive. By varying the magnetic parameters, we find that the transverse motion is strongly dependent on the skyrmion size. Smaller skyrmions are more rigid, which interferes with the size gradient and destroys the transverse motion. We also confirm the role of the size gradient by comparing our atomistic simulations with a particle-based model, where we find that the transverse motion is only transient. Our results are relevant for applications where skyrmions encounter repulsive magnetic walls, domain walls, or interfaces.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Melting, reentrant ordering and peak effect for Wigner crystals with quenched and thermal disorder

Abstract We consider simulations of Wigner crystals in solid state systems interacting with random quenched disorder in the presence of thermal fluctuations. When quenched disorder is absent, there is a well defined melting temperature determined by the proliferation of topological defects, while for zero temperature, there is a critical quenched disorder strength above which topological defects proliferate. When both thermal and quenched disorder are present, these effects compete, and the thermal fluctuations can reduce the effectiveness of the quenched disorder, leading to a reentrant ordered phase in agreement with the predictions of Nelson (1983 Phys. Rev. B 27 2902–14). There are two competing theories for the low temperature behavior, and our simulations show that both capture aspects of the actual response. The critical disorder strength separating ordered from disordered states remains finite as the temperature goes to zero, as predicted by Cha and Fertig (1995 Phys. Rev. Lett. 74 4867–70), instead of dropping to zero as predicted by Nelson. At the same time, the critical disorder strength decreases with decreasing temperature, as predicted by Nelson, instead of remaining constant, as predicted by Cha and Fertig. The onset of the reentrant phase can be deduced based on changes in the transport response, where the reentrant ordering appears as an increase in the mobility or the occurrence of a depinning transition. We also find that when the system is in the ordered state and thermally melts, there is an increase in the effective damping or pinning. This produces a drop in the electron mobility that is similar to the peak effect phenomenon found in superconducting vortices, where thermal effects soften the lattice or break down its elasticity, allowing the particles to better adjust their positions to take full advantage of the quenched disorder.

36 MATERIALS SCIENCE↗

Clogging, diode and collective effects of skyrmions in funnel geometries

Abstract Using a particle-based model, we examine the collective dynamics of skyrmions interacting with a funnel potential under dc driving as the skyrmion density and relative strength of the Magnus and damping terms are varied. For driving in the easy direction, we find that increasing the skyrmion density reduces the average skyrmion velocity due to jamming of skyrmions near the funnel opening, while the Magnus force causes skyrmions to accumulate on one side of the funnel array. For driving in the hard direction, there is a critical skyrmion density below which the skyrmions become trapped. Above this critical value, a clogging effect appears with multiple depinning and repinning states where the skyrmions can rearrange into different clogged configurations, while at higher drives, the velocity-force curves become continuous. When skyrmions pile up near the funnel opening, the effective size of the opening is reduced and the passage of other skyrmions is blocked by the repulsive skyrmion–skyrmion interactions. We observe a strong diode effect in which the critical depinning force is higher and the velocity response is smaller for hard direction driving. As the ratio of Magnus force to dissipative term is varied, the skyrmion velocity varies in a non-linear and non-monotonic way due to the pile up of skyrmions on one side of the funnels. At high Magnus forces, the clogging effect for hard direction driving is diminished.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Driven superconducting vortex dynamics in systems with twofold anisotropy in the presence of pinning

Abstract We examine the dynamics of superconducting vortices with twofold anisotropic interaction potentials driven over random pinning, and compare the behavior under drives applied along the hard and the soft anisotropy directions. As the driving force increases, the number of topological defects reaches a maximum near the depinning threshold, and then decreases as the vortices form one-dimensional (1D) chains. This coincides with a transition from a pinned nematic to a moving smectic aligned with the soft anisotropy direction. The system is generally more ordered when the drive is applied along the soft direction of the anisotropy. For driving along the hard direction, there is a critical value of the twofold anisotropy above which the system remains aligned with the soft direction. Hysteretic behavior appears upon cycling the driving force, with 1D vortex chains persisting during the decreasing leg below the threshold for chain formation for increasing drive. More anisotropic systems have a greater amount of structural disorder in the moving state. For lower anisotropy, the system forms a moving smectic-A state, while at higher anisotropy, a moving nematic state appears instead.

36 MATERIALS SCIENCE↗

Kibble-Zurek mechanism for nonequilibrium phase transitions in driven systems with quenched disorder

Abstract We show that the Kibble-Zurek mechanism applies to nonequilibrium phase transitions found in driven assemblies of superconducting vortices and colloids moving over quenched disorder where a transition occurs from a plastic disordered flowing state to a moving anisotropic crystal. We measure the density of topological defects as a function of quench rate through the nonequilibrium phase transition, and find that on the ordered side of the transition, the topological defect density ρ d scales as a power law, $${\rho }_{{{{{{{{\rm{d}}}}}}}}}\propto 1/{t}_{{{{{{{{\rm{q}}}}}}}}}^{\beta }$$ ρ d ∝ 1 / t q β , where t q is the quench time duration, consistent with the Kibble-Zurek mechanism. We show that scaling with the same exponent holds for varied strengths of quenched disorder and that the exponents fall in the directed percolation universality class. Our results suggest that the Kibble-Zurek mechanism can be applied to the broader class of systems that exhibit absorbing phase transitions.

97 MATHEMATICS AND COMPUTING↗

Skyrmion pinball and directed motion on obstacle arrays

Abstract We examine skyrmions interacting with a square array of obstacles under ac drives applied in one or two directions. For a single direction of ac driving, we find that the Magnus force in conjunction with the obstacle interactions can create elliptical skyrmion orbits of increasing size, leading to localized phases, chaotic phases, and translating or ratcheting orbits. Under two ac drives that are out of phase by 90° and applied in two directions, the skyrmions form localized commensurate orbits that encircle an integer number of obstacles, similar to the electron pinball effect observed for electrons in antidot lattices. As a function of ac amplitude, Magnus force strength, and obstacle size, we find that chaotic scattering regimes and directed motion can emerge even in the absence of asymmetry in the substrate. The directed motion follows different symmetry axes of the periodic substrate, and we observe a variety of reversed ratchet effects. The Magnus force in the skyrmion system produces a significantly larger number of directed motion regimes than are exhibited by overdamped systems. We discuss how these results could be used to move skyrmions in a controlled way for possible applications.

97 MATHEMATICS AND COMPUTING↗

Skyrmion dynamics and topological sorting on periodic obstacle arrays

Abstract We examine skyrmions under a dc drive interacting with a square array of obstacles for varied obstacle size and damping. When the drive is applied in a fixed direction, we find that the skyrmions are initially guided in the drive direction but also move transverse to the drive due to the Magnus force. The skyrmion Hall angle, which indicates the difference between the skyrmion direction of motion and the drive direction, increases with drive in a series of quantized steps as a result of the locking of the skyrmion motion to specific symmetry directions of the obstacle array. On these steps, the skyrmions collide with an integer number of obstacles to create a periodic motion. The transitions between the different locking steps are associated with jumps or dips in the velocity–force curves. In some regimes, the skyrmion Hall angle is actually higher than the intrinsic skyrmion Hall angle that would appear in the absence of obstacles. In the limit of zero damping, the skyrmion Hall angle is 90°, and we find that it decreases as the damping increases. For multiple interacting skyrmion species in the collective regime, we find jammed behavior at low drives where the different skyrmion species are strongly coupled and move in the same direction. As the drive increases, the species decouple and each can lock to a different symmetry direction of the obstacle lattice, making it possible to perform topological sorting in analogy to the particle sorting methods used to fractionate different species of colloidal particles moving over two-dimensional obstacle arrays.

36 MATERIALS SCIENCE↗