Search tips
Search criteria 


Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
Sci Rep. 2017; 7: 44370.
Published online 2017 March 15. doi:  10.1038/srep44370
PMCID: PMC5353626

Intrinsic optimization using stochastic nanomagnets


This paper draws attention to a hardware system which can be engineered so that its intrinsic physics is described by the generalized Ising model and can encode the solution to many important NP-hard problems as its ground state. The basic constituents are stochastic nanomagnets which switch randomly between the ±1 Ising states and can be monitored continuously with standard electronics. Their mutual interactions can be short or long range, and their strengths can be reconfigured as needed to solve specific problems and to anneal the system at room temperature. The natural laws of statistical mechanics guide the network of stochastic nanomagnets at GHz speeds through the collective states with an emphasis on the low energy states that represent optimal solutions. As proof-of-concept, we present simulation results for standard NP-complete examples including a 16-city traveling salesman problem using experimentally benchmarked models for spin-transfer torque driven stochastic nanomagnets.

The use of Ising computers to solve NP-hard problems has a rich heritage in both theory1 and practice. These computers seek to solve a wide range of optimization problems by encoding the solution to the problem as the ground-state of an Ising energy expression. Many diverse systems have been proposed to solve NP-hard optimization problems such as those based on simulated annealing2, DNA3,4, quantum annealing5,6, Cellular Neural Networks7,8,9, CMOS10, trapped ions11, electromechanics12, optics13,14,15,16,17,18,19,20, and magnets21,22,23. A common objective of many of the Ising-based approaches is the identification of hardware configurations that can efficiently solve optimization problems of interest.

In this letter, we demonstrate the possibility of a hardware implementation that does not just mimic the Ising model, but embodies it as a part of its natural physics21,22,23. It uses a network of N “soft” nanomagnets operating in a stochastic manner24, each with an energy barrier Δ comparable to kBT so that they switch between the two Ising states, ±1, on time scales An external file that holds a picture, illustration, etc.
Object name is srep44370-m1.jpg where τ0 ~ 0.1–1 ns. The natural laws of statistical mechanics guide the network through the 2N collective states at GHz rates, with an emphasis on low energy states. We show how an optimization problem of interest is solved by engineering the spin-mediated magnet-magnet interactions to encode the problem solution and to simulate annealing without any change in temperature simply by continuously adjusting their overall strength. As proof-of-concept for the potential applications of this natural Ising computer, we present detailed simulation results for standard NP-complete examples, including a 16-city traveling salesman problem. This involves using experimentally benchmarked modules to simulate a suitably designed network of 225 stochastic nanomagnets and letting the hardware itself rapidly identify solutions within the 2225 possibilities. It should be possible to integrate such hardware into standard solid state circuits, which will govern the scalability of the solution.

The Ising Hamiltonian for a collection of spins, Si, which can take on one of two values, ±1,

An external file that holds a picture, illustration, etc.
Object name is srep44370-m2.jpg

was originally developed to describe ferromagnetism where the Jij are positive numbers representing an exchange interaction between neighboring spins Si and Sj, while hi represents a local magnetic field for spin Si. Classically, different spin configurations σ{Si} have a probability proportional to An external file that holds a picture, illustration, etc.
Object name is srep44370-m3.jpg, T being the temperature, and kB, the Boltzmann constant. At low temperatures, the system should be in its ground state σG, the state with the lowest energy H(σ). With hi = 0, and positive Jij, it is easy to see that the ground state is the ferromagnetic configuration σF with all spins parallel.

Much of the interest in the Ising Hamiltonian arises from the demonstration of many direct mappings of NP-complete and NP-hard problems to the model1,25,26 such that the desired solution is represented by the spin configuration σ corresponding to the ground state. However, in general this mapping may require a large number of spins, and may require the parameters Jij and hi to take on a wide range of values, both positive and negative. Finding the ground state of this artificial spin glass is the essence of Ising computing, and broadly speaking it involves abstractly representing an array of spins, their coupling, and thermal noise through software and hardware that attempts to harness the efficiencies of physical equivalence27. These representations may take the form of abstract models of the spins, the use of random number generators to produce noise, and logical or digital adders for the weighted summing. If enough layers of abstraction can be eliminated, the underlying hardware will inherently solve a given problem as part of its natural, intrinsic operation and this should be reflected in increased speed and efficiency.

