Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spin models”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Artificial spin ice: Paths forward

Artificial spin ice systems are lithographically prepared arrays of interacting nanoscale magnetic moments with collective behavior resulting from the chosen array geometry. These many-body systems are unusual in that their simple constituent elements can be configured to design their interactions with exquisite control. They can then be probed experimentally over a remarkably large range of time scales and length scales, including imaging of the individual moments. The study of artificial spin ice has broadened well beyond the original focus on systems with connection to spin ice and ice models, and now the term is used in reference to many other structures that exhibit a much wider range of physical phenomena. Here in this perspective, we review progress in the field of artificial spin ice since its inception more than a decade ago. We then discuss prospects for future directions – continuing the exploration of nanomagnetism, statistical spin models, and the potential for technological devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Enhancing scalability of a matrix-free eigensolver for studying many-body localization

We propose several techniques to enhance the parallel scalability of a matrix-free eigensolver designed for studying many-body localization (MBL) of quantum spin chain models with nearest-neighbor interactions and on-site disorder. This type of problem is computationally challenging because the dimension of the associated Hamiltonian matrix grows exponentially with respect to the number of spins L, and we need to average over different realizations of the random disorder to obtain relevant statistical behavior. For each disorder realization, we need to compute eigenvalues from different regions of the spectrum and their corresponding eigenvectors. In previous work, the interior eigenstates for a single eigenvalue problem are computed via the shift-and-invert Lanczos algorithm. Due to the extremely high memory footprint of the LU factorizations, this technique is not well suited for large L’s. For example, we need thousands of compute nodes on modern high performance computing infrastructures to go beyond L = 24. The matrix-free approach does not suffer from this memory bottleneck, however, its scalability is limited by a computation and communication load imbalance. To reduce this imbalance and to significantly enhance the scalability of the matrix-free eigensolver, we reorder the matrix and leverage the consistent space runtime, CSPACER. We also show its efficiency in managing irregular communication patterns at scale compared to optimized MPI non-blocking two-sided and one-sided RMA implementation variants. This effort enables us to study MBL for spin chains with a larger number of spins. The efficiency and effectiveness of the proposed algorithm is demonstrated by computing eigenstates on a massively parallel many-core high performance computer.

METIS↗
DOE OSTI2025

Fluorescence quantum beats reveal tunable decoherence and the effect of transport in spin-entangled exciton pairs

Spin-decoherence in a population of entangled exciton pairs in an organic semiconductor can be tuned by the orientation of an applied magnetic field. Magnetic field dependent measurements of photoinduced fluorescence quantum beats and a model of spin coherence that accounts for transport-induced dephasing are used to derive the exciton hopping rate between inequivalent sites in a crystal lattice.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Spin Waves and Magnetic Exchange Hamiltonian in CrSBr

Abstract CrSBr is an air‐stable two‐dimensional (2D) van der Waals semiconducting magnet with great technological promise, but its atomic‐scale magnetic interactions—crucial information for high‐frequency switching—are poorly understood. An experimental study is presented to determine the CrSBr magnetic exchange Hamiltonian and bulk magnon spectrum. The A ‐type antiferromagnetic order using single crystal neutron diffraction is confirmed here. The magnon dispersions are also measured using inelastic neutron scattering and rigorously fit the excitation modes to a spin wave model. The magnon spectrum is well described by an intra‐plane ferromagnetic Heisenberg exchange model with seven nearest in‐plane exchanges. This fitted exchange Hamiltonian enables theoretical predictions of CrSBr behavior: as one example, the fitted Hamiltonian is used to predict the presence of chiral magnon edge modes with a spin‐orbit enhanced CrSBr heterostructure.

2D materials↗
DOE OSTI2020

Accelerating Scientific Computing in the Post-Moore’s Era

