Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Ising Model”

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 307 records · Page 17

The susceptibility critical exponent for a nonaqueous ionic binary mixture near a consolute point

We report turbidity measurements of a nonaqueous ionic solution of triethyl n-hexylammonium triethyl n-hexylboride in diphenyl ether. A classical susceptibility critical exponent gamma = 1.01 +/- 0.01 is obtained over the reduced temperature range t between values of 0.1 and 0.0001. The best fits of the sample transmission had a standard deviation of 0.39 percent over this range. Ising and spherical model critical exponents are firmly excluded. The correlation length amplitude xi sub 0 from fitting is 1.0 +/- 0.2 nm which is much larger than values found in neutral fluids and some aqueous binary mixtures.

Zhang, Kai C.↗

Defect-phase-dynamics approach to statistical domain-growth problem of clock models

The growth of statistical domains in quenched Ising-like p-state clock models with p = 3 or more is investigated theoretically, reformulating the analysis of Ohta et al. (1982) in terms of a phase variable and studying the dynamics of defects introduced into the phase field when the phase variable becomes multivalued. The resulting defect/phase domain-growth equation is applied to the interpretation of Monte Carlo simulations in two dimensions (Kaski and Gunton, 1983; Grest and Srolovitz, 1984), and problems encountered in the analysis of related Potts models are discussed. In the two-dimensional case, the problem is essentially that of a purely dissipative Coulomb gas, with a sq rt t growth law complicated by vertex-pinning effects at small t.

Kawasaki, K.↗

Kibble–Zurek mechanism of Ising domains

The formation of topological defects after a symmetry-breaking phase transition is an overarching phenomenon that encodes the underlying dynamics. The Kibble-Zurek mechanism (KZM) describes these non-equilibrium dynamics of second-order phase transitions and predicts a power-law relationship between the cooling rates and the density of topological defects. It has been verified as a successful model in a wide variety of physical systems, including structure formation in the early Universe and condensed-matter materials. However, it is uncertain if the KZM mechanism is also valid for topologically trivial Ising domains, one of the most common and fundamental types of domain in condensed-matter systems. Here we show that the cooling rate dependence of Ising domain density follows the KZM power law in two different three-dimensional structural Ising domains: ferro-rotation domains in NiTiO 3 and polar domains in BiTeI. However, although the KZM slope of NiTiO 3 agrees with the prediction of the 3D Ising model, the KZM slope of BiTeI exceeds the theoretical limit, providing an example of steepening KZM slope with long-range dipolar interactions. Finally, our results demonstrate the validity of KZM for Ising domains and reveal an enhancement of the power-law exponent for transitions of non-topological quantities with long-range interactions. The Kibble-Zurek mechanism is shown to apply to structural Ising domains in three-dimensional materials. Long-range interactions modify the critical exponents away from theoretical predictions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

An evaluation of three two-dimensional computational fluid dynamics codes including low Reynolds numbers and transonic Mach numbers

Full-potential, Euler, and Navier-Stokes computational fluid dynamics (CFD) codes were evaluated for use in analyzing the flow field about airfoils sections operating at Mach numbers from 0.20 to 0.60 and Reynolds numbers from 500,000 to 2,000,000. The potential code (LBAUER) includes weakly coupled integral boundary layer equations for laminar and turbulent flow with simple transition and separation models. The Navier-Stokes code (ARC2D) uses the thin-layer formulation of the Reynolds-averaged equations with an algebraic turbulence model. The Euler code (ISES) includes strongly coupled integral boundary layer equations and advanced transition and separation calculations with the capability to model laminar separation bubbles and limited zones of turbulent separation. The best experiment/CFD correlation was obtained with the Euler code because its boundary layer equations model the physics of the flow better than the other two codes. An unusual reversal of boundary layer separation with increasing angle of attack, following initial shock formation on the upper surface of the airfoil, was found in the experiment data. This phenomenon was not predicted by the CFD codes evaluated.

Hicks, Raymond M.↗

Adaptive variational ground state preparation for spin-1 models on qubit-based architectures