Engineering Correlations Through Spin Currents

Here we describe a natural hardware for an Ising computer based on the representation of an Ising spin Sk by the magnetization m of a stochastic nanomagnet (SNM), which we believe will compare well with other alternative representations. These SNMs are in the “telegraphic” switching regime24,28 requiring the existence of a small barrier in the magnetic energy (Δ  kBT), that gives a small, but definite preference for a given axis, with two preferred states ±1. In the absence of currents, these SNMs continually switch between +1 and −1 on the order of nanoseconds, and can be physically realized by a reduction of the magnetic grain volume29 or by designing weak perpendicular magnetic anisotropy (PMA) magnets30. Figure 1 shows the response of such a monodomain PMA magnet in the presence of an external spin current in the direction of the magnet’s easy axis.

Figure 1
Response of stochastic nanomagnet to spin current.

How do we couple the SNMs to implement the Ising Hamiltonian of Equation (1)? The usual forms of coupling involve dipolar or exchange interactions that are too limited in range and weightability. Instead, one possibility is an architecture23 that uses charge currents which can be readily converted locally into spin currents through the spin Hall effect (SHE). These charge currents can be arbitrarily long-range and the total number of cross-couplings is only limited by considerations of routing congestion and delay. The couplings may also be confined to nearest-neighbors, simplifying the hardware design complexity while promoting scalability and retaining universality26.

The Ising Hamiltonian of Equation (1) can be implemented by exposing each SNM mk to a spin current Ik

An external file that holds a picture, illustration, etc.
Object name is srep44370-m4.jpg

which has a constant bias determined by hk together with a term proportional to the magnetization of the jth SNM mj. The future state of magnet mk at time (t + Δt) is related to the state of the other magnets at time t through the current Ik. This expression is derived analytically in the following section using the Fokker-Planck equation for the system31.

The spin current Ik can be generated using well-established phenomena and the prospects for physical realization of such a system are discussed later in this paper. The distinguishing feature of the present proposal arises from the intrinsic stochasticity of SNMs and their biasing through the use of weighted spin currents (Fig. 1(a)). How the SNMs are interconnected to implement Equation (2) can evolve as the field progresses.

Getting a large system to reach its true ground state is non-trivial as it tends to get stuck in local minima32. It is common to guide the system towards the ground state through a process of “annealing”2 which is carried out differently in different hardware implementations. For example, systems based on superconducting flux qubits make use of quantum tunneling, which is referred to as quantum annealing33, whereas classical CMOS approaches make use of random number generators34 to produce random transitions out of local minima.

For our system of coupled SNMs, random noise is naturally present and can be easily controlled (Fig. 1(a)), causing the system of SNMs coupled according to Equation (2) to explore the configuration space of the problem on a nanosecond timescale. Annealing could be performed through a controlled lowering of the actual temperature, or equivalently through a controlled increase in the magnitude of the current Ik, even at room temperature. It has been noted that certain annealing schedules can guarantee convergence to the true ground state, but these schedules may be too slow to be used in practice35. This paper only presents a straightforward annealing process and does not seek out optimal annealing schedules. Consequently, as we show in one of our combinatorial optimization examples, we may find only an approximate solution which, however, may be adequate for many practical problems.

Steady-State Fokker-Planck Description

Our goal is to interconnect magnets such that their equilibrium state is governed by Boltzmann statistics with thermal noise as an inherent characteristic of the system. To see that this is possible, consider a system of N magnets where we want

An external file that holds a picture, illustration, etc.
Object name is srep44370-m5.jpg


