檢視 2D DDCC 的原始碼

2D DDCC is named from two dimensional Drift-diffusion Charge Control solver. This is 2D finite element based Poisson and drift-diffusion solver developed by Dr. Yuh-Renn Wu. This solver initially developed with the thermal solver. Then the Poisson and drift-diffusion solver was added into this project. This solver was initially developed to solves AlGaN/GaN HEMT structure. Therefore, the 1D Schrodinger cross section solver was added into the program for obtaining the confined state information. The electric field distribution was then used in Monte Carlo program for high field transport. After Dr. Wu returned NTU, the program was then modified to solve LED based current spreading problem. The mesh algorithm was then improve gradually in order to deal with certain problem. After years development, the 2D FEM based Schrodinger eigen solver was added. It also accept additional module to read in the optical field from 2D FD-TD program so that it can consider the solar cell problem. Then the 2D ray tracing program was added into this project to solve the light extraction problem. This solver now can solve many different problems such as trap problem, Gaussian shape tail state models, field dependent mobility, thermal, light extraction. Recently, localization landscape model was also added into this program so that it can calculate the effective quantum potential very efficiently. This code is written with Fortran language.

'''Command Manual'''
== parameters assigned command ==

[[$totalregion]]<br>
[[$ranges]]<br>
[[$gmshlineboundary]]<br>
[[$gmshsurfparameter]]<br>
[[$triregionnumber]]<br>
[[$triregiontype]]<br>
[[$triranges]]<br>
[[$assignparbyfunc]]<br>
[[$parameters]]<br>

[[$opertemp]]<br>
[[$Schottky]]<br>
[[$gatebias]]<br>
[[$drainbias]]<br>
[[$sourcebias]]<br>
[[$addrefbias]]<br>
[[$basebias]]<br>
[[$2ndgateaddbarrier]]<br>
[[$affinity]]<br>
[[$gatewithworkfunc]]<br>
[[$adddrainworkfunc]]<br>
[[$addsourceworkfunc]]<br>
[[$addbaseSB]]<br>
[[$addvrefSB]]<br>
[[$addbaseworkfunc]]<br>
[[$addvrefworkfunc]]<br>
[[$addsourceSB]]<br>
[[$adddrainSB]]<br>
[[$addbaseSB]]<br>
[[$addvrefSB]]<br>

[[$infile]]<br>
[[$checkmesh]]<br>
[[$outfile]]<br>
[[$maxiterations]]<br>
[[$errors]]<br>

[[$call2Dexciton]]<br>
[[$2Dexcitonpar]]<br>
[[$solveeigenvalue]]<br>
[[$solveeigenbias]]<br>
[[$eigensetting]]<br>
[[$usefeast]]<br>
[[$eigencalEl]]<br>
[[$assignsumrange]]<br>
[[$DoEalloyfluc 2D]]<br>
[[$polcharge]]<br>
[[$affinity]]<br>
[[$usefieldmufunc]]<br>
[[$gaussiantraps]]<br>
[[$traps (2D)]]<br>
[[$multitraps]] <br>

[[$recombine]]<br>
[[$UseAuger]]<br>

[[$landscape2D]]<br>
[[$landscapeDOS  (2D)]]<br>
[[$landshiftEcv]]<br>

[[$generation]]<br>
[[$usegenfunc]]<br>
[[$usegenspecfunc]]<br>
[[$usespecgen]]<br>
[[$ifimpact_ion]]<br>
[[$ifschockley]]<br>


[[$ifuseconstvb]]<br>
[[$iftrialmod]]<br>
[[$usepardiso]]<br>
[[$useevprofile]]<br>
[[$xend]]<br>
[[$yend]]<br>

[[$useevversusT]]<br>
[[$usedynataunrfunc]]<br>
[[$useemufunc]]<br>
[[$useemunfunc]]<br>
[[$useemupfunc]]<br>
[[$useemuTfunc]]<br>
[[$useconstantmu]]<br>
[[$useresistor]]<br>
[[$useconstIg]]<br>

[[$iffermiboltzman]]<br>


