Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Ising 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 127 records · Page 7

Magnetic properties of the Ising-like rare earth pyrosilicate: D-Er 2 Si 2 O 7

Ising-like spin-1/2 magnetic materials are of interest for their ready connection to theory, particularly in the context of quantum critical behavior. In this work we report detailed studies of the magnetic properties of a member of the rare earth pyrosilicate family, D-Er 2 Si 2 O 7 , which is known to display a highly anisotropic Ising-like g-tensor and effective spin-1/2 magnetic moments. We used powder neutron diffraction, powder inelastic neutron spectroscopy (INS), and single crystal AC susceptibility to characterize its magnetic properties. Neutron diffraction enabled us to determine the magnetic structure below the known transition temperature (T N = 1.9 K) in zero field, confirming that the magnetic state is a four-sublattice antiferromagnetic structure with two non-collinear Ising axes, as was previously hypothesized. Our powder INS data revealed a gapped excitation at zero field, consistent with anisotropic (possibly Ising) exchange. An applied field of 1 T produces a mode softening, which is consistent with a field-induced second order phase transition. To assess the relevance of D-Er 2 Si 2 O 7 to the transverse field Ising model, we performed AC susceptibility measurements on a single crystal with the magnetic field oriented in the direction transverse to the Ising axes. This revealed a transition at 2.65 T at 0.1 K, a field significantly higher than the mode-softening field observed by powder INS, showing that the field-induced phase transitions are highly field-direction dependent as expected. Furthermore, these measurements suggest that D-Er 2 Si 2 O 7 may be a candidate for further exploration related to the transverse field Ising model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Slowly decaying zero mode in a weakly nonintegrable boundary impurity model

The transverse field Ising model (TFIM) on the half-infinite chain possesses an edge zero mode. This work considers an impurity model: TFIM perturbed by a boundary integrability-breaking interaction. For sufficiently large transverse field, but in the ordered phase of the TFIM, the zero mode is observed to decay. The decay is qualitatively different from zero modes where the integrability-breaking interactions are nonzero all along the chain. It is shown that for the impurity model, the zero mode decays by relaxing to a nonlocal quasiconserved operator, the latter being exactly conserved when the opposite edge of the chain has no noncommuting perturbations so as to ensure perfect degeneracy of the spectrum. In the thermodynamic limit, the quasiconserved operator vanishes, and a regime is identified where the decay of the zero mode obeys Fermi's golden rule. Furthermore, a toy model for the decay is constructed in Krylov space and it is highlighted how Fermi's golden rule may be recovered from this toy model.

1-dimensional systems↗

Kinetic Monte Carlo Simulations of Rod Eutectics and the Surface Roughening Transition in Binary Alloys

In this paper we report on two related topics: Kinetic Monte Carlo simulations of the steady state growth of rod eutectics from the melt, and a study of the surface roughness of binary alloys. We have implemented a three dimensional kinetic Monte Carlo (kMC) simulation with diffusion by pair exchange only in the liquid phase. Entropies of fusion are first chosen to fit the surface roughness of the pure materials, and the bond energies are derived from the equilibrium phase diagram, by treating the solid and liquid as regular and ideal solutions respectively. A simple cubic lattice oriented in the {100} direction is used. Growth of the rods is initiated from columns of pure B material embedded in an A matrix, arranged in a close packed array with semi-periodic boundary conditions. The simulation cells typically have dimensions of 50 by 87 by 200 unit cells. Steady state growth is compliant with the Jackson-Hunt model. In the kMC simulations, using the spin-one Ising model, growth of each phase is faceted or nonfaceted phases depending on the entropy of fusion. There have been many studies of the surface roughening transition in single component systems, but none for binary alloy systems. The location of the surface roughening transition for the phases of a eutectic alloy determines whether the eutectic morphology will be regular or irregular. We have conducted a study of surface roughness on the spin-one Ising Model with diffusion using kMC. The surface roughness was found to scale with the melting temperature of the alloy as given by the liquidus line on the equilibrium phase diagram. The density of missing lateral bonds at the surface was used as a measure of surface roughness.

