# Computing with Dynamical Systems Based on Insulator-Metal-Transition Oscillators

## Abstract

In this paper we review recent work on novel computing paradigms using coupled oscillatory dynamical systems. We explore systems of relaxation oscillators based on linear state transitioning devices, which switch between two discrete states with hysteresis. By harnessing the dynamics of complex, connected systems we embrace the philosophy of “let physics do the computing” and demonstrate how complex phase and frequency dynamics of such systems can be controlled, programmed and observed to solve computationally hard problems. Although our discussion in this paper is limited to Insulator-to-Metallic (IMT) state transition devices, the general philosophy of such computing paradigms can be translated to other mediums including optical systems. We present the necessary mathematical treatments necessary to understand the time evolution of these systems and demonstrate through recent experimental results the potential of such computational primitives.

a]Abhinav Parihar b]Nikhil Shukla b]Matthew Jerry b]Suman Datta a]Arijit Raychowdhury a]Georgia Institute of technology b]University of Norte Dame

## 1 Introduction

Computing is the backbone of the modern society; from economics to security, scientific advancement to social welfare, each and every aspect of our lives has been enriched by the technology revolution. We have enjoyed the benefits of Moore’s Law over the last four decades as technology scaling brought the power of supercomputers to our smartphones. With increasing challenges in scaling, came ground-breaking innovations in the transistor technology. As we look ahead, limits of traditional scaling are already in sight. The demise of Dennard scaling and slowing down of Moore’s Law have further exposed the fundamental scaling limitations of the Von Neumann execution models of computing. This transition is accompanied by the realization that in a fast evolving, socially interconnected world, we are observing a seismic shift in the amount of unstructured data that need to be processed in real-time; and consequently future systems will be limited by the energy growth of data movement rather than compute. Therefore, we need fundamentally new approaches to sustain the exponential growth in performance beyond the end of the CMOS roadmap. This will require new execution models coupled with new devices to implement them. In particular, we observe that new models that deal with data analytics have compute and storage interleaved in a fine grained manner - not separated as in the Von Neumann world. Moving forward, computing technology will heavily penalize separation of data and compute and we need to marry them in better ways to handle emergent applications. This may necessitate that barriers of abstraction are broken. Next generation of computational models should map natively to the physics and dynamics of the physical devices, without a Boolean abstraction. Of course, it is a fool’s endeavor to assume that the Von Neumann architecture will perish; rather newer computing models and dedicated hardware accelerators will supplement traditional Von Neumann machines in “data centric” tasks.

In today’s computing landscape, ever harder problems are being encountered each day. From social networks to graph analytics, from weather prediction to scientific computation, computationally hard problems are present everywhere. Some of these problems are challenging just because of the sheer size of the data-sets. But many other problems like optimization problems, though small in size, are intractable because of their inherent complexity and often combinatorial nature of solution spaces[1]. The sequential Von Neumann machine has been useful for most problems we have encountered in the past, but as the complexity of problems increase, significant research is underway for alternative paradigms, architectures and hardware which can be used for solving complex problems more efficiently. Most of these efforts borrow essential ideas from natural computing processes. These include distributed computing, distributed memory, integrated processing and memory, shifting information representation from symbolic to physically meaningful quantities, and switching from sequential discrete time to continuous time dynamics.

Among the different problems of interest, associative computing, scientific computing (including solution of coupled Partial Differential Equations (PDEs)) and optimizations form important classes of data-centric as well as model-centric computations. Active research in all these areas suggest that analog or continuous time systems may offer alternative, faster and more energy-efficient solutions that their traditional digital counterparts. For example, in case of hard optimization problems alternative paradigms and architectures include, but are not limited to, cellular automata[2], quantum computing[3], Ising model based systems[4, 5], neural networks[6], stochastic searching architectures[7] and memcomputing[8]. Among architectures for solving PDEs cellular neural networks (CNN) have been studied[9, 10]. Some studies have also suggested using cellular automata for solving PDEs[11].

