# Modeling, Analysis, and Experimental Validation of Frequency Dividers with Direct Injection

Authors: Antonio Buonomo; Alessandro Lo Schiavo
Journal: Journal of Electrical and Computer Engineering
 (2013)

---

## Abstract

An analytical model and a methodology are presented for the analysis of CMOS injection-locked frequency dividers with direct injection, which are used in modern wireless communication systems. The amplitude and phase of the oscillation in the synchronous operation mode, as well as the locking range, are found in explicit form. The width of the locking range is determined by the condition of the existence and stability of synchronous oscillations. The accuracy of the model and of the presented formulas is validated through a comparison with experimental results, which are in good agreement with the analytical ones.


---

## Body

## 1. Introduction

The possibility offered by the modern RF-CMOS technology to realize on-chip injection-locked frequency dividers (ILFDs) has led to a renewed interest for the study of the behavior of LC oscillators under the action of an external signal, firstly performed in a pragmatic, and effective, manner in a pioneering paper by Adler [1] and then in [2, 3]. Many attempts have been devoted to seeking topologies able to broaden the locking range [4] and to increase the division ratio. A number of topologies of ILFDs have been devised to perform the frequency division with a power consumption lower than that of digital dividers and thus usable advantageously in PLLs.Investigations were also devoted to the issue of developing a model of the ILFDs and a methodology for their analysis, at first addressed in [5–7] and subsequently in many papers (see [8] and references therein), which is motivated by the need for a better understanding of their operation in order to provide useful design insights. Among these, mention should be made of papers [9, 10] where a comprehensive methodology is provided for the analysis of most commonly used topologies of ILFDs, namely, with injection through a tail device and with direct injection.However, the frequency dividers with direct injection have not been sufficiently treated in the literature from an analytical point of view, although they are widely used in applications for a very low input capacitance [11, 12] and a wide locking range [11, 13], which is one of the most important figures of merit of ILFDs. In the present work, we wish to demonstrate the accuracy of the model of ILFDs developed in [9], here made more general, as well as the accuracy of the methodology for their analysis under steady-state operating conditions. This allows us to determine the amplitude and the phase shift of the locked oscillation, as well as the locking range by analytical expressions. While the methodology was tested in previous papers through numerical simulations only, here the accuracy is verified for the first time through experimental tests on a circuit prototype that shows a good agreement between theory and experiments.


## 2. Modeling

