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 451 records · Page 25

Deconfined classical criticality in the anisotropic quantum spin- 1 2 XY model on the square lattice

The anisotropic quantum spin- 1 2 XY model on a linear chain was solved by Lieb, Schultz, and Mattis [] and shown to display a continuous quantum phase transition at the O(2) symmetric point separating two gapped phases with competing Ising long-range order. For the square lattice, the following is known. The two competing Ising ordered phases extend to finite temperatures, up to a boundary where a transition to the paramagnetic phase occurs, and meet at the O(2) symmetric critical line along the temperature axis that ends at a tricritical point at the Berezinskii-Kosterlitz-Thouless transition temperature where the two competing phases meet the paramagnetic phase. We show that the first-order zero-temperature (quantum) phase transition that separates the competing phases as a function of the anisotropy parameter is smoothed by thermal fluctuations into deconfined classical criticality. Published by the American Physical Society 2025

Mudry, Christopher (ORCID:0000000340746758)↗

CAIDA CUI Confirmation for OMRS & LCC SOCRRATES Emulator and iSEE Upgrade, Operations, & Analysis

This internship has focused on providing solutions for the Customer Avionics Interface Development and Analysis (CAIDA) subsystem. The main emulator that has been used during this internship is the Software-Only CEV (Crew Exploration Vehicle) Risk Reduction Analysis and Test Engineering Simulator (SOCRRATES). This emulator uses advanced math and physics methods to simulate specific points in a mission, such as ascent, entry, and orbit. Exploration Ground Systems (EGS) must ensure that all KSC based ground systems can be properly integrated with the flight vehicle software, therefore, the Modeling and Simulation Branch (NE-XM) of the KSC (Kennedy Space Center) Engineering Directorate supports a virtual environment that simulates the interface between ground systems and the flight vehicle. A primary component in assuring that is knowing whether the emulators have incorporated the correct CUIs (Compact Unique Identifiers) into their system, as well as understanding both the static and dynamic responses of the individual CUIs. Also, an effort for verification of Operational Maintenance Requirements Specifications (OMRS) and Launch Commit Criteria (LCC) that are supported by SOCRRATES for the Orion Crew Module were part of the project tasks this semester. This internship also focused on comparing two different simulation Commercial-off-the-shelf (COTS) products and to determine whether or not a COTS package was a viable replacement for the current software that the iSEE (Immersive Simulations and Engineering Environment) lab uses. Upon research and testing, I found that this software was not feasible for the lab. It could not easily load CAD (Computer-Aided Design) or CREO models, give live feedback while the user is in the environment, and was not compatible with a virtual reality headset, all of which are necessary for the lab.

Bundy, Caleb↗

Use of the thin sheath approximation for obtaining ion temperatures from the ISEE 1 limited aperture RPA

A procedure for analyzing low-energy (less than approximately 100 eV) ion data from the plasma composition experiment on ISEE 1 is set forth. The method is based on a derived analytic expression for particle flux to a limited aperture retarding potential analyzer (RPA) in the thin sheath approximation, which makes allowance for some effects of a charged spacecraft on plasma particle trajectories. Calculations using simulated data are employed in testing the efficacy and accuracy of the technique. On the basis of an analysis of these calculation results and the mathematical model, the method is seen as being able to provide accurate ion temperatures from all good plasmaspheric RPA data. It is noted that corresponding densities and spacecraft potentials should be accurate when spacecraft potentials are negative but that they are subject to error for positive spacecraft potentials, particularly when ion Mach numbers are much less than 1. An analysis of data from a representative ISEE 1 pass produces a plasmasphere temperature profile that is consistent in overall structure with previous observations.

Comfort, R. H.↗

Satellite observations of the spatial extent and structure of Pc 3, 4, 5 pulsations near the magnetospheric equator

Simultaneous observations of Pc 3, 4, 5 pulsations by five satellites in the pre-noon local time sector at and near synchronous orbit are examined. The periods of these simultaneous pulsations are not the same at the different observation points. This difference is attributed to site dependent resonant conditions. The spatial properties of the temporal phenomenon are demonstrated with observations by ISEE-1 and -2 as they pass through oscillations in a spatially limited region. Fundamental and second harmonic standing Alfven waves are observed simultaneously on the same field line. The periods are consistent with model predictions when the measured plasma composition, which by mass consists mainly of singly ionized oxygen, is taken into account.

Singer, H. J.↗

The observation of a coronal transient directed at earth