An external file that holds a picture, illustration, etc.
Object name is srep44370-m6.jpg

where mk represents the z-component of the magnets.

Suppose each magnet is driven by a spin current derived from the others. We start with the Fokker-Planck equation31 for the N-magnet system:

An external file that holds a picture, illustration, etc.
Object name is srep44370-m7.jpg

where An external file that holds a picture, illustration, etc.
Object name is srep44370-m8.jpg and ik = Ik/I0 with I0 as the critical switching spin current An external file that holds a picture, illustration, etc.
Object name is srep44370-m9.jpg. At equilibrium, An external file that holds a picture, illustration, etc.
Object name is srep44370-m10.jpg yielding from (3) and (5):

An external file that holds a picture, illustration, etc.
Object name is srep44370-m11.jpg
An external file that holds a picture, illustration, etc.
Object name is srep44370-m12.jpg

respectively. Comparing equations (6) and (7) while assuming symmetric coupling, An external file that holds a picture, illustration, etc.
Object name is srep44370-m13.jpg, for the system we find

An external file that holds a picture, illustration, etc.
Object name is srep44370-m14.jpg

and arrive at (2):

An external file that holds a picture, illustration, etc.
Object name is srep44370-m15.jpg

Stochastic Landau-Lifshitz-Gilbert (LLG) Model

In this section we briefly describe the simulation framework and stochastic LLG model used throughout this paper. We start with the LLG equation31 for a monodomain magnet with magnetization mi in the presence of a spin current An external file that holds a picture, illustration, etc.
Object name is srep44370-m16.jpg

An external file that holds a picture, illustration, etc.
Object name is srep44370-m17.jpg

The magnetic thermal noise enters the equation through the effective field of the magnet, Hi = H0 + Hn, as an uncorrelated external magnetic field in three dimensions with the following mean and variance:

An external file that holds a picture, illustration, etc.
Object name is srep44370-m18.jpg

The numerical model is implemented as an equivalent circuit for SPICE-like simulators and reproduces the equilibrium (Boltzmann) distribution from a Fokker-Planck Equation31.

A given system of magnets is simulated using a collection of independent, though current-coupled, stochastic LLG models. Delays associated with the communication from one magnet to the next are neglected assuming that the response time of the nanomagnets is much greater than associated wire-delays. Presently, the attempt time τ of experimental nanomagnets is on the order of ~μs to ~ms28,29,36. With additional scaling, the response times of these magnets will continue to improve37 and should approach the ~ns times discussed in this paper. With response times ~ns, our simulations show that even routing delays on the order of 100s of ps do not affect the results materially. Using nearest-neighbor Ising approaches or other constraining design decisions it should be possible to limit routing delays to shorter values. However, if the routing delay is comparable to the intrinsic response time of the nanomagnets then it would be important to include their effect in the simulation.

Many options exist, please see the final section, for physical realization of the proposed system of stochastic nanomagnets. For the simulations in this paper we simply use Equation (2) without assuming any specific hardware to implement it, since it is likely that better alternatives will emerge in the near future, given the rapid pace of discovery in the field of spintronics, see for example38,39,40,41.

Combinatorial Optimization

We will focus on two specific examples to demonstrate the ability of such an engineered spin glass to solve problems of interest42: an instructive example based on the satisfiability problem (SAT), and a representative example based on the traveling salesman problem (TSP). The first known NP-complete problem is the problem of Boolean satisfiability43, namely, deciding if some assignment of boolean variables {xi} exists that satisfies a given conjunctive normal form (CNF) expression. Finding the collection of inputs that makes the clauses of the CNF expression true is computationally difficult, but easy to verify.

It is known that any given CNF expression can be mapped to a collection of Ising contraints using the fundamental building blocks of NOT (An external file that holds a picture, illustration, etc.
Object name is srep44370-m19.jpg), AND (An external file that holds a picture, illustration, etc.
Object name is srep44370-m20.jpg), and OR (An external file that holds a picture, illustration, etc.
Object name is srep44370-m21.jpg) each subject to the Ising constraints given by44:

An external file that holds a picture, illustration, etc.
Object name is srep44370-m22.jpg
An external file that holds a picture, illustration, etc.
Object name is srep44370-m23.jpg
An external file that holds a picture, illustration, etc.
Object name is srep44370-m24.jpg

Using these building blocks, a network capable of finding the truth table for XOR An external file that holds a picture, illustration, etc.
Object name is srep44370-m25.jpg was prepared (Fig. 2). For simplicity, the solution uses a naive method to construct the network and leverages the use of ancillary spins to represent An external file that holds a picture, illustration, etc.
Object name is srep44370-m26.jpg and An external file that holds a picture, illustration, etc.
Object name is srep44370-m27.jpg respectively (note that four spins could have been used45). The array of spins from Fig. 2(b) are connected as specified by equations (11, 12, 13), driven by a reference current I0. As the magnets explore the configuration space, their outputs are digitized and used to compute the overall energy of the system (Fig. 2(c)). The regions of zero energy correspond to solutions of the problem. The digitized outputs are aggregated to determine their probability of occurrence. By looking at the first three bits of the most probable outputs, the solution to the problem can be directly found (Fig. 2(d,e)). While this problem helps convey the essence of the approach, a more demonstrative application is worth considering.

Figure 2
Boolean satisfiability with stochastic magnets.

The decision form of the TSP is NP-complete, that is, for a collection of N cities, does there exist a closed path for which each city is visited exactly once that has a tour length less than some value d? Finding tours that satisfy this problem is computationally challenging and also of great practical interest. There are well-known mappings that translate the TSP to the Ising model25,46. Here we adopt the following:

An external file that holds a picture, illustration, etc.
Object name is srep44370-m28.jpg

where xi,j is a Boolean variable that is TRUE when city i is stop number j and FALSE otherwise, and W(uv) are directed weights based on the distance between cities u and v. This Hamiltonian is mapped to a spin system by replacing each xij with 1/2(mij + 1) and weights W(uv) with i(uv) given by Equation (2).

If the interconnections between each city are symmetric, then a Boltzmann machine47 with each of the 2N×N states associated with an effective energy H is realized, and the probability of the system visiting a particular state is proportional to An external file that holds a picture, illustration, etc.
Object name is srep44370-m29.jpg. In order to find low-energy, optimized states, direct annealing of the glass can be performed. Using the ulysses16 reference dataset48, annealing of a problem specific magnetic array through control of the effective temperature was performed (Fig. 3). Two specific traits of interest arise, namely the energy decays in a sigmoidal relationship with the ln T, and the specific heat of the system, An external file that holds a picture, illustration, etc.
Object name is srep44370-m30.jpg, shows a defined peak about a critical temperature. At high temperatures, the system is disordered and corresponds to high energy states (Fig. 3(c)). As the temperature is reduced, the system continues to explore the energy landscape on a nanosecond timescale while gradually converging to a low-energy solution. For the given annealing profile and simulation duration, a low-energy, though not ideal, solution is found to the problem, highlighting the heuristic nature of the optimization46. Note that in principle these simulation results could be obtained directly from actual hardware. For example, Figs 2(d) and 3(d) could be obtained by continuously monitoring the states of the individual SNMs using spin valves.

Figure 3
Annealing stochastic nanomagnets for heuristic optimization of the traveling salesman problem.

Considerations for Physical Realization

Physical realization of these engineered spin glasses requires the integration of multiple functional elements as highlighted in Fig. 2(a). The magnetization of each magnet mi is first sensed with a read unit. The signal produced by this read unit is then propagated to all of the magnets with couplings dependent on the read magnetization mi. Each of these connections is independently weighted with weights W(ij) and provided as input, along with an on-site bias Bi to the write units. The write units in turn influence and control the state of magnet mj.