Bentz, Daniel N.↗

The Pioneer 11 1976 solar conjunction: A unique opportunity to explore the heliographic latitudinal variations of the solar corona

The 1976 Pioneer II Solar Conjunction provided the opportunity to accumulate a substantial quantity of doppler noise data over a dynamic range of signal closest approach point heliographic latitudes. The observed doppler noise data were fit to the doppler noise model ISED, and the deviations of the observed doppler noise data from the model were used to construct a (multiplicative) function to describe the effect of heliographic latitude. This expression was then incorporated into the ISED model to produce a new doppler noise model-ISEDB.

Berman, A. L.↗

Z 2 topological order and first-order quantum phase transitions in systems with combinatorial gauge symmetry

We study a generalization of the two-dimensional transverse-field Ising model, combining both ferromagnetic and antiferromagnetic two-body interactions, that hosts exact global and local Z 2 gauge symmetries. Using exact diagonalization and stochastic series expansion quantum Monte Carlo methods, we confirm the existence of the topological phase in line with previous theoretical predictions. Our simulation results show that the transition between the confined topological phase and the deconfined paramagnetic phase is of first order, in contrast to the conventional Z 2 lattice gauge model in which the transition maps onto that of the standard Ising model and is continuous. We further generalize the model by replacing the transverse field on the gauge spins with a ferromagnetic X X interaction while keeping the local gauge symmetry intact. We find that the Z 2 topological phase remains stable, while the paramagnetic phase is replaced by a ferromagnetic phase. The topological-ferromagnetic quantum phase transition is also of first order. For both models, we discuss the low-energy spinon and vison excitations of the topological phase and their avoided level crossings associated with the first-order quantum phase transitions.

2-dimensional systems↗

Effective and asymptotic criticality of structurally disordered magnets

Changes in magnetic critical behaviour of quenched structurally-disordered magnets are usually exemplified in experiments and in MC simulations by diluted systems consisting of magnetic and non-magnetic components. In our study we aim to show that similar effects can be observed not only for diluted magnets with non-magnetic impurities but may be implemented, e.g., by the presence of two (and more) chemically different magnetic components as well. Therefore we consider a model of the structurally-disordered quenched magnet where all lattice sites are occupied by Ising-like spins of different lengths L. In such a random spin length Ising model, the length L of each spin is a random variable governed by the distribution function p (L). We demonstrate that this model belongs to the universality class of the site-diluted Ising model. This proves that both models are described by the same values of asymptotic critical exponents. However, their effective critical behaviour differs. As a case study, we consider a quenched mixture of two different magnets with values of elementary magnetic moments L 1 = 1 and L 2 = s, and of concentration c and 1 - c, correspondingly. We apply field-theoretical renormalization group approach to analyse the renormalization group flow for different initial conditions, triggered by s and c, and to calculate effective critical exponents further away from the fixed points of the renormalization group transformation. We show how the effective exponents are governed by difference in properties of the magnetic components.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Susceptibility Critical Exponent of a 1D Ising Ring-Type Ferromagnetic

Using the Monte Carlo method, critical behavior of the one-dimensional ferromagnetic Ising model has been investigated with allowance for the interaction of the second and third neighbors and four-particle interaction. The obtained results on the critical temperature were compared with the critical temperature of the quasi-one-dimensional Ising magnetic [(CH{sub 3}){sub 3}NH] · FeCl{sub 3} · 2H{sub 2}O and with the magnitude of the exchange interaction J/k{sub B} = 17.4 K. Within the scope of the finite-dimensional scaling theory, the critical susceptibility exponent has been calculated. It has been shown that values of the susceptibility exponent for the one-dimensional Ising model with periodic boundary conditions are considerably less than the known values of the exponents for three-dimensional systems. The critical susceptibility exponent strongly depends on energy parameters; namely, it decreases with an increase in them.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Floquet operator dynamics and orthogonal polynomials on the unit circle

