Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “specific phase transitions”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Development of HPC based phase field simulations tool for modification of alloy morphology to enhance material properties during additive manufacturing (AM) process

Currently, additive manufacturing (AM) technology is extensively used to manufacture aerospace components because of the unique advantages of the technique in producing custom designs, complex geometries in a near-net fashion that significantly reduces the manufacturing cost, and improved energy efficiency. However, the thermodynamic properties of existing aerospace alloys and the thermal conditions that prevail in most of the powder-bed based additive manufacturing processes lead to undesirable columnar / columnar-dendritic microstructures that promote solidification cracking during the process, and the development of anisotropy in the mechanical properties of the as-built components. Time consuming, energy-intensive post AM treatments are required to recover the mechanical properties of the wrought components for which the existing alloy compositions were designed for. However, novel alloy compositions, that exploit the AM thermal conditions to produce fine, equiaxed grain structures in the as-built condition, with mechanical properties equal to or exceeding those of the wrought work-horse alloy, Ti-6Al-4V has recently been demonstrated in Ti-Cu alloys. Significant energy savings up to 66% can be obtained for a typical aerospace component through digitally designing the AM process, process parameters, and alloy composition. The current project focused on developing a fundamental understanding of the evolution of the columnar to equiaxed transition (CET) occurring during solidification under AM thermal conditions in Ti-Cu and Ti-Cu-X alloys using high-fidelity phase field (PF) simulations using the leadership class computing facilities that exist in the national laboratories. The simulations were able to capture the effect of alloy and process conditions on CET, and clearly showed the beneficial effect of a ternary solute addition to Ti-Cu binary alloys on CET. The simulations indicated that bulk nucleation in the liquid ahead of the solidifying epitaxial front was promoted by a large, transient, thermal undercooling promoted by the mismatch between the growth rate of the dendrite tips in the epitaxial front and the imposed solidification rate. AM experiments performed at Raytheon Technologies Research Center using the same thermal conditions used in the PF simulations indicated the formation of equiaxed grains, although the grain size was bigger than the ones reported in the literature. On the other hand, the equiaxed grain size predicted by PF simulations were smaller than the grain size of Ti-Cu alloys reported in the literature. The industry will follow up the existing work in a future program that would involve further optimization of the process for the Ti-Cu-X alloys to demonstrate the application of the process and the alloy to a specific aerospace component.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Performance Evaluation of a Microgrid System with Grid-Forming and Grid-Following Inverters with Diesel Generators: Insights From Hardware Experiments

This paper presents a comprehensive performance evaluation of a microgrid system integrating grid-forming (GFM) inverters, grid-following (GFL) inverters, and a diesel generator, focusing on their interactions and behavior under various dynamic events. The study is conducted using a pure hardware setup comprising two GFM inverters, one GFL inverter, a diesel generator, load banks, a point of common coupling (PCC), and an emulated main grid. The evaluation specifically examines dynamic scenarios, including voltage jumps, phase jumps, rate of change of frequency (ROCOF), synchronization, and islanding operations, which pose critical challenges to system stability. Among all the tests conducted, the overloading and phase jump tests proved to be the most challenging. The capacity of the DC side is crucial for withstanding overloading and grid disturbance tests; otherwise, GFM inverters frequently trip due to DC undervoltage. Throughout all grid disturbance tests, the diesel generator consistently stands out as the most robust and reliable GFM unit in the system. Overall, insights from these hardware experiments shed light on the response characteristics of different generation types during grid disturbances and identify potential stability concerns in such hybrid microgrid systems.

14 SOLAR ENERGY↗
DOE OSTI2024

Probabilistic measures for biological adaptation and resilience

