# Critical Delays in Clock Distribution Circuits

Authors: José R. C. Piqueira; Diego Paolo F. Correa
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

As the distribution of clock signals between the nodes of a network
became a critical operational requisite, the transmission delays have to
be studied because they affect the accurate recovering of the time basis.
In this work, the critical delay value is calculated considering simplifications compatible with practical situations. The main contribution is to
consider the dissipative terms in the node equations, expressing critical
delays, depending on the time constant of the dissipation.


---

## Body

## 1. Introduction

The phase control problem in a network of oscillators is a problem that appears in several engineering situations: world wide time signal distribution, synchronous communication networks, distributed computation, processes control, and synchronous operation of micro- and nanocircuits. Generally speaking, the distribution of clock signals can be performed controlling the phase of a node oscillation by using the information of the phases of the other nodes that compose the whole network.Modeling this kind of network presents two types of clock distribution strategies: master-slave (MS) and full-connected (FC), described by Lindsey et al., in a seminal paper that is an important reference for engineers working on design of networks [1].Concerning the networks conceived up to the end of the last century, Synchronous Digital Hierarchy (SDH), employing MS strategies, was the most important solution becoming a standard recommended by ITU (International Telecommunication Union) [2]. The mathematical analysis of this solution has been widely studied with a lot of different approaches: linear differential equations [3], nonlinear differential equations [4], and stochastic differential equations [5].Nowadays, it could be considered that clock signal distribution by using MS strategies is well studied and network designers have a lot of good references available for the several types of phase detection and signal filtering [6].However, the network architecture evolved to complex topologies with phase and frequency references depending, at each node, on the phase and frequency of almost all the other nodes [7] and, consequently, to obtain general solutions becomes feasible only for certain particular situations [8] as the resulting equations for the isolated nodes are second-order, including a dissipative term. Besides, the coupling terms between the nodes are composed of nonlinear delayed terms [1].Here, considering that operational frequencies of the node oscillators are normalized, the clock signal phase of each node is supposed to be described by [9]
(1)θ¨i(t)+1Tiθ˙i(t)=Ki[sin(ϕi(t)-θi(t))]+fi(t).
In (1), Ti represents a dissipative time constant term, and Ki represents a frequency node gain, both related to node i. The effect of weak forces is represented by fi(t) that considers noise and small perturbations. A wave propagation term appears as ϕi(t) and it is responsible for considering the dependence of the local phase on the phase of all nodes and is given by
(2)ϕi(t)=∑j=1n‍(ai,j)θj(t-τi,j).
In (2), ai,j represents the normalized connection matrix, that is, ∑j=1n‍(ai,j)=1, for each node, and τi,j is the propagation delay between the nodes i and j.The main goal of this paper is to study the several possible behaviors of a network, trying to express them to be related to the connection matrix, the nodes free-running frequencies, and the delays between the nodes. Considering the complexity of the problem, some simplifying assumptions are necessary in order to treat the problem.The first simplification is about the network structure that is supposed to be triangular [8] as the practical networks are trees, composed of cells that can be modeled as triangles [1, 10].  They  are supposed to have the same parameters. In spite of not being a general case, this assumption provides a simple reasoning with useful practical conclusions, as the circuit components are similar [10].In the next section, the equations for this type of configuration will be derived with the spatial phase errors being the measurable state variables. Then, in order to analyze the tracking mode of the system, a linear model is developed, resulting in two decoupled equations relating the spatial errors to the delayed spatial errors.Based on these equations, it is possible to calculate the critical delay that implies instability and to relate it to time constantTi of the dissipative terms. The calculation method for this critical value is presented, followed by numerical results and simulations confirming the conclusions.


## 2. Triangular PLL Network