The basic philosophy of most of these networks for optimization, e.g. artificial neural networks, is to first come up with an energy function which can be a penalty function or a rewarding function depending on how far the current solution is from the optimal solution[6]. The next step is to tune the parameters of the network such that as time evolves the dynamics of the system decreases the penalty function, or increases the reward. But even if such a massively parallel system is devised which can solve a NP-hard optimization problems, exponential resources of space, time or precision will be required[12, 13]. Another direction is often explored where instead of trying to solve the optimization exactly, an approximation is targeted which works well for most problems on average and allows less optimal solutions in harder instances. There can be other kinds of trade-offs like in the case of Hopfield networks where even though the optimal solution maps exactly to global minima, there can be too many local minima where the system can get trapped. Finally, the physical layer of computing, namely the semiconductor device platform needs to be able to support such systems and the CMOS transistor is not always an optmial device choice.

Recently, the development of novel phase transition materials like Vanadium Dioxide (VO\textsubscript2) and corresponding electronic devices, which show insulator-to-metal (IMT) transitions[14] have sparked interest in creating compact relaxation oscillators[15, 16]. These oscillators, when coupled to each other exhibit phase synchronization which can be used for phase based computing. Such new kinds of devices present interesting opportunities to create systems with novel synchornization dynamics. The impact of using such devices as basic units in circuits can break the abstraction between the physical and the algorithmic layers of computing. It should be noted that the synchronization dynamics of coupled oscillators not only have a wide variety of applications in engineering but they also explain many natural, chemical and biological synchronization phenomena like the synchronized flashing of fireflies, pacemaker cells in the human heart, chemical oscillations, neural oscillations, and laser arrays, to name a few. These novel computing primitives, of course, are neither drop-in replacements for CMOS transistors nor straightforward extensions of the existing computing architectures. It requires rethinking of the basic computational entities in new kinds of system architectures. In this paper we review some computational models using coupled relaxation oscillators based on VO\textsubscript2metal-insulator-transition devices focusing on how the system dynamics can be modeled and the applications they can enable. We limit our discussion on applications in image processing for the sake of brevity. However, the potential of dynamical systems extend far beyond image analytics and promises to be a competing computational model for post-CMOS technologies.

## 2 A Perspective on Coupled Oscillatory Networks

The area of coupled oscillators has been dominated mostly by theoretical models and numerical simulations, but very few successful physical implementations. The reason being the assumptions made for analytical simplification in those theoretical models are too difficult to realize in practice. Also an important limitation of such systems, which is true for any dynamical system, is that if a dynamical system is able to solve computationally hard (NP-class) problems exactly, then it necessarily have to be chaotic in nature which would require exponential precision in both simulations and physical implementations.

The most popular coupled oscillator models in this area are the Kuramoto oscillator models[17] which rely on sinusoidal oscillators coupled using “weak” and linear phase coupling. A Kuramoto system of N oscillators is described by

where and are the phase and frequency respectively of oscillator. Major challenges in this kind of coupled oscillator model is the notion of weak coupling, i.e. and the idea that the coupling effects phase only without disturbing the frequencies. Moreover, creating arrays of compact sinusoidal oscillators with many oscillators coupled to each other pose serious challenges given the requirements Similar models of weak linear phase coupling were also explored for Van Der Pol oscillators[18] which have an additional nonlinearity. But the implementation of such oscillators is also non-trivial and the coupling behavior becomes too complicated to tackle large connected networks[19, 20, 21, 22].

Nevertheless, the ability of coupled oscillatory systems to encode computing has long been realized. Associative computing, with applications in pattern detection and machine learning have been demonstrated in theory. Similarly oscillatory cellular neural networks have been show to possess extraordinary computing ability in solving problems as varied as template matching, PDEs and ODEs and so on. More recently, this effort has been augmented by advances in the development of compact oscillators in non-silicon technologies. One prominent effort is the use of Spin Torque Oscillators (STOs) coupled using spin diffusion currents and providing a computational platform for machine learning, spiking neural networks and others [23, 24]. However, the high current densities of STOs and the limited range of spin diffusion currents continue to pose serious challenges to technologists.

