Add files via upload

Untitled.m

@@ -0,0 +1,13 @@
+1310 & 1550
+
+c1= [0,1.61482486874064e-13,2.12799011750962e-13,2.41394344011323e-13,3.33808530344467e-13,3.33500010134680e-13,8.61242244085060e-13,1.28665583690878e-12,1.42636288947139e-12,2.07852687968410e-12,2.14801580716704e-12,2.14621692793392e-12,2.27564151241760e-12,2.29394064832009e-12,2.30026230579460e-12,2.64562416428374e-12,2.69690471147596e-12,3.15525540871802e-12,3.83953134804374e-12,3.89560667691090e-12,4.18723676270278e-12,4.26448560031409e-12,4.30719612441198e-12,4.39728133314245e-12,4.40743817102650e-12,4.41363194356490e-12,5.02351684719302e-12,5.23997384625533e-12,5.50701904882176e-12,6.20962343402480e-12,6.22987261099549e-12,6.26975102581784e-12,6.36513055314661e-12,6.42208529495365e-12,6.48148413392605e-12,6.67097789388911e-12,6.70090527965840e-12,7.37567055595332e-12,7.82127891366256e-12,7.91318228046313e-12,8.37327767189067e-12,8.38771011338391e-12,8.39369217317479e-12,8.50992207081992e-12,8.51366922193141e-12,8.55935626281336e-12,9.08930609397506e-12,9.19658385685518e-12,9.68527635735187e-12,1.02674920446932e-11,1.02787756259879e-11]
+c2= [0,1.94273347803102e-14,1.87500490470211e-14,2.33774920773370e-14,3.55044893444582e-14,3.54139127832416e-14,9.01807651597737e-14,1.35362027710150e-13,1.52535647846142e-13,2.18798022038262e-13,2.28366314926914e-13,2.27611561439293e-13,2.43707836404902e-13,2.41921040419164e-13,2.43265149939295e-13,2.80833955349065e-13,2.88284038463705e-13,3.36566737985195e-13,4.05774258570837e-13,4.13791669230902e-13,4.46175657167444e-13,4.56621302307953e-13,4.55591109082781e-13,4.67462503097529e-13,4.69454575125357e-13,4.70775839024229e-13,5.33150331579293e-13,5.58138466091541e-13,5.88124695953379e-13,6.59011448372637e-13,6.62508534875732e-13,6.68573916685589e-13,6.81025552106030e-13,6.80594286755156e-13,6.89000668652314e-13,7.11281889148646e-13,7.16083690982197e-13,7.84423757063775e-13,8.31511456156572e-13,8.43582784149797e-13,8.92032797890905e-13,8.94688325588769e-13,8.92846274264504e-13,9.07545985756750e-13,9.05763897035995e-13,9.12541732747340e-13,9.67519764862191e-13,9.81261003937951e-13,1.03239185104289e-12,1.09221811215822e-12,1.09438987318130e-12]
+p1= [0,0.0110171195518265,0.00812623424494349,0.00819990530973454,0.0105713771454885,0.0105713771434812,0.243552511679730,0.405801756317468,0.415816892095372,0.908472658174832,0.910423995749640,0.910423818113935,0.921576744415949,0.920802200164854,0.920803043144571,1.06619578068994,1.06707902639150,1.37629272763781,2.21605714658929,2.21806615963144,2.38939427209135,2.39480560034605,2.38729954609074,2.39562611483264,2.39563574245938,2.39563737858368,3.31164174507830,3.43080561293542,3.63760549930117,4.98082615891459,4.98105280625430,4.98289423888171,5.00518762020531,4.97511144644211,4.98101325175625,5.11540135966528,5.11620118596279,6.34482458251298,7.34691553256984,7.38086257977877,8.45862170439967,8.45878324876773,8.45734146355136,8.53442823533377,8.53357770831650,8.53949021968210,9.35466472764296,9.42685819196949,10.2334720417020,9.29420773994546,9.29428947243960]
+p2= [0,1.08876090163498e-05,1.46011884614790e-05,1.47994627842493e-05,1.79778842467287e-05,1.79778844389221e-05,0.000140231170879945,0.000219670257744334,0.000227564322301426,0.000400314908418301,0.000401963685864234,0.000401964674663927,0.000408592104300826,0.000409379965062980,0.000409383278614295,0.000460227994249685,0.000461040672790631,0.000547260816295970,0.000712058565349347,0.000713053981211099,0.000747401043532999,0.000749473878681900,0.000752185866094754,0.000755112384584509,0.000755125342020058,0.000755128459713983,0.000875689646448934,0.000894071654847695,0.000921578988539716,0.00103426313812744,0.00103434513415418,0.00103478468101492,0.00103800738197156,0.00104212507724180,0.00104324247170439,0.00105667997592564,0.00105689895671089,0.00113260811243862,0.00118209225859697,0.00118498391081048,0.00122892143974294,0.00122895507170018,0.00122915568645470,0.00123382588298495,0.00123395627875181,0.00123455283042778,0.00126696139004084,0.00127086594274557,0.00129896491911936,0.00129930546338922,0.00129932272184823]
+bias= [1.20687753883735,-0.916125631699111,0.438951507652081,0.882104601149073,-0.0124476489924152,2.28905602762798,2.04131578631506,1.11356648044584,2.55990176362186,0.748581736068777,-0.115571995253921,1.06701076140600,-0.619836252701880,0.157957370806974,1.82592171879253,0.626226756310116,2.12402673198734,2.62553233023885,0.660253593588904,1.66704765846693,0.795905676892591,-0.850972157772215,0.869348562703181,0.222014553169939,0.155586228670755,2.47106238691760,1.41788469248285,1.58951250311236,2.66235203656177,0.352617994975312,0.538546856255862,0.898141866775971,-0.955952326989885,0.683022415607648,1.31821794720619,0.451578523160691,2.60636052710749,2.09248226217062,0.879338622577838,2.12843034242249,0.282249392639314,-0.441262843285919,1.00453393035054,-0.395586147183687,0.584538731950554,2.29429205010941,0.960116651120547,2.19777798839787,2.41138894122015,0.238807615472093]
+o1= [0,8.83051567131176e-15,6.51339386045120e-15,6.57244318719213e-15,8.47324122348356e-15,8.47324122187464e-15,1.95213845239488e-13,3.25260949720124e-13,3.33288348626077e-13,7.28165107697419e-13,7.29729156898602e-13,7.29729014518838e-13,7.38668382928850e-13,7.38047565017667e-13,7.38048240688432e-13,8.54584404369836e-13,8.55292349397806e-13,1.10313539238150e-12,1.77622901062509e-12,1.77783928825429e-12,1.91516334785957e-12,1.91950067203318e-12,1.91348436900412e-12,1.92015833631633e-12,1.92016605311639e-12,1.92016736451373e-12,2.65436933765831e-12,2.74988236151326e-12,2.91563799562322e-12,3.99226524187830e-12,3.99244690577654e-12,3.99392286322599e-12,4.01179160398155e-12,3.98768474714845e-12,3.99241520178956e-12,4.10013086080061e-12,4.10077194294747e-12,5.08554642109843e-12,5.88874909099138e-12,5.91595855631152e-12,6.77981128978632e-12,6.77994077191340e-12,6.77878514254149e-12,6.84057224969534e-12,6.83989052956699e-12,6.84462956539438e-12,7.49801371299663e-12,7.55587870348736e-12,8.20240125479968e-12,7.44955582209405e-12]
+o2= [0,7.37547707559218e-18,9.89112766745397e-18,1.00254425312661e-17,1.21785667477846e-17,1.21785668779801e-17,9.49953093057737e-17,1.48808884278427e-16,1.54156476397748e-16,2.71181066993056e-16,2.72297980746752e-16,2.72298650578802e-16,2.76788199687670e-16,2.77321911816871e-16,2.77324156480665e-16,3.11767350943350e-16,3.12317875116249e-16,3.70725069103740e-16,4.82362253946355e-16,4.83036567917220e-16,5.06303932715927e-16,5.07708111365183e-16,5.09545264128729e-16,5.11527744395982e-16,5.11536522013613e-16,5.11538633999819e-16,5.93209115336403e-16,6.05661443606532e-16,6.24295443849515e-16,7.00629867763784e-16,7.00685413459320e-16,7.00983171010143e-16,7.03166291013029e-16,7.05955697486415e-16,7.06712642122364e-16,7.15815467562567e-16,7.15963809384830e-16,7.67250656813297e-16,8.00772175178634e-16,8.02731036355527e-16,8.32495168858160e-16,8.32517951796935e-16,8.32653852114514e-16,8.35817533635004e-16,8.35905866251267e-16,8.36309981902731e-16,8.58264167447065e-16,8.60909187021235e-16,8.79943977467997e-16,8.80174668747581e-16]
+
+for i=1:length(p1)
+np1(i)=sqrt((2*beta*nsp*p1(i)*rand(1))/(0.1*taup1));
+end