The paper reports the observation of a large coronal transient that can only be interpreted as a three-dimensional structure. Its form is one which has not been observed before: a gradually expanding, sun-centered disk of excess brightness, whose projected radius increased from 4 to 8 solar radii during 0832-0958 UT on November 27, 1979. This earth-directed transient was the source of an interplanetary shock wave that reached ISEE 3 at 0649 UT, November 30, and earth at 0738 UT, November 30. Fitting the shock speed at ISEE 3 and the average transit speed from the sun to ISEE 3 to a power law of the form V = (V sub 0)(r exp -n), it is found that V sub 0 = 1980 km/s and n = 0.294, in good agreement with shock wave models. The shock speed predicted by the power law at 10 solar radii is 1000 km/s, in good agreement with the estimated frontal speed of the transient.

Howard, R. A.↗

Investigation of the dayside boundary region of the magnetosphere

The region near the dayside boundary of the magnetosphere is of particular interest mainly because it is the most likely site for solar wind energy transfer. Exactly where and how this transfer occurs is one of the basic problems of magnetospheric physics. The global nature of the energy transfer is addressed in this research. In order to test whether energy transfer affects the large scale plasma flow in the magnetosheath, global properties of the magnetosheath were analyzed and compared to models with no energy transfer. Also, global models of merging sites on the magnetopause were modeled with the aid of the magnetosheath models. Since energy transfer and intrinsic magnetosheath properties depend upon the orientation of the interplanetary magnetic field (IMF), it was essential to have an IMF monitor for these studies. In the course of using ISEE 3 in this role, statistical studies were performed to test the accuracy of using data from an orbit so far upstream. Global properties of the magnetosheath energetic ion population, magnetic field strength and orientation, and flow pattern, determined with ISEE data are summarized.

Crooker, N.↗

Observations of Earthward and tailward propagating flux rope plasmoids: Expanding the plasmoid model of geomagnetic substorms

A survey of Interplanetary Monitoring Platform (IMP 8) magnetometer data for plasmoid signatures during magnetospheric intervals from 1981 through 1983 found 16 plasmoids and 37 traveling compression regions as well as two earthward propagating flux ropes and 19 south-north bipolar lobe signatures. The properties of these relatively near-Earth plasmoids, traveling compression regions, and earthward propagating flux ropes and a qualitative model for their formation are presented. The plasmoids have estimated sizes, durations, magnetic field signatures, downtail velocities, and substorm associations very similar to those of the plasmoids identified in International Sun-Earth Explorer (ISEE) 3 deep-tail observations. The occurrence frequency of these near-Earth plasma sheet plasmoids is significantly smaller than that of plasmoids found in the mid- and deep tail with ISEE 3. The earthward propagating flux ropes are characterized by a south-north bipolar turning in the Geocentric Solar Magnetospheric (GSM) B(sub z) component, are localized near the noon-midnight meridional plane, and are strongly correlated with interplanetary magnetic field B(sub z) north and small isolated high latitude geomagnetic substorms. These events are also apparently very rare and/or spatially localized. We propose that these structures are 'proto-plasmoids,' i.e., plasmoids for which near-Earth magnetic reconnection stopped before all the closed plasma sheet field lines were reconnected. The proto-plasmoids are then 'trapped' inside closed magnetic field lines and propagate earthward owing to the effect of the distant X-line's earthward plasma flow. We suggest that the two different 'types' of plasmoids are due to the different energy states of the magnetosphere during periods of southward and northward interplanetary magnetic field.

Moldwin, Mark B.↗

On the new universality class in structurally disordered n -vector model with long-range interactions

We study a stability boundary of a region where nontrivial critical behavior of an n-vector model with long-range power-law decaying interactions is induced by the presence of a structural disorder (e.g., weak quenched dilution). This boundary is given by the marginal dimension of the order parameter nc dependent on space dimension, d, and a control parameter of the interaction decay, σ, below which the model belongs to the new dilution-induced universality class. Exploiting the Harris criterion and recent field theoretical renormalization group results for the pure model with long-range interactions, we get nc as a three loop ɛ = 2σ – d-expansion. We provide numerical values for nc applying series resummation methods. Our results show that not only the Ising systems (n = 1) can belong to the new disorder-induced long-range universality class at d = 2 and 3.

Physics↗

Superdiffusion resilience in Heisenberg chains with two-dimensional interactions on a quantum processor

Superdiffusive spin transport in the one-dimensional (1D) Heisenberg model is a key theoretical discovery in nonequilibrium quantum many-body physics. Although extensively studied in 1D systems, the breakdown and sustenance of superdiffusion in two-dimensional (2D) lattices with integrability-breaking terms, as found in real materials, remains an open question. To address this, we develop a toy model that extends the 1D Heisenberg model with a representative set of 2D interaction types and tunable strengths. Our model exhibits varying degrees of superdiffusion breakdown depending on the interaction type, spanning ballistic to diffusive regimes. We establish and justify a hierarchy of 2D interactions based on their resilience against superdiffusion breakdown: Heisenberg >𝑋⁢𝑋 > Ising. This precise control over the superdiffusive behavior also enables rigorous benchmarking of quantum hardware, and our simulations on IBM's Heron devices confirm the hardware's ability to accurately capture these many-body nonequilibrium phenomena. Overall, our results are relevant not only to simulating superdiffusion in real materials, such as the 1D Heisenberg compound KCuF3, which contains modest nonintegrable 2D terms, but also to extending superdiffusive behavior to larger 2D qubit lattices and other 2D materials.