Novel uses of graphical processing units for accelerated computation revolutionized the field of high-performance scientific computing by providing specialized workflows tailored to algorithmic requirements. As the era of Moore’s law draws to a close, many new non–von Neumann processors are emerging as potential computational accelerators, including those based on the principles of neuromorphic computing, tensor algebra, and quantum information. While development of these new processors is continuing to mature, the potential impact on accelerated computing is anticipated to be profound. We discuss how different processing models can advance computing in key scientific paradigms: machine learning and constraint satisfaction. Significantly, each of these new processor types utilizes a fundamentally different model of computation, and this raises questions about how to best use such processors in the design and implementation of applications. While many processors are being developed with a specific domain target, the ubiquity of spin-glass models and neural networks provides an avenue for multi-functional applications. Furthermore, this also hints at the infrastructure needed to integrate next-generation processing units into future high-performance computing systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Spin–phonon interactions and magnetoelectric coupling in Co 4 B 2 O 9 ( B = Nb, Ta)

In order to explore the consequences of spin–orbit coupling on spin–phonon interactions in a set of chemically similar mixed metal oxides, we measured the infrared vibrational properties of Co 4 B 2 O 9 (B = Nb, Ta) as a function of temperature and compared our findings with lattice dynamics calculations and several different models of spin–phonon coupling. Frequency vs temperature trends for the Co 2+ shearing mode near 150 cm –1 reveal significant shifts across the magnetic ordering temperature that are especially large in relative terms. Bringing these results together and accounting for noncollinearity, we obtain spin–phonon coupling constants of –3.4 and –4.3 cm –1 for Co 4 Nb 2 O 9 and the Ta analog, respectively. Analysis reveals that these coupling constants are derived from interlayer (rather than intralayer) exchange interactions and that the interlayer interactions contain competing antiferromagnetic and ferromagnetic contributions. At the same time, beyond-Heisenberg terms are minimized due to fortuitous symmetry considerations, different from most other 4d- and 5d-containing oxides. Comparison with other contemporary oxides shows that spin–phonon coupling in this family of materials is among the strongest ever reported, suggesting an origin for magnetoelectric coupling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Spinning black hole scattering at $ \mathcal{O} $(G 3 S 2 )(G 3 S 2 ): Casimir terms, radial action and hidden symmetry

We resolve subtleties in calculating the post-Minksowskian dynamics of binary systems, as a spin expansion, from massive scattering amplitudes of fixed finite spin. In particular, the apparently ambiguous spin Casimir terms can be fully determined from the gradient of the spin-diagonal part of the amplitudes with respect to S 2 = −s(s+1)ħ 2 , using an interpolation between massive amplitudes with different spin representations. From two-loop amplitudes of spin-0 and spin-1 particles minimally coupled to gravity, we extract the spin Casimir terms in the conservative scattering angle between a spinless and a spinning black hole at $ \mathcal{O} $(G 3 S 2 ), finding agreement with known results in the literature. This completes an earlier study [Phys. Rev. Lett. 130 (2023), 021601] that calculated the non-Casimir terms from amplitudes. We also illustrate our methods using a model of spinning bodies in electrodynamics, finding agreement between scattering amplitude predictions and classical predictions in a root-Kerr electromagnetic background up to $ \mathcal{O} $(α 3 S 2 ). For both gravity and electrodynamics, the finite part of the amplitude coincides with the two-body radial action in the aligned spin limit, generalizing the amplitude-action relation beyond the spinless case. Surprisingly, the two-loop amplitude displays a hidden spin-shift symmetry in the probe limit, which was previously observed at one loop. We conjecture that the symmetry holds to all orders in the coupling constant and is a consequence of integrability of Kerr orbits in the probe limit at the first few orders in spin.

Classical Theories of Gravity↗
DOE OSTI2024

Emergence of the isotropic Kitaev honeycomb lattice $α$–RuCl 3 and its magnetic properties

Here, we present a comprehensive investigation of the crystal and magnetic structures of the van der Waals antiferromagnet $α$–RuCl 3 using single crystal x-ray and neutron diffraction. The crystal structure at room temperature is a monoclinic (C2/m). However, with decreasing temperature, a remarkable first-order structural phase transition is observed, leading to the emergence of a rhombohedral ($R\bar3$) structure characterized by three-fold rotational symmetry forming an isotropic honeycomb lattice. On further cooling, a zigzag-type antiferromagnetic order develops below $T_{\textrm{N}} = 6\sim6.6$ K. The critical exponent of the magnetic order parameter was determined to be $\beta = 0.11(1)$, which is close to the two-dimensional Ising model. Additionally, the angular dependence of the magnetic critical field of the zigzag antiferromagnetic order for the polarized ferromagnetic phase reveals a six-fold rotational symmetry within the ab–plane. These findings reflect the symmetry associated with the Ising-like bond-dependent Kitaev spin interactions and underscore the universality of the Kitaev interaction-dominated antiferromagnetic system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Beyond Kitaev physics in strong spin-orbit coupled magnets

