CMOS晶體震盪器模擬

monkeybad

之前有討論過有關CMOS Pierce crystal osc
3 t; I1 C2 u/ C1 c* L$ k. n4 H現在我的模擬好像有點問題 想請教一下
圖一為我簡化後的Model 基本上 就是一顆CMOS Inverter 然後並連上一個偏壓電阻
然後a b兩點再分別接到Crystal的兩端
" Q, R$ d6 l- S! f/ a% o8 HCrystal用的model是一個L, C, RS串聯在並連一個CO 震盪頻率為32.768K
6 N$ B& ^6 h8 T. U1 a( n這是要for時鐘用的
+ |2 b, F+ `- S: {8 O, v2 h# I以下是我模擬的值 這個是從網路上面抓的一組值
% Q! v4 e: x0 T3 g) \0 Z3 [9 sL=7000, C=0.003p, RS=200K, CO=1.7p
RBIAS=1E7,
C1=1.2p, C2=1.2p
  o! y0 p7 L, I7 {1 [% JCMOS INVERTER: MN W=0.22U L=10U
5 u' J/ I2 g$ J  _! k9 b0 X) l                         MP W=0.22U L=10U
6 e6 ~9 S9 \$ m0 f, E7 W/ iVDD=3.3V
根據理論 圖二看進去的負電組值 也就是電壓值實部的絕對值 要大於Crystal Osc裡面的寄生電阻RS
就是負電阻加正電阻要為負 才會震盪

+ G7 K+ W7 A2 z% B% _圖二是我模擬負電組的電路
& _1 [1 u& _6 j9 qaa和bb電壓偏在Inverter飽和區裡面 V(aa)=V(bb)=1.3V
我在aa bb這兩點加一個電流源 hspice的指令為
( C* t" U9 B! dIin aa bb ac=1
然後我跑.ac模擬
.ac DEC 1e5 1k 40k
再下
6 h0 u' v! g8 h.probe VR(aa,bb)
+ Z) T! o- m6 o3 r顯示aa,bb電壓的實部
看ac的圖 在Avanwaves裡面叫做Volts Real的圖 看VR(aa,bb)
然後發現在32.768K左右的時候 大約是-520K左右
  ^0 I0 c4 L4 a6 K- l" d% F, t; |請問一下這樣模擬負電阻的方法對嗎?

8 K; Q' V: d  o4 l假如是正確的 那負加正小於零應該會起震才對

/ `& a4 b2 y$ s; L5 g/ J& g在跑.TRAN模擬
- Q& ?- z5 r! k: A7 |* m我的初始值設定為
.ic V(A)=0 V(B)=1 V(C)=0 V(D)=0
可是後來發現一開始有震盪 但是越到後來振幅就越小 到最後整個衰減掉了 電壓收斂在1.3V
+ U+ K* }% C8 q7 _5 [# y請問是哪裡模擬出了問題呢?

finster

來補回答我的經驗
我之前是做0.18um的RTC的32.768KHz的Crystal Oscillator,條件規格是整體的Crystal Oscillator + Buffer + I/O PAD的current < 3uA,工作電壓要從1.2V ~ 1.95V

: d3 v/ V+ ?4 E' x/ x" d1. 首先,要做low current Crystal Oscillator時,要把I/O PAD給考慮進去,我之前在設計時,Crystal Oscillator設計一段時後,模擬及電流均在規格內可以動了,但在加入I/O PAD後,Crystal Oscillator卻不會起振了,後來追了快一個星期才發覺到因為整體current非常小,但因為I/O PAD本身的leaking current過大,導致Crystal Oscillator無法起振,所以,low current Crystal Oscillator無法使用Standard I/O PAD(leaking current會過大),必需使用額外的I/O PAD才行,這點非常重要,因為它會直接影響到起振與否的問題
+ K8 o$ M) n3 a) p
) \1 w# R7 L4 |( Z$ e- Z5 L2. 要作到low current Crystal Oscillator,本身的inverter不能使用單一個PMOS和一個NMOS,必需採用Bias voltage control的方式才行(架構有點類似1999年IEEE Conference paper那種,但需要再變化一下),同時,為了low current,本身的buffer以及送出clock的buffer需採用long channel的方式,而之所以要採用Bias voltage control的inverter架構,乃因要作到low current Crystal Oscillator,本身在XTAL IN和XTAL OUT兩端點的信號並無法達到rail-to-rail的幅度,故而若採用standard inverter方式,並不能在大的工作電壓中起振,再者,因為要low current,所以在XTAL IN端點的信號幅度會很小,所以,利用Bias voltage control方式來限制住inverter兩端點的voltage level,這個並不好作,必需不斷地試

