# Interaction of Self-Sustained Pulses in Tunnel-Diode Oscillator Lattices

Authors: Koichi Narahara
Journal: Mathematical Problems in Engineering
 (2021)

---

## Abstract

A one-dimensional lattice in tunnel-diode (TD) oscillators supports self-sustained solitary pulses resulting from the balance between gain and attenuation. By applying the reduction theory to the device’s model equation, it is found that two relatively distant pulses moving in the lattice are mutually affected by a repulsive interaction. This property can be efficiently utilized in equalizing pulse positions to achieve jitter elimination. In particular, when two pulses rotate in a small, closed lattice, they separate evenly at the asymptotic limit. As a result, the lattice loop can provide an efficient platform to obtain low-phase-noise multiphase oscillatory signals. In this work, the interaction between two self-sustained pulses in a TD-oscillator lattice is examined, and the properties of interpulse interaction are validated by conducting several measurements using a test breadboarded lattice.


---

## Body

## 1. Introduction

In this paper, the dynamics of interacting self-sustained pulses developed in a one-dimensional lattice in tunnel-diode (TD) oscillators are characterized. For this purpose, an oscillator, which is driven by a single DC voltagevB that biases the TD, is employed. As a result, the TD exhibits gain by its own negative differential resistance. Using bifurcational and experimental characterization, it is found that a self-sustained pulse in a TD-oscillator lattice has a rather long tail, and its velocity significantly depends on vB [1]. As a result, the pulse propagates faster and becomes larger as vB increases. Furthermore, the pulse tail exhibits a small oscillation if the inductors’ parasitic resistances are negligible. In contrast, the pulse oscillatory tail becomes exponentially damped for finite parasitic resistances. The measured solitary pulse is allowed to develop for a restricted range of vB values around the TD’s peak voltage. Below the lower bound of this range, the input pulse is monotonically attenuated and finally disappears. In contrast, above the upper bound of this range, the lattice exhibits synchronous oscillation, and the solitary pulse is replaced by phase waves.The long tail of the solitary pulse can affect the propagation properties of the subsequent pulse because the pulse velocity significantly depends on the baseline voltage. To examine the type of interaction between two successive pulses, the reduction theory [2] is applied to the governing equations of a TD-oscillator lattice. It is found that, for no parasitic resistances, this interaction becomes either repulsive or attractive, depending on the separation distance between the pulses. However, for finite parasitic resistances, the separation distance tends to approach a specific value. Similar wave interaction has been observed in liquid film coating a vertical fiber [3, 4], e.g., where a drop-like wave interacts with a nearby one. The existence of bound states has been clearly described using the generalized Kuramoto–Sivashinsky equation [5].The interacting self-sustained pulses can be utilized in several signal processing applications. Let us consider an infinite pulse train inserted into a straight lattice. When a repulsive interaction dominates the pulse dynamics, the separation between any two adjacent pulses is asymptotically equalized. In this case, the TD lattice can reduce the timing jitter, which is initially present in the pulse train. In addition, two rotary self-sustained pulses move in the lattice loop in such a way that one occupies the other’s antipodal point. This can be achieved by setting half of the loop size to be less than the asymptotically converging separation value. This mechanism is expected to be applied to phase-noise reduction. In this study, the pulse properties in the lattice are examined using TDs. Furthermore, validation was carried out with flexible low-frequency experiments. However, once resonant-tunneling diodes (RTDs) [6–9] are implemented in place of TDs, the device potentially achieves higher operation frequency than CMOS circuits.In Section2, the dynamics of relatively distant interacting nonlinear pulses are characterized in the reduction theory framework. A reduced model representing the interaction of self-sustained pulses in a TD-oscillator lattice is initially developed, after describing the fundamental properties of the lattice. The experimental results obtained using a test breadboarded lattice are presented in Section 3. In Section 4, the pulse dynamics in a lattice loop are predicted using the reduced model, and the equalization of interpulse spacing is achieved using the 4th-order Runge–Kutta calculation.


## 2. Reduced Description of the Pulse Interaction

