
2
risotrs, and make the connection between the statistical
properties of the circuit and ones of the Ising Hamilto-
nian. In the same section, we briefly discuss the nu-
merical approach used in our work. The results of our
simulations are presented in Sec. II with the emphasise
on the possibility of reaching FM and AFM interactions
in the circuit. The paper ends with a conclusion.
II. METHODS
Mathematically, we utilize an effective Ising-type
Hamiltonian to describe the probabilities observed in
the circuit simulations. For the circuit in Fig. 1(b), the
Hamiltonian has the form
H=−JX
i
σiσi+1 −J2X
i
σiσi+2 −hX
i
σi,(1)
where Jis the interaction coefficient for adjacent spins,
J2is the next-to-adjacent interaction, his the magnetic
field, and periodic coupling is assumed. Schematically,
these interactions are presented in Fig. 1(c). We consider
the electronic circuit as a physical system described by
the Boltzmann distribution
pi=1
Ze−Ei
kT .(2)
Here, Z=P
j
e−Ej
kT is the statistical sum, and Ej-s are the
“energies” of circuit states. We argue that for the circuit
in Fig. 1and similar circuits these “energies” correspond
to the Ising Hamiltonian (1).
To explain the co-existence of AFM and FM interac-
tions, consider a set of identical memristors in ROF F sub-
jected to a positive voltage pulse driving the OFF-to-ON
transition. Each memristor will have an equivalent prob-
ability of being the first to switch states. When one of
these memristors swaps states, it reduces the probability
to switch of its neighbors (reducing the voltage across
them). In this scenario, memristors with neighbors both
in the ROF F state will have the highest chance of switch-
ing. This leads to the tendency of antiferromagnetic or-
dering in the memristors under a positive voltage pulse.
However, under a negative voltage pulse the RON state
memristors with neighboring ROF F state memristors will
be favored to change states. Meaning, the configuration
will tend towards ferromagnetic ordering under a nega-
tive voltage pulse. The overall ordering of a memristive
circuit driven by an AC source will then be dependent
on the choice of model parameters for the memristors.
Based on the parameters, one type of ordering may be
dominant.
Next, we introduce the model of stochastic memris-
tors. According to experiments with certain electrochem-
ical metallization (ECM) cells [21,22] and valence change
memory (VCM) cells [23], the probability of switching be-
tween resistance states of these devices can be described
1000 1500 2000 2500 3000
0
2
4
6
8
1 0
1 2
1 4
1 6
S t a t e
T i m e ( a r b . u n i t s )
FIG. 2. Dynamics of the states in a circuit with N= 4
memristive spins. The circuit has 24= 16 states that are
labeled from 0 to 15. The 0000 state (all memristors are in
ROF F ) is labeled by 0, 0001 by 1, and so on. This plot was
obtained using the following set of parameters: R=r=
ROF F = 1 kΩ, RON = 100 Ω, τ01 = 3 ·105s, τ10 = 160 s,
V01 = 0.05 V, V10 = 0.5 V, Vpeak = 1 V, and T= 2 s.
by switching rates of the form
γ0→1(V) = τ01e−V /V01 −1, V > 0
0,otherwise ,(3)
γ1→0(V) = τ10e−|V|/V10 −1, V < 0
0,otherwise ,(4)
where Vis the voltage across the device, and τ01(10) and
V01(10) are device-specific parameters. Here, 0 and 1 cor-
respond to the high (ROF F ) and low (RON ) resistance
states, respectively. Under a constant voltage, the prob-
ability to switch follows the distribution [21,22]
P(t) = ∆t
τ(V)e−t/τ(V),(5)
where τ(V) is the inverse of the switching rate given by
Eq. (3) or (4) (depending on the sign of V). Previously,
we have developed a master equation approach for the
circuit of stochastic memristors [24] and designed its im-
plementation in SPICE [25].
Most of the results presented here are obtained through
numerical simulations of the circuit in Fig. 1(b) contain-
ing Nmemristive spins. The set of parameters defining
the circuit and the simulations such as the model con-
stants, voltage period, duration, resistances, etc. are first
set. The memristors are then initialized to their starting
states (typically all OFF). The voltages across each mem-
ristor are calculated for the current time-step through
Kirchhoff’s laws. The switching time is then generated
for each memristor randomly with Eq. (5) distribution.
The fastest switching time is extracted and compared to
the remaining time in the current voltage pulse. If there
is sufficient time remaining in the pulse, that memristor
switches states and the the time remaining in the period