3. low current Crystal Oscillator的外接兩端電容要非常小,我當時的外接電容方式是用可調的架構(類似1991年IEEE Conference paper那種,不過這兩個電容我是做在內部),我當時試過,外部電容不能超過10pf,否則會因為外部電容過大而無法起振,那時我只用3-4pf的電容(有考慮到封裝廠的寄生電容,I/O PAD的寄生電容,板子上的電容以及實際量測時probe的電容),另外,因為本身兩端的電容值很小,所流過的電流也異常地小,故而只要有稍大的noise或者leaking current,都會造成Crystal Oscillator無法起振

* g+ ~: b# J7 i" r4. 連接在XTAL IN和XTAL OUT兩端的回授電阻要很大,最好是50M ohm以上,而這是為了low current Crystal Oscillator而特別考量的,用MOS電阻即可
: w0 |* t  O8 t! u4 x2 a
5. 最後一點,這是經驗談,當初我負責的chip回來,量測function正常,量測到的current約在4uA左右,可是後來整合到digital counter時卻出現問題,追了快兩個星期,才發覺到因為是low current Crystal Oscillator,所以其送出外面的Buffer也是使用long channel,故而它的推力很weak,造成rising time/falling time會有過長的現象,而對digital counter而言,過長的rising time/falling time則會有重覆trigger的問題,以致兩個個別電路均沒有問題,但整合在一起時卻出現錯誤計數的問題,建議你,在送到digital counter時,請他們加一個schmitter buffer,以確保因為low current的緣故而產生過長的rising time/falling time的問題不會造成digital方面的誤判

finster

先給你幾篇我之前作RTC Crystal Oscillator的參考paper
( N: e  \) S$ y; w- A& b找之前也做過類似的circuit,而且current只有3uA
1. 1999, IEEE conference paper, An improved low power crystal oscillator
2. 1997, IEEE Journal paper, A 2.1MHz Crystal Oscillator Time Base with a current comsumption under 500nA
3. 1995, IEEE conference paper, Low Power and Very low EMI high efficiency high frequency crystal oscillator
4. 1991, IEEE conference paper, A Cell Based Suer low power CMOS crystal oscillator with "on chip" tuning capacitors
這幾篇的架構可參考看它們的架構和方式
晚些再分享一下當初我做的經驗以及遇到的相關問題和要注意的事項

monkeybad

圖二為負電阻模擬

圖三為負電阻模擬結果

圖四為b點電壓 .TRAN的模擬結果 發現振幅一直變小
& l4 A3 z* B1 i# _$ f0 [7 o
[ 本帖最後由 monkeybad 於 2007-8-1 05:03 PM 編輯 ]

adele

方法上，应该没有问题。
4 ^# N/ Y5 g% M8 b( u# |  V试一下：L=10.88KH, C=2.17fF, RS=30.9Kohm, CO=980fF.
9 x  c% i8 S  nC1=C2=15pF.  RBIAS=1E7,

建议可以先用一个10uA的电流源下面接一个nmos管的结构，nmos管gm设在40u.
你可以用spectre仿真吗？

m851055

在原本的invert處多串一個invert，並且把CO的aa端接到此invert的輸入端，再配合capacitor之調整就可以了。
一般invert之Length較小 width較大，你的設定值怪怪的，也可能因此造成驅動能力不佳而衰減。

wwm101

"我的初始值設定為
  K3 ?1 y- z. U  X7 w1 R        W8]@8n3FlU.ic V(A)=0 V(B)=1 V(C)=0 V(D)=0"
这样的设置没有多大的用处，直接给电感设置初始电流，而且初始电流也不能太小，一般在ua级比较合适！

monkeybad

