1. Introduction
The name “memristor” (memory resistor) was introduced in the early 1970s in the original work by L. Chua [
1]. The main feature of the memristor is the fact that its resistance depends on the charge passed through it. Therefore, this device retains the memory of past states. Chua’s work remained unnoticed until 2008 when D. Strukov and colleagues [
2] reported on an experimental implementation of the memristor, which aroused strong interest in the community. Later works challenged Chua’s conceptual line, questioning the memristor’s status as the fundamental fourth element [
3]. Despite the difficulties, D. Strukov’s study undoubtedly attracted the attention of researchers. It was quickly realized that memristors had the potential to revolutionize electronics, allowing powerless storage and logic operations [
4], as well as simulating the behavior of neuronal synapses [
5,
6]. These observations have opened up a wide field of applications in physical neural networks and neuromorphic architectures [
7,
8,
9,
10,
11,
12,
13,
14]. Moreover, the memory formalism [
15] is applicable to a wide range of physical and biological systems allowing to be extended far beyond the electronic realm.
In its most general formulation, a memristive device is determined by some input (
x) and output (
y) parameters (in the case of a memristor, these are current and voltage, respectively) and the so-called state parameter
, which changes according to a certain law over time [
16]. These parameters are related by the equations:
where
f and
g are the functions defining the system dynamics. Such memristive devices include, in particular, memcapacitors and meminductors. The input–output characteristic—the dependence of the output parameter
y on the input
x—is usually a hysteresis curve.
Mostly, in all memristors designed and studied so far, the input and output signals (x and y) are classical. A legitimate question arises about the possibility of creating a memristive device that exhibits a hysteresis curve in the input–output characteristic and at the same time has a truly quantum behavior that allows the manipulation of quantum information. According to the current trend, such a device is referred to as a quantum memristor. A quantum memristor has to demonstrate “memristive behavior” for mathematical expectations of quantum observables and must be able to coherently map a quantum input state to an output state. In this regard, a quantum memristor requires the development of an open quantum system in which a quantum device interacts with the environment, for example, via a measurement process. However, this is always associated with certain level of decoherence. This difficulty can be overcome by designing the interaction with the environment in such a way that it is strong enough to provide “memristive behavior”, but weak enough to sufficiently preserve quantum coherence.
To date, the development of quantum memristors is just beginning and their potential impact on neuromorphic computing is difficult to assess. The concept of a quantum memristor was introduced by Pfeiffer et al. [
17]. The motivation was to see if both the advantages of quantum computing and memristive behavior could be exploited simultaneously. The advantages of quantum systems provide certain advantages in computation power, while memristivity is a nonlinear feature that goes beyond purely unitary evolution and allows for the implementation of quantum neuromorphic computing.
Two different physical platforms have been studied so far: superconducting circuits [
18] and quantum photonics [
19,
20,
21]. Superconducting circuits use memory effects that naturally occur in Josephson junctions [
18]. In particular, such effects were used in the recently proposed design of the classical superconducting memristor [
22]. In [
18], a quantum memristor based on superconductors used quasiparticle tunneling, the control and measurement of which has recently advanced significantly [
23,
24].
For a quantum photonics platform, the key element of a quantum memristor is a beam splitter with a tunable reflectivity, which is modified depending on the results of measurements in one of the outgoing beams [
19]. This feedback mechanism provides memory and nonlinear behavior, while photonic degrees of freedom implement quantum coherence. A similar implementation of a quantum memristor can be achieved using frequency-entangled optical fields and a frequency mixer which creates superpositions of states operating similarly to a beam splitter [
20].
The first experimental implementation of a quantum memristor based on an integrated photonic circuit, which was fully reconfigurable using built-in phase shifters, was demonstrated in [
21]. This setup is able to exhibit “memristive behavior” of single-photon states using a measurement circuit and classical feedback. The phase shifter is controlled externally by connecting one of the output waveguides through an optical fiber to the micro-controller, which directly controls the phase shifter. The authors also analyzed the possibility of using quantum memristors to increase performance in machine learning problems. They performed numerical simulations of a neural network using purely classical resources and compared it with a network using only three quantum memristors [
21]. Numerical results showed that even with such a limited number of quantum memristors, a significant increase in performance can be achieved; an order of magnitude fewer operations were required to train the network using quantum memristors compared to the classical case. Moreover, a quantum network could classify images with up to 95% fidelity, while a classical analogue would only achieve 71% accuracy. Although it is difficult to compare neural networks with each other, whether classical or quantum, these results point to a possible advantage of using quantum memristors in machine learning.
An interesting proposal for the implementation of quantum memristors was made in [
25]. It has been shown that the nonlinear dynamics of a non-Markovian nature for quantum memristors can be obtained from the unitary dynamics of a larger system by tracking some of its degrees of freedom associated with auxiliary subsystems. This means that one can obtain efficient families of non-unitary and non-Markovian dynamics in quantum computers by considering a subset of quantum processor qubits. Using this approach, the authors [
25] simulated the dynamics of a quantum memristor on a quantum computer and showed that the memristive effect can be achieved via interaction between qubits.
At the moment, not all possible platforms for quantum memristors have been analyzed and demonstrated. In particular, trapped ions, being one of the most successful physical platforms for quantum computations, have not been considered yet. An ionic quantum memristor may have a number of advantages over other counterparts (photonic and superconducting) due to the presence of a sufficiently large number of ion levels with different lifetimes and transitions with different oscillator strengths. Furthermore, an ion-based memristor benefits from straightforward possibility to use common vibrational modes in the trap for ion–ion entanglement. An additional advantage of the ion-based memristive concept is the fact that it can be realized as a standard for ion-based quantum computing and optical atomic clock techniques [
26,
27].
In this paper, we propose and elaborate the idea of a quantum memristor based on ultracold ions in a Paul trap.
2. Trapped-Ion Quantum Memristor
For an ultracold trapped ion one can choose three electronic levels of interest: the ground state
, the excited long-living state
, and the excited short-living state
. Using various techniques (resonant laser or microwave excitation, Raman excitation), it is possible to initiate Rabi oscillations between levels
and
. The wave function of such a two-level system after the interaction with resonant to the
−
transition laser field is given by the superposition
where
. The excitation probability for
changes with time
t as [
28]
where
—Rabi frequency for
. The formula (
3) represents the wave function of the two-level system interacting with a resonant laser field.
It is also possible to electromagnetically couple the level
to a level
using a resonant field. The transition probability can be controlled by varying the laser pulse duration
and/or its amplitude. The transition
probability equals [
28]
with
being the corresponding Rabi frequency for
. Accordingly, the probabilities
and
can be changed from 0 to 1. It must be possible to detect the level
excitation with some technique (make a projection measurement).
Figure 1 schematically shows the level scheme.
If one excites
transition, the wave function of the system can be written as
As a result, the wave function describes the three-level system interacting with two resonant laser fields.
The input signal
x is defined by the mathematical expectation of the ion being in the state
before laser pulse 2 which excites the transition
,
As for the output signal
y, we set the mathematical expectation of the ion being in the state
after laser pulse 1 exciting the transition
,
In this case, the dependence of the output signal on the input is given by
As already noted,
can be varied over time by the laser 2 pulse shape, so it plays the role of the state parameter
s, i.e.,
. Equation (
9) takes a form similar to Equation (
1), which is typical for a classical memristor. It is possible to change the state parameter by selecting the duration of laser pulses and/or their amplitude in an arbitrary manner. For the experimental implementation of a quantum memristor, it is necessary to be able to measure
s during some time interval
and to organize feedback to change
s for the next interval
. The state parameter
s changes over time according to an expression similar to (
2):
where
c is an arbitrary constant. The dependence (
10) is not the only possible one, it can be chosen in an arbitrary way. It generally comes from the experimental observable behavior. We chose it as a closely corresponding to the equations that describe changes in resistance of classical memristors. By integrating (
10) we obtain
where
is the state parameter value at the initial time. The quantum memristor operation procedure is the following. At the initial time, a certain probability
is set. Parameters of the laser 2 pulse are set correspondiongly. The time interval
T is chosen, during which the transition
is excited and the population of
is detected subsequently. The probability of the ion being in the state
in this case equals
Hence, by experimental determination of
one can calculate the input signal
x on the first time interval
T. Let us denote this value as
. Using (
9) for
, one can obtain the output signal
. By substituting
into (
11) one can determine the new state parameter
s. For this value of
,
and
are calculated similarly to
and
, and then the value of
is found. After that the procedure is repeated. The result is a dependence of the output signal
y on the input
x.
The theoretical input–output characteristics calculated for
are shown in
Figure 2. It can be seen that with a change of
T, the form of
also changes and for some values of
T the input–output characteristic demonstrates memristive behavior.
A necessary condition for the implementation of the scheme described above is that the total registration time of the
state population
and the duration of the pulse that excites the transition
(
) should be much less than the period of Rabi oscillations between levels
and
:
Inequality (
13) limits the temporal characteristics of the laser field used in the proposed scheme: the temporal width of the second laser should be much less than both the temporal width of the first laser and the period of Rabi oscillations generated by the first laser.
3. Population Dynamics under the Action of Two Laser Fields
To study the influence of laser pulse parameters on the memristive characteristics, we simulated the population dynamics of a three-level system under the action of a sequence of resonant laser pulses. For this, we used the system of equations for the population amplitudes
a considering an atom interacting with two resonant laser fields as described in [
29]. In the absence of spontaneous relaxation, this system of equations allows us to calculate the population dynamics of atomic levels for all types of three-level schemes (
—scheme, cascade scheme, etc.):
where
is a time delay between laser pulses and
are laser pulse durations (
is a center of the pulse of the first laser,
is a center of the pulse of the second laser). By
, the complex conjugation is indicated. Calculations were carried out for the Rabi frequencies, the pulse durations and the time delays in non-dimensional units.
Figure 3a shows the dynamics of the level population amplitudes calculated for the following parameter values:
and
. Due to the impact of the first pulse, Rabi oscillations of the level populations
are observed. Under the action of the second resonant field pulse, the level
is excited. With given values of the field parameters and by changing the value
in accordance with Equations (
10) and (
11), we calculated
y and
x. As a result, a hysteresis dependence of
y on
x was obtained as shown in
Figure 3b.
Note that this dependence of
was calculated for a three-level system interacting with two resonant laser fields. It corresponds to one of the dependences which were calculated in accordance with (
9) and (
10) (please see
Figure 2). It demonstrates a principal opportunity to find values of y and x for a three-level system in the quantum case. It is important that the key control parameter is the time interval
T in Equation (
10).
Figure 3c shows the dependence of the
level population calculated for various time delays between pulses. It can be seen that the level population slightly changes with variation of the specified parameter (red curve in
Figure 3c). It is important to note that as the duration of the second pulse increases, the population of the
level can change significantly with the change of the time delay between pulses (black curve in
Figure 3c). This may be an additional possibility of obtaining a hysteresis dependence of
y on
x.
4. Three-Level Scheme for the Case of 171Yb+ Ion
For the experimental implementation of an ion-based quantum memristor, it is necessary to achieve a deeply cooled ion in a Paul trap. For single-ion manipulation, one can use three-dimensional traps, which are somewhat simpler in manufacturing than linear ones and have good optical access to the ion. There are a number of regular methods of ion cooling. For example, Doppler or sympathetic laser cooling [
30,
31] allows for the achievement of ion temperatures on the order of several mK or even lower. For cooling down to the ground vibrational state, side-band cooling [
32] and EIT cooling [
33,
34,
35,
36] are successfully used. The most difficult requirement in the experimental implementation of an ion-based memristor is the condition (
13). This difficulty is associated with the detection of the
state excitation. The observation of the spontaneous decay of this state on a single ion implies the high-fidelity detection of a single photon, which is extremely difficult. As a rule, the detection level population in an ion is determined by the well-known method of electron shelving [
37,
38,
39]. For this the lifetime of the state
must be much shorter than the lifetime of the state
to provide signal accumulation and reliable detection of the state
population. Typically, if one uses an electric dipole transition with a spectral linewidth of tens of MHz, the electron shelving requires several ms for robust measurement.
Accordingly, the lifetime of the level
should be much longer. For example, in the case of the
ion, the hyperfine sublevels of the
level can be assigned as
and
, respectively (
Figure 4). The energy difference between these states is 3.6 GHz. The level
can be initially populated by, e.g.,
-pulse from the ground state. To couple
and
(magnetic–dipole transition) one can use the radiofrequency pulses at 3.6 GHz (standing in for laser pulse 1 in
Figure 1). The only decay channel for the
state is the octupole transition to the ground
state. The rate of spontaneous emission at this transition is almost zero (the lifetime equals 6 years), which allows us to arbitrarily set the Rabi oscillation period
necessary to satisfy condition (
13). In turn, the ground state
can be chosen for the state
. Its population can be straightforwardly detected by electron shelving. A feature of the proposed scheme is that the level
has less energy compared to
; this is annoying, taking into account the characteristics of the structure of the
ion to fulfill the condition (
13).
Although in the discussion above we considered two laser fields and electric dipole transitions for simulation, the results remain true for the described case of magnetic dipole and electric octupole transitions in the ytterbium ion.
Experimental realization of the proposed ion-base quantum memristor concept is rather complicated: accurate pulse shaping, laser frequency tuning and high-fidelity state reading, in addition to a high-vacuum system with a radiofrequency trap must be under tight control. Furthermore, manipulating within the octupole transition using a narrow-line laser is required, but it can be achieved with the techniques used for ion-based quantum computing [
40] and atomic clocks. As for manipulating theoctupole transition, this was also already demonstrated [
41].
The necessary states for the implementation of a quantum memristor can be found in other ions as well.
5. Coupled Quantum Memristors
As we mentioned above, in the case of ions there is an advantage due to their strong coupling originating from the Coulomb interaction. It is possible to entangle two or more ions. For trapping of the ion chain, it is convenient to use a linear Paul trap. Laser cooling in all directions of motion results in chain crystallization along the trap axis, where the amplitude of the radiofrequency field vanishes. Due to the Coulomb repulsion between the ions, there is a number of ion oscillation normal frequencies. The lowest frequency corresponds to the common vibrational mode of the center of mass. Let us consider a system of two ions, which we designate as the “first” and the “second”. Similar to the case of a single ion, one can initiate Rabi oscillations with a frequency
between
and
of the first ion by means of resonant radiation at the frequency
. By setting the detuning
corresponding to the vibrational mode frequency, one can address adjacent common vibrational mode
and
. Here, the symbol
denotes the center-of-mass oscillation mode with
i quanta. The probability of this transition also can be changed by varying the laser pulse shape. Schematically, the level structure is shown in
Figure 5.
The transition indicates the birth of a phonon. The second ion, if it initially was in the state , will change the state to . An experimental proof of the phonon birth can be performed using electron shelving. Then as an input signal x, one can chose the mathematical expectation of the first ion being in the state , and as the output signal, the presence of the second ion in the state . The state parameter is the transition probability , which can be changed over time by implementing the feedback. The dependences will have a character similar to the corresponding dependences for the case of implementing a quantum memristor on a single ion.
The difficulties of implementing such a scheme include the need of cooling the ions to the ground vibrational state and controlling of the heating of the ions in the trap. The heating rates can be determined from the dephasing of the Rabi oscillations [
42] and should be much lower than the memristor coding procedure. The phonon interaction possesses a number of advantages. First, there is the possibility of encoding information on the second ion, i.e., setting the probability of finding the second ion in the state
by exposing it to an additional laser pulse. Second, one can straightforwardly scale up this system up to several tens of ions.