In our studies, we explore phase transition based relaxation oscillators with piecewise linear dynamics, which means the system is described by different linear dynamical systems in different “states” of the system. These states are basically charging and discharging of a capacitive element. The coupling is also electronic and is accomplished by linear capacitors and/or resistors. The repeated switching between these states gives rise to oscillatory behavior and the notion of phase. But because the switching is itself determined by the state variables (voltage thresholds), and not explicitly by time, the coupling dynamics and hence the overall system dynamics are often mathematically intractable and closed form approximations like the Kumamoto model are not possible. This further complicates the analysis when the goal is to perform computation in phase space. Hence model development coupled with numerical analysis and an intuitive understanding of how complex systems evolve over time become essential tools to engineer such systems. The relaxation oscillators we investigate, are built using phase change devices, which are devices that switch state between a metallic state with low resistance and and insulating state with high resistance depending on the voltage across them. In the next section, we will describe simple mathematical constructs that can assist in analyzing these oscillators and the results obtained once they are coupled electrically.

## 3 Relaxation oscillators based on phase change devices

### 3.1 State changing devices

The state changing devices we consider are essentially linear conductances (or resistances), but can transition between two conducting states - insulating state with conductance and metallic state with conductance . We assume that . They are also called Insulator-Metal-Transition (IMT) devices. We use Vanadium Dioxide (VO2) as the material choice. A state transition is triggered by the voltage applied across the devices as well as the history. When the voltage across the device increases above a higher threshold voltage the device changes its state to metallic with conductance and when the voltage decreases below a lower threshold it changes to insulating state with conductance . There is hysteresis in switching, i.e. , which means when the voltage applied is between and the device retains the last state it was in. An internal capacitance is associated with the device which ensures gradual buildup and decaying of the voltage (and energy) across the device. Rigorous mathematical analysis of such oscillator configurations can be found in [25].

### 3.2 Single oscillator configurations

There can be multiple configurations/circuits of a relaxation oscillator based on state changing device. On a simplistic level, two basic configurations exist - (a) two state changing devices in series, called D-D where D stands for device (figure 1a), and (b) a state changing device in series with a resistance, called D-R (figure 1b). In the former, the charging and discharging rates are equal, but they are different in the latter. The functioning of the two circuits is as follows.

In case of two devices in series (D-D), the two devices must be in opposite conduction states (one metallic and the other insulating) all the time for oscillations to occur. If the threshold voltages and are equal for the devices and the following condition holds:

(1) |

and at the devices are in different conduction states, then any time one device switches, the other will make the opposite transition as well. As , the devices can be considered as switches which are open in insulating state and closed in metallic with conductance . The mechanism of oscillations is essentially charging and discharging of the internal capacitances of the devices. The device in metallic state connects the circuit and charges (discharges) the lumped internal capacitance. The voltage at the output node increases (decreases) and eventually reaches the threshold voltage. Because of (1) both devices will switch at the same time causing their behavior to switch. The charging (discharging) becomes discharging (charging) and the cycle continues. The modeling of a DD oscillator is as follows. All the lowercase voltages referred in the paper are normalized voltages with respect to . Which means and . Also is used as normalized and hence .

The single D-D oscillator can be described by the following set of piecewise linear differential equations:

where and are metallic conductances of the two devices respectively. As there is no term involving in the equations. The equation can be re-written as:

where denotes the conduction state of the device (0 for metallic, and 1 for insulating) and and depend on the device conduction state as follows:

For D-R oscillators, the oscillations occur due to a lack of stable point as seen in the load line graph of figure 1c. Solid lines with slopes and are the regions of operation of the device in insulating and metallic states respectively. The system does not enter the dashed line region as a transition occurs to the other conduction state at the red points. The stable points, denoted by green points, are the points where the load line intersects the I-V curve of the device. These stable points in each state lie outside the region of operation the circuit and hence the circuit shows self sustained oscillations. This is a much more practical configuration from an electrical implementation point of view, as the conditions required for oscillations are not very strict. Following similar analysis as in the D-D oscillator case, the dynamics of the single D-R oscillator can be described as:

which can be re-written as:

where,

and denotes the conduction state of the system as before. Detailed analysis of configurations and modeling of such oscillators can be found in [25].

### 3.3 Pairwise coupling

Analysis of two coupled relaxation oscillators can give interesting insights into how such coupling dynamics can be used in various computing applications, and can also help understand and exploit dynamics from complex couplings. There can be many ways in which the oscillators can be coupled. We’ve looked at coupled oscillator circuits where the oscillators are coupled through their output nodes using a capacitance, a resistance, or a parallel RC combination (figure 2).

#### D-D oscillators

The D-D configuration, though difficult in an electrical implementation point of view, is very simple to analyse and gives interesting insights about dynamics of such piecewise linear systems. Two identical D-D oscillators coupled using a RC circuit can be modeled as follows. When coupled, the system has 4 conduction states corresponding to the 4 combinations of and . The coupled system can be described in matrix form as:

where is the state variable at any time instant . The matrices and , and vector are given by:

Here, is the ratio of the combined lumped
capacitance of oscillator to the coupling capacitance ,
and is the ratio of the metallic state
resistance of device of oscillator, where
and . The fixed point in a conduction state is
given by and the matrix determining the flow
(the *flow matrix* or the *velocity matrix*) is given by
.

When two identical D-D oscillators are coupled using a parallel RC circuit with coupling resistance and coupling capacitance , they can lock in-phase or anti-phase depending on the relative values of and . In the extreme case of purely capacitive coupling with the anti-phase locking orbit is stable and the in-phase locking orbit is unstable. In case of purely resistive coupling with the in-phase locking orbit is stable and the anti-phase locking orbit is unstable. For other values, the system always have in-phase locking periodic orbit as well as anti-phase periodic locking orbits, with the stable locking being the in-phase locking when the coupling is close to purely resistive, and when the coupling is close to purely capacitive the stable locking is the anti-phase locking. Interestingly, considering the parameter space of and there exist a region with bistable orbits, i.e. both the in-phase and anti-phase orbits are stable. In this region, the steady state locking depends on the initial starting voltages at of the oscillators.

#### D-R oscillators

From a computing application point of view pairwise coupled D-R oscillators have interesting applications. When the series resistances are replaced by transistors as shown in figure 3 and the coupling is a simple capacitive coupling, a pair of coupled D-R oscillators can be used as a analog comparator whose output has the form of a difference norm[26, 27, 16]. Typical steady state orbits of the coupled system plotted in a plane are shown in figure 4a. When increases, the steady state orbits of the oscillators gets deformed. Such deformation of the steady sate orbits can be measured using a simple averaged thresholding-and-XOR operation on the steady state outputs of the oscillators. This averaged XOR measure is defined as first thresholding the output to binary values, second applying XOR operation on these binary values at every time instant and finally averaging this XOR output over some time duration. The averaged XOR output as a function of and is shown in figure 4b. The XOR surface reaches minimum value along the line . Within the locking range, it rises as an even function of resembling a difference norm. Outside the locking range, it averages to about 0.5. These characteristics of the curve can be explained by realizing that the averaged XOR measure by construction is equal to the fraction of the time the system spends in the grey region (region where XOR output is 1—determined by the thresholds on and ) in steady state as shown in figure 4a. It can be seen that the XOR measure should have least value in the symmetric case when the system locks out-of-phase and should increase as values diverge. Such a system can be used as an analog comparator with output as a difference norm. Arrays of such comparators can be used for template matching applications where element-wise comparisons suffice to decide a match.

## 4 IMT Oscillators

### 4.1 IMT Basics