Operator spreading under stroboscopic time evolution under a unitary is studied. An operator Krylov space is constructed and mapped to orthogonal polynomials on a unit circle (OPUC), as well as to the Krylov space of the edge operator of the Floquet transverse field Ising model with inhomogeneous couplings (ITFIM). The Verblunsky coefficients in the OPUC representation are related to the Krylov angles parametrizing the ITFIM. The relations between the OPUC and spectral functions are summarized and several applications are presented. These include derivation of analytic expressions for the OPUC under persistent m-periodic dynamics, and the numerical construction of the OPUC for autocorrelations of the homogeneous Floquet-Ising model as well as the Z 3 clock model. The numerically obtained Krylov angles of the Z 3 clock model with long-lived period tripled autocorrelations show a spatial periodicity of six, and this observation is used to develop an analytically solvable model for the ITFIM that mimics this behavior.

Floquet systems↗

CFT and lattice correlators near an RG domain wall between minimal models

Conformal interfaces separating two conformal field theories (CFTs) provide maps between different CFTs, and naturally exist in nature as domain walls between different phases. One particularly interesting construction of a conformal interface is the renormalization group (RG) domain wall between CFTs. For a given Virasoro minimal model \mathcal{M}_{k+3,k+2} ℳ k + 3 , k + 2 , an RG domain wall can be generated by a specific deformation which triggers an RG flow towards its adjacent Virasoro minimal model \mathcal{M}_{k+2,k+1} ℳ k + 2 , k + 1 with the deformation turned on over part of the space. An algebraic construction of this domain wall was proposed by Gaiotto in [J. High Energy Phys. 12, 103 (2012)]. In this paper, we will provide a study of this RG domain wall for the minimal case k=2 k = 2 , which can be thought of as a nonperturbative check of the construction. In this case the wall is separating the Tricritical Ising Model (TIM) CFT and the Ising Model (IM) CFT. We will check the analytical results of correlation functions from the RG brane construction with the numerical density matrix renormalization group (DMRG) calculation using a lattice model proposed in [arXiv.1206.1332, Science 344, 280 (2014)], and find a perfect agreement. We comment on possible experimental realizations of this RG domain wall.

Cogburn, Cameron V.↗

Misanthropic entropy and renormalization as a communication channel

A central physical question is the extent to which infrared (IR) observations are sufficient to reconstruct a candidate ultraviolet (UV) completion. We recast this question as a problem of communication, with messages encoded in field configurations of the UV being transmitted to the IR degrees of freedom via a noisy channel specified by renormalization group (RG) flow, with noise generated by coarse graining/decimation. We present an explicit formulation of these considerations in terms of lattice field theory, where we show that the “misanthropic entropy” — the mutual information obtained from decimating neighbors — encodes the extent to which information is lost in marginalizing over/tracing out UV degrees of freedom. Our considerations apply both to statistical field theories as well as density matrix renormalization of quantum systems, where in the quantum case, the statistical field theory analysis amounts to a leading-order approximation. As a concrete example, we focus on the case of the 2D Ising model, where we show that the misanthropic entropy detects the onset of the phase transition at the Ising model critical point.

Physics↗

Diagnosis of information scrambling from Hamiltonian evolution under decoherence

We apply a quantum teleportation protocol based on the Hayden-Preskill thought experiment to quantify how scrambling a given quantum evolution is. It has an advantage over the direct measurement of out-of time ordered correlators when used to diagnose the information scrambling in the presence of decoherence effects stemming from a noisy quantum device. We demonstrate the protocol by applying it to two physical systems: Ising spin chain and SU(2) lattice Yang-Mills theory. To this end, we numerically simulate the time evolution of the two theories in the Hamiltonian formalism. The lattice Yang-Mills theory is implemented with a suitable truncation of Hilbert space on the basis of the Kogut-Susskind formalism. On a two-leg ladder geometry and with the lowest nontrivial spin representations, it can be mapped to a spin chain, which we call Yang-Mills-Ising model and is also directly applicable to future digital quantum simulations. Here, we find that the Yang-Mills-Ising model shows the signal of information scrambling at late times.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Permutation matrix representation quantum Monte Carlo