The unit-cell structure of the lattice under investigation is shown in Figure1. The lossy inductor represented by Lsh and Rsh is connected in parallel with a capacitance C. The nonlinear current source ID represents the TD. The Lsh and C pair establishes resonance resulting in voltage oscillation. By the application of the DC bias voltage vB the TD’s negative differential resistance compensates for the oscillation damping caused by Rsh. We consider the dynamics of self-sustained solitary waves developed in an inductively coupled system of unit-cell oscillators. The coupling is established by the lossy inductor represented by Lse and Rse. In this work, the Lse and Lsh values are commonly set as Lse=Lsh≡L (correspondingly, Rse=Rsh≡R). The resulting one-dimensional lattice supports a self-sustained pulse when the vB value is set in the restricted range around the TD’s peak voltage [1]. The near end of the lattice is connected to a pulse source to excite the self-sustained solitary pulses and the far end is connected to a load.Figure 1
A TD-oscillator lattice under investigation. The input voltage source excites a self-sustained solitary pulse in the TD-oscillator lattice. The resistanceRL is a load representing the subsequent circuit.The continuous spatial variablex is introduced instead of the cell address. The number of cells included in the unit length is defined as Nc. Then, in the long-wavelength approximation, the governing equations of the solitary pulses in the lattice are given as follows:(1)NcL∂i∂t+NcRi=−∂v∂x,(2)NcC∂v∂t+NcIDv+ish=−∂i∂x,(3)LshNc∂ish∂t+RNcish=vx−vB,where v, i, and ish represent the terminal voltage of ID, the series current, and the shunt current at x, respectively (see appendix for explicit derivations). Using an auxiliary variable w, equations (1)–(3) are simplified as follows:(4)∂∂twv=1LC−CR+LdIDdvw+1Nc2∂2v∂x2−v+vB−RIDw.The trigger wave [1, 10, 11] or dissipative solitons [3, 5, 12] including the pulse under investigation have no velocity flexibility. Then, according to the standard description, the solutions of equation (4) are assumed to be as follows:(5)w=Wx−ct−x1+Wx−ct−x2+bwx−ct,t,(6)v=vB+Vx−ct−x1+Vcx−ct−x2+bvx−ct,t,where V and W are the voltage and its temporal derivative of the self-sustained pulse, bv,w represents the small overlap function, and c is the pulse velocity. In the following, any function Xx−ct−xi is abbreviated as Xi (for i = 1, 2) for brevity.By substituting equations (5) and (6) into equation (4), the first-order dynamic equations for bw,v are derived as follows:(7)∂b∂t=Mb+g,where b=bw,bvT and the 2 × 2 operator M is defined as(8)M=−CR−LGvB+V1+cCL∂xCL−1−RGvB+V1−LKvB+V1W1+Nc−2∂x2CL1c∂x,where G=ID′ and K=G′ (′ represents the derivative with respect to the argument). In addition, the inhomogeneous term g is given as follows:(9)g=g1g2,g1=−RIDvB+V12−V12−CRW12−LGvB+V12W12+cCL∂xW12+Nc−2∂x2V12CL,+x˙1∂xW1+x˙2∂xW2,g2=W12+c+x˙1∂xV1+c+x˙2∂xV2,where X12 abbreviates X1+X2 for X=V,W. The second pulse is assumed to be far apart from the first. Then, the contributions of the second pulse in M can be neglected.The solvability condition is given byΨ,g=0 [2], where Ψ=ψ1,ψ2T is the zero eigenvector of the adjoint of M.(10)M†Ψ=−CR−LGvB+V1−cCL∂xCL1−1−RGvB+V1−LKvB+V1W1+Nc−2∂x2CL−c∂xψ1ψ2=0.The resulting equations are as follows:(11)ψ2=RL+GV1Cψ1+c∂xψ1,(12)cCL∂xψ2=−1−RGvB+V1−LKvB+V1W1ψ1+Nc−2∂x2ψ1.First, the solutionsW1 and V1 are numerically obtained by solving equation (4). These solutions are spatially shifted to obtain W2 and V2. Then, the zero adjoint eigenvector Ψ1=ψ11,ψ21T is obtained by solving(13)∂xϕψ11=cCR+LGvB+V1ϕ+1+RGvB+V1ψ11Nc−2−c2CLϕ.It is obtained by substituting equation (11) into (12) and introducing the auxiliary function ϕ.The solutionsV, W, and Ψ1 are used to evaluate the solvability condition Ψ,g=0. The dynamical equations of the pulse positions x1,2 are given as follows:(14)γ0x˙1+γ−lx˙2+Fl=0,(15)γlx˙1+γ0x˙2+F−l=0,where l≡x1−x2 and(16)γ0=∫−∞∞dxψ1x∂xWx+ψ2x∂xVx,(17)γl=∫−∞∞dxψ1x+l∂xWx+ψ2x+l∂xVx,(18)Fl=−1LC∫−∞∞dxψ1xRIDvB+Vx+Vx+l,−IDvB+Vx−IDvB+Vx+l,+LGvB+Vx+Vx+l−GvB+VxWx,+LGvB+Vx+Vx+l−GvB+Vx+lWx+l.Finally, equations (16) and (17) are reformulated as follows:(19)l˙=−γ0+γlFl−γ0+γ−lF−lγ02−γlγ−l≡Zl.Zl can be evaluated by setting Nc, L, C, and R to 30, 10 μ H, 500 pF, and 5 Ω, respectively. Furthermore, the following phenomenological model was employed for ID to simulate the TD’s fundamental properties including the presence of the peak and valley voltages and the relatively small negative differential resistance:(20)IDV=0⋅0050⋅175log1+exp−0⋅621+54V1+exp−0⋅621−54V×1⋅4+tan−11⋅85−19V+0.0018exp11.6V−1−0⋅187,where the numerical values were optimized to simulate the NEC 1S1763 Esaki diode.In Figure2Zl is plotted for l∈−2,2. Obviously, Zl is antisymmetric with respect to the vertical zero-crossing axis, and two stable equilibrium points P1,2 exist, corresponding to l=±0.2275≡l∗. If l is slightly smaller than P2, a positive l˙ means that the two pulses interact repulsively, whereas when l becomes larger than P2, a consequent attractive interaction diminishes l. Therefore, the two successive solitary pulses tend to move in the lattice with a fixed separation l∗.Figure 2
Reduced dynamics of two pulses. Sincel˙ is measured on the vertical axis, the pulses asymptotically arrive at either P1 or P2, regardless of their initial position.To examine this property in a practical discrete lattice, the original lattice equations given in equations (A.1)–(A.3) were solved in the time domain by using the 4th-order Runge–Kutta method. The results are presented in Figure 3. Equation (22) was used for calculating ID. A 200-cell lattice was calculated using L = 4 μ H and C = 500 pF. Two solitary pulses were inserted into the near end of the lattice, and the variation of their temporal separation between peaks was monitored. Initially, R was set to 4 Ω. The calculated results are shown in Figures 3(a) and 3(c). The waveform of the self-sustained pulse is shown in Figure 3(a). A long and exponential tail can be observed. Figure 3(c) shows the spatial variation of the pulse separation for nine different values of initial temporal separation, starting from 0.5 μ s and ending at 1.2 μ s with increments of 0.1 μ s. For initial separations of less than 1.0 μ s, the spatial separation increases with the moving distance. This is consistent with the reduced model prediction, which states that any two pulses with small separation are subject to repulsive interaction. On the other hand, Figures 3(b) and 3(d) correspond to the case of R being zero. The temporal waveform in Figure 3(b) exhibits an oscillatory tail, and correspondingly, several different converging separations exist, as shown in Figure 3(d). An effort was made to characterize this loss-free lattice using the reduced model. However, equation (15) fails to converge to the required solution.Figure 3
Interpulse separation predicted by Runge–Kutta analyses. The unit cell contains (a) 4-Ω resistors and (b) no resistors.
(a)(b)(c)(d)Note that no plots are shown for initial separations of less than 0.5μ s in either case. When the spacing becomes too small, the subsequent pulse is affected by the voltage drop caused by the tail of the leading pulse and decays and disappears gradually.