There are a number of design options available for each functional unit as shown in Table 1. Write-control of the magnets can be affected through a number of means including the spin Hall effect (SHE)39 or perhaps through voltage control40. The use of the SHE effect provides a convenient mechanism with which to sum several, independently weighted, input currents. Readout of the magnetization can be accomplished using well-established tunnel junctions49 which have been demonstrated for stochastic nanomagnets29. Alternatively, readout could perhaps be accomplished using the inverse SHE39. Assuming the use of a SHE material and tunnel junction stack, care must be given to accomodate the simultaneous use of write and read currents. One approach would be to introduce the use of a time-multiplexed scheme that disassociates the write and read operations50. Alternatively, structures that provide write and read isolation may be used51.

Table 1
Options for Physical Realization.

The ability to write and read the magnetization is of fundamental importance, however, once read, the likely weak signal must be amplified to satisfy the fanout requirements of the network. This transistor-like gain can be realized using all-spin based approaches23,51 or perhaps with the use of a hybrid-CMOS design52,53. These proposed approaches may introduce power dissipation challenges during the read operation, e.g. the short-circuit current produced with the use of amplifying inverters. Power dissipation considerations must be carefully evaluated to assess the viability of scaling the proposed system.

The output from the amplification stage can be selectively weighted so that a wide range of problems based on (1) can be encoded onto the network. The weighting of inputs can be based on an approach using re-programmable floating-gate voltages54 that would enable the use of analog weights for the circuit. While floating-gate regulation would enable convenient re-programmability, the design would be complicated with the requirement for peripheral drivers to control the floating-gate array. Others proposals have suggested the use of memristors24,53,55 or other programmable elements in a cross-bar like configuration50,56, though with constrained fanout. Note that one weighting scheme that still retains the ability to encode NP-hard problems onto the network is with the use of {−1, 0, 1} weights1. Using this simple approach removes the necessity for tunable weights and instead relegates the problem to one of routing, connectivity, and area.

All of the simulations used in this paper assume a fully connected network of magnets in which each magnet talks to all other magnets. For small networks this is reasonable, however, such an assumption is invalid for large networks as the number of routes grows rapidly. Instead, different topologies57 and routing considerations must be made to account for congestion and long-distance communication. By limiting the connections to local-neighbors9,10, the network may still be used to perform NP-hard optimization while also simplifying routing complexity. One design possibility is to leverage the lessons learned from the advances in the design of Field-Programmable Gate Array (FPGA) interconnects58. FPGAs are designed with routing topologies that facilitate both short and long-range interconnections while also providing re-programmability.

The fidelity of the programmed weights and number of high-fanout signals needed for robust solutions may impose challenges on the selected weighting and routing schemes. Additionally, the propagation delay of these high-fanout signals must be balanced with the response time of the magnets in order for the system to be governed by (1). While flexibility in the allowed weights and number of couplings is convenient for encoding problems onto the model25, it is important to note that discrete nearest-neighbor couplings still retain NP-hardness1 and may greatly simplify the hardware design, improving scalability at the expense of increased encoding complexity and area.

The main point of this paper is the remarkable high-speed search through Fock space enabled by the intrinsic physics of a network of stochastic nanomagnets interacting via spin-mediated interactions. We hope this work fosters an interest in the physical realization and exploration of stochastic nanomagnets as a viable Ising computer.


