currentz0.f

计算STM诱导金属表面等离激元发光的程序

	real*8 function current(z,w)
	implicit none
	real*8 PI,EL,MASS,H,Vb
	real*8 EFS,EFT,WS,WT,bat
	real*8 temp0,e0,e1,grunch,e,err 
	real*8 bohr,bbmin,fi,z,zz,w,ww
	integer nmax,ip
	common/om/zz,ww
	common/bias/Vb
	common/hyperb/bbmin,fi
	common/fermi/EFS,EFT,WS,WT
	common/int/nmax
	external bat 
	PI=3.14159265358979D0
	bohr=0.529177d0
	EL=1.60219D-19
	MASS=9.10953D-31
        H=1.05459D-34
	ww=w
	zz=z
	ip=0
	temp0=dsqrt(2.d0*MASS/H**2)*(2.d0*MASS/H**2)**3*EL/H/2.d0/pi
	e0=0.d0
	e1=EFS*EL+El*Vb-El*ww
	call qfint0(e0,e1,bat,grunch,err,e,nmax,ip)
	current=temp0*grunch
	return	
	end

	
	real*8 function bat(epar)
	implicit none
	real*8 EL,epar,ess,limit1,limit2 
	real*8 bi,gint,ee,err,zz
	real*8 Vb,EFS,EFT,WS,WT,omega
	integer ip,nmax
	common/om/zz,omega
	common/bias/Vb
	common/fermi/EFS,EFT,WS,WT
	common/int1/ess
	common/int/nmax
	external bi 
	ess=epar
	EL=1.60219D-19
	ip=0
	limit1=EFS*EL-epar
	limit2=EFS*EL+EL*Vb-EL*omega-epar
	call qfint1(limit1,limit2,bi,gint,err,ee,nmax,ip)
	bat=gint
	return
	end

	real*8 function bi(eper)
        implicit none
        real*8 kper,fact1,fact2,fact3,fact4,fact5
        real*8 epar,ess,eper,EFS,EFSEL,EFT,Vb,fis,fit,d,dd,EL,MASS,H,z
        real*8 kappas,kappat,a,zz,oomega,biapoyo
        real*8 bbmin,fi,b,theta,bohr,WS,WT,TEM
        real*8 omega,pi,r
	common/fermi/EFS,EFT,WS,WT
	common/bias/Vb
	common/sph/a(10),dd
	common/hyperb/bbmin,fi
        common/om/z,omega
        common/int1/epar
        common/int2/ess
        	if(eper.lt.0.d0) then
		biapoyo=0.d0
		else
	ess=eper
	zz=z
	oomega=omega
	EL=1.60219d-19
        MASS=9.10953d-31
        H=1.05459d-34
        PI=3.14159265358979D0
	bohr=0.529177d0
	b=(bbmin*bohr)*1.d-10
	theta=fi
	R=b*(dtan(theta))**2
	d=(dd+bbmin)*bohr*1.d-10
	fis=WS*EL
	fit=WT*EL
	EFSEL=EFS*EL
	kper=dsqrt(2.d0*MASS/H**2*eper)
	fact1=kper*dsqrt(epar+eper+El*omega)
	fact2=kper**2+(kappas(epar+eper))**2
	fact3=R**2*kappas(eper)/kappat(eper)
	fact4=dexp(2.d0*(kappat(eper+epar)-kappat(eper))*R)
	fact5=(TEM(zz,oomega))**2
	biapoyo=fact1/fact2*fact3*fact4*fact5
        endif
	bi=biapoyo
	return
        end