A divide-by-2 injection-locked frequency divider with complementary topology and direct injection (Figure1) can be modeled by the equivalent circuit in Figure 2 [9]. The model is formed by the parallel connection of an LC tank, whose losses are represented by an equivalent resistance R, a nonlinear two-terminal transconductor described by the current-voltage relationship inl(v), and a controlled current source described by the current-voltage relationship iin(v,vin), where v denotes the differential voltage between the drain nodes of the ILFD. The LC tank and the transconductor represent the basic oscillator that is perturbed by the injection current iin(v,vin), which depends on both the injection signal vin, and the tank voltage v.Figure 1
Circuit diagram of an injection-locked frequency divider with complementary topology and direct injection.Figure 2
Equivalent model of the ILFD in Figure1.In the circuit in Figure1, the complementary pairs of cross-coupled devices (M1-M2, M3-M4) are sized to have complementary characteristics that set the bias voltages at the drain nodes to half the supply voltage, in order to maximize the oscillation swing. These devices form the composite two-terminal shown in Figure 3(a), which is described by a nonlinear characteristic, inl(v) [14–17]. It is easy to realize that inl(0)=0 and, due to the circuit symmetry, inl(v) is an odd function of v. Since the two-terminal sub-circuit in Figure 3(a) has to compensate the tank losses, it exhibits a negative resistance, and thus, it is locally active with a current-voltage characteristic inl(v) lying into the second and fourth quadrants of the i-v plane (Figure 5).(a) Locally active two-terminal sub-circuit of the ILFD in Figure1; (b) injection circuit of the ILFD in Figure 1. V0 and VDD-V0 are bias voltages.
(a)
(b)The two complementary MOS switches (Ms1, Ms2) of the circuit in Figure 1 form the injection circuit shown in Figure 3(b), which is driven by the external periodic signal vin. Under the action of that signal, the injection circuit injects into the tank a current iin that depends in a nonlinear way on both the injection signal vin and the tank voltage v. Consequently, the injection circuit can be represented by a two-terminal element providing a controlled current described by the nonlinear function iin(v,vin). Since the injection circuit perturbs the basic LC oscillator in a passive way, iin(v,vin) is depicted by a family of current-voltage curves lying into the first and third quadrants of the i-v plane (Figure 6).From the previous observations, it follows that the behavior of the ILFD in Figure1 can be predicted through the equivalent circuit shown in Figure 2. The development of a behavioral model which allows us to analyze the divider in Figure 1 by an analytical approach relies on the possibility of identifying the functions inl(v) and iin(v,vin) through few parameters, which is possible, as we propose to show.First, we observe that usually the dividers operate with small amplitudes of the injection signal, that is, withvin small enough to make reasonable the following approximation:
(1)iin(v,vin)≈iin(v,0)+g(v)vin,
where
(2)g(v)=∂iin∂vin|vin=0≈iin(v,VinDC)-iin(v,-VinDC)2VinDC,VinDC being a constant voltage applied to the inputs, in addition to the bias voltages V0 and VDD-V0. Then, we observe that nonlinear functions inl(v), iin(v,0), and g(v) can be easily measured at the circuit terminals in Figures 3(a) and 3(b). To show this, a prototype of the circuit in Figure 1 was implemented using integrated circuits ST HCF4007UB (dual complementary pairs plus inverter) powered by a VDD=5V battery. While NMOS and PMOS in Figure 3(b) were implemented using a single device of the integrated circuit, each of the four cross-coupled active devices in Figure 3(a) was implemented with the parallel connection of two transistors of the integrated circuit to increase the transconductance parameter.In order to measure the nonlinear characteristics of the circuit under test, that is,inl(v), iin(v,0), and g(v), the auxiliary circuit in Figure 4 was connected to the nodes x and y of the circuit in Figure 3(a) and, then, of the circuit in Figure 3(b). Voltages at these nodes, formed by a bias component equal to VDD/2 and an AC component, were applied through 100Ω series resistors that allow us to measure the absorbed currents using voltage buffers and op-amp differential amplifiers. While in case of perfect matching of the transistors in Figure 3(a) we measure i1=i2, in case of transistor mismatch we measure slightly different currents. Therefore, we calculate the mean value of these currents; that is, inl=(i1+i2)/2 for the circuit in Figure 3(a) and iin=(i1+i2)/2 for the circuit in Figure 3(b). If the circuit under test is the locally active two-terminal sub-circuit in Figure 3(a), the measured current inl is a function of the applied voltage v, as shown in Figure 5. The same figure shows that the cubic approximant [18–20]
(3)inl(v)=anlv+cnlv3
is in good agreement with the actual characteristic.Figure 4
Auxiliary circuit used to apply voltagesvt1 and vt2 to the circuits in Figures 3(a) and 3(b) and to measure the absorbed currents i1 and i2.Figure 5
Nonlinear characteristic of the locally active two-terminal sub-circuit, shown in Figure3(a) with VDD=5V. Measurements are obtained using integrated circuits ST HCF4007UB and the testing circuit in Figure 4. Approximate model is obtained using a cubic polynomial with anl=-1.94mA/V and cnl=80uA/V3.(a) Nonlinear characteristics of the injection circuit shown in Figure3(b) for different values of the constant injection voltage vin=VinDC, obtained for VDD=5V and V0=1V. Measurements are obtained using integrated circuits ST HCF4007UB and the testing circuit in Figure 4. Approximate model for vin=0 is obtained using cubic polynomial (4) with ain=0.13mA/V, bin=40uA/V2, and cin=35uA/V3. (b) Voltage-dependent transconductance g(v). Measurement’s curve is obtained using (2) and curves shown in (a). Approximate model is obtained using the last term in (4) with k=-0.6mA/V2.
(a)
(b)On the other hand, if the circuit under test is the injection circuit shown in Figure3(b), the measured current is iin and it is plotted in Figure 6(a) as a function of the applied differential voltage, v, for different values of the applied input voltage vin=VinDC. Figure 6(a) shows that iin(v,vin=0) can be modeled by a cubic polynomial, including a quadratic odd function binv|v|; that is, iin(v,0)=ainv+binv|v|+cinv3, where |v| denotes the absolute value of v. As shown in Figure 6(a), the approximated function fits well the actual characteristic. Exploiting the measured curves iin(v,1) and iin(v,-1), shown in Figure 6(a), it is possible to estimate the function g(v) according to (2) as shown in Figure 6(b). Even if a fifth-order polynomial provides a very good approximation for g(v), a linear approximation; that is, g(v)≈kv, is sufficient to obtain simple and quite accurate analytical results, as we show in Section 4. Thus, the injection nonlinearity takes the following expression:
(4)iin(v,vin)=ainv+binv|v|+cinv3+kvvin.The schematization of the divider as in Figure2 and the identification of the functions (3) and (4) allow us to describe the circuit operation with the following model
(5)(D2+ω0QD+ω02)v=-Dω0RQ(inl(v)+iin(v,vin))
obtained by the Kirchhoff’s current law at the common node in Figure 2 and by using D to denote the derivative operator with respect to time t. ω0=1/LC is the resonant frequency of the tank and Q=ω0RC its quality factor.It should be observed that the values assumed by the parameters of the algebraic characteristics (3) and (4) and used to model the transconductor and the injection circuit, obviously depend on the technology used as well as on the transistor sizes, but this modeling approach can be advantageously applied independently on the device technology. The parasitics and the delay of devices, due to the high frequency, are included, as a rule [5, 6, 8], into the capacitance of the LC tank. This simplifies the analysis, making it possible to capture the essential aspects of the synchronization phenomenon in ILFDs.