This paper introduces an approach to quantifying ecological resilience in biological systems, particularly focusing on noisy systems responding to episodic disturbances with sudden adaptations. Incorporating concepts from nonequilibrium statistical mechanics, we propose a measure termed “ecological resilience through adaptation,” specifically tailored to noisy, forced systems that undergo physiological adaptation in the face of stressful environmental changes. Randomness plays a key role, accounting for model uncertainty and the inherent variability in the dynamical response among components of biological systems. Our measure of resilience is rooted in the probabilistic description of states within these systems and is defined in terms of the dynamics of the ensemble average of a model-specific observable quantifying success or well-being. Our approach utilizes stochastic linear response theory to compute how the expected success of a system, originally in statistical equilibrium, dynamically changes in response to a environmental perturbation and a subsequent adaptation. Importantly, the resulting mathematical derivations allow for the estimation of resilience in terms of ensemble averages of simulated or experimental data. Finally, through a simple but clear conceptual example, we illustrate how our resilience measure can be interpreted and compared to other existing frameworks in the literature. The methodology is general but inspired by applications in plant systems, with the potential for broader application to complex biological processes.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Time-Resolved X-ray Emission Spectroscopy and Resonant Inelastic X-ray Scattering Spectroscopy of Laser Irradiated Carbon

The existence of liquid carbon as an intermediate phase preceding the formation of novel carbon materials has been a point of contention for several decades. Experimental observation of such a liquid state requires nonthermal melting of solid carbon materials at various laser fluences and pulse properties. Reflectivity experiments performed in the mid-1980s reached opposing conclusions regarding the metallic or insulating properties of the purported liquid state. Time-resolved X-ray absorption studies showed shortening of C–C bonds and increasing diffraction densities, thought to evidence a liquid or glassy carbon state, respectively. Nevertheless, none of these experiments provided information on the electronic structure of the proposed liquid state. Herein, we report the results of time-resolved resonant inelastic X-ray scattering (RIXS) and time-resolved X-ray emission spectroscopy (XES) studies on amorphous carbon (a-C) and ultrananocrystalline diamond (UNCD) as a function of delay time between the irradiating pulse and X-ray probe. For both a-C and UNCD, we attribute decreases in RIXS or XES signals to transition blocking, relaxation, and finally, ablation. Increased signal at 20 ps following the irradiation of the UNCD is attributed to the probable formation of nanoscale structures in the ablation plume. Further, differences in the amount of signal observed between a-C and UNCD are explained by the difference in sample thickness and, specifically, incomplete melting of the UNCD film. Comparisons to spectral simulations based on MD trajectories at extreme conditions indicate that the carbon state in our experiments is crystalline. Normal mode analysis confirmed that symmetrical bending or stretching of the C–C bonds in the diamond lattice results in XES spectra with small intensity differences. Overall, we observed no evidence of melting to a liquid state, as determined by the lack of changes in the spectral properties for up to 100 ps delays following the melting pulses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

A General Model for the Ductility of Intermetallics Applied to Fe-Co Alloys

The mechanical properties, specifically ductility, of high performing soft magnets such as Fe-Co alloys are a limiting factor to their broader use in a number of systems. The understanding of the mechanical robustness in these materials is currently insufficient to be able to support the growing interest in applications such as magnetic shielding or electric motors. Fe-Co alloys provide the highest commercially available magnetic saturation and high magnetic permeability but have very poor ductility. The addition of vanadium to these alloys has allowed for significant commercialization and some ductility improvements, but the fundamental reasons for the observed improvements are not well understood. In most published work on mechanical properties in these alloys, the precise chemistry of the alloys investigated, often a critical aspect of intermetallics, is not reported or controlled and thermal history is unclear. This work creates a ductility model that is sensitive to changes in chemistry and can predict relative strain to failure as well as brittle fracture mode for intermetallics and is applied to Fe-Co alloys. Through the application of density functional theory (DFT), this model identifies defect stabilities, anti-phase boundary (APB) energies, the energy necessary to cross-slip, and cleavage energies and combines them through energetic competition to determine a relative failure strain. The model correctly predicts ductility improvements with the addition of vanadium as well as the transition from intergranular to transgranular cleavage, though more precise experiments are necessary to appropriately validate the various improvements observed.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Ice Phase Classification Made Easy with Score-Based Denoising

