NGSPICE TEMPLATES

==================Switchtest========================

1)Construct two voltages and a current PWL signals

2)Apply these signals to switches

3)Plot input and output waveforms and hysteresis loops

===================================================

Switchtest

*                          10

*                 R10         / s1

*             __/\  /\  /\___/   o______

*            |    \/  \/      _|_       |

*            |               / v1\      |

*            |               \___/      |

*            |                _|_       |

*            |                ///       |

*            |                          |

*            |   R20      20 /  s2      |

*            |_/\  /\  /\___/   o_______|

*            |   \/  \/      _|_        |

*   _________|              / v2\       |

*  _|_       |              \___/       |

* /v4 \      |               _|_        |

* \___/      |               ///        |

*  _|_       |                          |

*  ///       |    R30      30 / s3      |

*            |__/\  /\  /\___/ON o______|

*            |    \/  \/      _|_       |

*            |               / v3\      |

*            |               \___/      |

*            |                _|_       |

*            |                ///       |

*            |                          |

*            |                          |

*            |             40           |

*            |     R40         / w1     |

*            |___/\  /\  /\___/   o_____|

*                  \/  \/               |

*                   _  _       R40      |

*           _______/ \/ \____/\  /\  /\_|

*          _|_     \_/\_/      \/  \/  _|_

*          ///                         ///

*

.tran  2us 5ms

v1  1  0   DC  0.0 PWL(0 0   2e-3 2  4e-3 0)

v2  2  0   DC  0.0 PWL(0 0.9 2e-3 2  4e-3 0.4 )

i   3  3 0 DC  0.0 PWL(0 0   2e-3 2m 4e-3 0) 

v4  4  0   DC  2.0

r3  3  33  10k

vm3 33 0   dc  0 $ <--- measure t h e c u r r e n t

r10 4  10  10k

r20 4  20  10k

r30 4  30  10k

r40 4  40  10k

*

s1 10 0 1 0 switch1 OFF

s2 20 0 2 0 switch1 OFF

s3 30 0 2 0 switch1 ON

.model switch1 sw vt =1 vh =0.2 ron =1 roff =10k

* 

w1 40 0 vm3  wswitch1 off

.model wswitch1 csw it =1m ih =0.2m ron =1 roff =10k

*

.control

run

plot v ( 1 )     v ( 10 )

plot v ( 10 ) vs v ( 1 )    $ <-- get hysteresis loop

plot v ( 2 )     v ( 20 )   $ <--- differentinitial values

plot v ( 20 ) vs v ( 2 )    $ <-- gethysteresis loop

plot v ( 2 )     v ( 30 )   $ <--- different initial values

plot v ( 30 ) vs v ( 2 )    $ <-- get hysteresis loop

plot v ( 40 ) vs vm3#branch $ <--- current controlled sitch hysteresis

.endc

.end

* source /Users/don_sauer/Downloads/SI_Lib/Tests.cir

=====================END_OF_SPICE=======================

==============Non_linear_with_frequency_resistor===========

1)Define a B sources using the built in Hertz variable

2)Define a subcircuit nlres using the built in Hertz variable

3)Apply a parameter rb=rbase to nlres

4)Define function check using the vecmax function

5)Use the function check and the built in variable frequency with an if control

6)Plot AC output vres#branch

=================================================================

Non_linear_with_frequency_resistor

*

*                       Fx

*                ___   _  _   ___              _|_   _|_   _|_

*               |pos|_/ \/ \_|neg|            / B1\ / B2\ / B3\

*               |___| \_/\_/ |___|            \___/ \___/ \___/

*     R2               |  |                    _|    _|_   _|_

*  _/\  /\  /\___  2   Vx |_                   //    ///   ///

*  |  \/  \/    _|_         |

* _|_          / vx\        |                      Xres

* ///          \___/  1     |                33 _/\  /\  /\__ 10

*                |__________|                  |   \/  \/    |

*               _|_                           _|_           _|_

*              / Bx\                         / V1\         /Vre\

*              \___/                         \___/         \___/

*               _|_                           _|_           _|_

*               ///                           ///           ///

*

*      -1/{ rb }/sqrt(HERTZ)*v(pos ,neg)

*

B1    1   0     V = hertz*v(33)

B2    2   0     V = v(33)*hertz

b3    3   0     V = 6.283e3/(hertz+6.283e3)*v(33)

V1    33  0     DC      0   AC 1

Xres  33  10    nlres   rb =1k

Vres  10  0     DC      0

.Subckt   nlres pos neg  rb= rbase

Bx    1   0     v = -1/{ rb }/sqrt(HERTZ)*v(pos ,neg)

Vx    2   1     DC   0Volts

Rx    2   0     1

Fx    pos neg   Vx   1

.ends

.control

define check (a , b) vecmax(abs( a - b ))

ac     lin 10 100 1k

print  v(1) v(2) v(3)

if     check( v (1),frequency) < 1e-12

echo   " INFO: ok "

end

plot   vres#branch

.endc

.end

* ngspice 

* source  /Users/don_sauer/Downloads/SI_Lib/Non_linear_resistor.cir 

=====================END_OF_SPICE=======================

ngspice 3 -> source  /Users/don_sauer/Downloads/SI_Lib/Non_linear_resistor.cir 

Circuit: non_linear_resistor

Doing analysis at TEMP = 27.000000 and TNOM = 27.000000

No. of Data Rows : 10

                              non_linear_resistor

                              AC Analysis  Sat Mar 30 14:42:05  2013

--------------------------------------------------------------------------------

Index   frequency       v(1)                            

--------------------------------------------------------------------------------

0 1.000000e+02 1.000000e+02, 0.000000e+00

1 2.000000e+02 2.000000e+02, 0.000000e+00

2 3.000000e+02 3.000000e+02, 0.000000e+00

3 4.000000e+02 4.000000e+02, 0.000000e+00

4 5.000000e+02 5.000000e+02, 0.000000e+00

5 6.000000e+02 6.000000e+02, 0.000000e+00

6 7.000000e+02 7.000000e+02, 0.000000e+00

7 8.000000e+02 8.000000e+02, 0.000000e+00

8 9.000000e+02 9.000000e+02, 0.000000e+00

9 1.000000e+03 1.000000e+03, 0.000000e+00

                              non_linear_resistor

                              AC Analysis  Sat Mar 30 14:42:05  2013

--------------------------------------------------------------------------------

Index   frequency       v(2)                            

--------------------------------------------------------------------------------

0 1.000000e+02 1.000000e+02, 0.000000e+00

1 2.000000e+02 2.000000e+02, 0.000000e+00

2 3.000000e+02 3.000000e+02, 0.000000e+00

3 4.000000e+02 4.000000e+02, 0.000000e+00

4 5.000000e+02 5.000000e+02, 0.000000e+00

5 6.000000e+02 6.000000e+02, 0.000000e+00

6 7.000000e+02 7.000000e+02, 0.000000e+00

7 8.000000e+02 8.000000e+02, 0.000000e+00

8 9.000000e+02 9.000000e+02, 0.000000e+00

9 1.000000e+03 1.000000e+03, 0.000000e+00

                              non_linear_resistor

                              AC Analysis  Sat Mar 30 14:42:05  2013

--------------------------------------------------------------------------------

Index   frequency       v(3)                            

--------------------------------------------------------------------------------

0 1.000000e+02 9.843334e-01, 0.000000e+00

1 2.000000e+02 9.691501e-01, 0.000000e+00

2 3.000000e+02 9.544281e-01, 0.000000e+00

3 4.000000e+02 9.401466e-01, 0.000000e+00

4 5.000000e+02 9.262863e-01, 0.000000e+00

5 6.000000e+02 9.128287e-01, 0.000000e+00

6 7.000000e+02 8.997566e-01, 0.000000e+00

7 8.000000e+02 8.870535e-01, 0.000000e+00

8 9.000000e+02 8.747042e-01, 0.000000e+00

9 1.000000e+03 8.626939e-01, 0.000000e+00

 info: ok 

ngspice 4 -> 

===============Parameter_to_pwl_Bsource==========================

1)Define many x and y parameters

2)Apply the x and y parameters to B sources in PWL format

3)Do a DC sweep

4)Plot outputs with various offsets

=================================================================

Parameter_to_pwl_Bsource

.param x0=-4 y0=0 

.param x1=-2 y1=2

.param x2=2  y2=-2

.param x3=4  y3=1

.param xx0=  x0-1

.param xx3=  x3+1

Vin 1 0 DC=0V

R   1 0 2

Btest2 2 0 I = pwl( v(1),x0,y0,x1,y1,x2,y2,x3,y3)

Btest3 3 0 I = (    v(1)<x0) ? y0:(v(1)<x3) ?  pwl(v(1),x0,y0,x1,y1,x2,y2,x3,y3):y3

Btest4 4 0 I = pwl( v(1), xx0,y0,x0,y0,x1,y1,x2,y2,x3,y3,xx3,y3)

Btest5 5 0 I = pwl(-2*i(Vin),xx0,y0,x0,y0,x1,y1,x2,y2,x3,y3,xx3,y3)

Rint2  2 0 1

Rint3  3 0 1

Rint4  4 0 1

Rint5  5 0 1

.control

dc Vin -6 6 0.2