For adele:
$ h" l4 `) {$ ?7 C2 j$ A# ~有試過你給的值 而且用一顆NMOS gmN=44u 上面用定電流源PMOS來偏壓 I=37uA
可是好像還是不會震耶
另外我是用HSPICE模擬的 我這邊沒有spectre耶
" J6 ~. E: T0 H4 H" K% ]5 J8 e8 j& D
For m851055
' N9 s- P. k/ ~" u4 ]因為要做RTC(Real Time Clock)的OSC, 所以要求非常省電, Inverter的W/L值才會設這麼小

For wwm101 :
那要怎麼下指令去給電感初始電流? 還有初始條件要怎麼設才恰當?

4 [  Q5 n0 V1 X$ U$ D有沒有人有經驗的能指導一下 感溫感溫!

. V& q/ U* I0 `* {/ g[ 本帖最後由 monkeybad 於 2007-8-6 06:17 PM 編輯 ]

monkeybad

感謝finster大大提供這麼寶貴的經驗 真是造福人群 功德無量
請受我一拜!

monkeybad

昨天晚上試了一下以前在網路上找到的一個方法
* U7 i/ y% T6 [. N% U( j6 D$ u& i就是在一開始的時候 串聯一個起始電壓源來加速Hspice的模擬
起始電壓源為一個sine波 頻率設為crystal的震盪頻率 然後要加入damping factor 讓起始電壓慢慢衰減到零
! h. j' K4 x5 J+ @$ V結果神奇的事情發生了
Crystal Oscillator到最後居然開始震盪了
好感動喔!
趕快來跟大家分享一下!
! O' W. D, K9 P$ }# R1 j- P+ F
$ O5 y* L/ t/ P- {) y電路仍然是之前的Inverter電路
在d點加入起始電壓
VK d1 d2 sin ( 0V 800V 32K 0 1e5)
模擬時間約400ms
.TRAN 0.01u 400m
( h, f, T- Y, B& _' E; M5 W到最後就看到震盪器開始震盪了!

" b8 X+ x+ n8 o2 V5 w% j後來又試了另一個類似的方法 就是在電路另外並聯一個電流源
也是同樣的sin波然後慢慢衰減到零 結果也是能讓電路起震

; C* X* l$ |0 {9 G- d另外好像模擬時間的間隔不能設太大 設太大好像也不會起震
/ c4 {3 D- i: r. s
2 F; g! X1 Y( D0 d, Y附件為網路上找到的文件

, a; K+ r$ H* z. N! _" N2 W: t0 v[ 本帖最後由 monkeybad 於 2007-8-9 10:14 AM 編輯 ]

finster

建議你採用這個方式來試
VDD   VDD    0    PWL(0 0  10ms  3.3)

並加入initial voltage
.ic v(A)=3.3 v(B)=3.3
& q% @) n7 W2 v6 E
將電壓源採用ramp的方式來作,那是因應實際電壓源皆是從0慢慢爬升到穩定的電壓源
再者,將crystal model兩端設初使值為3.3V,最主要是因為inverter會將0-->1-->0-->1
inverter的兩端若皆為3.3,那在一開始時,inverter的output會因為input為3.3而讓output由3.3-->0,當inverter的input受到回授電阻的緣故而讓input由3.3-->0,又會讓output由0-->3.3,如此循環下去就可起振,當然,這要在crystal能夠起振才行
8 t6 E$ [% ~* U% d. m$ z4 y' u5 x如果一開始就是設v(A)=0, v(B)=3.3,或者v(A)=3.3, v(B)=0,那對inverter而言就會是一個穩定的狀態而不會進行起振的環境條件
$ [$ A/ @6 W9 k. j/ u
最後,實在不建議你用電壓源來作初使值的起振方式
crystal能不能夠起振,端看整個loop有否達到起振的條件
設初始值只在HSPICE一開始模擬才有效,真正在模擬時己經脫離初始值的狀能,設初始值是能夠幫助HSPICE在計算時有收斂的成效
* s: }/ _( F* V8 ~/ ^而使用電壓源或者電流源來作初始值的起振條件,相當於我們外加一個信號源來讓電路動作,但實際卻沒有這個信號,這對實際chip的環境而言是不對的,再者,用電壓源或者電流源作初始值要持續一段時間,而這種現象對crystal整體系統而言,無疑干擾了整個loop的運作

monkeybad

恩 這種加起始電壓的方法的確怪怪的 因為實際上晶片沒有這種起始電壓
應該是要用正統的加起始值以及VDD上升等條件一起模擬才符和實際的情況

/ x7 `1 h# z( h( \0 k1 c2 R( `看參考資料裡面是說 在晶體震盪器模擬都要花費相當長的時間 這種方法可以加速模擬的時間
' q; L' ]8 C! j/ t因為時間過很久後起始的電壓還是衰減到零了 對電路沒有影響 所以最後假如電路還是可以震盪的話
; R! m4 J8 _/ E4 V6 {4 z1 o. _我想OSC的設計應該是會震才對
6 K, f) V) `0 s9 K
我覺得這種方法只能表示最後會不會震盪 不能模擬剛開始的時候會不會起震 另外就是可以節省模擬的時間
% l# ^4 p; G% r
8 W- A: N! |+ Y* L[ 本帖最後由 monkeybad 於 2007-8-13 03:35 PM 編輯 ]

skyboy

原帖由 monkeybad 於 2007-8-1 04:50 PM 發表
之前有討論過有關CMOS Pierce crystal osc
, O9 d8 K6 S) \' Q: u& F1 E8 D+ {% G現在我的模擬好像有點問題 想請教一下
5 u9 M, }6 n/ F. Y圖一為我簡化後的Model 基本上 就是一顆CMOS Inverter 然後並連上一個偏壓電阻
7 b  a' K6 q6 Q. v% f& I然後a b兩點再分別接到Crystal的兩端
Cryst ...

; @! K& F- {: G* I2 K9 o: ?因為最近剛好也要做一個 1MHZ~10MHz 的 X'tal oscillator, 最近這兩篇討論對小弟
幫助很大..不過在模擬上有遇到問題, 怕是自己觀念有錯, 想請各位前輩指導...

+ a7 y) |5 W- K) P2 \3 A關於負電組的部份, 老實說是第一次看到這個觀念, 研讀了這兩篇討論串所提供的一些
  V2 n  l0 m) E' `, U! @文件, 總算有了一點點感覺...在推導的部份 monkeybad 大大提供的 Philips 那份文件