We present a quantum Monte Carlo algorithm for the simulation of general quantum and classical many-body models within a single unifying framework. The algorithm builds on a power series expansion of the quantum partition function in its off-diagonal terms and is both parameter-free and Trotter error-free. In our approach, the quantum dimension consists of products of elements of a permutation group. As such, it allows for the study of a very wide variety of models on an equal footing. To demonstrate the utility of our technique, we use it to clarify the emergence of the sign problem in the simulations of non-stoquastic physical models. We showcase the flexibility of our algorithm and the advantages it offers over existing state-of-the-art by simulating transverse- field Ising model Hamiltonians and comparing the performance of our technique against that of the stochastic series expansion algorithm. Furthermore, we also study a transverse-field Ising model augmented with randomly chosen two-body transverse-field interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Vanishing k-space fidelity and phase diagram’s bulk–edge–bulk correspondence

The fidelity between two infinitesimally close states and the fidelity susceptibility of a system are known to detect quantum phase transitions. Here we show that the k-space fidelity between two states far from each other and taken deep inside (bulk) of two phases, generically vanishes at the k-points where there are gapless points in the energy spectrum that give origin to the lines (edges) separating the phases in the phase diagram. We consider a general case of two-band models and present a sufficient condition for the existence of gapless points, given there are pairs of parameter points for which the fidelity between the corresponding states is zero. By presenting an explicit counter-example, we show that the sufficient condition is not necessary. Further, we show that, unless the set of parameter points is suitably constrained, the existence of gapless points generically implies accompanied pairs of parameter points with vanishing fidelity. Also, we show the connection between the vanishing fidelity and gapless points on a number of concrete examples (topological triplet superconductor, topological insulator, 1d Kitaev model of spinless fermions, BCS superconductor, Ising model in a transverse field, graphene and Haldane Chern insulator), as well as for the more general case of Dirac-like Hamiltonians. We also briefly discuss the relation between the vanishing fidelity and gapless points at finite temperatures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Bethe Approximation for a Two-Dimensional Spin-Pseudospin System

We have treated a two-dimensional spin-pseudospin model, which generalizes a diluted antiferromagnetic Ising model with charged nonmagnetic impurities in the case of two types of charges. The analytical results in the Bethe approximation are compared with the results of numerical simulation using the classical Monte Carlo method for various parameters.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Doppler noise considered as a function of the signal path integration of electron density

The hypothesis is advanced that observed Doppler noise during solar conjunctions is proportional to total columnar electron content along the signal path. This assumption leads directly to a geometrical model (ISED) for observed Doppler noise which is shown to be in very good agreement with Doppler noise data accumulated during the 1975 Pioneer 10, Pioneer 11 and Helios 1 solar conjunctions. An augmented model is constructed which quantitatively indicates correlation between earth observed sunspot activity and systematic, cyclical deviations from the ISED model.

Berman, A. L.↗

On the degeneracy of spin ice graphs, and its estimate via the Bethe permanent

The concept of spin ice can be extended to a general graph. We study the degeneracy of spin ice graph on arbitrary interaction structures via graph theory. We map spin ice graphs to the Ising model on a graph and clarify whether the inverse mapping is possible via a modified Krausz construction. From the gauge freedom of frustrated Ising systems, we derive exact, general results about frustration and degeneracy. We demonstrate for the first time that every spin ice graph, with the exception of the one-dimensional Ising model, is degenerate. We then study how degeneracy scales in size, using the mapping between Eulerian trails and spin ice manifolds, and a permanental identity for the number of Eulerian orientations. Furthermore, we show that the Bethe permanent technique provides both an estimate and a lower bound to the frustration of spin ices on arbitrary graphs of even degree. While such a technique can also be used to obtain an upper bound, we find that in all finite degree examples we studied, another upper bound based on Schrijver inequality is tighter.

97 MATHEMATICS AND COMPUTING↗

Discovering phases, phase transitions, and crossovers through unsupervised machine learning: A critical examination