## 3. Experimental Results

The nontrivial motion of two solitary pulses was examined by fabricating a test lattice on breadboards. TheL and C values were set to 5.0 μ H and 2.0 nF, respectively. No resistors were used for R. Instead, the parasitic L resistance values were used. An NEC 1S1763 Esaki diode was used to implement the TDs. The total number of connected cells NL was 50. A Keysight 1134B active probe was used to detect the time domain voltage variations in all cells of the lattice using an Agilent DS90254A oscilloscope.The measured variation of the pulse separation is shown in Figure4. The voltages were measured at nine different points (from the 1st cell to the 46th cell with increments of 5 cells) for an initial temporal separation between pulses varying from 0.0 to 6.0 μ s with increments of 0.055 μ s. The dependence of temporal waveforms based on the initially assigned temporal delays was measured in the first and 26th cells. The results are plotted in Figures 4(a) and 4(b), respectively. At the point of incidence, the leading pulse occupies a common temporal position, whereas the subsequent pulse delays linearly from the bottom (no separation) to the top (6.0-μ s separation). In contrast, in the 26th cell, the second pulse disappears for an initial delay of less than 0.5 μ s. The pulse delay exhibits only a small variation for an initial delay greater than 2.2 μ s. For other initial delay values, the second pulse survives, and its delay tends to converge to a common value (Figure 4(c)). Both the deceleration of the second pulse and the acceleration of the leading pulse contribute to a saturating delay. This observation indicates that mutual interaction is established between the pulses.Figure 4
The measured temporal evolution of interpulse temporal separation. The temporal waveforms for several different initial delays detected in (a) the first cell and (b) the 26th cell. (c) An enlarged view of the waveforms detected in the 26th cell for small interpulse delays.
(a)(b)(c)The pulse dynamics are shown in Figure5. This figure illustrates the variation of the pulse-to-peak separation, where the initial pulse separation ranges in 0.71–1.67 μ s with increments of 0.055 μ s. Beyond the 31st cell, the pulses with initial separation of less than 1.2 μ s exhibit a similar separation trend. However, the separation convergence value is not kept fixed. Instead, it keeps increasing with the distance traveled. Moreover, beyond the 46th cell, the second pulse gradually decays and finally disappears. This insufficiency may be due to fluctuation in device parameter values and breadboard’s bandwidth. A complete validation requires a thorough optimization in the selected parameter values and platform.Figure 5
Spatiotemporal variation of pulse separation.To simulate the measured results several Runge–Kutta calculations were carried out. The device parameter values were fluctuated with a 10% standard deviation. The average values ofL, C, and R were set to 5 μ H, 2 nF, and 1 Ω, respectively. Triangular pulses were applied to the near end with a 50 Ω input resistor, and the far end was terminated with a 50 Ω load resistor. The incident pulse had 100 ns fall and rise times and 0.5 V amplitude. Figure 6(a) compares the measured and calculated waveforms. The waveforms commonly have the short positive main pulse with a negative long tail. However, the measured pulse does not return to the zero level unlike the calculated one as indicated by arrows. This negative offset of the measured pulse may be caused by some technical problem. For example, the main pulse was steep sufficiently to exceed the breadboard’s bandwidth for reliable operation. Thus, we consider the numerical model under investigation as being suitable to simulate a well-composed TD-oscillator lattice. Figure 6(b) shows the calculated spatial waveform of the self-sustained pulse moving to the right. The pulse’s spatial width is estimated to be 17 cells with the 5-cell main pulse and 12-cell tail.Figure 6
Properties of calculated pulse. (a) The comparison between the measured and calculated waveforms and (b) the spatial waveform of the calculated pulse.
(a)(b)To examine the double-pulse response two triangular pulses are input to the near end of a 500-cell lattice. The initial pulse separationts was set to 0.63 μ s (the second pulse disappeared for ts≤ 0.62 μ s). The spatiotemporal behavior of two pulses is shown in Figure 7(a). The separation between the first and second pulses increases monotonically with decreasing growth rate (its spatial variation is shown in Figure 7(b)). The presence of the first pulse tail decelerates the second pulse. The almost asymptotic spatial waveform is shown in Figure 7(c), which was detected at the time indicated by the dashed line in Figure 7(a). The spatial width is estimated to be 34 cells, which is just twice as wide as that of the single pulse. Thus, the separation tends to increase until the pulses’ overlap is completely eliminated. The lattice size must be set more than 34 cells for the coexistence of two independent pulses.Figure 7
Calculated lattice’s response to the double-pulse incidence. (a) The spatiotemporal voltage map, (b) the spatial variation of the pulse separation, and (c) typical spatial waveform.
(a)(b)(c)Figure8 shows the calculated results when three evenly separated pulses are input to the 500-cell lattice (ts = 0.63 μs). The spatiotemporal voltage map shown in Figure 8(a) shows that three pulses move in the lattice and reach the far end without being eliminated on the way. Furthermore, the pulses 1 and 2 become more distant and the separation growth rate is independent of the presence of pulse 3. This is clearly observed in Figure 8(b), where the spatial variation of the interpulse separations is plotted. The red curves in Figures 7(b) and 8(b) are almost coincident. The separation between pulses 2 and 3 (blue curve) also has the increasing tendency. However, it initially decreases slightly and is always smaller than its counterpart up to the far end. The discrepancy in separations may start to decrease only after the overlap between pulses 1 and 2 becomes negligibly small.Figure 8
Calculated lattice’s response to the triple-pulse incidence. (a) The spatiotemporal voltage map and (b) the spatial variation of the pulse separation (red: the separation between pulses 1 and 2, blue: the separation between pulses 2 and 3).
(a)(b)Figure 9
Calculated lattice’s response to the multiple-pulse incidence. (a) The spatiotemporal voltage map, (b) the variation of the pulse separation with respect to the space and separation ID, and (c) the spatial waveform of the steadily moving periodic pulse train. The separation between thei th and (i + 1)th pulses has the separation ID of i in Figure 9(b).
(a)(b)(c)To clarify how the number of advancing pulses affects the evolution of the pulse separation we examined the lattice response to a pulse train. The lattice size was reduced to 400 for convenience. As a result, the device fluctuation differs from that used in the 500-cell lattice. A train of evenly separated triangular pulses (ts = 1.0 μ s) was input to the near end. Every pulse reaches the far end as shown in Figure 9(a). For brevity, the temporal separation between the m th and (m + 1)th pulses is simply denoted as the m th separation. Figure 9(b) shows the spatial dependence of the m th separation 1≤m≤31. At any position in the lattice the (m + 1)th separation becomes smaller than the m th one for m<20. Furthermore, the separation preserves its initial value even at the far end for m≥20. As a whole, the increasing separation is caused by the presence of the 1st pulse, i.e., the temporal edge of the pulse train. This edge effect gradually percolates through to the interior pulses at the rate of 20 cells/pulse (twenty pulses from the edge is subject to the edge effect by 400-cell propagation in the lattice). Without such an edge an input pulse train preserves its interpulse separation even in the lattice. This may be realized through the balance in multiple-pulse interactions. Figure 9(c) shows the spatial waveform of the pulse train preserving initial separation. The spatial period is estimated to be 11 cells, being much smaller than the width of the single pulse. It may be reduced to be in the neighborhood of the minimum value below which the following pulse disappears. It may also be expected that the pulse interaction equalizes the separation of each neighboring pulse pair. In the following, we consider two pulses in a lattice loop, which simulates the steadily moving pulse train in a straight lattice, to examine whether the equalization of pulse separation is feasible or not.