Alagarsamy Manikandan, Keerthi Kumaran [ORNL]↗

Problem-tailored Simulation of Energy Transport on Noisy Quantum Computers

The transport of conserved quantities like spin and charge is fundamental to characterizing the behavior of quantum many-body systems. Numerically simulating such dynamics is generically challenging, which motivates the consideration of quantum computing strategies. However, the relatively high gate errors and limited coherence times of today's quantum computers pose their own challenge, highlighting the need to be frugal with quantum resources. In this work we report simulations on quantum hardware of infinite-temperature energy transport in the mixed-field Ising chain, a paradigmatic many-body system that can exhibit a range of transport behaviors at intermediate times. We consider a chain with L = 12 sites and find results broadly consistent with those from ideal circuit simulators over 90 Trotter steps, containing up to 990 entangling gates. To obtain these results, we use two key problem-tailored insights. First, we identify a convenient basis – the Pauli Y basis – in which to sample the infinite-temperature trace and provide theoretical and numerical justifications for its efficiency relative to, e.g., the computational basis. Second, in addition to a variety of problem-agnostic error mitigation strategies, we employ a renormalization strategy that compensates for global nonconservation of energy due to device noise. We discuss the applicability of the proposed sampling approach beyond the mixed-field Ising chain and formulate a variational method to search for a sampling basis with small sample-to-sample fluctuations for an arbitrary Hamiltonian. This opens the door to applying these techniques in more general models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The isotopic composition of cosmic ray calcium

Data from the high energy cosmic ray experiment on the international sun earth explorer 3 (ISEE-3) spacecraft have been used to study the isotopic composition of cosmic ray calcium at an energy of approx. 260 MeV/amu. The arriving calcium is found to consist of (32 + or - 6)%. A propagation model consistent with both the light and the subiron secondary element abundances was used for the interpretation of the observed calcium composition. The measured 42Ca+43Ca+44Ca abundance is consistent with the calculated secondary production, while the 40Ca abundance implies a source ratio of 40Ca/Fe = (7.0 + or - 1.7)%.

Krombel, K. E.↗

The abundances of hydrogen, helium, oxygen, and iron accelerated in large solar particle events

Energy spectra measured in 10 large flares with the University of Maryland/Max-Planck-Institut sensors on ISEE I and Goddard Space Flight Center sensors on IMP 8 allowed us to determine the average H, He, O, and Fe abundances as functions of energy in the range of about 0.3-80 MeV/nucleon. Model fits to the spectra of individual events using the predictions of a steady state stochastic acceleration model with rigidity-dependent diffusion provided a means of interpolating small portions of the energy spectra not measured with the instrumentation. Particles with larger mass-to-charge ratios were relatively less abundant at higher energies in the flare-averaged composition. The Fe/O enhancement at low SEP energies was less than the Fe/O ratios observed in He-3-rich flares. Unlike the SEP composition averaged above 5 MeV/nucleon, the average SEP abundances above 0.3 MeV/nucleon were similar to the average solar wind.

Mazur, J. E.↗

A quantitative assessment of the role of the post-shock turbulent region in the formation of Forbush decreases

We present results of a survey of the relation between Forbush decreases, magnetic clouds, and interplanetary shocks during the period August 1978 to November 1982. We have used data from the ISEE-3 study of bidirectional ions associated with magnetic structures or clouds of Marsden et al. (1987), and ground-based observations of Forbush decreases from several neutron monitors. We use the two-step model of a Forbush decrease. We assume that the first step is due to the passage of the postshock turbulent region, and that the second is due to the passage of the magnetic cloud or structure which usually follows the postshock turbulent region. To determine the effectiveness of the postshock turbulent region in causing a Forbush decrease, we have evaluated the radial diffusion coefficient of the postshock turbulent region for the eight largest events during the above period using observations of the magnetic field. We have made a quantitative assessment of the relative importance of the postshock turbulent region in the formation of the Forbush decrease, concluding that the postshock turbulent region alone is not sufficient to cause a Forbush decrease.

Sanderson, T. R.↗

Magnetospherically trapped ions as a source of magnetosheath energetic ions