Accurate identification of ice phases is essential for understanding various physicochemical phenomena. However, such classification for structures simulated with molecular dynamics is complicated by the complex symmetries of ice polymorphs and thermal fluctuations. For this purpose, both traditional order parameters and data-driven machine learning approaches have been employed, but they often rely on expert intuition, specific geometric information, or large training data sets. In this work, we present an unsupervised phase classification framework that combines a score-based denoiser model with a subsequent model-free classification method to accurately identify ice phases. Further, the denoiser model is trained on perturbed synthetic data of ideal reference structures, eliminating the need for large data sets and labeling efforts. The classification step utilizes the smooth overlap of atomic position (SOAP) descriptors as the atomic fingerprint, ensuring Euclidean symmetries and transferability to various structural systems. Our approach achieves a remarkable 100% accuracy in distinguishing ice phases of test trajectories using only seven ideal reference structures of ice phases as model inputs. This demonstrates the generalizability of the score-based denoiser model in facilitating phase identification for complex molecular systems. The proposed classification strategy can be broadly applied to investigate structural evolution and phase identification for a wide range of materials, offering new insights into the fundamental understanding of water and other complex systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Unconventional fieldlike spin torques in CrPt 3

The topological semimetal CrPt 3 has potential for generating unconventional spin torques due to its ferrimagnetic ordering, topological band structure, and high anomalous Hall effect. CrPt 3 exhibits ferrimagnetic behavior only in its chemically ordered phase and is paramagnetic in its chemically disordered phase. By controlling the growth and annealing temperatures, we prepare epitaxial films of both chemically ordered and chemically disordered phases of CrPt 3 , allowing us to investigate the effect of magnetic ordering on unconventional-torque generation. We use angle-dependent spin-torque-ferromagnetic-resonance and second-harmonic Hall measurements to probe the spin torques generated from epitaxial CrPt 3 in CrPt 3 /Cu/Ni 81 ⁢Fe 19 heterostructures. With current applied along specific directions with respect to the crystal order, we reveal unconventional spin torques in both ordered and disordered films. When current flows parallel to the [$1\bar{1}1$] and [ $\bar{1}11$] directions, we observe an unconventional fieldlike torque that is opposite in sign for the two directions. Finally, our calculations reveal that this unconventional torque originates from an indirect nonlocal spin-orbit torque due to spin scattering at the CrPt 3 /Cu interface, in addition to symmetry breaking at this interface.

ferromagnetic resonance↗
DOE OSTI2025

Liquid Phase Modeling in Porous Media: Adsorption of Methanol and Ethanol in H-MFI in Condensed Water

Zeolites are used in the chemical and separation industries for their exceptional selectivity, adsorption capacity, regenerability, and stability in gas and liquid phase processing. Here, we developed an explicit solvation method for predicting solvent/condensed phase effects on adsorption free energies in microporous media such as zeolites based on the hybrid quantum mechanical/molecular mechanical free energy perturbation (QM/MM-FEP) technique. Our explicit solvation method for zeolite systems, called eSZS, aims to capture site-specific interactions during the adsorption process at the Brønsted acid sites of H-MFI zeolite while still considering the diverse configuration space of the solvent molecules. This strategy is ideal for chemical reactions or adsorbates that interact with the microporous medium in few distinct adsorbate/transition state configurations, i.e., the harmonic or similar approximations are acceptable for the adsorbate/transition state while such approximations break down for the solvent molecules that require extensive configuration space sampling. In this way, our approach effectively overcomes the limitations of implicit solvation models and classical force field methods for describing solvation effects on chemical reactions within porous materials such as zeolites. Specifically, in this study, we investigated various aspects of our hybrid QM/MM approach, including QM cluster size dependencies in a periodic electrostatically embedded cluster model (PEECM), rules for link atoms at the QM/MM boundary, and functional and basis set considerations for converged and reasonably accurate gas and aqueous phase methanol and ethanol adsorption free energy predictions in H-MFI. For gas phase adsorption of methanol and ethanol in H-MFI at a Brønsted acid site in T12 position, we compute adsorption free energies at 298 K of −0.61 and −0.75 eV, respectively, using a PEECM containing 50 Si and 1 Al atom with ωB97x-D/def2-TZVP level of theory. For solvent effect calculations, we sample the aqueous phase using grand canonical Monte Carlo (GCMC) simulations to (1) obtain a mean field of electrostatic interactions in the reaction system and (2) perform a rigorous free energy perturbation calculation. Similar to the experimentally and computationally observed endergonic solvation effects observed for hydrocarbon adsorption on metal surfaces, we also observe that a condensed aqueous environment destabilizes methanol and ethanol at these acid sites in H-MFI at 298 K. Specifically, the computed solvation free energies of adsorption (ΔΔG solv ) for methanol and ethanol are +0.44 and +0.54 eV, respectively. From this study, it is evident that adsorbates (methanol and ethanol) are competing with water for adsorption space inside the H-MFI zeolite, leading to an endergonic solvation effect. Here, we expect that the endergonic, aqueous solvent effect during adsorption in microporous zeolites is highly tunable by changing the pore size and hydrophobicity of the microporous material as this will affect the water density inside the pore structure.