plot v(2) v(3) v(4)-0.5 v(5)+0.5

.endc

.end

* ngspice 

* source  /Users/don_sauer/Downloads/SI_Lib/pwlfunctioninBsource.cir 

=====================END_OF_SPICE=======================

=================VCCS_examples============================

1)Create two parameters Vi and Offs

2)Apply parameters to devices

3)The set noaskquit will allow a quit without save

4)Transient simulate and plot outputs

=================================================================

VCCS_examples

.param Vi   = 1

.param Offs = 0.01*Vi

B21 int1 0   V = V( 3 )*V( 3 )

G1  21   22  int1 0  1

vm  22   0   dc 0

R21 21   0   1

G51 55   225 cur = V( 3 )*V(3)-Offs

vm5 225  0   dc 0

R51 55   0   1

B31 int2 0   V = V( 3 )*V( 3 )

E1  1    0   int2 0   1

R1  1    0   1

E41 4    0   vol = V( 3 )*V(3)-Offs

R4  4    0   1

V1  3    0   PWL(0 0 100u {Vi } )

.control

set   noaskquit

tran 10n     100u    uic

plot i( E1 ) i( E41 )

plot i( vm ) i( vm5 )

.endc

.end

* ngspice 

* source  /Users/don_sauer/Downloads/SI_Lib/VCCS_exmaples.cir 

=====================END_OF_SPICE=======================

===================Simple_PLL=================

1)Create VT as a PWL time voltage ramp

2)Create VIN as a sinwave source running off VT

3)Create LOSC like VIN, set to be phase modulated by voltage VPHASE

4)Multiply and low pass filter LOSC to create LPASS

5)Voltage LPASS is meant to by used as a frequency modulate signal

6)FM signal LPASS is translated to a PM signal VPHASE by integration

7)Transient simulation will show how LPASS settles to zero

8)Plot all the voltages at various points in time

9)Only eight lines of spice code needed to make a functional PLL 

=================================================================

Simple_PLL

*                                   ______________________________

*                                   |                             \

*                                   | V = sin(2*PI*Fre*T+IntFM*T)  \

*           VMULT    RLP      LPASS |                        ^      \

*              ____/\  /\  /\_______|     RINT     VPHAS    /|\      \

*            _|_     \/  \/   |     | __/\  /\  /\_____ _ _ _|        \_ LOSC

*   VIN     /   \            _|_    |     \/  \/  _|_                 / |

*    ____|\/  \/ \     CLP   ___    | FM          ___      SYSTEM    /  |

*  _|_   |/\  /\ /           _|_    |       CINT  _|_      MODELED  /   |

* /_  \     \___/ BMULT      ///    |             ///        VCO   /    |

*// \  \      ^                     |_____________________________/     |

*\   \//     /_\                                                        |

* \___/       |_________________________________________________________|

*  _|_

*  ///   WinSpice/MacSpice 2.19.11_6.58PM dsauersanjose@aol.com  

*

*=========All_System_Modeled_PLL==================

VT        VT     0      DC     0      PWL(   0      0      10      10      )

BVIN      VIN    0      V =           sin( 6.28*1000*V(VT) ) 

BLOSC     LOSC   0      V =           sin( 6.28*1000*V(VT)  + 10*V(VPHAS)*V(VT) )

BMULT     VMULT  0      V =           V(VIN)*V(LOSC)

RLP       VMULT  LPASS  1k

CLP       LPASS  0      1.6u

RINT      LPASS  VPHAS  1meg

CINT      VPHAS  0      1u

.control

*set       pensize =     2

tran      .01m   5    0  .01m

plot      lpass 

plot      vin lpass losc vmult xlimit 9m 10m

plot      vin lpass losc vmult xlimit 498m 500m

plot      vin lpass losc vmult xlimit 4998m 5000m

.endc

.end

*display

.endc

.end

* source /Users/don_sauer/Downloads/SI_Lib/Tests.cir

=====================END_OF_SPICE=======================

=================PINCH_RESISTOR========================

1)Create Vtest as a PWL voltage ramp

2)Create Rpinch as the voltage node subcircuit

3)Use VRset to set the zero volt resistance.

4)Assume the tub is connected to Vtest as in drawing

5)Use B sources to define the V to I relationship.

6)Plot the current versus voltage to see the pinchoff effect.  

=================================================================

PINCH_RESISTOR

*

*                XRpinch

*     Vtest   ________

*       _____|_/\  /\  /\_

*     _|__       \/  \/  _|_

*    /_   \      :       ///

*   // \   \     :

*   \   \_//     :        Rset

*    \____/      :./\  /\  /\_

*      _|_           \/  \/  _|_

*      ///     VRset         ///

*    Vramp

*                    Attach Pin 2 to highest voltage

*

*         ________________________________________________________

*        /               _____________________      Nepi         /

*       /  ____________ /..................../_______________   /

*      / / ________    /                    /     ________  /  /

*     / / /____ 1 /   /         2          /     /____ 3 / /  /

*    / / //___/| /   /                    /     //___/| / /  /

*   / / / |___|//   /                    /     / |___|// /  /

* _/_/_/_______/___/____________________/_____/_______/_/__/___

*  | | | p+imp |  |       n+imp         |     | p+imp | |  |

*  |  \ \_____/    \___________________/ base  \_____/ /   |

*  |   \______________________________________________/    | Psub

*  |  _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _  _ _ _ _ _ _ _ |

*  |/   N BURIED LAYER                                    \|

*   \ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ /

*          P-Substrate

*

Vramp     Vtest   0    PWL(   0     0.1  5  5.0)

XRpinch      Vtest   0    VRset  Rpinch

Rset         VRset   0    1k

.control

tran         100m    5     0     100m

plot        -vramp#branch

.endc

.SUBCKT   Rpinch  INP   INN  VRS

Isen         VRS     0     -1u

Bpinc        INP     INN    I= 1u*(v(INP)-v(INN))/(v(VRS) + v(VRS)*(v(INP)-v(INN)) )

.ENDS       Rpinch

.end

* source /Users/don_sauer/Downloads/SI_Lib/Tests.cir

=====================END_OF_SPICE=======================

===================test_pz_sanity===============================

1)Apply an ac current source to a unity LRC circuit

2)Suppress  printing of statistics  of the Initial Transient Solution using .options noacct

3)Do PZ in terms of current and poles

4)print out the two poles

5)Do for various values of R

6)Note the frequency is in terms of 1 rad/sec which is 0.159Hz.

7)Sanity check using the web page at http://www.idea2ic.com/PlayWithJavascript/L_C_R_F.html   

=================================================================

test_pz_sanity

*                        Iin

*                   _  _

*           _______/ \/ \_____

*          _|_     \_/\_/     |   _  _  _

*          ///                |  / \/ \/ \

*                          1  |  | () () |

*              _/\  /\  /\____|__|       |__

*            _|_  \/  \/    |        L     _|_

*            ///     R     _|_             ///

*                          ___ C

*                           |

*                          _|_

*                          ///

iin    1  0   ac 

r      1  0   1

c      1  0   1

l      1  0   1

.print pz all

.options noacct

.control

pz 1 0 2  0 cur pol

print pole(1) pole(2)

.endc

.end

* source /Users/don_sauer/Downloads/SI_Lib/Tests.cir

=====================END_OF_SPICE=======================

r     1  0   .100

c     1  0   1

l     1  0   1

pole(1) = -9.89898e+00,0.000000e+00

pole(2) = -1.01021e-01,0.000000e+00

r      1  0   .3

c      1  0   1

l      1  0   1

pole(1) = -3.00000e+00,0.000000e+00

pole(2) = -3.33333e-01,0.000000e+00

r      1  0   .501

c      1  0   1

l      1  0   1

pole(1) = -9.98004e-01,6.315087e-02

pole(2) = -9.98004e-01,-6.31509e-02

r      1  0   1

c      1  0   1

l      1  0   1

pole(1) = -5.00000e-01,8.660254e-01

pole(2) = -5.00000e-01,-8.66025e-01

r     1  0   100

c     1  0   1

l     1  0   1

pole(1) = -5.00000e-03,9.999875e-01

pole(2) = -5.00000e-03,-9.99987e-01

r     1  0   100k

c     1  0   1

l     1  0   1

pole(1) = -5.00000e-06,1.000000e+00

pole(2) = -5.00000e-06,-1.00000e+00

====================NPN_Ftau_15uu_TF==========================

1)Using an ac current source applied to the emitter of an NPN

2)At ftau, the magnitude of AC current in the base equals that of the collector

3)Do a foreach for five levels of emitter current myic

3)Add a foreach mytf for three levels of TF for the NPNV model

4)Do AC from 10KHz to 100Ghz each time

5)Compose five element arrays ie, ft1p, ft2p, ft4p to use for waveforms

6)Fill array ie with the five levels of emitter current myic

7)Fill arrays ft1p, ft2p, ft4p the found ftau levels

8)Plot ft1p, ft2p, ft4p versus ie to see the result

9)Use The destroy all clears all the plots.

=================================================================

NPN_Ftau_15uu_TF 

*                    

*           _________

*      b   _| c     _|_

*     ___|'QN1     /vcc\