It has been suggested that energetic ions observed in the magnetosheath may be due to the direct leakage of trapped magnetospheric ions. To test this hypothesis, three-dimensional ion spectra from the energetic particle experiment on ISEE 1 for a magnetopause crossing on Nov. 10, 1977 are utilized to construct three-dimensional distribution functions in the magnetosphere and in the sheath. Using the observed magnetic field, a simple one-dimensional model of the magnetopause is developed. Ions are then followed in the model, starting in the magnetosphere, through the magnetopause and ending up in the sheath. Using Liouville's Theorem a model sheath distribution function is then built up by following the magnetospheric distribution function through the model fields. The model distribution is then compared with the observed sheath distribution. For this case it is found that the main features of the observed ions in the sheath are consistent with direct leakage and with no energization or de-energization processes, and an inward-pointing normal component is required. The energetic particles mapped in this case apparently follow a flux tube which does not penetrate the magnetopause where local tangential electric fields have been reported.

Speiser, T. W.↗

Dynamics During Thrust Maneuvers of Flexible Spinning Satellites with Axial and Radial Booms

The dynamic response to operational maneuvers of spinning symmetric spacecraft with radial and axial booms was analyzed as part of the prelaunch dynamic analysis of the ISEE-3 spacecraft placed in a halo orbit around an Earth-Sun libration point, and later renamed ICE when it was directed to fly-by comet Giacobini-Zinner. The results presented use simple spacecraft models, and frequently give predictions that are good and easily obtained when the results from using a general purpose multibody dynamics program were very time consuming to obtain. Deployment of radial booms, spin-up after partial deployment, stationkeeping, and trajectory changes are analyzed. The latter two can involve both axial thrusting and pulsed radial thrusting once per revolution.

Longman, R. W.↗

Focused transport of energetic particles along magnetic field lines draped around a coronal mass ejection

Evidence is presented for focused transport of energetic particles along magnetic field lines draped around a coronal mass ejection. This evidence was obtained with the University of Maryland/Max-Planck-Institute experiment on the ISEE-3 spacecraft during the decay phase of the June 6, 1979, solar particle event. During the early portion of the decay phase of this event, interplanetary magnetic field lines were apparently draped around a coronal mass ejection, leading to a small focusing length on the western flank where ISEE 3 was located. A period of very slow decrease of particle intensity was observed, along with large sunward anisotropy in the solar wind frame, which is inconsistent with predictions of the standard Fokker-Planck equation models for diffusive transport. It was found possible to fit the observations, assuming that focused transport dominates and that the particle pitch angle scattering is isotropic.

Tan, L. C.↗

The matrix product approximation for the dynamic cavity method

Stochastic dynamics of classical degrees of freedom, defined on vertices of locally tree-like graphs, can be studied in the framework of the dynamic cavity method which is exact for tree graphs. Such models correspond for example to spin-glass systems, Boolean networks, neural networks, and other technical, biological, and social networks. The central objects in the cavity method are edge messages—conditional probabilities of two vertex variable trajectories. In this paper, we discuss a rather pedagogical derivation for the dynamic cavity method, give a detailed account of the novel matrix product edge message (MPEM) algorithm for the solution of the dynamic cavity equation as introduced, and present optimizations and extensions. Matrix product approximations of the edge messages are constructed recursively in an iteration over time. Computation costs and precision can be tuned by controlling the matrix dimensions of the MPEM in truncations. Without truncations, the dynamics is exact. Data for Glauber–Ising dynamics shows a linear growth of computation costs in time. In contrast to Monte Carlo simulations, the approach has a much better error scaling. Hence, it gives for example access to low probability events and decaying observables like temporal correlations. Here, we discuss optimized truncation schemes and an extension that allows to capture models which have a continuum time limit.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Observation of Unconventional Dynamics of Domain Walls in Uniaxial Ferroelectric Lead Germanate

Application of scanning probe microscopy techniques such as piezoresponse force microscopy (PFM) opens the possibility to re-visit the ferroelectrics previously studied by the macroscopic electrical testing methods and establish a link between their local nanoscale characteristics and integral response. It is reported that the nanoscale PFM studies and phase-field modeling of the static and dynamic behavior of the domain structure in the well-known ferroelectric material lead germanate Pb 5 Ge 3 O 11 . Here, several unusual phenomena are revealed: 1) domain formation during the paraelectric-to-ferroelectric phase transition, which exhibits an atypical cooling rate dependence; 2) unexpected electrically induced formation of the oblate domains due to the preferential domain walls motion in the directions perpendicular to the polar axis, contrary to the typical domain growth behavior observed so far; 3) absence of the bound charges at the 180 degrees head-to-head (H-H) and tail-totail (T-T) domain walls, which typically exhibit a significant charge density in other ferroelectrics due to the polarization discontinuity. This strikingly different behavior is rationalized by the phase field modeling of the dynamics of uncharged H-H and T-T domain walls. The results provide a new insight into the emergent physics of the ferroelectric domain boundaries, revealing unusual properties not exhibited by conventional Ising-type walls.

36 MATERIALS SCIENCE↗