Adsorption↗
DOE OSTI2025

Phenomena Identification and Ranking Table (PIRT) for heat pipes

Heat pipes are advanced passive thermal management devices that utilize phase change and capillary action to achieve efficient heat transfer. However, due to the complexity of the phenomena coupled in heat pipes, including capillary, phase change, turbulence, and compressibility effects, there are high uncertainties in the predictability of their operational regimes and performance. This PIRT exercise, conducted as a collaborative effort involving the Department of Energy (DOE) Microreactor Program (MRP), the Nuclear Regulatory Commission (NRC), and university partners systematically identifies, reviews, and prioritizes critical phenomena affecting the operation of heat pipes based on their importance and knowledge levels. Additional analyses and discussion are provided for phenomena with high importance and low knowledge, such as wick de-wetting, critical heat flux, contact angles, and pressure dynamics. The discussions included the recognizing challenges and proposing future research directions for both modeling and simulation and experimental efforts. Additionally, the report addresses phenomena with medium importance and low knowledge that could impact heat pipe operation during non-normal or transient operation, including frozen startup, laminar to turbulent transition, geyser boiling, wick priming, underfilling conditions, surface roughness of the wick, NCGs trapped in the wick, and the timescales of startup and shutdown. In conclusion, this comprehensive evaluation serves as a valuable resource for guiding future research and development efforts, supporting the successful integration of heat pipes into critical applications such as nuclear reactors, and contributing to the advancement of heat pipe technologies in safety-critical industries.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2023

Thermodynamic investigation of the NaCl-KCl salt system from 25 to 950 °C

NaCl-KCl molten salt system has been proposed as a primary component of several promising heat transfer eutectics and fuel host for Generation IV molten salt reactors (MSR). In this work, several key thermodynamic parameters, including high temperature enthalpies, isobaric heat capacities (C p ) and molar enthalpies of mixing (ΔH mix ), of the NaCl -KCl system (NaCl, 75mol% NaCl – 25mol% KCl, 51mol% NaCl – 49mol% KCl, 25mol% NaCl – 75mol% KCl, and KCl) were measured by high temperature drop calorimetry (HTDC) utilizing laser sealed nickel and aluminum crucibles. Salts were checked for phase purity by X-ray diffraction (XRD), and water content by Karl Fischer coulometric titrimetry and thermogravimetric analysis (TGA). Thermo-mechanical analysis (TMA) utilizing custom-built boron nitride (BN) crucibles was employed to determine temperatures of phase transitions (solid $\rightarrow$ liquid) for the salts and plotted against the pseudobinary phase diagram for the NaCl-KCl system. Small deviations (3-10 %) were found among C p values obtained from HTDC, differential scanning calorimetry (DSC), and NIST, comparable to the average uncertainty based on two standard deviations of measured data that are ~6% for DSC and ~3% for HTDC. Additionally, the measured enthalpies and C p values of mixed NaCl-KCl salts (e.g., 51mol% NaCl – 49mol% KCl eutectic) were shown to behave as statistical mixtures of its pure endmembers, consistent with ΔH mix = ~ 0 kJ/mol at 800 °C, and a regular interaction parameter of Ω = 1.20 ± 0.76 kJ/mol. In conclusion, this study represents a step towards improved accuracy and precision in determinations of thermodynamic parameters of molten salts by using DSC, TMA, and HTDC, in completion of the Molten Salt Thermal Properties Database (MSTDB).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Two-step electronic response to magnetic ordering in a van der Waals ferromagnet