*   _|_  |`->e     \___/

*  /vb \    |        |

*  \___/   _|_      _|_

*    |    / _ \i1   ///

*   _|_   \/ \/

*   ///   /\_/\

*         \___/

*          _|_

*          ///

*

.OPTIONS GMIN=1e-15 METHOD=gear ABSTOL=1e-15 TEMP=27 srcsteps = 1 gminsteps = 1

VCC C  0    DC 5V

QN1 C  B E  NPNV

VB  B  0 0

I1  E  0    DC 100u AC .01u

.control

destroy all

foreach mytf 1p  2p   4p

foreach myic 10u 100u 1e-3 10e-3 50e-3 100e-3

 alter  i1 = $myic

 altermod   npnv   tf = $mytf

 ac dec 10 10k 100g

end

end

echo $plots

plot abs(ac1.vb#branch) abs(ac1.vcc#branch)

plot abs(ac2.vb#branch) abs(ac2.vcc#branch)

compose ie    start = 0 stop = 5 step =1

compose ft1p  start = 0 stop = 5 step =1

compose ft2p  start = 0 stop = 5 step =1

compose ft4p  start = 0 stop = 5 step =1

let index = 0 

while (  index < length(abs(ac1.vb#branch)) )

if ( abs(ac1.vcc#branch[index]) >  abs(ac1.vb#branch[index]) )  

let ft1p[0] = mag(frequency[index] )

let ie[0] = 10e-6   

end

if ( abs(ac2.vcc#branch[index]) >  abs(ac2.vb#branch[index]) )  

let ft1p[1] = mag(frequency[index] )

let ie[1] = 100e-6   

end

if ( abs(ac3.vcc#branch[index]) >  abs(ac3.vb#branch[index]) )  

let ft1p[2] = mag(frequency[index] )

let ie[2] = 1e-3     

end

if ( abs(ac4.vcc#branch[index]) >  abs(ac4.vb#branch[index]) )  

let ft1p[3] = mag(frequency[index] )

let ie[3] = 10e-3     

end

if ( abs(ac5.vcc#branch[index]) >  abs(ac5.vb#branch[index]) )  

let ft1p[4] = mag(frequency[index] )

let ie[4] = 50e-3    

end

if ( abs(ac6.vcc#branch[index]) >  abs(ac6.vb#branch[index]) )  

let ft1p[5] = mag(frequency[index] )

let ie[5] = 100e-3    

end

if ( abs(ac7.vcc#branch[index]) >  abs(ac7.vb#branch[index]) )  

let ft2p[0] = mag(frequency[index] )

let ie[0] = 10e-6   

end

if ( abs(ac8.vcc#branch[index]) >  abs(ac8.vb#branch[index]) )  

let ft2p[1] = mag(frequency[index] )

end

if ( abs(ac9.vcc#branch[index]) >  abs(ac9.vb#branch[index]) )  

let ft2p[2] = mag(frequency[index] )

end

if ( abs(ac10.vcc#branch[index]) >  abs(ac10.vb#branch[index]) )  

let ft2p[3] = mag(frequency[index] )

end

if ( abs(ac11.vcc#branch[index]) >  abs(ac11.vb#branch[index]) )  

let ft2p[4] = mag(frequency[index] )

end

if ( abs(ac12.vcc#branch[index]) >  abs(ac12.vb#branch[index]) )  

let ft2p[5] = mag(frequency[index] )

end

if ( abs(ac13.vcc#branch[index]) >  abs(ac13.vb#branch[index]) )  

let ft4p[0] = mag(frequency[index] )

let ie[0] = 10e-6   

end

if ( abs(ac14.vcc#branch[index]) >  abs(ac14.vb#branch[index]) )  

let ft4p[1] = mag(frequency[index] )

end

if ( abs(ac15.vcc#branch[index]) >  abs(ac15.vb#branch[index]) )  

let ft4p[2] = mag(frequency[index] )

end

if ( abs(ac16.vcc#branch[index]) >  abs(ac16.vb#branch[index]) )  

let ft4p[3] = mag(frequency[index] )

end

if ( abs(ac17.vcc#branch[index]) >  abs(ac17.vb#branch[index]) )  

let ft4p[4] = mag(frequency[index] )

end

if ( abs(ac18.vcc#branch[index]) >  abs(ac18.vb#branch[index]) )  

let ft4p[5] = mag(frequency[index] )

end

let index = index +1

end

*print ft1p[0] ft1p[1] ft1p[2] ft1p[3] ft1p[4] ft1p[5]

plot ft1p ft2p ft4p vs ie loglog 

.endc

.MODEL NPNV NPN(

+ IS=10.15e-18 NF=1.0 BF=89 VAF=107 IKF=.4

+ NR=1.006 BR=0.4822 VAR=4.286 IKR=0.0002472

+ ISE=9.15E-17 NE=2

+ ISC=1E-21 NC=2

+ RB=12 RBM=3 IRB=.1

+ RE=2 RC=33

+ CJE=37e-15 VJE=0.75 MJE=0.35

+ CJC=45E-15 VJC=0.6399 MJC=0.3531

+ CJS=2.939E-15 VJS=0.3488 MJS=0.1813 XCJC=0.4201

+ TF=1e-12 XTF=1 VTF=1.5 ITF=.8

+ TR=5.9e-9 FC=0.5 PTF=30

+ KF=1.000E-16 AF=1

+ XTB=2 EG=1.11 XTI=5 TNOM=25 )

.end

* source /Users/don_sauer/Downloads/SI_Lib/Tests.cir

=====================END_OF_SPICE=======================

==========================Model_internals===================================

1)Do a normal curve tracer DC sweep on an NMOS transistor m1

2)The vss#branch has the normal output

3)A print @m1  will print out the parameter available for viewing

4)The simulated parameters can be plotted using the format of @m1[id]

=================================================================

Model_internals

*            _______

*         __| d1   _|_

*       ||____    /Vdd\

*  g1 __||    |   \___/

*   _|_ ||->  |     |

*  /Vsg\    |_|    _|_

*  \___/   _|_     ///

*    |    /Vss\

*   _|_   \___/

*   ///     |

*          _|_

*          ///

*

vdd    vdd  0    2.0  

vss    vss  0    0

vsig   g1   vss  1

m1     vdd  g1   vss vss  nmosc w=1e-5 l=3.5e-007

.control

dc     Vdd       0     5    .1    Vsig 0 5 1

plot   vss#branch 

save   @m1[id]   @m1[vds]  @m1[vdsat]  @m1[vth]  @m1[gm] 

dc     Vdd       0     5    .1    Vsig 0 5 1

plot   @m1[id] 

plot   @m1[vds]

*plot  @m1[vgs]

plot   @m1[vdsat]

plot   @m1[vth]

plot   @m1[gm]

*print @m1      $ to view model parameters

.endc

*  source /Users/don_sauer/Downloads/SI_Lib/Tests.cir

.model NMOSC NMOS

+ Level= 8 Tnom=27.0 version =3.3.2

*------------------Process-----------------------------------------------

+ tox=160e-10 xj=0.25e-06 nch=0.5e+17

*------------------V_threshold-------------------------------------------

+ vth0=0.72 nlx=0.12e-06

*------------------Bulk--------------------------------------------------

+ k1=1.04 k2=-1.209E-01

*------------------mobility----------------------------------------------

+ u0=678 ua=8.964e-10

+ ub=1.472e-18 uc=-4.441E-17 vsat=86000

*------------------Subthresshold-----------------------------------------

+ nfactor=1.8

+ cit=-5.0E-04 voff=-7.862E-02

+ eta0=4.441e-16 etab=-2.E-01 dsub=0.7

*------------------Hot electrons-----------------------------------------

+ alpha0=1.61e-05 beta0=36.68

*------------------VAF---------------------------------------------------

+ lint=.12e-06 pclm=.19 pscbe1=3.79e+08 pscbe2=9.4e-05

+ delta=0.01655 pvag=0.4484

*------------------Bulk_diode--------------------------------------------

+ js=5.858e-08

*------------------Resistance--------------------------------------------

+ rsh=70 rdsw=375

+ wr=0.7586 prwb=0 prwg=-4.441E-17

*------------------Capacitance-------------------------------------------

+ cj=0.0002424 cjsw=2.73e-10 mj=0.3551 mjsw=0.3873

+ cgso=9e-13 cgdo=9e-13 cgbo=7e-10

+ cdsc=-2.4E-4 cdscd=-1.506E-01 cdscb=-2.219E-04

+ pb=0.5614 pbsw=0.8 xpart=0

+ dlc=5e-08 dwc=1.5e-07

*------------------BulkChargeEffect--------------------------------------

* a0=0.7 a1=0 a2=1 ags=0.05583

* b0=6.305e-08 b1=6.579e-08 keta=-1.531E-02

*------------------ShortChannel------------------------------------------

+ dvt0=2.2 dvt1=0.53 dvt2=-1.521E-01 drout=0.76

+ pdiblcb=.4 pdiblc1=0.00886 pdiblc2=0.00029

*------------------NarrowChannel-----------------------------------------

+ w0=2.6e-04 wint=0.16e-06

+ ww=-9.525E-14 wwn=1.0

+ dvt0w=0 dvt1w=5.3e6 dvt2w=-1.E-01

+ k3=2.53 k3b=-5 dwg=0 dwb=0

*------------------Noise-------------------------------------------------

* af=1 kf=1e-28 ef=0.95

*------------------Temperature-------------------------------------------

* pvsat=0 ute=-1.258E+00 kt1=-3.85E-01

* kt1l=0 kt2=-3.098E-02 ua1=5.705e-09

* ub1=-1.147E-17 uc1=-1.302E-01 at=20380

* prt=-3.287E+02 lk1=0 lk2=0

* lvsat=0 la0=0 lags=0 lute=0

.end

*  source /Users/don_sauer/Downloads/SI_Lib/Tests.cir

=====================END_OF_SPICE=======================

ngspice 9 -> @m1  $ to view model parameters

@m1: no such command available in ngspice

ngspice 10 -> print  @m1  $ to view model parameters

Vds [vds]=4.11659

Vgs [vgs]=4.55899

Vbs [vbs]=-0.441014

Ids [id]=0.0063002

Vth [vth]=-0.00673526

Vdsat [vdsat]=1.03402

Gds [gds]=0.000302297

Gm [gm]=0.00139937

Gmb [gmbs]=-0.00150947

AC NQS model selector [acnqsmod]=0

Non-quasi-static model selector [nqsmod]=0

ERROR: enumeration value `CP_BOOL' or `CP_LIST' not handled in vec_get

