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 37 records · Page 2

Monte Carlo Simulation of the Rapid Crystallization of Bismuth-Doped Silicon

In this Letter we report Ising model simulations of the growth of alloys which predict quite different behavior near and far from equilibrium. Our simulations reproduce the phenomenon which has been termed 'solute trapping,' where concentrations of solute, which are far in excess of the equilibrium concentrations, are observed in the crystal after rapid crystallization. This phenomenon plays an important role in many processes which involve first order phase changes which take place under conditions far from equilibrium. The underlying physical basis for it has not been understood, but these Monte Carlo simulations provide a powerful means for investigating it.

Jackson, Kenneth A.↗

A Multilevel Approach For SolvingLarge-Scale QUBO Problems With Noisy Hybrid Quantum Approximate Optimization

Quantum approximate optimization is one ofthe promising candidates for useful quantum computation,particularly in the context of finding approximate solutionsto Quadratic Unconstrained Binary Optimization (QUBO)problems. However, the existing quantum processing units(QPUs) are of relatively small size, and canonical mappingsof QUBO via the Ising model require one qubit per vari-able, rendering direct large-scale optimization infeasible.In classical optimization, a general strategy for addressingmany large-scale problems is via multilevel/multigrid meth-ods, where the large target problem is iteratively coarsenedand the global solution is constructed from multiple small-scale optimization runs. In this work, we experimentallytest how existing QPUs perform when used as a sub-solverwithin such a multilevel strategy. To this aim, we com-bine and extend (via additional classical processing steps)the recently proposed Noise-Directed Adaptive Remapping(NDAR) and Quantum Relax&Round (QRR) algorithms.We first demonstrate the effectiveness of our heuristicextensions on Rigetti’s superconducting transmon deviceAnkaa-2. We find approximate solutions to10instances offully connected82-qubit Sherrington-Kirkpatrick graphswith random integer-valued coefficients obtaining normal-ized approximation ratios (ARs) in the range∼0.98−1.0,and the same class with real-valued coefficients (ARs∼0.94−1.0). Then, we implement the extended NDAR andQRR algorithms as subsolvers in the multilevel algorithmfor6large-scale graphs with at most∼27,000variables.In practice, the QPU (with classical post-processing steps)is used to find approximate solutions to dozens of at most82-qubit problems, which are iteratively used to constructthe global solution. We observe that quantum optimizationresults are competitive in terms of the quality of solutionswhen compared to classical heuristics used as subsolverswithin the multilevel approach.Reproducibility: source code and data are available at[TBA upon acceptance]

quantum computing↗

Three-dimensional analytical model for the spatial variation of the foreshock electron distribution function - Systematics and comparisons with ISEE observations

A model which is consistent with the solar wind and shock surface boundary conditions for the foreshock electron distribution in the absence of wave-particle effects is formulated for an arbitrary location behind the magnetic tangent to the earth's bow shock. Variations of the gyrophase-averaged velocity distribution are compared and contrasted with in situ ISEE observations. It is found that magnetic mirroring of solar wind electrons is the most important process by which nonmonotonic reduced electron distributions in the foreshock are produced. Leakage of particles from the magnetosheath is shown to be relatively unimportant in determining reduced distributions that are nonmonotonic. The two-dimensional distribution function off the magnetic field direction is the crucial contribution in producing reduced distributions which have beams. The time scale for modification of the electron velocity distribution in velocity space can be significantly influenced by steady state spatial gradients in the background imposed by the curved shock geometry.

Fitzenreiter, R. J.↗

Detailed observations of the plasma sheet during a substorm on April 24, 1979