Here, we apply the adaptive variational quantum imaginary time evolution (AVQITE) method to prepare ground states of one-dimensional spin S = 1 models. We compare different spin-to-qubit encodings (standard binary, Gray, unary, and multiplet) with regard to the performance and quantum resource cost of the algorithm. Using state-vector simulations, we study two well-known spin-1 models: the Blume-Capel model of transverse-field Ising spins with single-ion anisotropy, and the XXZ model with single-ion anisotropy. We consider system sizes of up to 20 qubits, which corresponds to spin-1 chains up to length 10. We determine the dependence of the number of CNOT gates in the AVQITE state preparation circuit on the encoding, the initial state, and the choice of operator pool in the adaptive method. Independent of the choice of encoding, we find that the CNOT gate count scales cubically with the number of spins for the Blume-Capel model and quartically for the anisotropic XXZ model. However, the multiplet and Gray encodings present smaller prefactors in the scaling relations. These results provide useful insights for the implementation of AVQITE on quantum hardware.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Unified theory of hastatic order and antiferromagnetism in URu 2 Si 2

The hidden-order phase of URu 2 Si 2 has eluded identification for over 35 years. A compelling proposal that explains the Ising heavy-fermion nature of the material is hastatic order: a symmetry-breaking heavy Fermi liquid arising from a spinorial microscopic hybridization. The original hastatic proposal cannot microscopically model the pressure-induced antiferromagnetic phase; and while it predicts a spinorial order parameter, it does not provide any detectable signatures of the spinorial nature. Here, we present a more realistic microscopic model of hastatic order in URu 2 Si 2 based on two conducting electron bands and explore its phase diagram in detail. Our model nontrivially preserves the Ising heavy-fermion physics of the original, while allowing us to tune between the antiferromagnet and hidden order using pressure analogs and magnetic field. Furthermore, our model is also consistent with recent phenomenological predictions of not one, but two, vector order parameters associated with the spinorial order that provide the first microscopic predictions for detecting the spinorial nature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NeuroSymbolic Approaches as a Vector for Assured Artificial Intelligence

The deployment of artificial intelligence systems in critical applications requires higher levels of assurance for safety, security, and interpretability. While neurosymbolic (NESY) approaches combining neural networks with symbolic reasoning offer potential advantages for assured AI, existing differentiable neurosymbolic frameworks face significant limitations including computational overhead and performance constraints. This report investigates the ISED (InferSampleEstimateDescend) framework as an alternative approach that enables neurosymbolic learning without requiring endtoend differentiability. We evaluate ISED’s utility for geointelligence applications by comparing neurosymbolic models against standard neural networks on aircraft classification tasks using the RarePlanes and MTARSI imagery datasets. Our results demonstrate that while ISEDbased models achieve slightly lower accuracy (89.7% vs 92.1% on RarePlanes; 91.1% vs 92.5% on MTARSI), they provide critical explainability capabilities that enable tracing incorrect predictions back to specific attribute misclassifications. We also present an automated pipeline that generates both attributeclass mappings and neurosymbolic model architectures from natural language descriptions, significantly reducing the manual effort required for NESY model deployment. These findings suggest that ISED offers a promising direction for developing assured AI systems where interpretability and reasoning transparency are prioritized alongside performance.

97 MATHEMATICS AND COMPUTING↗

Evidence for slow-mode shocks in the deep geomagnetic tail

Plasma and field data from the ISEE-3 space probe provide evidence that the lobe-plasma sheet boundary is the central part of the distant geomagnetic tail is often a slow-mode shock. Such shocks are predicted by Petchek's (1964) model of the reconnection. If this model applies, then the ISEE-3 observations place the general location of the reconnection diffusion region closer to the earth than x of about -100 earth radii.

Feldman, W. C.↗

Dynamic scaling of order parameter fluctuations in model B