Following the reasoning and notation exposed in the former section, a three-node PLL network can be modeled by the following equations:(3)θ¨(1)(t)+1T(1)θ˙(1)(t)=K(1)sin[a11θ(1)(t)+a12θ(2)(t-τ21)+a13θ(3)(t-τ31)-θ(1)(t)],θ¨(2)(t)+1T(2)θ˙(2)(t)=K(2)sin[a22θ(2)(t)+a21θ(1)(t-τ12)+a23θ(3)(t-τ32)-θ(2)(t)];θ¨(3)(t)+1T(3)θ˙(3)(t)=K(3)sin[a33θ(3)(t)+a31θ(1)(t-τ13)+a32θ(2)(t-τ23)-θ(3)(t)].At this point, it is assumed that the nodes have the same constitutive parameters, that is,T=T(1)=T(2)=T(3) and K=K(1)=K(2)=K(3). These assumptions are reasonable considering a practical communication point of view, as the nodes are composed of standard blocks [6]. Besides, due to the same motive, it is considered that all nodes have the same accuracy and, consequently, the phases of the three nodes are used to adjust the phase of the local node with the same weight, resulting that ai,j=1/3,∀(i,j) [1].In order to be more accurate from an engineering point of view, the delays between the nodes ought to be different. But in real modern networks, they assume low values, being fractions about10% of the period of the VCO oscillations, that is, 10% of the node free-running period. Therefore, here all the delays are considered to be equal, that is, τi,j=τ,∀(i≠j), and the results can be related to maximum permitted delays for real networks [2].Taking the former conditions into consideration, the equations of the phase nodes become(4)θ¨(1)(t)+1Tθ˙(1)(t)=Ksin[13[θ(2)(t-τ)+θ(3)(t-τ)]-23θ(1)(t)],θ¨(2)(t)+1Tθ˙(2)(t)=Ksin[13[θ(1)(t-τ)+θ(3)(t-τ)]-23θ(2)(t)]θ¨(3)(t)+1Tθ˙(3)(t)=Ksin[13[θ(1)(t-τ)+θ(2)(t-τ)]-23θ(3)(t)].To proceed with the analysis of the network, it must be considered that to obtain the node clocks in a correct sequence, the phase and frequency errors between the nodes have to vanish after any kind of perturbation, with as minimum overshoot and tracking time as possible. Consequently, the important state variables are the spatial phase (Ψi,j) and frequency (Ψ˙i,j) errors, defined as
(5)Ψi,j=θ(i)-θ(j);(6)Ψ˙i,j=θ˙(i)-θ˙(j).In this case, a three-node network, there are six index combinations(i,j) defining the spatial errors (Ψi,j) and (Ψ˙i,j). However, if (Ψ1,2) and (Ψ1,3), combined with (Ψ˙1,2) and (Ψ˙1,3), are chosen, all the other spatial errors can be expressed as linear combinations of them. Therefore, combining (4), the system can be described by
(7)Ψ¨1,2+1TΨ˙1,2=K{sin[13[θ(2)(t-τ)+θ(3)(t-τ)]-23θ(1)(t)]-sin[13[θ(1)(t-τ)+θ(3)(t-τ)]-23θ(2)(t)]},(8)Ψ¨1,3+1TΨ˙1,3=K{sin[13[θ(2)(t-τ)+θ(3)(t-τ)]-23θ(1)(t)]-sin[13[θ(1)(t-τ)+θ(2)(t-τ)]-23θ(3)(t)]}.


## 3. Analysis of the Tracking Mode

An important indicator of the performance of a clock distribution systems is  regarding  its behavior in the tracking mode; that is, how the system recovers the synchronism when a perturbation around the synchronous state occurs [3].In this case, the phase deviation can be considered to be small and the argument of thesin functions in (7) and (8) can replace the functions [1]. Therefore, the equations describing the system are
(9)Ψ¨1,2+1TΨ˙1,2+23Ψ1,2-13Ψ1,2(t-τ)=0;(10)Ψ¨1,3+1TΨ˙1,3+23Ψ1,3-13Ψ1,3(t-τ)=0.Consequently, the model for the tracking mode of a triangular PLL network is given by two decoupled equations, which are ordinary, homogeneous, and linear, containing a delay term. This equation admits an analytical solution. According to [11, 12], this type of delay differential equation corresponds to a transcendental characteristic equation, decomposable as
(11)∑j=0N‍ajλj+e-λτ∑j=0M‍bjλj=0,
which permits to calculate the critical stability value for the delay τ.


## 4. Critical Delay Value