Plasma, magnetic field, and energetic particle data obtained by ISEE 1 and 2 satellites for the April 24, 1979 substorm are studied in relation to the neutral line model and the boundary layer model. The ISEE 1 and 2 instruments and experiments utilized to collect the data are discussed. The major reconfiguration of the tail plasma and magnetic field plasma region, and the plasma ion flows observed support the neutral line model for interpreting substorms. The plasma ion distribution function and plasma flow are examined. The region of lobe-plasma sheet interface referred to as the separatrix layer is identified. The differences in times of plasma sheet dropout and recovery, and absence or presence of flux anisotropies are investigated. Energetic particle measurements are analyzed to study the relationship between energetic ions and plasma ions, and the velocity distributions. The data support the application of the neutral line model to the evaluation of substorms; however, the data are inconsistent with the boundary layer dynamics model.

Hones, E. W., Jr.↗

Evolution and folding of repeat proteins

Repeat proteins are made with tandem copies of similar amino acid stretches that fold into elongated architectures. These proteins constitute excellent model systems to investigate how evolution relates to structure, folding, and function. Here, we propose a scheme to map evolutionary information at the sequence level to a coarse-grained model for repeat-protein folding and use it to investigate the folding of thousands of repeat proteins. We model the energetics by a combination of an inverse Potts-model scheme with an explicit mechanistic model of duplications and deletions of repeats to calculate the evolutionary parameters of the system at the single-residue level. These parameters are used to inform an Ising-like model that allows for the generation of folding curves, apparent domain emergence, and occupation of intermediate states that are highly compatible with experimental data in specific case studies. We analyzed the folding of thousands of natural Ankyrin repeat proteins and found that a multiplicity of folding mechanisms are possible. Fully cooperative all-or-none transitions are obtained for arrays with enough sequence-similar elements and strong interactions between them, while noncooperative element-by-element intermittent folding arose if the elements are dissimilar and the interactions between them are energetically weak. Additionally, we characterized nucleation-propagation and multidomain folding mechanisms. We show that the global stability and cooperativity of the repeating arrays can be predicted from simple sequence scores.

Ezequiel A. Galpern↗

Droplet model for autocorrelation functions in an Ising ferromagnet

The autocorrelation function of Ising spins in an ordered phase is studied via a droplet model. Only noninteracting spherical droplets are considered. The Langevin equation which describes fluctuations in the radius of a single droplet is studied in detail. A general description of the transformation to a Fokker-Planck equations and the ways in which a spectral analysis of that equation can be used to compute the autocorrelation function is given. It is shown that the eigenvalues of the Fokker-Planck operator form (1) a continuous spectrum of relaxation rates starting from zero for d = 2, (2) a continuous spectrum with a finite gap for d = 3, and (3) a discrete spectrum for d greater than 4, where d is the spatial dimensionality. Detailed solutions for various cases are presented.

Tang, Chao↗

Magnetopause modeling - Flux transfer events and magnetosheath quasi-trapped distributions

Three-dimensional distribution functions for energetic ions are studied numerically in the magnetosphere, through the magnetopause, and in the magnetosheath using a simple one-dimensional quasi-static model and ISEE 1 magnetopause crossing data for November 10, 1977. Quasi-trapped populations in the magnetosheath observed near flux transfer events (FTEs) are investigated, and it is shown that the population in the sheath appears to sandwich the FTE distributions. These quasi-trapped distributions are due to slow, large pitch angle, outward moving particles left behind by the outward rush of the ions more field-aligned at the time the flux was opened. It is found that sheath convective flows can map along the connected flux tube without drastically changing the distribution function, and results suggest that localized tangential fields above the upper limit may exist.

Speiser, T. W.↗

Accelerating Continuous Variable Coherent Ising Machines Via Momentum