The two-dimensional material Cr 2 ⁢Ge 2 ⁢Te 6 is a member of the class of insulating van der Waals (vdW) magnets. Here, using high resolution angle-resolved photoemission spectroscopy in a detailed temperature dependence study, we identify a clear response of the electronic structure to a dimensional crossover in the form of two distinct temperature scales marking onsets of modifications in the electronic structure. Specifically, we observe Te 𝑝-orbital-dominated bands to undergo changes at the Curie transition temperature 𝑇 𝐶 while the Cr 𝑑-orbital-dominated bands begin evolving at a higher temperature scale. Combined with neutron scattering, density functional theory calculations, and Monte Carlo simulations, we find that the electronic system can be consistently understood to respond sequentially to the distinct temperatures at which in-plane and out-of-plane spin correlations exceed a characteristic length scale. Finally, our findings reveal the sensitivity of the orbital-selective electronic structure for probing the dynamical evolution of local moment correlations in vdW insulating magnets.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Universal finite-time thermodynamics of many-body quantum machines from Kibble-Zurek scaling

We demonstrate the existence of universal features in the finite-time thermodynamics of quantum machines by considering a many-body quantum Otto cycle in which the working medium is driven across quantum critical points during the unitary strokes. Specifically, we consider a quantum engine powered by dissipative energizing and relaxing baths. We show that under very generic conditions, the output work is governed by the Kibble-Zurek mechanism; i.e., it exhibits a universal power-law scaling with the driving speed through the critical points. We also optimize the finite-time thermodynamics as a function of the driving speed. The maximum power and the corresponding efficiency take a universal form, and are reached for an optimal speed that is governed by the critical exponents. We exemplify our results by considering a transverse-field Ising spin chain as the working medium. For this model, we also show how the efficiency and power vary as the engine becomes critical.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Nanoscale Engineering of Polymorphism in Cu 2 Se-Based Composites

Crystal polymorphism selection during synthesis is extremely challenging. However, promoting the formation of a specific metastable polymorph enables modulation of the functional properties of phase-change materials through alteration of the relative abundance of various polymorphs. Here, we demonstrate the stabilization of the superionic β-Cu 2 Se phase under ambient conditions and the direct control over the relative ratio between the α-Cu 2 Se and β-Cu 2 Se polymorphs in (x)CuGaSe 2 /(1-x)Cu 2 Se composites using CuGaSe 2 nanoseeds. Here, we found that the small lattice mismatch between β-Cu 2 Se (cubic) and the ab plane of tetragonal CuGaSe 2 nanoseeds promotes the formation of low-energy coherent CuGaSe 2 /β-Cu 2 Se interfaces leading to preferential stabilization of β-Cu 2 Se. Astonishingly, the hierarchical microstructure of the resulting composites enables a remarkable decoupling of charge and heat transport, which is manifested by a breakdown of the Wiedemann-Franz law

36 MATERIALS SCIENCE↗
DOE OSTI2020

A practical guide for crystal growth of van der Waals layered materials