Aborting...

Fatal error in NGSPICE - Return

================RC_3dB=====================================

1)Finding the 3dB point of a AC plot of a RC lowpass

2)Find size of the frequency array after simulation

3)Use a repeat 39 to step through all the array values 

4)Use if statement find the frequency where the dB value is -3dB 

5)Print of the found value using the  $&f format

6)Sanity check the meas AC statement

=================================================================

RC_3dB

*

*

*              R1

*  VIN                          ___

*       _____/\  /\  /\________|OUT|

*      |       \/  \/    |     |___|

*     _|__               |

*    /_   \             _|_

*   // \   \            ___ C1  Find 3dB freq

*   \   \_//             |

*    \____/             _|_

*     _|_               ///

*     ///

*

*

.OPTIONS GMIN=1e-15 METHOD=gear ABSTOL=1e-15 temp=27

VIN  IN   0   DC 0V AC 1V

R1   IN   OUT 1k

C1   OUT  0   1u

.control

run

ac   dec  10  1   10khz

plot db(v(out)) ylimit -3 0 title Find_R_C_3dB

let  n = 0

let  f = 0

repeat 39

let  n = n+1

if   ( db( out[n]) >-3)

let  f = mag(frequency[n])

endif

endrepeat

echo "3dB point = $&f Hz"

meas ac freq_3d_at WHEN  vdb(out) =-3

.endc

.end

*  source /Users/don_sauer/Downloads/SI_Lib/Tests.cir

=====================END_OF_SPICE=======================

3db point           = 158.489 hz

freq_3d_at          = 1.587800e+02

======================BSIM_NMOS_Caps=================================

1)Sanity checking the stray capacitance to NMOSC device.

2)Find size of the frequency array or length(g)

3)Use a while (ii < length(g) )to step through all the array values 

4)Use if statement find the frequency where the dB value is -3dB 

5)Print of the found value using the  $&f format

6)Note schematic expected model values included with  spice file

=================================================================

BSIM_NMOS_Caps

*

*      IN

*      ___________/\  /\  /\___________________________________

*      |     RG 10K \/  \/                                     |

*      |                                                       |

*      |                                                       |

*      |__________/\  /\  /\_____________                      |

*      |    RDS 10K \/  \/               |                     |

*      |                                 |                     |

*      |                                 |                     |

*      |__________/\  /\  /\__           | 9.1fF               |

*      |    RB  10K \/  \/    |          |                     |

*      |                      |B         |DS                 G | 3.3fF

*      |                      |12fF      |                     |

*     _|__ VIN         MB     |      MDS |            MG       |

*    /_   \          _________|       ___|            ___      |

*   // \   \        |   |_    |      |   |_          |   |_    |

*   \   \_//        |     ||__|      |     ||_       |     ||__|

*    \____/      __/|\____||      __/|\____|| |    _/|\____||

*     _|_       |   |   <-||     |   |   <-|| |   |  |   <-||

*     ///       |   |___|        |   |___|    |   |  |___|

*               |                |            |   |      |

*               |________________|____________|___|______|

*                            _|_

*                            ///

.OPTIONS GMIN=1e-15 METHOD=gear ABSTOL=1e-15 temp=27

VIN   IN 0  DC  0V  AC  1V

RG    G  IN 10k

RDS   DS IN 10k

RB    B  IN 10k

MB    B  B  B   0 NMOSC W=3u L=1u AD=7p AS=7p PD=10u PS=10u

MDS   DS 0  DS  0 NMOSC W=3u L=1u AD=7p AS=7p PD=10u PS=10u

MG    0  G  0   0 NMOSC W=3u L=1u AD=7p AS=7p PD=10u PS=10u

.control

op

ac    dec      50       100Meg    10ghz

plot  db(v(g)) db(v(b)) db(v(ds)) ylimit -3 0

set   stopp =  length(g)

*compose fg    start = 0 stop = $&stopp step =1

let   fg  = 0

let   fb  = 0

let   fds = 0

let   ii  = 0

while (ii < length(g) )

if   (db(g)[ii] > -3 )

let   fg = mag( frequency[ii] )

end

if (  db(b)[ii] > -3 )

let   fb = mag( frequency[ii] )

end

if    (db(ds)[ii] > -3 )

let   fds = mag( frequency[ii] )

end

let   ii = ii + 1

end

print fg fb fds

let   cg = 1/(6.28e4*fg)

echo "3dB point f = $&fg Hz cg = $&cg"

let   cb = 1/(6.28e4*fb)

echo "3dB point f = $&fb Hz cb = $&cb"

let   cds = 1/(6.28e4*fds)

echo "3dB point f = $&fds Hz cds = $&cds"

let   cbs =         cds/2

let   cgb =         cb - cds

let   cgs =        (cg - cgb)/2

echo               "cbs = cbd =    $&cbs"

echo               "cgs = cgd =    $&cgs"

echo               "cgb =          $&cgb" 

.endc

* ngspice

* source /Users/don_sauer/Downloads/SI_Lib/Tests.cir

*model = bsim3v3

*Berkeley Spice Compatibility 

* Lmin= .35 Lmax= 20 Wmin= .6 Wmax= 20

.model NMOSB NMOS

+Level=        8 version=3.2.2

+Tnom=27.0

+Nch= 2.498E+17  Tox=9E-09 Xj=1.00000E-07

+Lint=9.36e-8 Wint=1.47e-7

+Vth0= .6322    K1= .756  K2= -3.83e-2  K3= -2.612 

+Dvt0= 2.812  Dvt1= 0.462  Dvt2=-9.17e-2 

+Nlx= 3.52291E-08  W0= 1.163e-6 

+K3b= 2.233 

+Vsat= 86301.58  Ua= 6.47e-9  Ub= 4.23e-18  Uc=-4.706281E-11 

+Rdsw= 650  U0= 388.3203 wr=1

+A0= .3496967 Ags=.1    B0=0.546    B1= 1   

+ Dwg = -6.0E-09 Dwb = -3.56E-09 Prwb = -.213

+Keta=-3.605872E-02  A1= 2.778747E-02  A2= .9 

+Voff=-6.735529E-02  NFactor= 1.139926  Cit= 1.622527E-04 

+Cdsc=-2.147181E-05   

+Cdscb= 0  Dvt0w =  0 Dvt1w =  0 Dvt2w =  0 

+ Cdscd =  0 Prwg =  0 

+Eta0= 1.0281729E-02  Etab=-5.042203E-03 

+Dsub= .31871233 

+Pclm= 1.114846  Pdiblc1= 2.45357E-03  Pdiblc2= 6.406289E-03 

+Drout= .31871233  Pscbe1= 5000000  Pscbe2= 5E-09 Pdiblcb = -.234

+Pvag= 0 delta=0.01

+ Wl =  0 Ww = -1.420242E-09 Wwl =  0 

+ Wln =  0 Wwn =  .2613948 Ll =  1.300902E-10 

+ Lw =  0 Lwl =  0 Lln =  .316394 

+ Lwn =  0

+kt1=-.3  kt2=-.051 

+At= 22400 

+Ute=-1.48 

+Ua1= 3.31E-10  Ub1= 2.61E-19 Uc1= -3.42e-10 

+Kt1l=0 Prt=764.3

.model NMOSC NMOS

+ Level= 8 Tnom=27.0 version=3.2.2

*------------------Process-----------------------------------------------

+ tox=160e-10 xj=0.25e-06 nch=0.5e+17

*------------------V_threshold-------------------------------------------

+ vth0=0.72 nlx=0.12e-06

*------------------Bulk--------------------------------------------------

+ k1=1.04 k2=-1.209E-01

*------------------mobility----------------------------------------------

+ u0=678 ua=8.964e-10

+ ub=1.472e-18 uc=-4.441E-17 vsat=86000

*------------------Subthresshold-----------------------------------------

+ nfactor=1.8

+ cit=-5.0E-04 voff=-7.862E-02