Simulations based on the modular framework for spintronics38 were used to produce the results in this work. Within the framework, a stochastic Landau-Lifshitz-Gilbert (LLG) model was used to simulate each nanomagnet. The magnetic parameters for the telegraphic PMA magnets used in the simulations are: effective anisotropy of PMA, An external file that holds a picture, illustration, etc.
Object name is srep44370-m31.jpg, saturation magnetization, An external file that holds a picture, illustration, etc.
Object name is srep44370-m32.jpg, damping coefficient, α = 0.01, and PMA diameter, Φ = 45 nm, amounting to a barrier height of Δ = 1 kT. In all simulations, the initial state of the magnetic array was randomly selected. Figure 1 was produced using a modular stochastic LLG simulation element with the input current swept from −2 μA to 2 μA in increments of 800 nA. At each current, the response of the magnet is observed for 10 μs. Figure 2 was simulated for 100 μs using the coupling depicted in the Figure and a reference current I0 of 2 μA. Figure 3 used an annealing schedule of Ti+1 = 0.9 Ti and Lagrange multiplier of An external file that holds a picture, illustration, etc.
Object name is srep44370-m33.jpg. At each temperature the magnets were allowed to randomly walk for 1 μs and were measured every 200 ps. The SAT and TSP magnetic networks were simulated using coupled stochastic LLG models with the intermagnet-coupling and on-site biases produced via the spin current term of the LLG equation. The magnetization of each magnet was digitized using Schmitt trigger based thresholds. HSPICE was used to solve the simultaneous coupled differential equations of the magnetic network.

Additional Information

How to cite this article: Sutton, B. et al. Intrinsic optimization using stochastic nanomagnets. Sci. Rep. 7, 44370; doi: 10.1038/srep44370 (2017).

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


This work was supported in part by C-SPIN, one of six centers of STARnet, a Semiconductor Research Corporation program, sponsored by MARCO and DARPA and in part by the National Science Foundation through the NCN-NEEDS program, contract 1227020-EEC.


The authors declare no competing financial interests.