Brookhaven National Laboratory has performed a dedicated set of experiments, called the Beam Energy Scan (BES), with the goal of location a possible critical point in the phase diagram of nuclear matter, This point is analogous to the endpoint of the liquid-gas phase transition in water, where the liquid phase corresponds to a phase of hadrons, and the gas phase to a plasma of quarks and gluons. In order to describe the data generated by these experiments we have studied numerical simulations of the stochastic diffusion equation with a conserved charge. Here, we focus on the dynamics in the vicinity of a critical point in the Ising universality class. The model we consider is expected to describe the critical dynamics near a possible QCD critical point if the coupling of the order parameter to the momentum density of the fluid can be neglected. The simulations are performed on a spatial lattice, and the time evolution is performed using a Metropolis algorithm. We determine the dynamical critical exponent z ≃ 3.972(2), which agrees with predictions of the epsilon expansion. We also study non-equilibrium sweeps of the reduced temperature and observe approximate Kibble-Zurek scaling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Magnetic field draping against the dayside magnetopause

Interplanetary magnetic fields observed upstream of Earth's magnetosphere at ISEE 3 form input for a gasdynamic model of magnetic field draping in the dayside magnetosheath. Model results near the magnetopause are compared with appropriately lagged observations at ISEE 1. In 16 to 24 cases, the angle between the transverse component of the model and observed fields is less than 20 deg. The agreement is surprisingly good in view of the uncertainty introduced by the large distances between ISEE 1 and ISEE 3. The results indicate that magnetohydrodynamic and energy transfer processes at the magnetopause do not cause large distortions of the magnetosheath magnetic field. In addition, a comparison between observed and model field magnitudes indicates that immediately outside the magnetopause the observed field behaves like the model field at a distance of approx. 0.5 R sub E from the magnetopause, outside the region where magnetohydrodynamic effects make the gasdynamic model inapplicable. Patterns of model magnetic field orientation at the magnetopause are presented for practical application.

Crooker, N. U.↗

Magnetic field draping against the dayside magnetopause

Interplanetary magnetic fields observed upstream of earth's magnetosphere at ISEE 3 form input for a gasdynamic model of magnetic field draping in the dayside magnetosheath. Model results near the magnetopause are compared with appropriately lagged observations at ISEE 1. In 16 to 24 cases, the angle between the transverse component of the model and observed fields is less than 20 deg. The agreement is surprisingly good in view of the uncertainty introduced by the large distances between ISEE 1 and ISEE 3. The results indicate that magnetohydrodynamic and energy transfer processes at the magnetopause do not cause large distortions of the magnetosheath magnetic field. In addition, a comparison between observed and model field magnitudes indicates that immediately outside the magnetopause the observed field behaves like the model field at a distance of approx. 0.5 R sub E from the magnetopause, outside the region where magnetohydrodynamic effects make the gasdynamic model inapplicable. Patterns of model magnetic field orientation at the magnetopause are presented for practical application.

Crooker, N. U.↗

Heliomagnetic latitude dependence of the heliospheric magnetic field

Previous studies have revealed systematic variations of the interplanetary magnetic field with heliographic latitude. Luhmann et al. (1987) modeled Pioneer Venus (PVO) and ISEE-3 observations by assuming an asymmetric dependence on heliolatitude with stronger fields in the northern hemisphere. In a subsequent study, using data from ISEE-3/ICE and IMP-8, Burton et al. (1990) found evidence for a similar asymmetry. However, neither model has been completely successful. The model derived from PVO/ICE observations agrees quite well near solar maximum but shows significant discrepancies during the descending phase of the solar cycle. The model derived from the ICE/IMP-8 comparison suffers from significant phase delays between the difference in field magnitude at the two spacecraft and their latitude difference. In an attempt to account for these phase shifts, the IMP-8 and ICE data have been reexamined in heliomagnetic coordinates which are defined by the orientation of the solar magnetic dipole. The latitude and longitude of the dipole inferred from the data have then been compared with those implicit in source surface calculations. The IMP/ICE correlations have been extended into the recent solar maximum and descending phase. Comparisons have also been carried out between IMP-8 and Ulysses as it traveled to -30 deg south heliographic latitude.

Burton, M. E.↗

Diffusive equilibrium distributions of He(+) in the plasmasphere

Recent satellite observations of thermal ion composition in the near-equatorial plasmasphere have shown that He(+) comprises 5-10 percent typically and occasionally 25 percent or more of the total thermal ion density. A steady state diffusive equillibrium model for the distributions of H(+), He(+) and O(+) along a plasmaspheric flux tube is used to elicit effects that may help explain these observed high He(+) fractional concentrations. The model indicates that both the ionospheric composition and the temperature distribution along the flux tubes are important factors controlling the equatorial He(+) composition, through the phlasma scale height and thermal diffusion effects. Direct comparison of the model results with thermal ion observations by ISEE-1 indicates that the effects incorporated into the model may explain some of the elevated He(+) concentrations. In some instances, however, effects not included in the model may be of importance.