+ eta0=4.441e-16 etab=-2.E-01 dsub=0.7

*------------------Hot electrons-----------------------------------------

+ alpha0=1.61e-05 beta0=36.68

*------------------VAF---------------------------------------------------

+ lint=.12e-06 pclm=.19 pscbe1=3.79e+08 pscbe2=9.4e-05

+ delta=0.01655 pvag=0.4484

*------------------Bulk_diode--------------------------------------------

+ js=5.858e-08

*------------------Resistance--------------------------------------------

+ rsh=70 rdsw=375

+ wr=0.7586 prwb=0 prwg=-4.441E-17

*------------------Capacitance-------------------------------------------

+ cj=0.0002424 cjsw=2.73e-10 mj=0.3551 mjsw=0.3873

+ cgso=9e-13 cgdo=9e-13 cgbo=7e-10

+ cdsc=-2.4E-4 cdscd=-1.506E-01 cdscb=-2.219E-04

+ pb=0.5614 pbsw=0.8 xpart=0

+ dlc=5e-08 dwc=1.5e-07

*------------------BulkChargeEffect--------------------------------------

* a0=0.7 a1=0 a2=1 ags=0.05583

* b0=6.305e-08 b1=6.579e-08 keta=-1.531E-02

*------------------ShortChannel------------------------------------------

+ dvt0=2.2 dvt1=0.53 dvt2=-1.521E-01 drout=0.76

+ pdiblcb=.4 pdiblc1=0.00886 pdiblc2=0.00029

*------------------NarrowChannel-----------------------------------------

+ w0=2.6e-04 wint=0.16e-06

+ ww=-9.525E-14 wwn=1.0

+ dvt0w=0 dvt1w=5.3e6 dvt2w=-1.E-01

+ k3=2.53 k3b=-5 dwg=0 dwb=0

*------------------Noise-------------------------------------------------

* af=1 kf=1e-28 ef=0.95

*------------------Temperature-------------------------------------------

* pvsat=0 ute=-1.258E+00 kt1=-3.85E-01

* kt1l=0 kt2=-3.098E-02 ua1=5.705e-09

* ub1=-1.147E-17 uc1=-1.302E-01 at=20380

* prt=-3.287E+02 lk1=0 lk2=0

* lvsat=0 la0=0 lags=0 lute=0

*      CG  = CGS + CGD + Cox      CGS = (CG-Cgb)/2

*      CB  = CBS + CBD + Cox      Cgb = CB -CDS

*      CDS = CBS + CBD            CBS = CDS/2

**

**       SPICE MODEL                      ___

**                           __/\  ______|   | Source

**          CGS             |    \/ RS   |___|

**             _____________|__________

**           _|_     _|_        _|_   _|_

**     Gate  ___    / _ \ I_ds  \ /   ___ CBS  6fF<- (Caps for L/W 1u/2.6u)

**      ___   |     \/ \/       _V_    |      

**     |   |__|     /\_/\  |     |_____|     ___   

**     |___|  |___  \___/  V  ___|     |____|   | Bulk   10K10F 100p

**            |   |   |      |  _|_   _|_   |___|

**           _|_  |__/|\__||_|   ^    ___     |__..no..__

**           ___      |2fF||    /_\    |  CBD           _|_

**            |_______|__________|_____|  6fF            ^ 

**           CGD        CGD |             ___           /_\

**                          |__/\  ______|   | Drain    _|_

**                               \/ RD   |___|         \sub/

**                                                      \ /

**                         <Bulk_hidden = gnd!>          V

**   

**                  (G)  

**           (D)   ___|__   (S)        (B)         (Sub)

**           _|_  |______|  _|_        _|_          _|_ 

**    ______|___|_|______|_|___|______|___|________|___|__ _

**     |  |   P+  | <++++|   P+  |  |   N+  | |  |   P+  |   Xj

**     |   \_____/        \_____/    \_____/  |   \_____/  _

**     |                NWELL                 |

**     \______________________________________/ 

**   

**   

.end

* source /Users/don_sauer/Downloads/SI_Lib/Tests.cir

=====================END_OF_SPICE=======================

No. of Data Rows : 101

fg  = 2.754229e+09

fb  = 1.000000e+09

fds = 1.148154e+09

3db point f = 2.75423E+09 hz cg = 5.7815E-15

3db point f = 1E+09 hz cb = 1.59236E-14

3db point f = 1.14815E+09 hz cds = 1.38688E-14

cbs = cbd =    6.93442E-15

cgs = cgd =    1.86339E-15

cgb =          2.05472E-15

======================NPN_GUMMEL=========================================

1)Sanity checking NPN device.

2)DC sweep the base voltage

3)Use current of VC and VB to measure collector and base current

4)Plot i(VC) and  i(VB) versus base voltage V(b)

5)Plot i(VC) / i(VB) versus collector current  

=================================================================

NPN_GUMMEL

* MEASURE IS,ISE,NF,NE,RE,IKF,BETA 

*

*                 ________

*                |        | VC

*                |C      _|_

*         B     _|      /2v \

*         ____|'  npnv  \___/

