^{1}

^{*}

^{1}

^{*}

^{1}

^{*}

^{2}

^{*}

A Colpitts oscillator, working around a 3 GHz frequency, contains a double gate Metal Oxide Semiconductor transistor (DGMOS). A mixed-mode analysis is involved, applying a quantum model to the device, whereas the rest of the considered circuit is governed by Kirchhoff’s laws. The Linear Time Variant (LTV) model of phase noise is based on the Impulse Sensitivity Function of the Colpitts Oscillator which describes carefully the sensitivity of an oscillator to any impulse current injection in any node of the circuit. Finally, we improve the phase noise modeling, confronting some analytical developments to mixed-mode simulations.

The development of wireless communications imposes more and more stringent requirements on the phase noise. The adequate phase noise computer modeling and simulation allow the investigation, prediction and minimization of this important characteristic at the design stage.

In wireless transceivers, a significant part of the power is consumed by the VCO (Voltage Controlled Oscillator). VCOs are implemented in silicon-based technologies which are mainly realized as differential cross-coupled topologies due to the poor substrate isolation properties of silicon based technologies. Their efficiency is low because in addition to the oscillator core, an output buffer is needed [1-3]. Phase and frequency fluctuations have therefore been the subject of numerous studies [4-9].

Since the (low frequency) noise sources present in the circuit play a significant role in generating the phase noise, MOS oscillators are especially prone to higher phase noise because of the notoriously poor 1/f noise performance. Therefore, it is essential to design them carefully in order to implement better RF CMOS integrated circuit (ICs) and to push the performance limits of CMOS implementations.