Waite, J. H., Jr.↗

Surface magnetism in Fe 3 GeTe 2 van der Waals ferromagnet

The surface magnetization of Fe 3 GeTe 2 was examined by low-energy electron microscopy (LEEM) using an off-normal incidence electron beam. We found that the 180o domain walls are of Bloch type. Temperature-dependent LEEM measurements yield a surface magnetization with a surface critical exponent β1 = 0.79±0.02. This result is consistent with surface magnetism in the 3D semi-infinite Heisenberg (β1 = 0.84±0.01) or Ising (β1 = 0.78±0.02) models, which is distinctly different from the bulk exponent (β = 0.34±0.07). Further, the measurements reveal the power of LEEM with a tilted beam to determine magnetic domain structure in quantum materials, without need for the use of spin-polarized electrons. Single crystal diffraction measurements reveal inversion symmetry-breaking weak peaks and yield space group P-6m2. This Fe site defect-derived loss of inversion symmetry enables the formation of skyrmions in this Fe 3 GeTe 2 crystal.

36 MATERIALS SCIENCE↗

Phase diagram of generalized XY model using the tensor renormalization group

We use the higher-order tensor renormalization group method to study the two-dimensional generalized XY model that admits integer and half-integer vortices. This model is the deformation of the classical XY model and has a rich phase structure consisting of nematic, ferromagnetic, and disordered phases and three transition lines belonging to the Berezinskii-Kosterlitz-Thouless and Ising class. We explore the model for a wide range of temperatures, T , and the deformation parameter, Δ , and compute specific heat along with integer and half-integer magnetic susceptibility, finding both Berezinskii-Kosterlitz-Thouless-like and Ising-like transitions and the region where they meet. Published by the American Physical Society 2024

2-dimensional systems↗

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↗

Two-fermion negativity and confinement in the Schwinger model

We consider the fermionic (logarithmic) negativity between two fermionic modes in the Schwinger model. Recent results pointed out that fermionic systems can exhibit stronger entanglement than bosonic systems, exhibiting a negativity that decays only algebraically. The Schwinger model is described by fermionic excitations at short distances, while its asymptotic spectrum is the one of a bosonic theory. We show that the two-mode negativity detects this confining, fermion-to-boson transition, shifting from an algebraic decay to an exponential decay at distances of the order of the de Broglie wavelength of the first excited state. We derive analytical expressions in the massless Schwinger model and confront them with tensor network simulations. We also perform tensor network simulations in the massive model, which is not solvable analytically, and close to the Ising quantum critical point of the Schwinger model, where we show that the negativity behaves as its bosonic counterpart. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Tail Lobe Revisited: Magnetic Field Modeling Based on Plasma Data

Plasma data from the ISEE-1 and -2 spacecraft during 1977-1980 have been used to determine the distribution of data points in the magnetotail in the range of distances -20 < XGSM < --15, i.e. which of the records that were located in the current sheet, in the tail lobe, in the magnetosheath and in the boundary layers respectively. The ISEE-1 and -2 magnetic field data for the records in the tail lobe were then used to model the tail lobe magnetic field dependence on the solar wind dynamic pressure, on the Interplanetary Magnetic Field (IMF) and on the Dst index. The tail lobe magnetic field was assumed to be dependent on the square root of the dynamic pressure based on the balance between the total magnetic pressure in the tail lobes and the dynamic pressure of the solar wind. The IMF dependent terms, added to the pressure term, were sought in many different forms while the Dst dependence of the tail lobe magnetic field was assumed to be linear. The field shows a strong dependence on the square root of the dynamic pressure and the different IMF dependent terms all constitute a significant contribution to the total field. However, the dependence on the Dst index turned out to be very weak at those down-tail distances. The results of this study are intended to be used for parameterizing future versions of the data-based models of the global magnetospheric magnetic field.

Karlsson, S. B. P.↗