7 m$ z5 w7 A; v$ q' q5 w有詳細說明, 裡頭推論出滿足負電阻條件時 inverter 最小的 gm 值, 不過就是在這邊
遇到一點問題...
2 B! Q  q- r3 o* f
我模擬的條件為:
; [6 @; C; q  ~' o* b; ~10M X'tal model 為: RS=50, CS=0.004pF, L=60039.21315uH, CP=5pF (前人留下的..)
; `+ Q3 t, P2 W% u: AInverter size 為: WN=3u, WP=10u, L=0.5u, 偏壓在約 1.65V (VDD/2)
+ k1 q! {6 _# C& hRP 為 2E6
CX1=CX2=10p (X'tal兩邊的負載電容)
% |7 N' [- c0 a+ I+ @0 x
假設先不考慮 Cpar 的效應, 根據公式算出來的 gmmin 約為 115u 左右, 用本篇一開始
monkeybad 大大提出的方式去模擬, HSPICE 模擬出來 NMOS 的 gm 約為 521u, 是理論
# f3 ^4 ^) Z3 @" x' u最小值將近五倍, 然而 AC 分析看到的阻值卻是正的(請參考附圖), 我試著調整參數, 發現
將 CX1, CX2 以及 RP 縮小才會得到類似 monkeybad 大提供的圖 (整個從1MEG到10MEG
的頻域都是負電阻..) 不過阻值也不大..
" M; v+ p! i2 H- K
- [4 H9 T% z7 `0 i2 K另外一個問題是, 根據公式, RP 如果愈大, gmmin 會愈小, 所以我自己的推論是如果 RP 愈大,
( `( u2 C4 a( o, G" \" j, `模擬結果應該會往"更負"的那邊移動, 不過我用 50E6 模擬結果卻與預期完全相反...

研究了兩天實在不知道是哪邊弄錯了, 煩請前輩們指導一下..謝謝..

monkeybad

根據推導的公式 RP應該是越大 所需要的gmmin會越小  這應該是沒錯
8 d5 Z6 R  H' m: c' `. Y" y  l! L
另外說實話 我模擬負電阻時也有遇到全部變成正的情形 但是我也不知道為什麼 我猜可能是相位的問題吧 Hspice把他反過來了
所以我後來就直接根據推導的公式去計算 就是之前那份Philip的文件 移項一下就可以得到正電阻跟負電阻的計算公式
) ?, t8 a4 f/ L1 H& B4 \1 b然後至少負電阻絕對值要大於正電阻 根據之前網友的說法是說要大五倍比較保險

然後在時域上面去模擬 有一個可以參考的方法來驗證你OSC會不會震 就是一開始起始值的電壓用很小
7 B/ T' m8 ]1 C然後你去看震盪器裡面的波形振幅 如果振幅慢慢被一直放大 則表示你的OSC是會震盪的
* f, `* x) _; W$ L! A6 O) P9 o  a& w因為表示負電阻的力量是比較強的
: D3 L/ ^3 [7 t. {
若是要看steady state時的波形 就把起始值加大 加速他transient的時間
: T2 K# W* _" l9 t
另外注意好像.TRAN時間間隔也不能設太小的樣子 不然也不會震 目前我模擬是設0.01U 但是你頻率比較快 可能在設小一點吧

skyboy

哎哎...負電阻的部分還是參不透玄機...模擬不出來...

不過我按照推導的結果確認 gm 值足夠後把整個線路丟進 ADiT 去模擬, 把 X'tal 兩端設成
! b" n* b: h+ k3 Q; l# x0 kSPICE_NODE, 結果是會順利起震的, 可是又遇到另一個問題, 就是用 HSPICE 模擬無法起震,
& O: F" H" @& T3 ^  a4 R0 ?: v: h我覺得應該是模擬條件設定的問題, 我有參考了前面 finster 大大提過的技巧, 不過兩個端點
最後還是會停在偏壓點 (剛好約 VDD/2) 的地方, 我知道 time step 的設定似乎會影響到模擬
$ S# K* P4 E& U  e1 H4 o會不會起震, 不知有沒有前輩們能分享一下除了 time step 之外還有什麼地方要注意的??
1 g7 D2 p1 s: p5 K6 s8 M( K
另外再請教 monkeybad 大大, 您提到的 "一開始起始值的電壓用很小" 指的就是 X'tal 兩端的
初始值嗎??

monkeybad

請問一下ADiT是什麼呀

/ E/ b, j, ?* L9 B另外Hspice模擬中 我是參考之前的附件裡面設的起始電壓 就是原來附圖中的d點加入一個起始電壓
VK d1 d2 sin ( 0V 800V 32K 0 1e4) 其中d拆成了d1 d2
& d1 d1 I- a: l+ |) ~那起始電壓用小一點的話 就是sin的震幅用小一點 例如設成50V這樣
+ P& P, W! K" c7 A  @7 @0 @8 S用這樣設起始條件模擬時間會比較快一點

5 v. A8 k0 h$ r, B6 F另外再模擬起震時 參考finster的設法 VDD用ramp模擬上升 另外V(A) V(B)就設成3.3V
那因為crystal從零一直到起震需要很久的時間 所以一開始你會看到XTALIN XTALOUT兩端都會被拉到VDD/2
4 C0 r% m" D: g/ r但是你去看crystal model裡面的c點 應該會開始震盪 而且會慢慢變大
2 u) M& d' n5 |% D3 p根據穩態模擬的結果 震盪到穩態時電壓振幅都會到上百V至上千V
* Y1 v4 z$ ^/ ~
另外我看你crystal model最好是去找看看 實際上廠商在用的model參數下去設
) ]+ }% \: u  X5 a$ X
' Z2 \4 i0 Z. d$ ]2 S2 p9 P, D2 S[ 本帖最後由 monkeybad 於 2007-8-24 05:09 PM 編輯 ]