*        |    |`->        |

*    VB _|_      |  0     |

*      /.7v\     |________|

*      \___/      _|_

*        |        ///

*       _|_

*       ///  0      1u

.OPTIONS GMIN=1e-15 METHOD=gear ABSTOL=1e-15

*======== ====== ====== ====== ====== ====== ====== ====== ======

VC        C      0      DC     5V

VB        B      0      0V

Q1        C      B      0      NPNV

.MODEL    NPNV   NPN(

*==========================================================

+IS=1.1E-18 NF=1.005 BF=220 VAF=130 IKF=5e-03

+ISE=9.15E-15 NE=2

*==========================================================

+CJE=2E-12 CJC=2E-13 CJS=3E-12 TF=.6E-9 )

.control

*DC SOURC1 VSTART VSTOP VSTEP SOURC2 START2 STOP2 STEP2

dc  vb     .4     1.2   .1

plot mag(-i(vc)) mag(-i(vb)) vs mag(V(b))      ylog   title Gummel

plot mag(i(vc)/i(vb))        vs mag(vc#branch) loglog title Beta_vs_IC

.endc

.end

* source /Users/don_sauer/Downloads/SI_Lib/Tests.cir

=====================END_OF_SPICE=======================

=============NMOS_bsim3_SubSanity=================

1)Sanity checking subthresshold behavior of NMOS.

2)DC sweep the gate voltage

3)Plot i(V1)  loglog 

4)Note lab data included at end of spice file. 

=================================================================

NMOS_bsim3_SubSanity

    *            _______

    *       3  _|_     _|_

    *         /V1 \   /Vds\

    *         \___/   \___/

    *         __| 1     |

    *       ||____     _|_

    *   2 __||    |    ///

    *   _|_ ||->  |

    *  /Vgs\    |_|

    *  \___/   _|_

    *    |     ///

    *   _|_

    *   ///

    *

Vds 3 0 dc 5v

Vgs 2 0 dc 1.2v

v1  3 1 dc 0v

m1  1 2 0  0 N1 W=3u L=1u AD=7p AS=7p PD=10u PS=10u

.OPTIONS GMIN=1e-15 METHOD=gear ABSTOL=1e-15

.dc Vds 1m 5 0.01 vgs 0 1 .1

.control

run

plot mag(i(v1)) loglog

.endc

***************SILICON_DATA****************************

* see reference curve below

*Lmin= .35 Lmax= 20 Wmin= .6 Wmax= 20

.model N1 NMOS

+ Level= 8 Tnom=27.0

*------------------Process-----------------------------------------------

+ tox=160e-10 xj=0.25e-06 nch=0.5e+17

*------------------V_threshold-------------------------------------------

+ vth0=0.72 nlx=0.12e-06

*------------------Bulk--------------------------------------------------

+ k1=1.04 k2=-1.209E-01

+ cdsc=-2.4E-4 cdscd=-1.506E-04 cdscb=-2.219E-04

*------------------mobility----------------------------------------------

+ u0=678 ua=8.964e-10

+ ub=1.472e-18 uc=-4.441E-17 vsat=86000

*------------------Subthresshold-----------------------------------------

+ nfactor=1.8

+ cit=-5.0E-04 voff=-7.862E-02

+ eta0=4.441e-16 etab=-2.E-01 dsub=0.7

*------------------Hot electrons-----------------------------------------

+ alpha0=1.61e-05 beta0=36.68

*------------------VAF---------------------------------------------------

+ lint=.12e-06 pclm=.19 pscbe1=3.79e+08 pscbe2=9.4e-05

+ delta=0.01655 pvag=0.4484

*------------------Bulk_diode--------------------------------------------

+ js=5.858e-08

*------------------Resistance--------------------------------------------

+ rsh=70 rdsw=375

+ wr=0.7586 prwb=0 prwg=-4.441E-17

*------------------Capacitance-------------------------------------------

+ cj=0.0002424 cjsw=2.73e-10 mj=0.3551 mjsw=0.3873

+ cgso=9e-13 cgdo=9e-13 cgbo=7e-10

+ pb=0.5614 pbsw=0.8 xpart=0

+ dlc=5e-08 dwc=1.5e-07

*------------------BulkChargeEffect--------------------------------------

+ a0=0.7 a1=0 a2=1 ags=0.05583

+ b0=6.305e-08 b1=6.579e-08 keta=-1.531E-02

*------------------ShortChannel------------------------------------------

+ dvt0=2.2 dvt1=0.53 dvt2=-1.521E-01 drout=0.76

+ pdiblcb=.4 pdiblc1=0.00886 pdiblc2=0.00029

*------------------NarrowChannel-----------------------------------------

+ w0=2.6e-04 wint=0.16e-06

+ ww=-9.525E-14 wwn=1.0

+ dvt0w=0 dvt1w=5.3e6 dvt2w=-1.E-01

+ k3=2.53 k3b=-5 dwg=0 dwb=0

*------------------Noise-------------------------------------------------

+ af=1 kf=1e-28 ef=0.95

*------------------Temperature-------------------------------------------

+ pvsat=0 ute=-1.258E+00 kt1=-3.85E-01

+ kt1l=0 kt2=-3.098E-02 ua1=5.705e-09

+ ub1=-1.147E-17 uc1=-1.302E-01 at=20380

* prt=-3.287E+02 lk1=0 lk2=0

+ lvsat=0 la0=0 lags=0 lute=0

+ luc=0

*Lmin= .35 Lmax= 20 Wmin= .6 Wmax= 20

.model N2 NMOS

+Level= 8 Tnom=27.0

****************Process***********************

+Tox=16e-9 Xj=250E-09 Nch=.5E+17

****************V_thresshold***********************

+Vth0=.72 Nlx=0.12E-06

****************Bulk***********************

+K1=1.04 K2=-1.2e-1 K3=-2.612

****************Mobility***********************

+U0=478 Ua=6.47e-9

+Ub=4.23e-18 Uc=-4.706281E-11 Vsat=86301

****************Subthresshold***********************

+NFactor=1.14

+Cit= 1.6E-04 Voff=-6.74E-02

+Eta0=1.03E-02 Etab=-5.04E-03 Dsub= .32

****************VAF***********************

+Lint=.12e-6 Pclm=.19 Pscbe1=4e8 Pscbe2=5E-09

+delta=0.01 Pvag=.44

****************Resistance***********************

+Rdsw=650 Prwg =0 Prwb = -.213 wr=1

****************Capacitance***********************

+Cdsc=-2.15E-05 Cdscb=0 Cdscd=0

****************BulkChargeEfffect***********************

+A0=.35 A1= 2.8E-02 A2= .9 Ags=.1

+B0=0.546 B1=1 Keta=-3.6E-02

****************ShortChannel***********************

+Dvt0=2.812 Dvt1= 0.462 Dvt2=-9.2e-2 Drout= .31871233

+Pdiblcb= -.234 Pdiblc1= 2.5E-03 Pdiblc2=6.4E-03

****************NarrowChannel***********************

+W0=1.163e-6 Wint=1.47e-7 Ww=-1.42E-09 Wwn=.2613948

+Wl=0 Wwl=0 Wln=0

+Dvt0w=0 Dvt1w=0 Dvt2w=0

+K3b= 2.233 Dwg=-6.0E-09 Dwb= -3.56E-09

****************Process***********************

+Ll=1.30E-10 Lw=0 Lwl=0

+Lln=.316394 Lwn=0

+Ute=-1.48

****************Temperature***********************

+Kt1l=0 kt1=-.3 kt2=-.051

+Ua1= 3.31E-10 Ub1= 2.61E-19 Uc1= -3.42e-10

+At= 22400 Prt=764.3

.model NMOSC NMOS

.model NMOSC NMOS

+ Level= 8 Tnom=27.0

*------------------Process-----------------------------------------------

+ tox=160e-10 xj=0.25e-06 nch=0.5e+17

*------------------V_threshold-------------------------------------------

+ vth0=0.72 nlx=0.12e-06

*------------------Bulk--------------------------------------------------

+ k1=1.04 k2=-1.209E-01

*------------------mobility----------------------------------------------

+ u0=678 ua=8.964e-10

+ ub=1.472e-18 uc=-4.441E-17 vsat=86000

*------------------Subthresshold-----------------------------------------

+ nfactor=1.8

+ cit=-5.0E-04 voff=-7.862E-02

+ eta0=4.441e-16 etab=-2.E-01 dsub=0.7

*------------------Hot electrons-----------------------------------------

* alpha0=1.61e-05 beta0=36.68

*------------------VAF---------------------------------------------------

+ lint=.12e-06 pclm=.19 pscbe1=3.79e+08 pscbe2=9.4e-05

+ delta=0.01655 pvag=0.4484

*------------------Bulk_diode--------------------------------------------

* + js=5.858e-08

*------------------Resistance--------------------------------------------

+ rsh=70 rdsw=375

+ wr=0.7586 prwb=0 prwg=-4.441E-17

*------------------Capacitance-------------------------------------------

+ cj=0.0002424 cjsw=2.73e-10 mj=0.3551 mjsw=0.3873

+ cgso=9e-13 cgdo=9e-13 cgbo=7e-10

+ cdsc=-2.4E-4 cdscd=-1.506E-01 cdscb=-2.219E-04

+ pb=0.5614 pbsw=0.8 xpart=0

+ dlc=5e-08 dwc=1.5e-07

*------------------BulkChargeEffect--------------------------------------

* a0=0.7 a1=0 a2=1 ags=0.05583

* b0=6.305e-08 b1=6.579e-08 keta=-1.531E-02

*------------------ShortChannel------------------------------------------

+ dvt0=2.2 dvt1=0.53 dvt2=-1.521E-01 drout=0.76

+ pdiblcb=.4 pdiblc1=0.00886 pdiblc2=0.00029

*------------------NarrowChannel-----------------------------------------

+ w0=2.6e-04 wint=0.16e-06

+ ww=-9.525E-14 wwn=1.0

+ dvt0w=0 dvt1w=5.3e6 dvt2w=-1.E-01

+ k3=2.53 k3b=-5 dwg=0 dwb=0

*------------------Noise-------------------------------------------------

* af=1 kf=1e-28 ef=0.95

*------------------Temperature-------------------------------------------

* pvsat=0 ute=-1.258E+00 kt1=-3.85E-01

* kt1l=0 kt2=-3.098E-02 ua1=5.705e-09

* ub1=-1.147E-17 uc1=-1.302E-01 at=20380

* prt=-3.287E+02 lk1=0 lk2=0

* lvsat=0 la0=0 lags=0 lute=0

* luc=0

.end

* source /Users/don_sauer/Downloads/SI_Lib/Tests.cir

***************SILICON_DATA*****************************

*

*         NMOS_Subthresshold_Drain_Current and Gate_Voltage

*    1mA  _______________________________________________

*        |           .           .           .           .

*        |  w=20um   . l = 1u    .           .           .

*        |           .           .           .           .

*  100uA |.......................................960mV....

*        |           .           .        960mV          .

*        |           .         960mV         .   840mV   .

*        |           .           .        840mV          .

*   10uA |.....................840mV......................

*        |          960mV        .           .           .

*        |           .           .           .   720mV   .

*        |          840mV      720mV       720mV         .

*    1uA |................................................

*        |           .           .           .           .

*        |960mV     720mV        .           .           .

*        |840mV      .           .           .   600mV   .

*  100nA |.....................600mV......600mV...........

*        |           .           .           .           .

*        |720mV     600mV        .           .           .

*        |           .           .           .   480mV   .

*   10nA |................................480mV...........

*        |           .         480mV         .           .

*        |600mV     480mV        .           .           .

*        |           .           .           .           .

*    1nA |.......................................360mV....

*        |           .           .        360mV          .

*        |           .         360mV         .           .

*        |           .           .           .           .

*  100pA |480mV.....360mV........................240mV....

*        |           .           .           .           .

*        |           .           .        240mV          .

*        |           .         240mV         .           .

*   10pA |360mV..................................120mV....

*        |           .           .           .           .

*        |          240mV        .           .           .

*        |           .           .           .           .

*    1pA |..240mV..............120mV......120mV...0mV.....

*        |           .           .           .           .

*        |          120mV        .          0mV          .

*        |           .           .           .           .

*  100fA |__120mV__________________0mV___________________.

*        1mV        10mV        100mV       1V         10V

*

*                         Drain Voltage

=====================END_OF_SPICE=======================

==================Match_NPN_TC_to_Silicon======================

1)Sanity checking the temperature behavior of npnv

2)Do a foreach tempval

3)Set temp =  $tempval each time

4)Simulate and measure for each temp simulation

5)echo Temp = $temp C_deg Vbe = $&vbe Beta = $&beta 

6)Note the actual silicon data is included at end of spice file 

=================================================================

Test Match NPN TC to Silicon

*   _____

*  |    0|   Use this to simulate Vbe and Beta over temp.

*  |    _|_  Tweek the model terms IS, XTB, and XTI.

*  |   / _ \ Repeat until simulations match actual silicon.

* _|_  \/ \/ Different simulators often give different results!

* ///  /\_/\ I1

*      \___/ V2 R2

*        |                     __/\  /\  /\__

* _______|_________ C         _|_  \/  \/   _|_

* |                 |C  ____ /   \          ///

* |        R1      _|       |EGbas|

* |___/\  /\  /\_|' npn ____|     |

*       \/  \/   |`-> B      \___/