Vanadium dioxide undergoes a first order metal-insulator phase transition marked by an abrupt change in conductivity up to five orders in magnitude. This abrupt change in conductivity which is also accompanied by changes in optical properties and can be triggered externally using thermal [28], optical [29], electronic [30] or strain [31] stimuli can pave the way for novel electronic devices including steep slope switching devices, memory elements and ultra-fast optical switches. One of the characteristic features of the electrically driven first order metal-insulator transition in VO2 is its inherently hysteretic nature as shown in figure 5(a). In previous publications, we demonstrated an electric field driven non-hysteretic phase transition in VO2 and showed novel device functionalities like coupled oscillations that may enable efficient implementation of novel, non-Boolean computing models. The electrically driven first-order phase transition in results in abrupt switching in conductivity but always comes at the cost of an intrinsic hysteresis, because the electrical field at which insulator to metal (IMT) transition occurs is always higher than that at which the metal to insulator transition happens (MIT). HIN our experiments, the VO2 is epitaxially grown on (001) TiO2 (-0.86% compressive strain) using Molecular Beam Epitaxy (MBE), then patterned to form channels and followed by deposition of Au/Pd contacts to electrically access the VO2 channel (details of growth and fabrication are given the supplementary section). First, the device is electrically driven across the phase transition boundary with zero external series resistance. A current compliance is set to limit the current in the metallic state (to prevent excess joule heating resulting in permanent damage). The insulator to metal transition (IMT) and the reverse metal to insulator transition (MIT) occur at two critical fields, and , respectively. The critical field is the field required to attain the Mott criteria and thereby trigger the formation of the metallic phase. Since transport in the insulating phase is dominated by hopping transport, we use the field dependent hopping conductance to understand its electrical properties. . The Negative Differential Resistance (NDR) region, where the VO2 is characterized by a conductivity intermediate between the metallic and insulating states is referred to as the phase coexistence region. The net conductivity in this region is due to contributions from the metallic and insulating phase4. In-situ nano XRD characterization performed simultaneously with the transient waveform measurement confirms that the nature of the insulating phase is Monoclinic M1 and that of the metallic phase is Rutile which is expected as the films are -0.86% compressively strained. When such a phase transition device is operated in the hystertic region, it breaks into spontaneous oscillations as shown in figure 5(b). This is a relaxation oscillator with piecewise linear dynamics as has been discussed in the previous section.

### 4.2 Coupling IMT Oscillators: Electrical and Optical

In our experiments and model development, we focus on electrical coupling using resistors and capacitors. The coupling dynamics has been discussed in the previous section and in more details in our earlier publications [25, 27]. For example, two oscillators with slightly different native frequencies when coupled show frequency locking as observed in experiments and simulations 3(c). A major challenge in scaling up chip-scale coupled oscillator circuits to solve practical computational problems lies in the difficulty in implementing an all-to-all oscillator coupling schemes. Recently, it has been proposed that cavity field coupled oscillators can overcome this challenge [32]. Lipson et. al. [33] has demonstrated experimentally the synchronization of a pair of silicon nitride optomechanical oscillators (OMOs) that are optically coupled through the radiation pressure field as opposed to mechanically coupled. The researchers were able to turn the optical coupling on and off using a heating laser via thermooptic effect. The ability to manipulate the coupling strength using well known nanophotonics techniques such as colossal electro-optic or thermo-optic effects, increase their future potential to realize large scale on-chip nonlinear dynamical systems. This is an area of research which lies unexplored from a computational point of view, and success in creating complex systems with many to many connections will be key to achieving hardware platforms capable to truly delivering the promise of coupled dynamical systems as computational elements.

## 5 Oscillator Networks and Applications

### 5.1 Image Data Processing and Analytics

Arrays of such IMT oscillator based comparators can be used for template matching applications where element-wise comparisons suffice to decide a match. Figure 6(a,b) illustrate the XOR-ed output of the phases of two coupled oscillators configured as figure 3. We observe from figure 6(c,d) a close match between experiments and simulations and it demonstrates that the XOR measure between the outputs of the two oscillators is a measure of distance between two inputs as given by .