skyboy

ADiT 是本土公司 (華凱, EverCAD) 開發的一套 simulation tool, 這套 tool 標榜模擬的速度
比 HSPICE 快上許多, 也有提供準確度較高的 spice mode...Mentor 為了加強他在電路模擬
這個領域的 tool, 在去年(前年?? 忘了..)購併了 EverCAD...所以這套 tool 現在也是爆貴...-_-
7 O7 r3 R9 |" {) j4 i3 n" ]
以上是題外話....

X'tal 的 model 似乎不容易要到, 我 12MHz 的 model 是參考 Philips 文件上的...那應該算是
' m. W. `9 i: q) c0 \. W: v廠商的 model 吧??....
" |# {  L0 X0 e7 ?  |4 _% U
我去看了不會起震時 C 點的波形, 初始值是參考 finster 大大的建議, 震幅與操作點一開始從
% l, l' ~+ w" C接近 full swing  衰減, 然後隨著 VDD 逐漸上升, 當 VDD 到 3.3V 時 C 點會上升到偏壓點附近,
接著就是維持著約 350mV 的振幅到模擬結束...:f33 試過調整 time step 與 VDD ramp up
% s& G2 ]+ e' v8 H) a的時間, 看來都沒有什麼關係...
0 I' p" r) J) s+ k/ c2 C
! R7 S2 I0 L$ n我還沒有試過 monkeybad 大提供的外加 power 的作法, 不過在網路上找資料時有看到大大
& j& P: m; P+ r* X提供的那份文件, 只是如同 finster 大所說, 這是一個不存在的電源, 總覺得這樣模擬可能不是
很正確...