The Coherent Ising Machine (CIM) is a non-conventional architecture that takes inspiration from physical annealing processes to solve Ising problems heuristically. Its dynamics are naturally continuous and described by a set of ordinary differential equations that have been proven to be useful for the optimization of continuous variables nonconvex quadratic optimization problems. The dynamics of such Continuous Variable CIMs (CV-CIM) encourage optimization via optical pulses whose amplitudes are determined by the negative gradient of the objective; however, standard gradient descent is known to be trapped by local minima and hampered by poor problem conditioning. In this work, we propose to modify the CV-CIM dynamics using more sophisticated pulse injections based on tried-and-true optimization techniques such as momentum and Adam. Through numerical experiments, we show that the momentum and Adam updates can significantly speed up the CV-CIM’s convergence and improve sample diversity over the original CV-CIM dynamics. We also find that the Adam-CV-CIM’s performance is more stable as a function of feedback strength, especially on poorly conditioned instances, resulting in an algorithm that is more robust, reliable, and easily tunable. More broadly, we identify the CIM dynamical framework as a fertile opportunity for exploring the intersection of classical optimization and modern analog computing.

Ising Model↗

A model of the near-earth plasma environment and application to the ISEE-A and -B orbit

A model of the near-earth environment to obtain a best estimate of the average flux of protons and electrons in the energy range from 0.1 to 100 keV for the International Sun-Earth Explorer (ISEE)-A and -B spacecraft. The possible radiation damage to the thermal coating on these spinning spacecraft is also studied. Applications of the model to other high-altitude satellites can be obtained with the appropriate orbit averaging. This study is the first attempt to synthesize an overall quantitative environment of low-energy particles for high altitude spacecraft, using data from in situ measurements.

Chan, K. W.↗

Cosmic ray composition investigations using ICE/ISEE-3

The analysis of data from the high energy cosmic experiment on ISEE-3 and associated modeling and interpretation activities are discussed. The ISEE-3 payload included two instruments capable of measuring the composition of heavy cosmic rays. The designs of these two instruments incorporated innovations which made it possible, for the first time, to measure isotopic as well as the chemical composition for a wide range of elements. As the result of the demonstrations by these two instruments of the capability to resolve individual cosmic ray isotopes, a new generation of detectors was developed using very similar designs, but having improved reliability and increased sensitive area. The composition measurements which were obtained from the ISEE-3 experiment are summarized.

Wiedenbeck, Mark E.↗

Energetic ion beams at the edge of the plasma sheet - ISEE 1 observations plus a simple explanatory model

It has been found that energetic ions with energies in the range from 10 keV to 2 MeV are an important part of the geomagnetic tail particle population. A description is presented of energetic ion, three-dimensional distribution observations from ISEE 1 during 0400-2000 UT, May 3, 1978, a time when several encounters with the geomagnetic tail plasma sheet occurred. Pronounced energetic ion streaming is observed at each entry to and exit from the plasma sheet. Azimuthal gradients observed during these times indicate that this streaming occurs within two gyroradii of the plasma sheet edge. Streaming characteristics show both narrow and broad angular flows at arbitrary angles to B. These flows are observed simultaneously and individually. The time, energy, and pitch angle evolution of the ion fluxes show that ISEE 1 is encountering a dynamic and not a static equilibrium situation at the edge of the plasma sheet.

Williams, D. J.↗

Quantum-Assisted Variational Segmentation for Image-to-Image Wildfire Detection Using Satellite Data