Here, we apply unsupervised machine learning techniques, mainly principal component analysis (PCA), to compare and contrast the phase behavior and phase transitions in several classical spin models - the square and triangular-lattice Ising models, the Blume-Capel model, a highly degenerate biquadratic-exchange spin-one Ising (BSI) model, and the 2D XY model, and examine critically what machine learning is teaching us. We find that quantified principal components from PCA not only allow exploration of different phases and symmetry-breaking, but can distinguish phase transition types and locate critical points. We show that the corresponding weight vectors have a clear physical interpretation, which is particularly interesting in the frustrated models such as the triangular antiferromagnet, where they can point to incipient orders. Unlike the other well-studied models, the properties of the BSI model are less well known. Using both PCA and conventional Monte Carlo analysis, we demonstrate that the BSI model shows an absence of phase transition and macroscopic ground-state degeneracy. The failure to capture the 'charge' correlations (vorticity) in the BSI model (XY model) from raw spin configurations points to some of the limitations of PCA. Finally, we employ a nonlinear unsupervised machine learning procedure, the 'antoencoder method', and demonstrate that it too can be trained to capture phase transitions and critical points.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Using the D-Wave 2X Quantum Computer to Explore the Formation of Global Terrorist Networks

Social networks with signed edges (+/-) play an important role in an area of social network theory called structural balance. In these networks, edges represent relationships that are labeled as either friendly (+) or hostile (-). A signed social network is balanced only if all cycles of three or more nodes in the graph have an odd number of hostile edges. A fundamental property of a balanced network is that it can be cleanly divided into 2 factions, where all relationships within each faction are friendly, and all relationships between members of different factions are hostile. The more unbalanced a network is, the more edges will fail to adhere to this rule, making factions more ambiguous. Social theory suggests unbalanced networks should be unstable, a finding that has been supported by research on gangs, which shows that unbalanced relationships are associated with greater violence, possibly due to this increased ambiguity about factional allegiances (Nakamura et al). One way to estimate the imbalance in a network, if only edge relationships are known, is to assign nodes to factions that minimize the number of violations of the edge rule described above. This problem is known to be computationally NP-hard. However, Facchetti et al. have pointed out that it is equivalent to an Ising model with a Hamiltonian that effectively counts the number of edge rule violations. Therefore, finding the assignment of factions that minimizes energy of the equivalent Ising system yields an estimate of the imbalance in the network. Based on the Ising model equivalence of the signed-social network balance problem, we have used the D-Wave 2X quantum annealing computer to explore some aspects of signed social networks. Because connectivity in the D-Wave computer is limited to its particular native topology, arbitrary networks cannot be represented directly. Rather, they must be “embedded” using a technique in which multiple qubits are chained together with special weights to simulate a collection of nodes with the required connectivity. This limits the size of a fully connected network in the D-Wave to about 50 simulated nodes, using all of the approximately 1150 qubits in the machine. In order to keep within this limitation, while exploring a problem of potential social relevance, we constructed time series of historical network snapshots from Stanford’s Mapping Militants Project, where nodes represent militant organizations, and edges represent either alliances or rivalries between organizations. We constructed two series from different theaters – Iraq and Syria – spanning timelines from about 2000 to 2016, each with networks whose maximum size was in the 20-30 node range. Computationally, our experience suggests D-Wave technology is promising, providing fast, nearly constant scaling of computational effort in the main part of the calculation that relies on the quantum annealing cycle. However, the cost of embedding an arbitrary network of interest in the D-Wave native topology scales poorly. If the embedding cost can be amortized relative to the annealing cycle, it may be possible to gain a substantial advantage over classical computing methods, provided a large enough network can be accommodated by partitioning into subnetworks or some similar strategy. In terms of our application to networks of militant organizations, we found a rise in network imbalance in the Syrian theater that appears to correspond roughly with the entrance of the Islamic State into a milieu already populated with other groups, a phenomenon we plan to explore in more detail. In these very preliminary results, we also noticed that during at least one period where both the size and imbalance of the network increased substantially, the imbalance per edge seemed to remain fairly steady. This may suggest some adaptive behavior among the participating factions, which may also warrant further exploration.

97 MATHEMATICS AND COMPUTING↗