不知道是哪邊設錯了, 請各位前輩能再指導一下...

finster

如果說Crystal Model沒有問題的話,有用Power ramp和加入初使值的設定而依然不能夠起盪的話
2 G- J: E1 m7 ?- J  J1 Z+ p) O那我只能說應該是你電路本身的問題

5 u* q' h( B* C, _- l+ w$ M假如是用hspice在作模擬的話,32.768khz的要跑上一個星期是常有的事
- E" [  J6 b0 S2 E! {: e4 Y( Q如果是使用spectre來作模擬的話,大概半天到一天應該就會有結果出來
建議你兩個方向試試看
2 g: j/ [) {% i. L一是先測試一下反相器的輸入電壓對輸出電壓的關係,你說你的偏壓大約在VDD/2,但我看了一下你反相器的size再對應映到你的反相器的轉換電壓卻似乎並不在這個區間,先單獨對反相器作模擬,如果你的輸入電壓並不在這個範圍的話,那你的反相器根本無法作出反應,也根本無法起振
, t* V9 ~8 P$ _9 X" }% H先讓反相器的輸入電壓可以在接近VDD/2作出反轉,然後再試著加大size試試看
7 K/ D8 D4 ]) Y% T1 c另外一個作法則是加大回授電阻,建議你觀察一下回授電阻和兩端電容的電流值及Crystal Model兩端各為多少以及流向,這點純粹是個人經驗,當然,這必須是在上面那個建議改了之後仍然無法起振才用這個,而這點也是當初我在作極低low current crystal oscillator時所觀察到的現象

1 F5 Q) b1 _% u; H* I
% p. |: J$ y, {: k  }$ }4 x; r

原帖由 skyboy 於 2007-8-27 12:32 PM 發表
ADiT 是本土公司 (華凱, EverCAD) 開發的一套 simulation tool, 這套 tool 標榜模擬的速度
比 HSPICE 快上許多, 也有提供準確度較高的 spice mode...Mentor 為了加強他在電路模擬
& F- W4 I* Q! {$ [5 s這個領域的 tool, 在去年(前 ...

monkeybad

哇 32768的crystal OSC模擬要一個禮拜喔!
原來要這麼久 難怪我跑模擬過半個小時去看 怎麼感覺都一樣勒

monkeybad

另外附帶提一下 RTC晶片規格裡面 要求非常省電
8 Q8 s* R* g8 Q# y# ]因為可以延長電池的壽命 或是外掛的精電容可以撐的比較久
1 i( w3 C# D% g" ]( Z目前我找到比較厲害的是seiko 還有一家瑞士做的 兩個差不多強
7 F6 t: U9 g3 m: ~, w, g4 Wtime keeping current=0.25uA, VDD=5V-1.3V

' ?& ^; p& ]9 s/ @& I  ]4 f真的能做到這麼省電 而且VDD可以吃這麼大的範圍的電壓嗎?
  u: H# D, u: q4 \9 f# s5 }! a* T
這種規格實在很嚇人

猜想是不是他先內建一個LDO把電壓降下來 然後再用這個LDO Vout去做類比和數位的電路的power嗎