We investigate the application of pairwise coupled oscillators for visual saliency approximations (detecting parts of the image that visually standout). Oscillator-based edge detection is performed using an array of pairwise oscillators to approximate the degree of dissimilarity between a given image pixel and its immediate neighbors. Different edges, vertical, horizontal, diagonal, are detected based on the selection of neighboring pixels for comparison. As this concept is expanded to include the comparison of pixels within a larger neighborhood (pixels surrounding reference pixel; a 3x3 neighborhood is used here), the output approximates the visual saliency. We also note that pairwise coupling of oscillators lead to the inherent calculation of fractional norms between templates and inputs, a task which is notoriously painstaking on a digitally abstracted computer. It is evident from figure 6(e) that the coupled oscillators show higher sensitivity to image contrast in comparison to a CMOS ASIC accelerator that uses a linear norm.

We further demonstrate this by comparing images of faces and hand-written numbers as shown in figure 6(f). We first use the XOR measure for each pixel and calculate the number of pixels with XOR output below a threshold value . Figure 6(f) shows the results of comparing faces with a relaxation comparator, where the grey shade corresponds to the fraction of pixels with positive match, white being the highest. Such system followed by a winner-take-all (WTA), i.e. a threshold on the number of pixels that give a positive match, can be used to decide if the input image matches a stored template pattern. The value of is chosen around 0.2 considering the minimum values of the XOR surface in the operating range of values. The two thresholds described above depend on different factors. The threshold of the number of pixels for WTA would depend on the database and the error statistics required or estimated. On the other hand, would be decided more by the nature of the XOR surface (Figure 10) and its minimum values. A coupled VO2-MOSFET configuration cascaded with a XOR provides a way of measuring a form of fractional distance using FSK. Such associative networks can be used in more complex pattern matching and classification problems with potentially large benefits in energy efficiency. The advantage of such oscillator based computing systems can be truly harnessed is they are miniaturized, made compact and integrated. We expect such systems to provide large improvements in energy-efficiency opening up possibilities in areas as varied as surveillance, consumer electronics and in-sensor processing. For comparison, a digital baseline design is designed and simulated. All digital circuits are implemented with 11nm node transistor models. We observe that the coupled oscillators provide a power reduction of 20X over CMOS reflecting the advantage of ‘let physics do the computing’ approach and potentially removing the Boolean bottleneck. For further details interested readers are pointed to the previous work by the authors [16].

### 5.2 Complex Global Connections and Possibilities for Computation

The full computational power of such coupled dynamical systems can be further harnessed using complex global interactions among oscillators instead of just pairwise interactions. Such coupled oscillators provide a Hopfield type networks, but with piecewise linear dynamics and hysteretic switching of oscillators. When connected in complex networks with both global and local connectivity, instead of pairwise coupling, these networks of oscillators can be much more powerful and can possibly compute approximate solutions of many hard optimization problems. In one implementation of globally coupled D-R oscillators (with a series resistance), we have recently demonstrated using theory and experimental implementations that such networks are capable of approximating the solution of NP-hard minimum graph coloring problems. When identical D-R oscillators are coupled using identical capacitors, such network settles to a steady state wherein the relative phases of the oscillators get ordered in a way that corresponds to the solution of graph coloring. The time evolution of the piecewise linear dynamical system of such coupled D-R oscillators inherits the properties of and hence mimics approximate graph coloring algorithms because of the construction of such networks. Basic functioning of such a system is shown in figure 7.

## 6 Concluding Remarks

In this article, we have demonstrated that the dynamics that evolve from the complex interactions among oscillators can be a powerful computing paradigm. However, this requires innovations in fabrication of compact and coupled networks of oscillators; and the current implementation using phase transition devices is one promising candidate. Even with advances in post-silicon devices and technology, the true potential of integrated dynamical systems can only be harnessed when controllable and programmable coupling can be realized and the phase and frequency dynamics carefully measured and read out. A plethora of challenges remain; however the opportunities offered by dynamical systems to make a significant impact in a post-CMOS world in undeniable.