[[$outputwave]]<br>
[[$ifgainwithpol]]<br>
[[$ifcalculategain]]<br>


[[$thermaldiff]]<br>
[[$heatiterative]]<br>
[[$heatbndistance]]<br>
[[$maxTlimit]]<br>
[[$heatsolvebytime]]<br>
[[$heatsource]]<br>
[[$timeevolve]]<br>
[[$usechannelaveT]]<br>

[[$usemuTfactor]]<br>


$[[solvetimestep2D]]<br> 
$[[savetimestep2D]]<br>
[[$usetaunrbyfunc]]<br>
[[$usemubydopeT]]<br>

[[$usegenfile]]<br>

[[$DOEreftoEc]]<br>

[[$callpoissonexternal]] <br>
[[$call2DNEGF]] <br>
[[$Trapmemorymode]]<br>
[[$RandomizedR]]<br>

[[$usenonRAugRasheat]]<br>
[[$$activateStimEL]]<br>

== time dependent ion move ==
[[$Withionmove]]

[[$IonMovewithPoisson]]

[[$IonMovewithPatDt]]

== Ray Tracing parameters assigned command ==

[[$element_total/$photon_element_pass]]<br>
[[$trifilename]]<br>
[[$fetfilename]]<br>

[[$totalparameterregion]]<br>
[[$ranges (RayTracing)]]<br>
[[$parameters (RayTracing)]] <br>
[[$pattern_tri_region]]<br>
[[$pattern_tri_period]]<br>
[[$Metalregion]]<br>
[[$setradseparately]]<br>
[[$tracephase]] <br>
[[$Metaltri]]<br>
[[$diffusortype]]<br>

'''gmsh setting''' <br>

[[$setradseparate]]<br>
[[$gmshfilename]]<br>
[[$gmshdiffusor]]<br>
[[$gmshparameterregion]]<br>

'''Output setting''' <br>

[[$Light_intensity_region_TBLR]]<br>
[[$photon_total]]<br>
[[$initialanglerange]]<br>
[[$consider3dipoles]]<br>
[[$TETMratio]]<br>
[[$setcountedexitlength]]<br>
[[$traceraypath]]<br>

== 2D RCWA commands ==

[[$RCWAtotallayer]]<br>
[[$RCWAfourier]]<br>
[[$RCWAanglelength]]<br>
[[$RCWAangle2]]<br>
[[$RCWAwidth]]<br>
[[$RCWAoxide2]]<br>
[[$RCWAperiod]]<br>
[[$RCWAlamdalength]]<br>
[[$RCWAsourcelength ]]<br>
[[$RCWAsourceposition ]]<br>
[[$RCWAsourceintensereal ]]<br>
[[$RCWAsourceintenseimage ]]<br>
[[$RCWAlamda2]]<br>
[[$RCWArefractiveindexreal]]<br>
[[$RCWAabsorption2]]<br>
[[$RCWAdielectric]]<br>
[[$RCWAdielecabsorption]]<br>
[[$RCWAoxidedielectric]]<br>
[[$RCWAoxideabsorption]]<br>
[[$RCWAthickness]]<br>
[[$RCWAcavitystart]]<br>
[[$RCWAcavityend]]<br>
[[$$RCWArcaa]]<br>

== phasing out commands ==

[[$other]]<br>
[[$useschrodn]]<br>

[[$usedvprofile]]<br>
[[$fermitable]]<br>
[[$ifpulsemode]]<br>
[[$ifdynarecombine]]<br>
[[$initial]]<br>
[[$ifshiftef]]<br>
[[$ifshiftgateef]]<br>
[[$ifshiftallef]]<br>

== Mesh setting input command ==

[[$xsecysecmaxob]]<br>
[[$xnode]]<br>
[[$ynode]]<br>
[[$xdiv]]<br>
[[$ydiv]]<br>
[[$ygradualdiv]]<br>
[[$xgradualdiv]]<br>
[[$xd]]<br>
[[$contact]]<br>
[[$bottomcontactnum]]<br>
[[$areacontact]]<br>
[[$emptysection]]<br>
[[$emptyperiodshape]]<br>
[[$emptyperiodtriangle]]<br>