We review the recent advances and current challenges in the field of strong spin-orbit coupled Kitaev materials, with a particular emphasis on the physics beyond the exactly-solvable Kitaev spin liquid point. To this end, we present a comprehensive overview of the key exchange interactions in candidate materials with a specific focus on systems featuring effective J eff = 1/2 magnetic moments. This includes, but not limited to, 5d 5 iridates, 4d 5 ruthenates and 3d 7 cobaltates. Our exploration covers the microscopic origins of these interactions, along with a systematic attempt to map out the most intriguing correlated regimes of the multi-dimensional parameter space. Our approach is guided by robust symmetry and duality transformations as well as insights from a wide spectrum of analytical and numerical studies. We also survey higher spin Kitaev models and recent exciting results on quasi-one-dimensional models and discuss their relevance to higher-dimensional models. Lastly, we highlight some of the key questions in the field as well as future directions.

Kitaev materials↗
DOE OSTI2025

Spin Waves in Dirac Semimetal Ca 0.6⁢ Sr 0.4⁢ MnSb 2 Investigated with Neutrons by the Diffraction Method

The tunability of Dirac semimetals with antiferromagnetic Mn layers is of great interest. The observed sign change of interlayer magnetic coupling between CaMnBi 2 and SrMnBi 2 suggests ionic substitution as a potential tuning mechanism. If so, novel behavior near the compensation point could be expected. Here, to explore this, we study a mixed-cation analog, Ca 0.6⁢ Sr 0.4⁢ MnSb 2 , where Bi is replaced by Sb. Conventional inelastic neutron scattering is impractical due to the small crystal size (m ≈ 0.28 g) available for compositional studies; however, we find that using a neutron diffractometer with a wide-angle area detector we can obtain a good quality spin-wave signal, which is shaped by energy-momentum conservation and retains spectroscopic information even without direct energy analysis. Spin-wave modeling reveals an interlayer coupling quantitatively similar to SrMnSb 2 , indicating it is not directly tuned by ionic size and that the sign change in Bi-based compounds likely arises from the observed change in lattice symmetry. Beyond this key insight, our results present an efficient method for parametric and compositional studies of spin dynamics in small crystals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Angular Momentum Removal by Neutron and γ-Ray Emissions during Fission Fragment Decays

In this work, we investigate the angular momentum removal from fission fragments (FFs) through neutron and γ-ray emission, finding that about half the neutrons are emitted with angular momenta ≥ 1.5 ℏ and that the change in angular momentum after the emission of neutrons and statistical γ rays is significant, contradicting usual assumptions. Per fission event, in our simulations, the neutron and statistical γ-ray emissions change the spin of the fragment by 3.5 – 5 ℏ , with a large standard deviation comparable to the average value. Such wide angular momentum removal distributions can hide any underlying correlations in the fission fragment initial spin values. Within our model, we reproduce data on spin measurements from discrete transitions after neutron emissions, especially in the case of light FFs. The agreement further improves for the heavy fragments if one removes from the analysis the events that would produce isomeric states. Finally, we show that while in our model the initial FF spins do not follow a sawtoothlike behavior observed in recent measurements, the average FF spin computed after neutron and statistical γ emissions exhibits a shape that resembles a sawtooth. This suggests that the average FF spin measured after statistical emissions is not necessarily connected with the scission mechanism as previously implied.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Spin dynamics of the generalized quantum spin compass chain