The quantum computing community has been searching for suitable applications to demonstrate the potential of near-term quantum devices. Quantum machine learning is a potential candidate, particularly using models that cannot be efficiently simulated with classical computers [1, 2]. This work focuses on a transition phase of quantum computers where the quantum machine learning model is still simulable classically but projected not to be simulable as the size of the model grows. Ultimately quantum computers may have advantages for high-dimensional real-world problems. Due to the limited number of qubits in current noisy intermediate-scale quantum (NISQ) devices, the direct application of quantum computers in high dimensional data is not feasible. To remedy this problem, an encoder-decoder architecture can be utilized. The encoder model would transform the high-dimensional data into a compact representation, to a level that small quantum computers can be used today (or in the near future), and the decoder would take the quantum processed outputs back to the high-dimensional space. Addressing the two challenges of quantum machine learning, this work investigates a hybrid supervised generative model with a quantum Ising Born machine embedded as the latent distribution. The model contains four main parts (Figure 1.a.): (1) a U-NET architecture responsible for learning segmentation flow, (2) a Prior network responsible for learning an encoded latent distribution of the input data, (3) a Born machine which represents the latent distribution, and (4) a Posterior network in charge of learning the joint encoded latent distribution of inputs and target data. The initial model, proposed by [3], is optimized by (1) maximizing the overlap of the prior and posterior latent distributions, and (2) minimizing the segmentation loss. The proposed model is designed to be investigated in a simulation environment applied to the real-world application of wildfire segmentation. Specifically, the model is designed to solve the patchy wildfire segmentations of Moderate Resolution Imaging Spectroradiometer (MODIS) by taking the MODIS observations and using Visible Infrared Imaging Radiometer Suite’s (VIIRS) consistent wildfire product as the target. The model solves patchy wildfire segmentations and provides insight into the epistemic errors sourced from model variation. The model utilizes the Born machine as a QUBO solver to represent the latent space as a Bernoulli distribution. The proposed configuration allows the variational segmentation model to leverage the true quantum probabilistic nature and derive a more expressive latent configuration, increasing the model performance in describing wildfire segmentations. The quantum probabilistic information of the Born machine is directly incorporated in the Kullback-Leibler divergence loss in the prior and posterior distributions, forcing the Bernoulli latent distribution to maximize the overlap of input and joint input-target distributions. The proposed model is then trained and compared with a baseline only consisting of direct Bernoulli latent distribution with no Born machine representing the latent space. The models are evaluated based on the segmentation metrics, such as precision, recall, intersect of union, with uncertainty boundaries accounting for the stochastic nature of the model. Our findings show that even in low latent-dimensional space (due to the limit in computational power of the classical quantum simulator), we are able to effectively capture the latent representation and hence the model performs better than the baseline. The findings are a projection for scaling the model into higher dimensional latent space with the Born machine surpassing the baseline performance. Figure 1. Sub-figure (a) demonstrates the architecture for the training phase. The model consists of a Prior and Posterior network that encode inputs and joint input-target data into compact representations, respectively. The Born machine represents the latent distribution, and the U-NET branch learns the segmentation patterns of the data. The stochasticity is introduced to the U-NET through its last layer to create meaningful but stochastic segmentations. Sub-figure (b) represents the inference phase where the model takes the stochastic behavior from the prior network and injects that into the U-NET. Each attempt of inference will generate different but similar segmentations from the same distribution of the wildfire event. REFERENCES [1] Coyle, B., Mills, D., Danos, V., & Kashefi, E. (2020). The Born supremacy: quantum advantage and training of an Ising Born machine. npj Quantum Information, 6(1), 1-11. [2] Liu, J. G., & Wang, L. (2018). Differentiable learning of quantum circuit born machines. Physical Review A, 98(6), 062324. [3] Kohl, S., Romera-Paredes, B., Meyer, C., De Fauw, J., Ledsam, J. R., Maier-Hein, K., ... & Ronneberger, O. (2018). A probabilistic u-net for segmentation of ambiguous images. Advances in neural information processing systems, 31.

quantum machine learning↗

The TAME Project: Towards improvement-oriented software environments

Experience from a dozen years of analyzing software engineering processes and products is summarized as a set of software engineering and measurement principles that argue for software engineering process models that integrate sound planning and analysis into the construction process. In the TAME (Tailoring A Measurement Environment) project at the University of Maryland, such an improvement-oriented software engineering process model was developed that uses the goal/question/metric paradigm to integrate the constructive and analytic aspects of software development. The model provides a mechanism for formalizing the characterization and planning tasks, controlling and improving projects based on quantitative analysis, learning in a deeper and more systematic way about the software process and product, and feeding the appropriate experience back into the current and future projects. The TAME system is an instantiation of the TAME software engineering process model as an ISEE (integrated software engineering environment). The first in a series of TAME system prototypes has been developed. An assessment of experience with this first limited prototype is presented including a reassessment of its initial architecture.

Basili, Victor R.↗

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.↗

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.↗

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.↗

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.↗