carrier1.m

@@ -0,0 +1,33 @@
+function c1=carrier1(t, bias,c1,p1,p2)
+q       = 1.6e-10;          %1.6e-19;            %electron charge in A/ns
+vg      = 8.36e4;           %8.36e9;             %group velocity
+d       = 24;               %24e-9;              %thickness of active layer in cm
+l       = 1000;             %1e-6;               %length of cavity in VCSEL cm 
+ao      = 3.5e-2;           %3.5e-16;            %differential gain in cm2
+beta    = 1;                %1;                  %spontaneous emission rate
+ntr     = 1.33e-3;             %1.33e18;            %transparency carrier density
+nsp     = 2;                %2;                  %inversion factor
+gamma_1 = 0.63;             %0.75;               %confinement factor for mode 1
+gamma_2 = 0.35 ;            %0.35;               %confinement factor for mode 2
+ks      = 8.6e-7;           %8.6e-7;             %gain compression coefficient
+r1      = 4000;             %4e-6;               %radius of mode 1
+r2      = 7500;             %7.5e-6;             %radius of mode 2
+alpha   = 1e-6;             %1000;               %material loss of the active layer
+rf      = 0.9991;           %0.9991;             %reflectivity of the front fece
+rb      = 0.9998;           %0.9998;             %reflectivity of the back face
+taup1   = 2;                %2e-12;              %photon lifetime of mode 1
+taup2   = 1.88;             %1.88e-12;           %photon lifetime of mode 2
+taue    = 3000;             %3e-9;               %carrier lifetime
+
+
+v1=pi*(r1^2)*d;
+v2=pi*(r2^2)*d;
+
+g1l=(vg*(d/l)*ao*(c1-ntr));              % modal gain
+% g2l=(vg*(d/l)*ao*(c2-ntr));              % modal gain
+
+g1=((vg*(d/l)*ao*(c1-ntr))*(1-ks*(gamma_1*p1+gamma_2*p2)));             %linear modal gain
+% g2=((vg*(d/l)*ao*(c2-ntr))*(1-ks*((p1*(1-gamma_1))+(p2*(1-gamma_2))))); %linear modal gain
+
+c1=(bias/(v1*q)) - (c1/taue) - ((gamma_1*p1 + gamma_2*p2)*(g1/v1));
+end