[[$deformboundary]]<br>
[[$deformafter]]<br>
[[$deformbeforemod]]<br>

[[$emptytriangle]]<br>
[[$emptysineregion]]<br>
[[$emptysinesum]]<br>

== External function ==
[[function getmunbyET(ind,Ex,Ey,T,xp,yp)]]<br>
[[function getmupbyET(ind,Ex,Ey,T,xp,yp)]]<br>
[[function usefielddependent_n(indEl, matind ,Ex,Ey,T,xp,yp)]]<br>
[[function usefielddependent_p(indEl, matind, Ex, Ey, T, xp, yp)]]<br>
[[subroutine parameter_mod]]<br>
[[subroutine updatetaunp]]<br>
[[subroutine callgeneration(x,y,z,generation,ind) ]]<br>
[[subroutine callspecgeneration(nnode,x,y,eg,generation) ]]<br>
[[subroutine callspecgen(Annode,x,y,Eg,generation)]]<br>
[[subroutine exciton2Dquench(  )]]<br>


== Execution files ==

[[2D-ddcc-win64.exe & 2D-ddcc-win32.exe]]<br>
[[2D-ddcc.exe (Linux version)]]<br>
[[2D-ddcc-dyna.exe (Linux version with all dynamic link modules]]<br>
[[raytrace2D-win32.exe | raytrace2D-win64.exe | raytrace2D.exe ]] <br>
[[ddcc_2d.m & ddcc_2d.fig]] <br>
[[resultopen.m resultopen.fig ]] <br>
[[Material.mat]] <br>
[[load_gmsh.m]] <br>
[[createTable.m]] <br>
[[libiomp5md.dll and libiomp5md.lib]] <br>

== output file format ==

[[*.tri]]<br>
[[*.fet]]<br>
[[*.e]]<br>
[[*.jn]]<br>
[[*.iv (2D)]]<br>
[[*.QE]]<br>
[[*.overlap]]<br>
[[*.wv & *.wvp]]<br>
[[*.timeiv]]<br>
[[*.IQEALL]]<br>
[[*.Nsum]]<br>
[[*.time_ion]] <br>

== output file format (Ray Tracing)==

[[*.photon.out]]<br>
[[*.window.out]]<br>
[[*.raytrace.out]]<br>
[[*.angle.out]]<br>
[[*.absorption.out]]<br>
[[*.info.out]]<br>
[[*raypathsave.way]]<br>
[[*.phasesaved]]<br>
[[*.angle.4surf.out]]<br>
[[*.E-field-phasesum.out]]<br>

== output file format (RCWA)==

[[*.rcwa2d.dat]]<br>
[[*.rcwatopreflectivity.dat]]<br>
[[*.rcwatoptransmission.dat]]<br>
[[*.rcwabottomreflectivity.dat]]<br>
[[*.rcwabottomtransmission.dat]]<br>
[[*.rcwatopemmision.dat]]<br>
[[*.rcwabottomemmision.dat]]<br>
[[*.rcwaef2d.dat]]<br>
[[*.x.dat]]<br>
[[*.z.dat]]<br>
[[*.degree.dat]]<br>
[[*.farfield.dat]]<br>

== Frequently used global variables used in the program ==

Region based variables <br>

[[par%xxxx]] <br>
affinitypar(:) <br>

Node based variables <br>

[[me%xxxxx]]<br>

[[sprsvc(:)  ]]       <br> 
[[sprst(:)    ]]       <br>        
[[sprsefn(:) ]]        <br>
[[sprsefp(:)      ]]   <br>
[[gaussiantrapNt(:,:)]]<br>
[[gaussiantrapEt(:,:)]]<br>
[[gaussiantrapdegenercy(:,:)]]<br>
[[gaussiantraptaun(:,:)]]<br>
[[gaussiantraptaup(:,:)]]<br>
[[gaussiantrapsigma(:,:)]]<br>

[[affinity(:)]]<br>

Element based variables <br>

[[affinityEl(:)]]<br>