## 3. Analysis

If a sinusoidal synchronization signalvin=Vincos(2ωt) with frequency ω close to the tank resonant frequency ω0=1/LC is applied to the circuit in Figure 1, it operates as a divider-by-2 because, by virtue of the last term in (4), the injection circuit essentially performs a mixing between the input signal vin=Vincos(2ωt) and the tank voltage v. Consequently, the current injected into the tank has a component at fundamental frequency ω, which determines the tank voltage, while the harmonics are filtered out by the LC tank. We will therefore seek a solution to (5) in the form of a quasiharmonic function of time, namely, v=V(t)cos[ωt+θ(t)]. Here, V(t) and θ(t) denote the envelope of the oscillation and the phase difference between the forced oscillation and the forcing signal, respectively, which do not vary in time when synchronization occurs. Thus, in order to reveal the locking conditions, we need to formulate the explicit equations describing the evolution of V(t) and θ(t). Locking mean that there are stable fixed points in these equations, so we have to find the conditions for these points to exist and to be stable.It should be observed that the exact periodic solution of (5) for arbitrary values of the parameters ω0, Q, Vin, and ω cannot be found analytically with the mathematical tools available. However, we can analytically find the approximate solutions that describe the unknown solutions with a certain degree of accuracy. This is possible resorting to the methods of asymptotic analysis for weakly nonlinear systems. Applying the well-known solution method of slowly varying amplitude and phase (SVA) to (5), as shown in detail in [9], we obtain a couple of equations for the amplitude and the phase, that is, the so-called averaging equations
(6)V˙(t)=-ω0aR2QV(1aR-1+8b3πV+34cV2-kcVin2cos2θ),θ˙(t)=Δ-ω0QaRkcVin4sin2θ,
where the dot denotes the derivative with respect to t, Δ=ω0-ω is the frequency detuning of the output voltage v, a=-(anl+ain), b=bin/a, c=(cnl+cin)/a, and kc=-k/a.The formulation (6) of the described equations of the ILFD in terms of amplitude and phase is a useful tool for analysis. These equations have a fundamental importance in the theory of synchronization and their significance is due to the fact that the analysis of (5) that is non autonomous, that is, depending explicitly on time, is reduced to the analysis of the autonomoussystem of (6). In terms of bifurcation theory, instead of analyzing the periodic orbits of (5), we can analyze the fixed points in (6), which is obviously much easier. This makes it possible to determine the following: (i) the amplitude and the phase of the output voltage in steady state, (ii) the locking range; (iii) the stability characteristics of locked oscillations, (iv) the amplitude and the phase of the output voltage during the transient preceding the steady-state. Note that previous approaches for analyzing ILFDs are limited only to the calculation of the locking range with the simplifying assumption that the amplitude of the output voltage is not dependent on the frequency [1, 5, 6]. The fixed points of (6) that correspond to V˙=θ˙=0 can be found by solving the resulting nonlinear algebraic system for the amplitude V and the phase θ. It can be shown that this system can be solved in a closed form [9], leading to the following physically meaningful solution:
(7)V=43caR-1aR+2kcVin3c1-(4CΔakcVin)2+(16b9πc)2-16b9πc,θ=12arcsin(4CΔakcVin)+nπ.
As solution (7) is defined only if |Δ|<akcVin/(4C), it exists only within the following frequency range:
(8)2|Δ|MAX=akcVin2C
centred around ω=ω0. In order to investigate the stability of solution (7) in the interval of existence (8), we linearize the nonlinear system (6) around the solution (7) and calculate the eigenvalues of the corresponding Jacobian matrix. Defining m=kcVin/4 and Δn=CΔ/a, we obtain the following eigenvalues:
(9)λa=-2m1-(Δnm)2,λb=λa-(1-1aR-4b3πV).
Taking into account the start-up condition of the basic unforced LC differential oscillator, that k<0, and hence, m>0, both of the eigenvalues are negative within the interval of existence of the solution. Thus, we conclude that (8) provides the frequency interval wherein the injection-locked frequency divider is able to lock the input signal, that is, the so-called locking range. Expression (8) clearly shows that the locking range of direct-injection frequency dividers, differently from the locking range of frequency dividers with injection via tail device [6, 9], does not depend on the circuit losses, and thus, it can be widened without lowering the quality factor and the phase noise of the divider.Numerical integration of system (6) confirms the validity of formulas (7)–(9), as shown in Figure 7, where the phase portrait of the system is represented for two set of parameters, the first of which corresponds to the second set of experimental measurements shown in Figures 8(a) and 8(b).Phase portrait of system (6) with bin=40μA/V2, cnl+cin=115μA/V3, k=-0.6mA/V2, R=1950Ω, and Vin=0.5V. (a) for anl+ain=-1.81mA/V and Δn=0, (V,θ)=(3.9V,0°) is a stable node with λa=-0.083 and λb=-0.763; (b) for anl+ain=0.53mA/V and Δn=0.1, (V,θ)=(1.02V,22.2°) is a stable node with λa=λb=-0.202.
(a)
(b)(a) Steady-state amplitude and (b) phase for the ILFD in Figure1 with VDD=5V, V0=1V, L=4.759 μH, and C=110 pF (f0=6.956MHz). Results are obtained by measurements (lines with circles) using integrated circuits ST HCF4007UB and by (7) (lines without circles) using parameter values deduced by Figures 5 and 6.
(a)
(b)Finally, note that the presented formulas provide several interesting and useful design insights to circuit designers, which are usually interested in design variables such as oscillation amplitude, lock range, and stability of lock, which are here given in a simple explicit form in terms of circuit parameters.