*                B |          _|_

*                 _|_         ///

*                 ///

I1     0   C DC 1e-6

R1     C   B 1

R2     V2  0 1k

Q1     C   B 0 npnv 1

EGbas  V2  0 C B 1

.control

foreach tempval -55 -35 -15 5 25 45 65 85 105 125

set temp = $tempval

op

let beta = mag(1e-6/v2)

let vbe = mag(b)

echo Temp = $temp C_deg Vbe = $&vbe Beta = $&beta

end

.endc

.model npnv npn (

*==========================================================

+IS=1.1E-17 NF=1.0 BF=120 VAF=30 IKF=6e-03

+ NR=1.0 BR=0.5 VAR=4 IKR=3e-04

+ISE=9E-17 NE=2

+ISC=1E-21 NC=2

+RB=750 RBM=450 IRB=8E-04

+RE=17 RC=110

*==========================================================

+CJE=2E-14 VJE=0.65 MJE=0.3

+CJC=2E-14 VJC=0.65 MJC=0.3

+CJS=3E-14 VJS=0.35 MJS=0.19 XCJC=0.42

+TF=2E-11 XTF=1.25 VTF=1 ITF=0.0035

+TR=6E-09 FC=0.9 PTF=210

*==========================================================

+KF=1.0E-16 AF=1

+XTB=1.4 EG=1.11 XTI=8 TNOM=25

.end

*********Actual Silicon DATA***********************

** Temp NPN_1uA Beta

** -55 8.23E-01 67%

** -35 7.87E-01

** -15 7.43E-01 75%

**   5 6.96E-01

**  25 6.53E-01 100%

**  45 6.09E-01

**  65 5.64E-01 125%

**  85 5.19E-01

** 105 4.74E-01

** 125 4.17E-01 160%

** tweek XTI & IS XTB

* source /Users/don_sauer/Downloads/SI_Lib/Tests.cir

=====================END_OF_SPICE=======================

temp = -55 c_deg vbe = 0.823219 beta = 62.3451

temp = -35 c_deg vbe = 0.782302 beta = 71.2528

temp = -15 c_deg vbe = 0.740177 beta = 80.6774

temp =   5 c_deg vbe = 0.696933 beta = 90.6232

temp =  25 c_deg vbe = 0.652636 beta = 101.095

temp =  45 c_deg vbe = 0.607372 beta = 112.1

temp =  65 c_deg vbe = 0.561192 beta = 123.644

temp =  85 c_deg vbe = 0.514146 beta = 135.734

temp = 105 c_deg vbe = 0.466281 beta = 148.378

temp = 125 c_deg vbe = 0.417638 beta = 161.583

==============X-Y Plot example (CDHW Version)=====================

X-Y Plot example (CDHW Version)

vin 1 0 pulse 0 1 10u

r1  1 2 10

c1  2 3 1u

l1  3 0 1m

.control

tran 1u 1000u

plot v(3)

plot i(vin)

plot i(vin) vs v(3)

plot i(vin) vs v(3) pointplot

.endc

.end

=====================END_OF_SPICE=======================

===================RMS_Noise_Tests======================

1)Create a Vsin sine wave and Vsq square wave having 1V peak magnitude

2)Create Vn TRNOISE noise wave having magnitude 1 and period 1/2 msec

3)The sqrt(mean(vsin*vsin)) will do an RMS of vsin 

4)The RMS values for vsin is printed in $&RmsValsin format

5)The RMS values for vsin, Vsq and Vn are 0.707, 1, and 0.87

6)Sanity test the measure trans statements

7)The noise Vn should have a bandwidth of 1KHz for a 1/2msec sample rate

8)The total sample time is 1 sec for a 1Hz resolution

9)Set the spectrum bandwidth and resolution at 4KHz and 1Hz

10)The RMS of Vn should be spread over 1000 frequency bins 

11)The average value of Vn over 1000 frequency bins  should be 1/sqrt(1k)

12)The bandwidth of Vn should be 1KHz. 

=================================================================

RMS_Tests

VT   Vtime  0  PWL(0 0 1 1 )

Bsin Vsin   0  V = sin(6.283185307179586232*2*v(Vtime) )

Bsq  Vsq    0  V = 2*u(V(Vsin)) -1

Vn   Vn     0  TRNOISE( 1 .5m 0 0)           $  1Vrms,  0.5m period, 1/fslope=0, 1/fmag=0

.control

*destroy    all

tran        .5m   1 

plot        Vsin  Vsq 

plot        Vsq   Vn

let         RmsValsin = sqrt(mean(vsin*vsin))

let         RmsValsq  = sqrt(mean(vsq*vsq))

let         RmsValn   = sqrt(mean(vn*vn))

echo       "RMS level sin $&RmsValsin"

echo       "RMS level sq   $&RmsValsq"

echo       "RMS level nois $&RmsValn"

meas tran   vsin_rms  RMS vsin from=0m to=1

meas tran   vsq_rms   RMS vsq  from=0m to=1

meas tran   vn_rms    RMS vn   from=0m to=1

let         FFT_BandWidth_Hz =  4k

let         FFT_resolution_Hz = 1

echo       "FFT_BandWidth_Hz= $&FFT_BandWidth_Hz"

echo       "FFT_resolution_Hz= $&FFT_resolution_Hz"

set         specwindow= "rectangular"

spec        $&FFT_resolution_Hz $&FFT_BandWidth_Hz $&FFT_resolution_Hz v(vn)

let         Vexpect = 1/ (sqrt(1000)*(1+sqrt(frequency/1000)*sqrt(frequency/1000)*sqrt(frequency/1000)))

plot        Vexpect mag(vn) loglog

let         freq1 = frequency[100]

echo       "Freq at index 100 = $&freq1"

*display

.endc

.end

* source /Users/don_sauer/Downloads/SI_Lib/Tests.cir

=====================END_OF_SPICE=======================

rms level sin          0.706959

rms level sq           0.999943

rms level nois     0.885789

vsin_rms            =  7.07108e-01 from=  0.00000e+00 to=  1.00000e+00

vsq_rms             =  9.99999e-01 from=  0.00000e+00 to=  1.00000e+00

vn_rms              =  8.50116e-01 from=  0.00000e+00 to=  1.00000e+00

fft_bandwidth_hz=      4000

fft_resolution_hz=     1

freq at index 100 =    101

===============1_F_Tests=====================

1)Create a Vsin sine wave and Vsq square wave and Vn noise wave at 1V peak magnitude

2)Create Vn TRNOISE noise wave having magnitude 1 and period 1/2 msec

2)Create Vn2 TRNOISE noise wave having two 1/f values of unity

3)The RMS values for vsin, Vsq,  Vn and Vn2 are 0.707, 1, , 0.87, and 1.938

4)The noise Vn2 should have a bandwidth of 1KHz for a 1/2msec sample rate

5)The total sample time is 1 sec for a 1Hz resolution

6)Set the spectrum bandwidth and resolution at 4KHz and 1Hz

7)The spectrum for Vn2 is plotted with next to a 1/sqrt(frequency) reference 

=================================================================

1_F_Tests

VT   Vtime 0  PWL(0 0 1 1 )

Bsin Vsin  0  V = sin(6.283185307179586232*2*v(Vtime) )

Bsq  Vsq   0  V = 2*u(V(Vsin)) -1

Vn   Vn    0  TRNOISE( 1 .5m 0 0)    $  1Vrms,  0.05m period, 1/fslope=0, 1/fmag=0

Vn2  Vn2   0  TRNOISE( 1 .5m 1 1)    $  1Vrms,  0.05m period, 1/fslope=1, 1/fmag=1

.control

*destroy all

tran  .05m   1 

plot   Vsin  Vsq 

plot   Vsq   Vn

plot   Vn    Vn2

let    RmsValsin = sqrt(mean(vsin*vsin))

let    RmsValsq  = sqrt(mean(vsq*vsq))

let    RmsValn   = sqrt(mean(vn*vn))

let    RmsValn2  = sqrt(mean(vn2*vn2))

echo  "RMS level sin $&RmsValsin"

echo  "RMS level sq    $&RmsValsq"

echo  "RMS level nois $&RmsValn"

echo  "RMS level oneF $&RmsValn2"

let    FFT_BandWidth_Hz = 10k

let    FFT_resolution_Hz = 1

echo  "FFT_BandWidth_Hz= $&FFT_BandWidth_Hz"

echo  "FFT_resolution_Hz= $&FFT_resolution_Hz"

set    specwindow= "rectangular"

spec   $&FFT_resolution_Hz $&FFT_BandWidth_Hz $&FFT_resolution_Hz v(vn2)

let    Vexpect = .7/( sqrt(frequency))

plot   Vexpect mag (vn2) loglog

let    freq1 = frequency[100]

echo  "Freq at index 100 = $&freq1"

*display

.endc

.end

* source /Users/don_sauer/Downloads/SI_Lib/Tests.cir

=====================END_OF_SPICE=======================

No. of Data Rows :  26005

rms level sin       0.707039

rms level sq        0.999981

rms level nois   0.848564

rms level onef   1.93869

fft_bandwidth_hz=   10000

fft_resolution_hz=  1