## 4. Application

An application utilizing pulse interaction is the phase-noise reduction. The case where two or more solitary pulses rotate in a lattice loop is considered. By implementing the multiple output ports, the device can operate as a multiphase oscillator [13–17], which contributes to the relaxation of the base frequency to drive the ultrafast systems operating in submillimeter wave or terahertz frequencies using time-interleaving [18] or output combining [19].The oscillation frequency increases with the number of rotary pulses. However, the fluctuation in the position of each pulse may result in degradation in the oscillator’s phase-noise performance. The repulsive interaction observed can force rotary pulses to be equally spaced, contributing to phase-noise reduction.In a loop with sizeL0, the same equation for l2=L0−l holds as l in equation (19). Thus the following equation is obtained:(21)−l˙=−γ0+γL0−lFL0−l−γ0+γl−L0Fl−L0γ02−γL0−lγl−L0.By subtracting equation (21) from equation (19), l satisfies the following equation:(22)l˙=12−γ0+γlFl−γ0+γ−lF−lγ02−γlγ−l+γ0+γL0−lFL0−l−γ0+γl−L0Fl−L0γ02−γL0−lγl−L0≡ZL0l.Notice that the equilibrium point of equation (22) occurs for l=±L0/2. To quantify the pulse dynamics, the same lattice parameter values and ID model presented in Figure 2 were employed. The results shown in Figure 10(a) correspond to the case of L0 being set to 1.0. Two stable equilibrium points P1,2 and an unstable one Q corresponding to the half-loop are observed. Since the half-loop size exceeds l∗, two pulses tend to move in the loop with an l∗ separation as that in the straight lattice. On the other hand, the case of L0 being set to 0.5 is shown in Figure 10(b), where Q becomes stable. If the loop size is designed to be smaller than l∗ a novel strategy for reducing phase noise in coupled oscillator systems may emerge.Figure 10
Predictions of pulse reduction in a loop. The loop size is set to (a) 1.0 and (b) 0.5.
(a)(b)The equalization of pulse separation in a TD lattice loop was demonstrated by carrying out a Runge–Kutta calculation for the structure shown in Figure11(a). The L, R, and C values were set to 10.0 μ H, 5 Ω, and 2.0 nF, respectively. Two 10-cell straight lattices were prepared for the signal input (orange) and output (green) signals. In addition, two coupled TD lattices (yellow) were implemented to introduce and extract self-sustained pulses in the loop. One of the coupled lattices was connected to the input/output lattices, and the other partially formed the lattice loop with the supplementary 26-cell lattice (blue). The coupling capacitance Ccpl was set to 470 pF. When two pulses were inserted into the input lattice, they started to rotate in the loop in an anticlockwise fashion. The repeated rotation of these two pulses in the loop results in the pulse train at the OUT port. The separation variation of the two pulses with the rotation was observed. Figure 11(b) shows the trajectories of the two rotary pulses in the lattice loop. Time progresses horizontally, and the spatial position is measured on the vertical axis. The vertical axis bottom and top correspond to the first and 36th cells, respectively. Due to the closed structure, the pulse reaches the first cell after passing the 36th cell, so that the tilted stripes in Figure 11(b) represent one or the other rotary pulses, alternatively. Initially, the pulse separation set was rather unbalanced. At the beginning, the first pulse leads the second one, as shown by P1, which is relatively larger than P2=36−P1. After a six-time rotation, the separation becomes Q1, which is equal to Q2=36−Q1. Thus, the pulse separation is successfully and promptly equalized.Figure 11
Numerically obtained pulse dynamics in a tunnel-diode (TD) lattice loop. (a) Structure of the TD lattice loop and interface lattices and (b) calculated spatiotemporal voltage map.
(a)(b)To quantify the phase-noise reduction achieved using repulsive interactions, we obtained the perturbation projection vector (PPV) [20, 21] or phase sensitivity for the oscillation on the TD lattice loop. The lattice equations shown in eqs. (A.1)–(A.3) were solved for the TD-oscillator lattice loop with Nc=30. The L, C, and R values were set to 10 μ H, 2 nF, and 5 Ω, respectively. No parameter fluctuations were assumed. Equations (A.1)–(A.3) are conveniently written as follows:(23)dXdt=FX,where the first, second, and last Nc components of the 3 Nc vector X correspond to in, vn, and jsh,n, respectively. The components of F are given by the right-hand side of equations (A.1)–(A.3) with proper variable replacement. The phase-noise performance is significantly dependent on which noise source mainly influences it. In this study, the white noise voltage source was connected in series with Rse,sh in each cell; thus, the number of sources becomes 2 Nc. As a result, the governing equation is given as follows:(24)dXdt=FX+Bξ,where the 2 Nc vector ξ enumerates the white noise voltage sources and B represents the 3Nc×2Nc transfer matrix between ξ and X⋅. The deterministic part of equation (24) has a limit cycle solution denoted by Xs. We then linearized F for X=Xs and solved the resulting linear equation to obtain the 3Nc×3Nc monodromy matrix Φ. The PPV, denoted by Z, can be defined as the eigenvector of the adjoint of Φ corresponding to a zero eigenvalue with an appropriate normalization [20, 21]. The phase noise is characterized by a unique parameter c, defined as follows:(25)c=1T0∫0T0dtZ†tBXstB†XstZt,where T0 represents the period of limit cycle oscillation. For sufficiently small c, the phase noise of the oscillator can then be calculated using the following equation:(26)ℒfm=10log10f02cfm2,where f0 and fm are the oscillation and offset frequencies, respectively [21].First, the single, double, and triple rotary pulse solutions were obtained by solving equation (23) in the time domain by using the Runge–Kutta method. Figure 12(a) shows the properties of the single pulse rotation. The upper panel shows the spatiotemporal voltage map. A single pulse is shown to rotate in the loop to the direction the cell address increases. The rotation frequency was 241.6 kHz. The lower panel shows the locations of the eigenvalues of Φ. Because the rotary pulse is stable, all eigenvalues are included in the unit circle except for λ = 1 (indicated by black arrow). Figures 12(b) and 12(c) show the properties of the double and triple rotary pulse solutions, respectively. The mismatch in rotation direction does not matter because the lattice was homogeneous. The rotation frequency decreases as the pulse number increases: 235.1 and 215.1 kHz for the double and triple cases, respectively (for comparison, two- and three-pulse durations were employed as the oscillation period for the double and triple cases, respectively). The second largest eigenvalue is indicated by green arrows in each lower panel. As suggested in [22], when this eigenvalue approaches the unit circle more, the oscillator’s original Q factor increases. The presence of two or three pulses in the loop contributes to the enhanced Q factor.Figure 12
Eigenvalues of monodromy matrix. (a) Single, (b) double, and (c) triple rotary pulse(s) in the loop. In each figure, the spatiotemporal voltage map is shown at the top and the eigenvalues of the corresponding monodromy matrix are shown at the bottom. For clarity, the vertical and horizontal scales are shown only for the left voltage map, although the other two are shown in the common scale. Red curve shows the unit circle.
(a)(b)(c)To evaluate PPV the white noise source intensity for resistanceR was assumed to be 4kBTRΔf, where the temperature T and the noise frequency bandwidth Δf were set to 300 K and 100 kHz, respectively (kB: Boltzmann constant). The c value was calculated to be 3.53 ×10−13, 1.46 ×10−13, and 1.58 ×10−13 s for the single, double, and triple cases, respectively. By using equation (26), the phase noise is calculated to be −76.86, −80.92, and −81.37 dBc/Hz at 1 kHz offset for the single, double, and triple cases, respectively. By the presence of the second pulse, the phase noise was reduced by 4 dBc/Hz. Further investigations are required to optimize the impact of the interpulse repulsion on phase-noise reduction, including the loop size, device parameters, and bias voltage.As mentioned above the velocity of the self-sustained pulse significantly depends onvB [1]; therefore, the rotation frequency also exhibits significant dependence on vB. This means that the oscillator under investigation operates as a voltage-controlled oscillator, which can be efficiently implemented in a phase-locked loop (PLL), where the low-frequency phase noise inside the loopback bandwidth can be reduced to the same level as the reference oscillator. The interpulse repulsion promptly equalizes separation; therefore, it contributes to the reduction of the high-frequency phase noise. Furthermore, the present scheme does not disturb the PLL operation and has the synergetic potential to achieve broadband phase-noise reduction.


## 5. Conclusions

The interactions between two successive self-sustained pulses were characterized in a TD-oscillator lattice. The reduction theory, which converts the original circuit equations into the dynamic equations with respect to the pulse positions, predicts that the lossy TD-oscillator lattice has a unique stable point in the pulse separation, if half of the loop size becomes less than the asymptotically converging stable spacing. Furthermore, the antipodal points exhibit unique stable positions for two pulses rotating in the lattice loop. The results obtained from several measurements conducted using a test breadboarded lattice clearly demonstrated that the delay between the two comoving pulses tends to converge to a fixed value. In addition, the interpulse delay equalization was successfully calculated in the time domain. We believe that the lattice loop investigated provides an efficient platform for generating low-phase-noise multiphase oscillatory signals.



---

*Source: 5027127-2021-12-15.xml*