carrier11.m

@@ -0,0 +1,36 @@
+function c1=carrier11(t, bias,c1,p1,p2)
+q       = 1.6e-10;          %1.6e-19;            %electron charge in A/ns
+vg      = 8.36e4;           %8.36e9;             %group velocity
+d       = 24;               %24e-9;              %thickness of active layer in cm
+l       = 1000;             %1e-6;               %length of cavity in VCSEL cm 
+ao      = 3.5e-2;           %3.5e-16;            %differential gain in cm2
+beta    = 1;                %1;                  %spontaneous emission rate
+ntr     = 1.33e-3;             %1.33e18;            %transparency carrier density
+nsp     = 2;                %2;                  %inversion factor
+gamma_1 = 0.63;             %0.75;               %confinement factor for mode 1
+gamma_2 = 0.35 ;            %0.35;               %confinement factor for mode 2
+ks      = 8.6e-7;           %8.6e-7;             %gain compression coefficient
+r1      = 4000;             %4e-6;               %radius of mode 1
+r2      = 7500;             %7.5e-6;             %radius of mode 2
+alpha   = 1e-6;             %1000;               %material loss of the active layer
+rf      = 0.9991;           %0.9991;             %reflectivity of the front fece
+rb      = 0.9998;           %0.9998;             %reflectivity of the back face
+taup1   = 2;                %2e-12;              %photon lifetime of mode 1
+taup2   = 1.88;             %1.88e-12;           %photon lifetime of mode 2
+taue    = 3000;             %3e-9;               %carrier lifetime
+eps11   =2e-7;              %self-gain saturation coefficient mode 1
+eps22   =2e-7;              %self-gain saturation coefficient mode 2
+eps12   =0.5e-7;            %cross- gain saturation coefficient mode 1
+eps21   =0.5e-7;            %cross- gain saturation coefficient mode 2
+
+v1=pi*(r1^2)*d;
+v2=pi*(r2^2)*d;
+
+g1l=(vg*(d/l)*ao*(c1-ntr));              % modal gain
+% g2l=(vg*(d/l)*ao*(c2-ntr));              % modal gain
+
+g1=((vg*(d/l)*ao*(c1-ntr))*(1-(p1*eps11)-(eps12*p2))*(1-ks*(gamma_1*p1+gamma_2*p2)));             %linear modal gain
+% g2=((vg*(d/l)*ao*(c2-ntr))*(1-ks*((p1*(1-gamma_1))+(p2*(1-gamma_2))))); %linear modal gain
+
+c1=(bias/(v1*q)) - (c1/taue) - ((gamma_1*p1 + gamma_2*p2)*(g1/v1));
+end

carrier2.m