freq at index 100 = 101

ngspice 25 -> 

================CREATE/SAVE_PWL_WAVEFORM_THD=====================

1)Create a DC voltage source V1 to receive array data into it’s PWL port

2)Create pwl_pars1 array using vector(2*129)*tstep*0.5 

3)Install 129 point sine wave data points into pwl_pars1 in PWL format 

4)Write pwl_pars1 to output file "~output.txt" for viewing

5)Install pwl_pars1 into V1 PWL’s port 

6)Run transient for 128msec and linearize

7)Create a 129 point square wave squar

8)Multiply output V(1)by square wave squar and average to find magnitude magg

9)Subtract a perfect sinewave fund of magnitude magg from output V(1) 

10)Plot this error signal thdd

11)RMS both fundamental fund and distortion thdd and find ration in percent format

12)Print both input_magnitude and THD_%

=================================================================

NON_FFT_THD_VIN

*                          ^  VC

*                         /_\

*                          |

*             __/\  /\  /\_|__/\  /\  /\__

*            |    \/  \/        \/  \/    |

*            |         R1        R2       |

*   VIN     _|VCN1                   VCN2 |_

*    _____|'QN1                          QN2`|__

*   |     |`->1X                        1X<-'|  |

*  _|_       |VEN1                        |     |

* /V1 \      |____________________________|    _|_

* \___/                _|_            VEN1     ///

*   |                 / _ \1I

*  _|_                \/ \/

*  ///                /\_/\

*                     \___/

*                      _|_

*                      ///

.OPTIONS GMIN=1e-15 METHOD=gear ABSTOL=1e-15 temp=27

*====== ====== ====== ====== ====== ====== ======

VCC     VC      0      DC     10

V1      V1      0      0      dc 

VMAG    VMAG    0      DC     1m 

BIN     BIN     0      V=     V(VMAG)*V(V1)

QN1     VCN1    BIN    VEN1   NPN1 1.00

QN2     VCN2    0      VEN1   NPN1 1.00

I1      VEN1    0      100u

R1      VCN1    VC     1

R2      VCN2    VC     1

BOUT    1       0      V = V(VCN1)-V(VCN2)

.control

*================Want_129_1ms_steps=======

let     n =      129

let     tstep =  1ms

*================The_PWL_needs_258_points_.5ms_each=======

let pwl_pars1 =  vector(2*n)*tstep*0.5

*================View_The_PWL_array=======

let ii=          vector(2*129)

*plot pwl_pars1  vs ii

*================Create_PWL_array=======

let index =      0

repeat           $&n

  let pwl_pars1[ 2*index+1] =sin(2*3.1415926*index/128)

  let index =    index + 1

end

*================Write_PWL_array=======

set outfile =   "/Users/don_sauer/Downloads/SI_Lib/output.txt" 

let index =      0

repeat           255

let val =        pwl_pars1[index]

echo            "$&index $&val" >> $outfile

let index =      index + 1

end

*================Install_the_PWL_arrays=======

alter @v1[pwl] = pwl_pars1

*================Run_and_Plot=======

tran 1ms 128m 

linearize

*================Create_SquareWave=======

let squar =      vector(129)

let index =      0

repeat           129

let squar[index] = 1

if (index >63 )

let squar[index] = -1

end

let index =      index + 1

end

let FF =         1.0075

let magg =       FF*3.1415926*mean( v(1)*squar)/2

*plot v(1)       pointplot

*print magg

*================Create_Fundamental=======

let fund=        vector( 129)

let index =      0

repeat 129

let fund[index] =magg*sin(2*3.1415926*index/128)

let index =      index + 1

end

*================Subtract_Fundamental=======

let thdd =       v(1) -fund

*================View=======

plot thdd        pointplot

plot thdd        fund pointplot

let RmsValfund = sqrt(mean(fund*fund))

let RmsValthdd = sqrt(mean(thdd*thdd))

let thdPC =      100*RmsValthdd/RmsValfund

let Vin =        mean(VMAG)

echo            "THD = $&thdPC %   At Vmag = $&Vin"

.endc

.model NPN1 NPN( BF=2100 VAF=216 )

.model PNP1 PNP( BF=2100 VAF=21)

.end

*  source /Users/don_sauer/Downloads/SI_Lib/Tests.cir

=====================END_OF_SPICE=======================

No. of Data Rows : 570

thd = 0.0105797 %   at vmag = 0.001

thd = 16.7357 %     at vmag = 0.1

=============ALTERNATIVE_THD_VS_VIN=============

1)Create a DC voltage source V1 to receive array data into it’s PWL port

2)Create pwl_pars1 array using vector(2*129)*tstep*0.5 

3)Install 129 point sine wave data points into pwl_pars1 in PWL format 

4)Write pwl_pars1 to output file "~output.txt" for viewing

5)Install pwl_pars1 into V1 PWL’s port 

6)Run transient for 128msec and linearize

7)Create a 129 point square wave squar

8)Multply output V(1)by square wave squar and average to find magnitude magg

9)Subtract a perfect sinewave fund of magnitude magg from output V(1) 

10)Plot this error signal thdd

11)RMS both fundamental fund and distortion thdd and find ration in percent format

12)Print both input_magnitude and THD_%

13)Do foreach for VinMag  at levels of 1m 10m 30m 100m 300m 1000m

14)Alter VMAG = $VinMag and transient simulate 

15)Extract values from each run as tran2.fund etc

16)Use vector(6) to create two arrays Vthd and Vinpk

17)Plot Vthd vs Vinpk

18) Use destroy all to clear the plots

========================================================================

THD_VIN_magnitude

*                          ^  VC

*                         /_\

*                          |

*             __/\  /\  /\_|__/\  /\  /\__

*            |    \/  \/        \/  \/    |

*            |         R1        R2       |

*   VIN     _|VCN1                   VCN2 |_

*    _____|'QN1                          QN2`|__

*   |     |`->1X                        1X<-'|  |

*  _|_       |VEN1                        |     |

* /VIN\      |____________________________|    _|_

* \___/                _|_            VEN1     ///

*   |                 / _ \1I

*  _|_                \/ \/

*  ///                /\_/\

*                     \___/

*                      _|_

*                      ///

.OPTIONS GMIN=1e-15 METHOD=gear ABSTOL=1e-15 temp=27

*====== ====== ====== ====== ====== ====== ======

VCC    VC     0      DC     10

V1     V1     0      0      dc 

VMAG   VMAG   0      DC     1m 

BIN    BIN    0      V=     V(VMAG)*V(V1)

QN1    VCN1   BIN    VEN1 NPN1 1.00

QN2    VCN2   0      VEN1 NPN1 1.00

I1     VEN1   0      100u

R1     VCN1   VC     1

R2     VCN2   VC     1

BOUT   1      0      V = V(VCN1)-V(VCN2)

.control

destroy         all

*===============Want_129_1ms_steps=======

let n =         129

let tstep =     1ms

*===============The_PWL_needs_258_points_.5ms_each=======

let pwl_pars1 = vector(2*n)*tstep*0.5

*===============View_The_PWL_array=======

let ii=         vector(2*129)

*plot pwl_pars1 vs ii

*===============Create_PWL_array=======

let index =     0

repeat          $&n

  let           pwl_pars1[2*index+1] =sin(2*3.1415926*index/128)

  let index =   index + 1

end

*===============Write_PWL_array=======

set outfile = "/Users/don_sauer/Downloads/SI_Lib/OpAmpCornerData.txt" 

let index =     0

repeat          255

let val =       pwl_pars1[index]

echo           "$&index $&val" >> $outfile

let index =     index + 1

end

*===============Install_the_PWL_arrays=======

alter @v1[pwl] = pwl_pars1

*===============Loop_VIN=======

foreach VinMag  1m 10m 30m 100m 300m 1000m

alter  VMAG =   $VinMag

*===============Run_and_Plot=======

tran 1ms 128m 

linearize

*===============Create_SquareWave=======

let squar=      vector(128)

let index =     0

repeat 129

let squar[index] = 1

if (index >63 )

let squar[index] = -1

end

let index =     index + 1

end

let FF =        1.008

let magg =      FF*3.1415926*mean( v(1)*squar)/2

*plot v(1)      pointplot

*print magg

*===============Create_Fundamental=======

let fund=       vector(129)

let index =     0

repeat          129

let fund[index] =magg*sin(2*3.1415926*index/128)

let index =     index + 1

end

*===============Subtract_Fundamental=======

let thdd =      v(1) -fund

*===============View=======

*plot thdd      pointplot

*plot thdd      fund pointplot

let RmsValfund = sqrt(mean(fund*fund))

let RmsValthdd = sqrt(mean(thdd*thdd))

let thdPC =     100*RmsValthdd/RmsValfund

let Vin =       mean(VMAG)

echo  "THD =    $&thdPC %   At Vmag = $&Vin"

end

echo            $plots

*set   plot tran3

*plot  tran2.thdd 

*plot  tran4.thdd

*plot  tran2.fund 

*plot  tran4.fund

let   Vthd = vector(6)

let   Vinpk= vector(6)

let   RmsValfund = sqrt(mean(tran2.fund*tran2.fund))

let   RmsValthdd = sqrt(mean(tran2.thdd*tran2.thdd))