Here we calculate the dynamical spin structure factor of the generalized spin-1/2 compass spin chain using the density matrix renormalization group. The model, also known as the twisted Kitaev spin chain, was recently proposed to be relevant for the description of the spin chain compound CoNb 2 O 6 . It features bond-dependent interactions and interpolates between an Ising chain and a one-dimensional variant of Kitaev's honeycomb spin model. The structure factor, in turn, is found to interpolate from gapped and nondispersive in the Ising limit to gapless with nontrivial continua in the Kitaev limit. In particular, the component of the structure factor perpendicular to the Ising directions changes abruptly at the Kitaev point into a dispersionless continuum due to the emergence of an extensive ground-state degeneracy. We show this continuum is consistent with analytical Jordan-Wigner results. We also discuss implications for future inelastic scattering experiments and applications to materials, particularly CoNb 2 O 6 .

1-dimensional spin chains↗
DOE OSTI2021

Scaling functions for eigenstate entanglement crossovers in harmonic lattices

For quantum matter, eigenstate entanglement entropies obey an area law or log-area law at low energies and small subsystem sizes and cross over to volume laws for high energies and large subsystems. This transition is captured by crossover functions, which assume a universal scaling form in quantum critical regimes. We demonstrate this for the harmonic lattice model, which describes quantized lattice vibrations and is a regularization for free scalar field theories, modeling, e.g., spin-0 bosonic particles. In one dimension, the ground-state entanglement obeys a log-area law. For dimensions d ≥ 2, it displays area laws, even at criticality. The distribution of excited-state entanglement entropies is found to be sharply peaked around subsystem entropies of corresponding thermodynamic ensembles in accordance with the eigenstate thermalization hypothesis. Numerically, we determine crossover scaling functions for the quantum critical regime of the model and do a large-deviation analysis. We show how infrared singularities of the system can be handled and how to access the thermodynamic limit using a perturbative trick for the covariance matrix. Eigenstates for quasi-free bosonic systems are not Gaussian. Here, we resolve this problem by considering appropriate squeezed states instead. For these, entanglement entropies can be evaluated efficiently.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Three-dimensional checkerboard spin structure on a breathing pyrochlore lattice

The standard approach to realize a spin-liquid state is through magnetically frustrated states, relying on ingredients such as the lattice geometry, dimensionality, and magnetic interaction type of the spins. While Heisenberg spins on a pyrochlore lattice with only antiferromagnetic nearest-neighbor interactions are theoretically proven disordered, spins in real systems generally include longer-range interactions. The spatial correlations at longer distances typically stabilize a long-range order rather than enhancing a spin-liquid state. Both states can, however, be destroyed by short-range static correlations introduced by chemical disorder. Here, using disorder-free specimens with a clear long-range antiferromagnetic order, we refine the spin structure of the Heisenberg spinel ZnFe 2 O 4 through neutron magnetic diffraction. The unique wave vector (1, 0, $\frac{1}{2}$) leads to a spin structure that can be viewed as alternatively stacked ferromagnetic and antiferromagnetic tetrahedra in a three-dimensional checkerboard form. Stable coexistence of these opposing types of clusters is enabled by the bipartite breathing pyrochlore crystal structure, leading to a second-order phase transition at 10 K. The diffraction intensity of ZnFe 2 O 4 is an exact complement to the inelastic scattering intensity of several chromate spinel systems which are regarded as model classical spin liquids. Our results challenge this attribution, and suggest instead of the six-spin ring mode, spin excitations in chromate spinels are closely related to the (1, 0, $\frac{1}{2}$) type of spin order and the four-spin ferromagnetic cluster locally at one tetrahedron.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Optimizing temperature distributions for training neural quantum states using parallel tempering

Parametrized artificial neural networks (ANNs) can be very expressive ansatzes for variational algorithms, reaching state-of-the-art energies on many quantum many-body Hamiltonians. Nevertheless, the training of the ANN can be slow and stymied by the presence of local minima in the parameter landscape. One approach to mitigate this issue is to use parallel tempering methods, and in this work, we focus on the role played by the temperature distribution of the parallel tempering replicas. Using an adaptive method that adjusts the temperatures in order to equate the exchange probability between neighboring replicas, we show that this temperature optimization can significantly increase the success rate of the variational algorithm with negligible computational cost by eliminating bottlenecks in the replicas' random walk. Furthermore, we demonstrate this using two different neural networks, a restricted Boltzmann machine and a feedforward network, which we use to study a toy problem based on a permutation invariant Hamiltonian with a pernicious local minimum and the 𝐽 1 −𝐽 2 model on a rectangular lattice.

