Utilisateur:Mathieu Perrin/Travail en cours

%--------------------------------------------------------------------------------------------------%
%            paramètres physiques des binaires (recommandés par Pryor et. al)                      %
%--------------------------------------------------------------------------------------------------%
% Les paramètres sont :  %
% g1 g2 et g3 : paramètres gamma --> calcul de la masse effective des électons et des trous        %
% a0 :paramètre de maille                                                                          %
% Cdilat coefficient de dilatation thermique                                                       %
% Eg0 : gap à 0K                                                                                   %
% alpha et beta : coefficients de dépendance en température du gap : Eg(T)=Ego-alpha*T^2/(beat +T) %
% Evav :offset de la band de valence                                                               %
% Delta0 : spin-orbit splitting                                                                    %
% EKane: énergie de Kane                                                                           %
% C11 et C12 : coefficients d'élasticité                                                           %
% ac et av : potentiels de deformation hydrostatique                                               %
% Cb=ac-av coefficient de deformation                                                              %
% b : potentiel de deformation de cisaillement (expérimentales)                                    %
% bow et Cquat : coefficients d'ajustement des ternaires et quaternaires (empiriques)              %
%--------------------------------------------------------------------------------------------------%




%constantes physiques
m0=9.109e-31;
q=1.602e-19;
h=1.054e-34;
c=3e8;



%---------------------------------------------------------------------------------------------------
% GaAs

g1GaAs=6.98;
g2GaAs=2.06;
g3GaAs=2.93;

a0GaAs=5.653;
CdilatGaAs=6.86E-6; 

Eg0GaAs=1.519;
alphaGaAs=0.5405E-3;
betaGaAs=204;

EvavGaAs=-0.80;
Delta0GaAs=.341; 

EKaneGaAs=28.8;

C11GaAs=1.221;      
C12GaAs=.566;     

acGaAs=-7.17;  
avGaAs=1.16; 
CbGaAs=acGaAs-avGaAs;

bGaAs=-2;


%---------------------------------------------------------------------------------------------------
% GaSb

g1GaSb=13.4;
g2GaSb=4.7;
g3GaSb=6;

a0GaSb=6.09593;    
CdilatGaSb=7.75E-6;

Eg0GaSb=0.812; 
alphaGaSb=0.417E-3;
betaGaSb=140;

EvavGaSb=-0.03;
Delta0GaSb=.76;

EKaneGaSb=27; 

C11GaSb=0.8842;        
C12GaSb=.4026; 

acGaSb=-7.5; 
avGaSb=0.8;  
CbGaSb=acGaSb-avGaSb;

bGaSb=-2;


%---------------------------------------------------------------------------------------------------
% AlAs

g1AlAs=3.76;
g2AlAs=0.82;
g3AlAs=1.42;

a0AlAs=5.6611; 
CdilatAlAs=5.2E-6;

Eg0AlAs=3.099;
alphaAlAs=0.885E-3;
betaAlAs=530;

EvavAlAs=-1.33;
Delta0AlAs=.28; 

EKaneAlAs=21.1; 

C11AlAs=1.250;      
C12AlAs=.534;            

acAlAs=-5.64; 
avAlAs=2.47;
CbAlAs=acAlAs-avAlAs;

bAlAs=-2.3;


%---------------------------------------------------------------------------------------------------
% AlSb

g1AlSb =5.18;
g2AlSb =1.19;
g3AlSb =1.97;

a0AlSb=6.1355;
CdilatAlSb=4.2E-6;

Eg0AlSb=2.386;   
alphaAlSb=0.42E-3;
betaAlSb=140;

EvavAlSb=-0.41;
Delta0AlSb=.676;

EKaneAlSb=18.7;

C11AlSb=.877;        
C12AlSb=0.434;

acAlSb=-4.5;
avAlSb=1.4; 
CbAlAs=acAlSb-avAlSb;

bAlSb=-1.35;