Phase noise in RF oscillators is an important item in the design of integrated transceivers. As a result, different phase noise theories have emerged. Several papers addressing phase noise have been published [5,10-14]. Recently, Hajimiri and Lee [

In this work, we want to investigate some phase noise phenomena. We will apply the theory described by Hajimiri and Lee in 1998, which makes use of the ISF function.

The paper is organized as follows. In Section 2, we present the double gate transistor structure and the Poisson-Schrodinger model is considered. In Section 3, we present the static characteristics (DC) for the DGMOS and we point the reduction of short channel effect (SCE) in this structure. Section 4 presents a Colpitts oscillator using a mixed-mode simulation, which is carried out using the DESSIS tool in the ISE-TCAD [

The geometry that we consider is bidimensional. _{G} = 10 nm, oxide thickness T_{OX} = 1.5 nm, silicon thickness T_{Si} = 1.5 nm; source/drain length L_{SD} = 5 nm. The source and drain regions are doped N_{D} = 10^{20} cm^{−3}, the channel is undoped N_{A} = 10^{−10} cm^{−3} (_{GS} applied to the two gates are identical. All calculations were done at the temperature 300 K.

This device was simulated with a 2D solution to Poisson equation coupled to a 1D Schrödinger equation. Usually, 2D Poisson equation describes electrostatic transport and 1D Schrodinger equation is used to handle the quantum transport. The coupling of Poisson and Schrödinger equations is necessary when the oxide thickness is decrease and the thickness of the channel takes values close to the wavelength of the electrons. This approach is justified by several results presented in the literature [16,17].

The Poisson (1) and Schrodinger (2) equations are given by the following equations:

where: V(x,y) is the electrostatic potential; ρ(x,y) is the charge density; ε is the permittivity of the dielectric material constant; m* is the electron effective mass, q is the electron charge, h is the Planck’s constant and ψ(y) is the wave function corresponding to the eigenvalue E.

We clearly see that the Equations (1) and (2) are coupled. It is therefore self consistency in their resolution. We can illustrate the self consistent system by the following one:

where, the functions S[V(y)] and P[ρ(y)] represent the Schrödinger and Poisson equations.

The Poisson and Schrodinger equations are discretized using the finite difference method:

For the Poisson’s equation:

The discrete Poisson equation is nonlinear, and therefore written as:

where: F(V_{i}) is the Poisson equation functional.

The potential V, solution of Equation (4) can be written as: V = V_{i} + ΔV_{i}.

For a limited development in Taylor series of a first order, we can write:

where: J is the Jacobian system matrix (4), which equals to the derivative of the Poisson equation functional.

From the expression (5) we obtain:

For the Schrödinger equation:

Discretized Schrödinger equation can be written in a matrix form:

where: H is the Hamiltonian system defining by:

For solving the Schrödinger equation, one must know the electrostatic potential V(x,y) and to solve the Poisson equation, we must know the carrier concentration or charge density ρ(x,y). There is therefore self consistency between the resolutions of the 2D Poisson and the 1D Schrödinger equations.

The static characteristic of the considered DGMOS is obtained by self-consistent Poisson-Schrodinger system. _{DS} = 0.6 V, V_{GS} = 0.6 V bias voltages. The potential energy barrier between source and channel is pointed.

The evolution of the potential energy of double gate structure for different gate voltages V_{G} is presented

drain I_{DS} is then obtained.

_{DS} (V_{GS}) characteristics for different channel lengths. It can be observed that when channel length (L_{G}) decrease, drain current (I_{DS}) increase also I_{OFF} (off-current) increase too; as well as lead to a decrease of the threshold voltage dramatically. In order to have an acceptable threshold voltage, we study the influence of the gate work function (

On conclude that for a nanometric length of channel, and in order to optimise the threshold voltage, the solution is to choose a gate metal with a high value of work function.

A mixed-mode analysis involve, applying a 2D Poisson-

Schrodinger model to the device, whereas the rest of the circuit used is governed by Kirchhoff’s laws (

Although all components of a circuit create noise, the major source of noise in an oscillator is the active device. The noise sources in any oscillator circuit ultimately combine to form amplitude noise and phase noise. The phase noise is usually characterized in the frequency

domain. For an ideal oscillator, the output can be expressed as V_{out}(t) = V_{0}cos[ω_{0}t + Φ_{0}], where the amplitude V_{0}, the frequency ω0, and the phase reference Φ_{0} are constants. The phase noise is given by the following [

where P_{sideband}(ω_{0} + Δω, 1 Hz) represents the single sideband power of the Δω frequency offset, and P_{carrier} is the total power. Therefore, the unit for phase noise by definition is dB/Hz, which reads decibel with respect to carrier per Hertz (

The phase noise models describe the phase noise generation mechanism in oscillators. There are, at present, two models of phase noise oscillator output: Leeson’s Model and Hajimiri’s one. In order to increase the accuracy of the phase noise predictions a more sophisticated theory is needed Hajimiri’s model.

Hajimiri’s model is based on the time-varying properties of the oscillator current waveform and the resultant implications on the phase noise production.

The main idea of Hajimiri’s phase noise model is to characterize the phase noise of an oscillator by its phase impulse response. It implies that a given voltage or current impulse cause phase shift. So the phase noise can be characterized by an impulse response [

where u(t – τ) is the unit step function, q_{max}, the maximum charge displacement, Γ(ω_{0}τ) is the impulse sensitivity function (ISF). This later is periodic; it can be expressed as a Fourier series:

where C_{0} is the flicker noise conversion factor and C_{n} is the white noise conversion factor. The phase noise predicted using Hajimiri’s model is [

where Γ_{rms} is the root mean square value of the ISF and the power spectral density of the white noise source. The LTV phase noise model can predict the phase noise, but it is difficult to define the ISF function by simulation.

In this section, we use the linear time-variant model described by Hajimiri to analyze the Colpitts oscillator sensitivity. The method is based on the ISF calculation; this function is used to study the sensitivity of the elements (L, C_{1}, and C_{2}) of the Colpitts oscillator to a current disturbance.

For determining the ISF, we inject a pulse i(t) at time τ of period T and simulate the oscillator response for a few periods after the injection. By measuring the deviation ∆t of the moment of zero crossing, we can measure the phase difference ΔΦ = 2πΔt/T, produced by the injection (

The source i(t) is injected into the node of the oscillator output. In our measurement, the setting time of injection is increased in steps of 0.00674 ns to cover a period T_{0} = 1/f_{0} = 1/2.97 GHz 0.3358 ns. The total number of points of the ISF is 50 points per period.

The main interest of this study is to obtain the sensitivity of the elements (L, C_{1 }and C_{2}) of the Colpitts oscillator (

Now, we would like to extend the discussion to any value

of the Φ phase (Φ Î [0, 2π]), not only for a small charge injection.

If we consider a parasitic pulse injection (a “dirac”) at τ time, the global output voltage of an oscillator can be considered as the superposition of the permanent voltage and the response to the parasitic injection, because any time or phase shifted solution remains a solution of the harmonic oscillator:

With:

V_{T}: total output voltage V_{s}: output voltage of the free oscillator (without parasitic pulse)

C: global capacity of the Colpitts tank, C_{1}C_{2}/(C_{1} + C_{2}), (see

This equation can be reformulated as:

With

*note: we get directly the ISF, by deriving (16bis):

Now, from these general formulas, and assuming a limited expansion of φ, we obtain:

If we consider only the first order, we get:

If we compare this latter formula to [10, appendices], we find the well known result:

Γ = –sin(ω_{0}t)-see Equation (11), with qmax = CV_{0}; this result is valid, but for q very small. Note that a small signal analysis cannot yield any information on the amplitude of the oscillation, as this is determined by the nonlinearities of the system. In a real circuit, the frequency of oscillation depends on the transistor junction capacitances (that can be handled by our mixed mode simulations, anymore).

We represent on the

By. In this typical case, we have from the “mixed-mode” CV_{0} = 100 fc, in very good accordance with the “arctan” fit of equation 16 bis, considering also the graph asymmetry (shift) along the y axis. We have chosen, for this characteristic calculation, an angle very near of 90˚, for a high φ magnitude (not exactly 90˚, to “observe” the influence of the denominator of equation (22 bis), i.e.: 81˚ (90˚ - 10% of 90˚).

Up to our knowledge, it is the first time that such comparisons are presented. Their interest is twofold: our “microscopic” simulations seem validated by our analytical formulations, and we do not need to stay in the linear domain, like in the standard ISF theory–see, for instance: ([

According to the Equations (16) and (16bis), phase or/and amplitude can change after the perturbation. We represent on _{DD} .Our observations are:

- At zero crossing For (, or τ = T/4 or 3T/4):

For zero crossing falling (τ = T/4): φ = –φ_{max}; For zero crossing rising (τ = 3T/4): φ = φ_{max}.

This is quite negligible if q is small (q/C = V_{0}: our case).

- At the peaks For φ = φ_{min} = 0: ω_{0}τ = π or 2π, or τ = T/2 or T):

That mean for ω_{0}τ = 2π

and ω_{0}τ = π: (22 bis)

Our “raw” simulations (q < 0-current pulse negative) show these results (cf.

To go further, an interesting question is: how to calculate phase noise, when Dirac pulses are applied?

The φ(τ)/ω fluctuations are to the phase noise origin, depending on the noise sources (e.g. process induced traps for the devices, jitter noise).

If we are interested by time fluctuations, we can study the following expression:

Taking the square of its RMS value:

If now, we say that this process is a second order stationary one, the RMS value does not depend on time, i.e.

.

The second term in the right side-Equation (24)-is the so called correlation function: C_{φ}(τ).

Then Equation (24) is reduced to:

Take back the parasitic waveform which comes from a capacitive charge injection, at f_{0} frequency, as treated above; we have (for weak values of φ; otherwise: take Equation (18).

And (27bis)

Or (28)

This gives the associated noise power density:

The last goal is to calculate the oscillator power spectral density output (PSD), via the autocorrelation output voltage:

(FT: Fourier Transform)

It is easy to verify that the variance satisfy the relations:

With these conditions and the fact that φ(t) varies slowly over time, Equation (30) can be approximated as:

With φ^{2} is given by (28).

Then, the Equation (32) can be rewritten:

Find the Fourier Transform in Equation (33) comes down to find the Fourier transform of i.e.:

This general formula (numerical solution), can be reduced, for the specific case where q = CV_{0}; Equation (33) can be rewritten:

That gives:

In this paper, the analysis of a Colpitts oscillator designed in a nanometer CMOS technology is presented. Firstly, we have analyzed the drain current in DGMOS devices; we see that if the gate lengths or work functions are increased, the drain current is decreasing. The future trends will analyse the Double Gate transistor with Highk gate dielectric.

Secondly, we go further on phase noise theory of such oscillator, starting with mixed mode analysis of a LCtype oscillator; we derive some analytical solutions for the phase noise, through a new lecture of the impulse sensitivity function (ISF). Our first numerical simulation raw results are very consistent with such a numerical model. So, we have been invited to study more in-depth of the ISF, and try to extend its efficiency.

The future trends will lead to apply these later theoretical insights to a resonant tunnelling nano-device, not only on a first order approach, but also on a noise point of view.

This work is supported by UPM (Union pour la méditéranée)