The critical delay value for the three-node clock distribution network can be determined considering that the spatial errors are described by(12)Ψ¨(t)+1TΨ˙(t)+23Ψ-13Ψ(t-τ)=0,
with parameter T∈ℜ+ representing the time constant of the dissipative term and the bifurcation parameter τ∈ℜ+ representing the signal processing and propagation delay.In order to determine how the bifurcation parameterτ determines the system stability it is necessary, for a given T, to know the set Γ⊂ℜ+ such that, for a given τ=τc∈Γ, solutions for (12) change stability.For the dynamical system described by (12), the characteristic equation is given by
(13)P(λ;T,τ)=λ2+λT+23-13e-λτ=0
and can be written shortly as
(14)P(λ;T,τ)=a(λ;T)-e-λτ=0,
where a(λ;T)=3λ2+(3/T)λ+2.It is supposed that there existsτ=τc, called critical τ, that leads roots of the characteristic equation up to the imaginary axis, such that λ=±jω with (ω∈ℜ+), implying the possible onset of a Hopf bifurcation [13].Consequently, the module condition(15)|a(jω;T)|=|e-jωτc|=1
must hold. Therefore, considering only positive values of ω, from (15), it is possible to find
(16)ω=ω(T).The number of values forω depends on the degree of the polynomial a(λ;T). Considering only the positive and real ones, the roots ω have to also satisfy the angle condition for (14), given by
(17)tan-1[ℑ(a(λ;T))ℜ(a(λ;T))]λ=jω±2nπ=-ωτc,
for some n∈Z.Then, from (16) and (17), it is possible to find a relationship between a given T and its associated critical τc:
(18)τc(T)=-1ω(T)[tan-1(ℑ(a(λ;T))ℜ(a(λ;T)))λ=jω(T)±2nπ].By computing the module condition given by (15), for a(jω;T), there are four possible roots:
(19)ω(T)=±[16T2(4T2-3±4T4-24T2+9)]1/2,
and replacing (16) into (18) results in
(20)τc(T)=-1ω(T)[tan-1(3ω(T)T(2-3ω(T)2))±2nπ].Equation (20) shows that, for a given T, there are infinite periodical critical values for τ, with period:
(21)2nπω(T).It can be noticed that there are some situations with (19) having no positive roots, meaning that, for some values of T, there is no τc that leads the roots up to the imaginary axis. In this case, system stability will depend on the position of the roots for τ=0, for the given T.


## 5. Numerical Results and Critical Delays

For a system given by (12), it is possible, by using (19), to calculate which values of T allow critical time delays τc, responsible for bifurcations characterized by changing stability.Figure1 shows the real positive values of ω, given by (19). It can be seen that these roots suddenly appear for values of T greater than 3π/4. Analyzing characteristic equation (13), for τ=0, it can be seen that the system described by (12) is stable for every T∈ℜ+. Then, for 0<T<3π/4, the system described by (12) remains stable for every value of time delay.Figure 1
Real positive values forω(T) given in (19). These values appear for T>3π/4.However, considering a practical point of view, it is important to know how the critical delay values (τc) are related to the dissipative time constant (T). Equations (19) and (20) show that, for a set of values of T, there exist families of possible τc, depending on n (whenever T>3π/4), corresponding in (19) to two real positive roots, for each given T, each one associated with a critical τc, as shown in Figure 2.Figure 2
τ
c
X
T  following (19) and (20). One curve for each n.


## 6. Engineering Results

The results presented here enable the design of a clock distribution network with second-order dissipative nodes [1] considering delays between nodes as a complement of the designing technics derived in [7] for delay-free networks.The general idea is to maintain the delays of the network smaller than the critical value derived here. Under this condition, the reachability of the synchronous state follows the criteria shown in [7].As the equations considered here are for normalized angular frequency and the time constantT represents the node filter, it must be equal to 2π, in order to avoid double-frequency jitter [14, 15].Consequently, takingT=2π and using (19) result in ω1=0.9804 rad/seg and ω2=0.5889 rad/seg. From (20), the corresponding critical delay values are τc1=3.7013 and τc2=10.1853, for n=1.Simulations of (12) were performed by using DDE23 Matlab’s subroutine, with integration step Δt=0.1 and simulation time t=1000s, for three different situations: Figure 3(a): τ<τc1, Figure 3(b): τ=τc1, and Figure 3(c): τ>τc1.Simulations for (12) with parameter T=2π for three different values of τ. (a) τ<τc1; initial condition:  (1.0x10-1,0);  attractor: equilibrium point. (b) τ=τc1; initial condition: (1.0x10-1,0);  attractor: limit-cycle. (c) τ>τc1; initial condition:  (1.0x10-3,0);  attractor: infinity.
(a)
(b)
(c)Finally, it can be concluded thatτc1 is a bifurcation parameter. Figure 4 shows a periodic solution, related to Figure 3(b). The angular frequency experimentally measured for ψ(t) is ω1=0.98 rad/seg.Figure 4
ψ
(
t
)
×
t under Figure 3(b) conditions.



---

*Source: 145936-2013-11-11.xml*