This Tutorial provides an overview of the techniques that are most commonly utilized to grow bulk van der Waals crystals. The materials discussed were selected to highlight various challenges that are often encountered during crystal growth. Additionally, in relatively equal parts, the text covers melt-based techniques, vapor transport growths, and the characterization of crystal quality with an emphasis on structural and chemical homogeneities. Furthermore, pertinent details are given regarding the growth and characterization of many specific compounds, with examples mostly drawn from our own research, and an effort is made to highlight cases where the growths offer a particular lesson or the conditions have a significant impact on the crystal’s physical properties. A primary goal is to motivate more researchers to grow crystals by providing general descriptions and considerations for different growth techniques and equipment while sharing some of our own lessons learned and best practices for the growth and characterization of layered van der Waals crystals. The Tutorial is not written solely for aspiring crystal growers, however, because any researcher who collaborates with a crystal grower can benefit from having a greater understanding and appreciation of the processes of crystal growth and materials development.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Disorder recovers the Wiedemann-Franz law in the metallic phase of VO 2

At temperatures higher than 341 K, vanadium dioxide (VO 2 ) is a strongly correlated metal with resistivity exceeding the Mott-Ioffe-Regel limit. Its electronic thermal conductivity is lower than that predicted by the Wiedemann-Franz (WF) law, and can be explained by nonquasiparticle transport where heat and charge currents follow separate diffusive modes. In contradiction, the Wiedemann-Franz law is a direct consequence of quasiparticle transport where charge carriers are elastically scattered. In this work, we enhance elastic electron scattering in VO 2 by introducing atomic disorder with ion irradiation. Furthermore, a gradual and eventually full recovery of the WF law is observed at high defect densities. This observation provides an example that connects hydrodynamic quasiparticle transport to nonquasiparticle transport in metallic systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Witnessing Nonequilibrium Entanglement Dynamics in a Strongly Correlated Fermionic Chain

Many-body entanglement in condensed matter systems can be diagnosed from equilibrium response functions through the use of entanglement witnesses and operator-specific quantum bounds. Here, we investigate the applicability of this approach for detecting entangled states in quantum systems driven out of equilibrium. We use a multipartite entanglement witness, the quantum Fisher information, to study the dynamics of a paradigmatic fermion chain undergoing a time-dependent change of the Coulomb interaction. Our results show that the quantum Fisher information is able to witness distinct signatures of multipartite entanglement both near and far from equilibrium that are robust against decoherence. Here, we discuss implications of these findings for probing entanglement in light-driven quantum materials with time-resolved optical and x-ray scattering methods.

1-dimensional spin chains↗
DOE OSTI2022

Role of electromagnetic gauge-field fluctuations in the selection between chiral and nematic superconductivity

Motivated by the observation of nematic superconductivity in several systems, we revisit the problem of the leading pairing instability of two-component unconventional superconductors on the triangular lattice—such as (p x , p y )-wave and (d x 2 –y 2 , d xy )-wave superconductors. Such a system has two possible superconducting states: the chiral state (e.g., p + ip or d + id), which breaks time-reversal symmetry, and the nematic state (e.g., p + p or d + d), which breaks the threefold rotational symmetry of the lattice. Weak-coupling calculations generally favor the chiral over the nematic superconducting state, raising the question of what mechanism can stabilize the latter. Here, we show that the electromagnetic field fluctuations can play a crucial role in selecting between these two states. Specifically, we derive and analyze the effective free energy for the two-component superconducting order parameter after integrating out the gauge-field fluctuations, which is formally justified if the spatial order parameter fluctuations can be neglected. A nonanalytic cubic term arises, as in the case of a conventional s -wave superconductor. However, unlike the latter, the cubic term depends on the relative phase and on the relative amplitudes between the two order parameter components, in such a way that it generally favors the nematic state. This result is a direct consequence of the fact that the stiffness of the superconducting order parameter is not isotropic. Competition with the quartic term, which favors the chiral state, leads to a renormalized phase diagram in which the nematic state displaces the chiral state over a wide region in the parameter space. Here, we analyze the stability of the fluctuation-induced nematic phase, generalize our results to tetragonal lattices, and discuss their applicability to candidate nematic superconductors, including twisted bilayer graphene.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

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↗