Neural network simulations↗
DOE OSTI2025

Strong Hilbert space fragmentation and fractons from subsystem and higher-form symmetries

Here, we introduce a route to Hilbert space fragmentation in high dimensions leveraging the group-word formalism. We show that taking strongly fragmented models in one dimension and “lifting” to higher dimensions using subsystem symmetries can yield strongly fragmented dynamics in higher dimensions, with subdimensional (e.g., lineonic) excitations. This provides a route to higher-dimensional strong fragmentation, and also a route to fractonic behavior. Meanwhile, lifting one-dimensional fragmented models to higher dimensions using higher-form symmetries yields models with topologically robust fragmentation. In three or more spatial dimensions, one can also “mix and match” subsystem and higher-form symmetries, leading to canonical fracton models such as X cube. We speculate that this approach could also yield a route to non-Abelian fractons. These constructions unify a number of phenomena that have been discussed in the literature, as well as furnishing models with unique properties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Magnetic field-temperature phase diagram of the spin-$\frac{1}{2}$ triangular lattice antiferromagnet KYbSe 2

A quantum spin liquid (QSL) is a state of matter characterized by fractionalized quasiparticle excitations, quantum entanglement, and a lack of long-range magnetic order. However, QSLs have evaded definitive experimental observation. Several Yb 3+ -based triangular lattice antiferromagnets with effective 𝑆 = $\frac{1}{2}$ have been suggested to stabilize the QSL state as the ground state. Here, in this work, we build a comprehensive magnetic temperature phase diagram of a high-quality single crystalline KYbSe 2 via heat capacity and magnetocaloric effect down to 30 mK with magnetic field applied along the 𝑎 axis. At zero magnetic field, we observe the magnetic long-range order at 𝑇 N =0.29 K entering 120 degrees ordered state in heat capacity, consistent with neutron scattering studies. Analysis of the low-temperature (𝑇) specific heat (𝐶) at zero magnetic field indicates linear 𝑇 dependence of 𝐶/𝑇 and a broad hump of 𝐶/𝑇 in the proximate QSL region above 𝑇 N . By applying magnetic field, we observe the up-up-down phase with 1/3 magnetization plateau and oblique phases, in addition to two new phases. These observations strongly indicate that while KYbSe 2 closely exhibits characteristics resembling an ideal triangular lattice, deviations may exist, such as the effect of the next-nearest-neighbor exchange interaction, calling for careful consideration for spin Hamiltonian modeling. Further investigations into tuning parameters, such as chemical pressure, could potentially induce an intriguing QSL phase in the material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Toward Accelerating Discovery via Physics-Driven and Interactive Multifidelity Bayesian Optimization

Both computational and experimental material discovery bring forth the challenge of exploring multidimensional and often nondifferentiable parameter spaces, such as phase diagrams of Hamiltonians with multiple interactions, composition spaces of combinatorial libraries, processing spaces, and molecular embedding spaces. Often these systems are expensive or time consuming to evaluate a single instance, and hence classical approaches based on exhaustive grid or random search are too data intensive. This resulted in strong interest toward active learning methods such as Bayesian optimization (BO) where the adaptive exploration occurs based on human learning (discovery) objective. However, classical BO is based on a predefined optimization target, and policies balancing exploration and exploitation are purely data driven. In practical settings, the domain expert can pose prior knowledge of the system in the form of partially known physics laws and exploration policies often vary during the experiment. Here, we propose an interactive workflow building on multifidelity BO (MFBO), starting with classical (data-driven) MFBO, then expand to a proposed structured (physics-driven) structured MFBO (sMFBO), and finally extend it to allow human-in-the-loop interactive interactive MFBO (iMFBO) workflows for adaptive and domain expert aligned exploration. These approaches are demonstrated over highly nonsmooth multifidelity simulation data generated from an Ising model, considering spin–spin interaction as parameter space, lattice sizes as fidelity spaces, and the objective as maximizing heat capacity. Detailed analysis and comparison show the impact of physics knowledge injection and real-time human decisions for improved exploration with increased alignment to ground truth. Here, the associated notebooks allow to reproduce the reported analyses and apply them to other systems.

97 MATHEMATICS AND COMPUTING↗