%---------------------------------------------------------------------------------------------------
%InAs

g1InAs=19.67;
g2InAs=8.37;
g3InAs=9.2;

a0InAs=6.0583; 
CdilatInAs=4.52E-6;

Eg0InAs=0.415;   
alphaInAs=0.276E-3;
betaInAs=93;

EvavInAs=-0.59;
Delta0InAs=.39;   

EKaneInAs=21.5; 

C11InAs= 0.833;      
C12InAs=.453;     

acInAs=-5.08; 
avInAs=1.00; 
CbInAs=acInAs-avInAs;

bInAs=-1.8;


%---------------------------------------------------------------------------------------------------
%InSb

g1InSb=34.8;
g2InSb=15.5;
g3InSb=16.5;
                
a0InSb=6.4794;
CdilatInSb=5.37E-6;

Eg0InSb=0.235;    
alphaInSb=0.21E-3;
betaInSb=170;

EvavInSb=-0.0;
Delta0InSb=0.81;    

EKaneInSb=23.3;

C11InSb=.6847;        
C12InSb=.3735;

acInSb=-6.94;
avInSb=.36; 
CbInSb=acInSb-avInSb;

bInSb=-2;

%---------------------------------------------------------------------------------------------------
%InP

g1InP=5.08;
g2InP=1.60;
g3InP=2.10;

a0InP=5.8688;
CdilatInP=4.75E-6;

Eg0InP=1.4236;
alphaInP=0.363E-3;
betaInP=162;

EvavInP=-.94;
Delta0InP=.108;

EKaneInP=20.7;

C11InP=1.011;
C12InP=0.561;

acInP=-6;
avInP=0.6;
CbInP=acInP-avInP;

bInP=-2;


%---------------------------------------------------------------------------------------------------
%GaP

g1GaP=4.05;
g2GaP=0.49;
g3GaP=2.93;

a0GaP=5.4417;
CdilatGaP=4.65E-6;

Eg0GaP=2.9;
alphaGaP=0.5771E-3;
betaGaP=372;

EvavGaP=-1.27;
Delta0GaP=.08;

EKaneGaP=31.4;

C11GaP=1.405;
C12GaP=0.6203;

acGaP=-8.2;
avGaP=1.7;
CbGaP=acGaP-avGaP;

bGaP=-1.6;

%---------------------------------------------------------------------------------------------------
%bowing parameters :

%band gap bowing :

bowGaAsSb=1.3;
bowAlAsSb=0.6;
bowInGaAs=-0.12;
bowInAsSb=0.67;
bowInGaSb=0.415;
%bowAlGaSb=0.368;
%bowAlGaAs=-0.127;
%bowAlGaAs=-0.127+1.31*(1-xGaA);
%bowAlGaSb=-0.44+1.22*(1-xGaA);
bowInGaP=0.0;
bowInAsP=0.0;
bowGaAsP=0.0;


%splitt-off bowing :
bowdeltaAlGaAs=0;
bowdeltaGaAsSb=0.6;
bowdeltaAlAsSb=0.18;
bowdeltaAlGaSb=0.3;
bowdeltaInGaAs=0.15;
bowdeltaInAsSb=1.2;
bowdeltaInGaSb=0.1;
bowdeltaInGaP=0;
bowdeltaInAsP=0.16;
bowdeltaGaAsP=0;

%valence band offset bowing
bowVBOAlGaAs=0;
bowVBOGaAsSb=-1;
bowVBOAlAsSb=-0.93;
bowVBOAlGaSb=0;
bowVBOInGaAs=-0.35;
bowVBOInAsSb=0;
bowVBOInGaSb=0;
bowVBOInGaP=0;
bowVBOInAsP=0;
bowVBOGaAsP=0;
%---------------------------------------------------------------------------------------------------
%paramètre d'ajustement des quaternaires Deltaquat = Cquat*x*(1-x)*y*(1-y)

Cquat=0;