@@ -0,0 +1,33 @@
+function c2=carrier2(t, bias,c2,p1,p2)
+q       = 1.6e-10;          %1.6e-19;            %electron charge in A/ns
+vg      = 8.36e4;           %8.36e9;             %group velocity
+d       = 24;               %24e-9;              %thickness of active layer in cm
+l       = 1000;             %1e-6;               %length of cavity in VCSEL cm 
+ao      = 3.5e-2;           %3.5e-16;            %differential gain in cm2
+beta    = 1;                %1;                  %spontaneous emission rate
+ntr     = 1.33e-3;             %1.33e18;            %transparency carrier density
+nsp     = 2;                %2;                  %inversion factor
+gamma_1 = 0.63;             %0.75;               %confinement factor for mode 1
+gamma_2 = 0.35 ;            %0.35;               %confinement factor for mode 2
+ks      = 8.6e-7;           %8.6e-7;             %gain compression coefficient
+r1      = 4000;             %4e-6;               %radius of mode 1
+r2      = 7500;             %7.5e-6;             %radius of mode 2
+alpha   = 1e-6;             %1000;               %material loss of the active layer
+rf      = 0.9991;           %0.9991;             %reflectivity of the front fece
+rb      = 0.9998;           %0.9998;             %reflectivity of the back face
+taup1   = 2;                %2e-12;              %photon lifetime of mode 1
+taup2   = 1.88;             %1.88e-12;           %photon lifetime of mode 2
+taue    = 3000;             %3e-9;               %carrier lifetime
+
+
+v1=pi*(r1^2)*d;
+v2=pi*(r2^2)*d;
+
+% g1l=(vg*(d/l)*ao*(c1-ntr));              % modal gain
+g2l=(vg*(d/l)*ao*(c2-ntr));              % modal gain
+
+% g1=((vg*(d/l)*ao*(c1-ntr))*(1-ks*(gamma_1*p1+gamma_2*p2)));             %linear modal gain
+g2=((vg*(d/l)*ao*(c2-ntr))*(1-ks*((p1*(1-gamma_1))+(p2*(1-gamma_2))))); %linear modal gain
+
+c2=(bias/(v2*q)) - (c2/taue) - (((1-gamma_1)*p1 + (1-gamma_2)*p2)*(g2/v2));
+end

carrier22.m

@@ -0,0 +1,37 @@
+function c2=carrier22(t, bias,c2,p1,p2)
+q       = 1.6e-10;          %1.6e-19;            %electron charge in A/ns
+vg      = 8.36e4;           %8.36e9;             %group velocity
+d       = 24;               %24e-9;              %thickness of active layer in cm
+l       = 1000;             %1e-6;               %length of cavity in VCSEL cm 
+ao      = 3.5e-2;           %3.5e-16;            %differential gain in cm2
+beta    = 1;                %1;                  %spontaneous emission rate
+ntr     = 1.33e-3;             %1.33e18;            %transparency carrier density
+nsp     = 2;                %2;                  %inversion factor
+gamma_1 = 0.63;             %0.75;               %confinement factor for mode 1
+gamma_2 = 0.35 ;            %0.35;               %confinement factor for mode 2
+ks      = 8.6e-7;           %8.6e-7;             %gain compression coefficient
+r1      = 4000;             %4e-6;               %radius of mode 1
+r2      = 7500;             %7.5e-6;             %radius of mode 2
+alpha   = 1e-6;             %1000;               %material loss of the active layer
+rf      = 0.9991;           %0.9991;             %reflectivity of the front fece
+rb      = 0.9998;           %0.9998;             %reflectivity of the back face
+taup1   = 2;                %2e-12;              %photon lifetime of mode 1
+taup2   = 1.88;             %1.88e-12;           %photon lifetime of mode 2
+taue    = 3000;             %3e-9;               %carrier lifetime
+eps11   =2e-7;              %self-gain saturation coefficient mode 1
+eps22   =2e-7;              %self-gain saturation coefficient mode 2
+eps12   =0.5e-7;            %cross- gain saturation coefficient mode 1
+eps21   =0.5e-7;            %cross- gain saturation coefficient mode 2
+
+
+v1=pi*(r1^2)*d;
+v2=pi*(r2^2)*d;
+
+% g1l=(vg*(d/l)*ao*(c1-ntr));              % modal gain
+g2l=(vg*(d/l)*ao*(c2-ntr));              % modal gain
+
+% g1=((vg*(d/l)*ao*(c1-ntr))*(1-ks*(gamma_1*p1+gamma_2*p2)));             %linear modal gain
+g2=((vg*(d/l)*ao*(c2-ntr))*(1-(p2*eps22)-(eps21*p1))*(1-ks*((p1*(1-gamma_1))+(p2*(1-gamma_2))))); %linear modal gain
+
+c2=(bias/(v2*q)) - (c2/taue) - (((1-gamma_1)*p1 + (1-gamma_2)*p2)*(g2/v2));
+end