### References

- Michael R. Garey and David S. Johnson. Computers and Intractability: A Guide to the Theory of NP-Completeness. W. H. Freeman & Co., New York, NY, USA, 1979.
- Stephen Wolfram, editor. Theory and applications of cellular automata: including selected papers, 1983-1986. Number v. 1 in Advanced series on complex systems. World Scientific, Singapore, 1986.
- P. Shor. Polynomial-Time Algorithms for Prime Factorization and Discrete Logarithms on a Quantum Computer. SIAM Journal on Computing, 26(5):1484–1509, October 1997.
- Andrew Lucas. Ising formulations of many NP problems. Interdisciplinary Physics, 2:5, 2014.
- Zhe Wang, Alireza Marandi, Kai Wen, Robert L. Byer, and Yoshihisa Yamamoto. Coherent Ising machine based on degenerate optical parametric oscillators. Physical Review A, 88(6):063853, December 2013.
- J. J. Hopfield and D. W. Tank. “Neural” computation of decisions in optimization problems. Biological Cybernetics, 52(3):141–152, July 1985.
- Hesham Mostafa, Lorenz K. Müller, and Giacomo Indiveri. An event-based architecture for solving constraint satisfaction problems. Nature Communications, 6:8941, December 2015.
- F. L. Traversa, C. Ramella, F. Bonani, and M. Di Ventra. Memcomputing NP-complete problems in polynomial time using polynomial resources and collective states. Science Advances, 1(6):e1500031–e1500031, July 2015.
- L. O. Chua and L. Yang. Cellular neural networks: applications. IEEE Transactions on Circuits and Systems, 35(10):1273–1290, October 1988.
- L. O. Chua and L. Yang. Cellular neural networks: theory. IEEE Transactions on Circuits and Systems, 35(10):1257–1272, October 1988.
- Tommaso Toffoli. Cellular automata as an alternative to (rather than an approximation of) differential equations in modeling physics. Physica D: Nonlinear Phenomena, 10(1):117–127, January 1984.
- Hava T. Siegelmann. Computation Beyond the Turing Limit. Science, 268(5210):545–548, April 1995.
- Anastasios Vergis, Kenneth Steiglitz, and Bradley Dickinson. The complexity of analog computation. Mathematics and Computers in Simulation, 28(2):91–113, April 1986.
- Masatoshi Imada, Atsushi Fujimori, and Yoshinori Tokura. Metal-insulator transitions. Reviews of Modern Physics, 70(4):1039–1263, October 1998.
- Nikhil Shukla, Abhinav Parihar, Eugene Freeman, Hanjong Paik, Greg Stone, Vijaykrishnan Narayanan, Haidan Wen, Zhonghou Cai, Venkatraman Gopalan, Roman Engel-Herbert, Darrell G. Schlom, Arijit Raychowdhury, and Suman Datta. Synchronized charge oscillations in correlated electron systems. Scientific Reports, 4, May 2014.
- N. Shukla, A. Parihar, M. Cotter, M. Barth, X. Li, N. Chandramoorthy, H. Paik, D.G. Schlom, V. Narayanan, A. Raychowdhury, and S. Datta. Pairwise coupled hybrid vanadium dioxide-MOSFET (HVFET) oscillators for non-boolean associative computing. In Electron Devices Meeting (IEDM), 2014 IEEE International, pages 28.7.1–28.7.4, December 2014.
- Steven H. Strogatz. From Kuramoto to Crawford: exploring the onset of synchronization in populations of coupled oscillators. Physica D: Nonlinear Phenomena, 143(1):1–20, 2000. read.
- Balth Van Der Pol. The Nonlinear Theory of Electric Oscillations. Proceedings of the Institute of Radio Engineers, 22(9):1051–1086, 1934.
- R. H. Rand and P. J. Holmes. Bifurcation of periodic motions in two weakly coupled van der Pol oscillators. International Journal of Non-Linear Mechanics, 15(4–5):387–399, 1980.
- D. W. Storti and R. H. Rand. Dynamics of two strongly coupled van der pol oscillators. International Journal of Non-Linear Mechanics, 17(3):143–152, 1982.
- Nancy Kopell and David Somers. Anti-phase solutions in relaxation oscillators coupled through excitatory interactions. Journal of Mathematical Biology, 33(3):261–280, December 1995.
- Norimichi Hirano and Sławomir Rybicki. Existence of limit cycles for coupled van der Pol equations. Journal of Differential Equations, 195(1):194–209, November 2003.
- K. Yogendra, D. Fan, and K. Roy. Coupled spin torque nano oscillators for low power neural computation. IEEE Transactions on Magnetics, 51(10):1–9, Oct 2015.
- Abhronil Sengupta, Priyadarshini Panda, Parami Wijesinghe, Yusung Kim, and Kaushik Roy. Magnetic tunnel junction mimics stochastic cortical spiking neurons. Scientific Reports, 6:30039 EP –, 07 2016.
- Abhinav Parihar, Nikhil Shukla, Suman Datta, and Arijit Raychowdhury. Synchronization of pairwise-coupled, identical, relaxation oscillators based on metal-insulator phase transition devices: A model study. Journal of Applied Physics, 117(5):054902, February 2015.
- Suman Datta, Nikhil Shukla, Matthew Cotter, Abhinav Parihar, and Arijit Raychowdhury. Neuro Inspired Computing with Coupled Relaxation Oscillators. In Proceedings of the The 51st Annual Design Automation Conference on Design Automation Conference, DAC ’14, pages 74:1–74:6, New York, NY, USA, 2014. ACM.
- A. Parihar, N. Shukla, S. Datta, and A. Raychowdhury. Exploiting Synchronization Properties of Correlated Electron Devices in a Non-Boolean Computing Fabric for Template Matching. IEEE Journal on Emerging and Selected Topics in Circuits and Systems, PP(99):1–10, 2014.
- M. M. Qazilbash, M. Brehm, Byung-Gyu Chae, P.-C. Ho, G. O. Andreev, Bong-Jun Kim, Sun Jin Yun, A. V. Balatsky, M. B. Maple, F. Keilmann, Hyun-Tak Kim, and D. N. Basov. Mott transition in vo2 revealed by infrared spectroscopy and nano-imaging. Science, 318(5857):1750–1753, 2007.
- A. Cavalleri, Cs. Tóth, C. W. Siders, J. A. Squier, F. Ráksi, P. Forget, and J. C. Kieffer. Femtosecond structural dynamics in during an ultrafast solid-solid phase transition. Phys. Rev. Lett., 87:237401, Nov 2001.
- Bong-Jun Kim, Yong Wook Lee, Sungyeoul Choi, Jung-Wook Lim, Sun Jin Yun, Hyun-Tak Kim, Tae-Ju Shin, and Hwa-Sick Yun. Micrometer x-ray diffraction study of films: Separation between metal-insulator transition and structural phase transition. Phys. Rev. B, 77:235401, Jun 2008.
- CaoJ., ErtekinE., SrinivasanV., FanW., HuangS., ZhengH., YimJ. W. L., KhanalD. R., OgletreeD. F., GrossmanJ. C., and WuJ. Strain engineering and one-dimensional organization of metal-insulator domains in single-crystal vanadium dioxide beams. Nat Nano, 4(11):732–737, 11 2009.
- C. A. Holmes, C. P. Meaney, and G. J. Milburn. Synchronization of many nanomechanical resonators coupled via a common cavity field. Phys. Rev. E, 85:066203, Jun 2012.
- Mian Zhang, Gustavo S. Wiederhecker, Sasikanth Manipatruni, Arthur Barnard, Paul McEuen, and Michal Lipson. Synchronization of micromechanical oscillators using light. Phys. Rev. Lett., 109:233906, Dec 2012.