Author Contributions B.S. designed the networks, performed the simulations, and wrote the paper. K.C. developed the modules describing the biased superparamagnets, B.B. provided insights on PMA magnets and S.D. developed the quantitative connection between spin-torque coupling and the Ising model. B.B. and S.D. conceived the use of stochastic magnets for solving problems of this type, while K.C. and S.D. helped write the paper. All authors discussed the results and helped refine the manuscript.


  • Barahona F. On the computational complexity of Ising spin glass models. J. Phys. Math. Gen. 15, 3241–3253 (1982).
  • Kirkpatrick S., Gelatt C. D. & Vecchi M. P. Optimization by simulated annealing. Science 220, 671–680 (1983). [PubMed]
  • Adleman L. M. Molecular computation of solutions to combinatorial problems. Science 266, 1021–1024 (1994). [PubMed]
  • Ouyang Q., Kaplan P. D., Liu S. & Libchaber A. DNA solution of the maximal clique problem. Science 278, 446–449 (1997). [PubMed]
  • Johnson M. W. et al. . Quantum annealing with manufactured spins. Nature 473, 194–198 (2011). [PubMed]
  • Perdomo-Ortiz A., Dickson N., Drew-Brook M., Rose G. & Aspuru-Guzik A. Finding low-energy conformations of lattice protein models by quantum annealing. Sci. Rep. 2, 571 (2012). [PMC free article] [PubMed]
  • Chua L. O. & Yang L. Cellular neural networks: applications. IEEE Trans. Circuits Syst. 35, 1273–1290 (1988).
  • Chua L. O. & Roska T. The CNN paradigm. IEEE Trans. Circuits System. I, Fundam. Theory Appl. 40, 147–156 (1993).
  • Ercsey-Ravasz M., Roska T. & Néda Z. Cellular Neural Networks for NP-hard optimization. EURASIP J. Adv. Signal Process 2009, 2:1–2:7 (2009).
  • Yamaoka M. et al. . A 20 k-spin ising chip to solve combinatorial optimization problems with CMOS annealing. IEEE J. Solid-State Circuits 51, 303–309 (2016).
  • Kim K. et al. . Quantum simulation of frustrated Ising spins with trapped ions. Nature 465, 590–593 (2010). [PubMed]
  • Mahboob I., Okamoto H. & Yamaguchi H. An electromechanical Ising Hamiltonian. Sci. Adv. 2, e1600236 (2016).
  • Shaked N. T., Messika S., Dolev S. & Rosen J. Optical solution for bounded NP-complete problems. Appl. Opt. 46, 711–724 (2007). [PubMed]
  • Oltean M. Solving the Hamiltonian path problem with a light-based computer. Nat. Comput. 7, 57–70 (2008).
  • Utsunomiya S., Takata K. & Yamamoto Y. Mapping of Ising models onto injection-locked laser systems. OOpt. Express 19, 18091–18108 (2011). [PubMed]
  • Wu K., García de Abajo J., Soci C., Ping Shum P. & Zheludev N. I. An optical fiber network oracle for NP-complete problems. Light Sci. Appl. 3, e147 (2014).
  • Wang Z., Marandi A., Wen K., Byer R. L. & Yamamoto Y. Coherent Ising machine based on degenerate optical parametric oscillators. Phys. Rev. A. 88, 063853 (2013).
  • Marandi A., Wang Z., Takata K., Byer R. L. & Yamamoto Y. Network of time-multiplexed optical parametric oscillators as a coherent Ising machine. Nat. Photon. 8, 937–942 (2014).
  • McMahon P. L. et al. . A fully-programmable 100-spin coherent Ising machine with all-to-all connections. Science aah5178 (2016). [PubMed]
  • Inagaki T. et al. . A coherent Ising machine for 2000-node optimization problems. Science 354, 603–606 (2016). [PubMed]
  • Bhanja S., Karunaratne D. K., Panchumarthy R., Rajaram S. & Sarkar S. Non-Boolean computing with nanomagnets for computer vision applications. Nature Nano. 11, 177–183 (2015). [PubMed]
  • Arnalds U. B. et al. . A new look on the two-dimensional Ising model: thermal artificial spins. New J. Phys. 18, 023008 (2016).
  • Behin-Aein B., Diep V. & Datta S. A building block for hardware belief networks. Sci. Rep. 6 29893 (2016). [PMC free article] [PubMed]
  • Locatelli N., Cros V. & Grollier J. Spin-torque building blocks. Nature Mater. 13, 11–20 (2014). [PubMed]
  • Lucas A. Ising formulations of many NP problems. Front. Physics. 2, 5 (2014).
  • De las Cuevas G. & Cubitt T. S. Simple universal models capture all classical spin physics. Science 351, 1180–1183 (2016). [PubMed]
  • Khasanvis S. et al. . Physically equivalent magneto-electric nanoarchitecture for probabilistic reasoning. Proceedings of the International Symposium on Nanoscale Architectures (NANOARCH), pp. 25–26 (2015).
  • Bapna M. et al. . Magnetostatic effects on switching in small magnetic tunnel junctions. Appl. Phys. Lett. 108, 022406 (2016).
  • Locatelli N. et al. . Noise-enhanced synchronization of stochastic magnetic oscillators. Phys. Rev. Applied 2, 034009 (2014).
  • Cowburn R. P., Koltsov D. K., Adeyeye A. O., Welland M. E. & Tricker D. M. Single-domain circular nanomagnets. Phys. Rev. Lett. 83, 1042–1045 (1999).
  • Butler W. H. et al. . Switching distributions for perpendicular spin-torque devices within the macrospin approximation. IEEE Trans. Magn. 48, 4684–4700 (2012).
  • Aaronson S. Guest column: NP-complete problems and physical reality. SIGACT News 36, 30–52 (2005).
  • Kadowaki T. & Nishimori H. Quantum annealing in the transverse Ising model. Phys. Rev. E 58, 5355–5363 (1998).
  • Cheemalavagu S., Korkmaz P., Palem K. V., Akgul. B. E. S. & Chakrapani L. N. A probabilistic CMOS switch and its realization by exploiting noise. Proceedings of the IFIP international conference on very large scale integration (2005).
  • Geman S. & Geman D. Stochastic relaxation, Gibbs distributions, and the Bayesian restoration of images. IEEE Trans. Pattern Anal. Mach. Intell. 6, 721–741 (1984). [PubMed]
  • Koch R. H. et al. . Thermally assisted magnetization reversal in submicron-sized magnetic thin films. Phys. Rev. Lett. 84, 23 (2000). [PubMed]
  • Urazhdin S., Birge N. O., Pratt W. P. & Bass J. Current-driven magnetic excitations in permalloy-based multilayer nanopillars. Phys. Rev. Lett. 91, 14 (2003). [PubMed]
  • Camsari K. Y., Ganguly S. & Datta S. Modular approach to spintronics. Sci. Rep. 5, 10571 (2015). [PMC free article] [PubMed]
  • Liu L. et al. . Spin-torque switching with the giant spin Hall effect of tantalum. Science 336, 555–558 (2012). [PubMed]
  • Heron J. T. et al. . Deterministic switching of ferromagnetism at room temperature using an electric field. Nature 516, 370–373 (2014). [PubMed]
  • Rojas Sánchez J. C. et al. . Spin-to-charge conversion using Rashba coupling at the interface between non-magnetic materials. Nature Commun. 4, 2944 (2013). [PubMed]
  • Karp R. M. Reducibility among combinatorial problems. In Complexity of Computer Computations (eds Miller R. E. & Thatcher J. W.) pp. 85–103 (Plenum Press, New York, 1972).
  • Cook S. A. The complexity of theorem-proving procedures. In Proc. 3rd Ann. Symp. on Theory of Computing 151–158 (ACM, 1971).
  • Bian Z., Chudak F., Macready W. G. & Rose G. The Ising model: teaching an old problem new tricks. D-Wave Systems 2 (2010).
  • Biamonte J. D. Nonperturbative k-body to two-body commuting conversion Hamiltonians and embedding problem instances into Ising spins. Phys. Rev. A 77, 052331 (2008).
  • Schneider J. J. & Kirkpatrick S. Stochastic Optimization (Springer, 2006).
  • Ackley D. H., Hinton G. E. & Sejnowski T. J. A learning algorithm for Boltzmann machines. Cognitive Sci. 9, 147–169 (1985).
  • Reinelt G. TSPLIB–A traveling salesman problem library. ORSA Journal on Computing 3, 376 (1991).
  • Parkin S. et al. . Magnetically engineered spintronic sensors and memory. Proceedings of the IEEE. 91, 661–680 (2003).
  • Sengupta A., Parsa M., Han B. & Roy K. Probabilistic deep spiking neural systems enabled by magnetic tunnel junction. IEEE Trans. Electron Dev. 63, 2963–2970 (2016).
  • Datta S., Salahuddin S. & Behin-Aein B. Non-volatile spin switch for Boolean and non-Boolean logic. Appl. Phys. Lett. 101, 252411 (2012).
  • Tangel A. & Choi K. “The CMOS Inverter” as a comparator in ADC designs. Analog Integr. Circuits Signal Process 39, 147–155 (2004).
  • Sengupta A., Choday S. H., Kim Y. & Roy K. Spin orbit torque based electronic neuron. App. Phys. Lett. 106, 143701 (2015).
  • Diep V. Q., Sutton B., Behin-Aein B. & Datta S. Spin switches for compact implementation of neuron and synapse. Appl. Phys. Lett. 104, 222405 (2014).
  • Yang J. J., Strukov D. B. & Stewart D. R. Memristive devices for computing. Nature Nanotech. 8, 13–24 (2013). [PubMed]
  • Sengupta A., Shim Y. & Roy K. Proposal for an all-spin artificial neural network: emulating neural and synaptic functionalities through domain wall motion in ferromagnets. IEEE Trans. Biomed. Circuits Syst. 99, 1–9 (2016). [PubMed]
  • Bunyk P. I. et al. . Architectural considerations in the design of a superconducting quantum annealing processor. IEEE Trans. Appl. Supercond. 24, 1–10 (2014).
  • Lemieux G. & Lewis D. Design of interconnection networks for programmable logic (Springer, Boston, 2004).

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group