## 4. Experimental Results

In order to experimentally verify the circuit model and the presented formulas, a prototype of the circuit in Figure1 was implemented using integrated circuits ST HCF4007UB (dual complementary pairs plus inverter) powered by a VDD=5V battery. The RLC tank was realized using an RF choke of BC series by Epcos, a ceramic disc capacitor of 140–500 N5 series by Xicon and a metal film fixed resistor of MF-RC series by Xicon.The presented methodology for predicting the behavior of the injection-locked frequency divider provides the analytical expressions for the output amplitudeV, the phase θ, and the locking range LR=2|Δ|MAX as a function of the active and passive circuit parameters. In particular, the active parameters of the ILFD appear through the parameters of the nonlinearities, that is, anl, cnl, ain, bin, cin, and k. Thus, we first determined the parameters of the nonlinearities using the auxiliary circuit in Figure 4, obtaining anl=-1.94mA/V and cnl=80μA/V3 for the nonlinear characteristic of the locally active two-terminal sub-circuit, and ain=0.13mA/V, bin=40μA/V2, cin=35μA/V3, and k=-0.6mA/V2 for the nonlinear characteristic of the injection circuit. The equivalent capacitance due to the combined effect of the tank capacitor, the active device parasitics, and the layout parasitics was determined by measuring the actual free-running oscillation frequency f0=ω0/(2π), that is, the output divider frequency when vin=0. The overall tank resistance was estimated by measurements at the tank resonant frequency.Then, we compared the results obtained by formulas (7) with the measurements of the output amplitude and phase of the circuit in Figure 1. The results are shown in Figure 8, for different values of the amplitude of the input signal, Vin, the frequency of the input signal, fin=2f, and circuit losses, R. Figures 8(a) and 8(b) show that formulas (7) proved an overall acceptable estimation. Moreover, it is interesting to observe that a quite good approximation to measurement results is obtained also by formula (8) for the locking ranges, which are represented in Figures 8(a) and 8(b), by the width of the frequency intervals, over which the analytical and numerical results are available. However, formulas are not able to predict a constant phase offset θ0=5°, that is, due to the reactive behavior of nonlinearities, which were here assumed to be memoryless. Finally, it should be observed that the accuracy of amplitude formula is better for small injection values as it is derived assuming small injection values and neglecting output harmonics [21–26]. Thus, we can conclude that, even if small discrepancies are present between experimental and analytical results, formulas are able to capture the main characteristics of the injection-locking phenomenon, providing a quite good approximation to experimental results.Finally, it should be highlighted that even if our experiments were performed in MHz range as in [8, 16], similar results are expected at GHz range because it was shown in [5] that approximate models based on algebraic characteristics for the active parts of the divider can be effectively used at GHz frequencies.


## 5. Conclusions

We have experimentally verified the possibility of predicting accurately the behavior of injection-locked frequency dividers with direct injection through simple formulas that provide useful design insights to circuit designers. This is made possible by the dynamical model (6) of frequency dividers that describes their nonlinear dynamic to first-order, according to the well-known averaging method. This method is a powerful tool since it allows us to predict not only the oscillations in the synchronized states and their stability, which have been analyzed in this investigation, but also their transient behavior before reaching the steady states.



---

*Source: 365692-2013-03-31.xml*