隧穿效应的设计.docx

上传人:b****8 文档编号:11068233 上传时间:2023-02-24 格式:DOCX 页数:22 大小:116.66KB
下载 相关 举报
隧穿效应的设计.docx_第1页
第1页 / 共22页
隧穿效应的设计.docx_第2页
第2页 / 共22页
隧穿效应的设计.docx_第3页
第3页 / 共22页
隧穿效应的设计.docx_第4页
第4页 / 共22页
隧穿效应的设计.docx_第5页
第5页 / 共22页
点击查看更多>>
下载资源
资源描述

隧穿效应的设计.docx

《隧穿效应的设计.docx》由会员分享,可在线阅读,更多相关《隧穿效应的设计.docx(22页珍藏版)》请在冰豆网上搜索。

隧穿效应的设计.docx

隧穿效应的设计

量子隧穿效应的可视化

1技术指标

此为设计一个用户界面,通过输入粒子质量,势垒高度等参数,进行量子隧穿规律的演示。

它包括以下要求:

有效户任意输入可调参数的界面。

可调参数包括输入的粒子质量,粒子能量,势垒高度,势垒宽度,且反射率和透射率能够直观的显示出来。

依照量子力学中对一维方势垒的求解,设计出各个模块的参数(例如透射系数大小、反射系数大小等);

用Matlab来进行模拟。

编写隧穿效应程序及利用GUI(图形用户接口)设定用户界面。

通过给定的四个参数,能够给出该条件下,量子隧穿规律的动画演示。

而且在有参数改变的情形下,能够观看动画演示的改变。

2大体原理

在量子力学里,量子隧穿效应为一种量子特性,是一种衰减波耦合效应,是如电子等微观粒子能够穿过它们本来无法通过的“墙壁”的现象。

其量子行为遵守薛定谔波动方程。

假假设条件适当,任何波动方程都会显示出出衰减波耦合效应。

依照量子力学,微观粒子具有波的性质,而有不为零的概率穿过位势障壁。

假设要使隧穿效应发生,必需有一个2型介质的薄区域,像三明治一样,夹在两个1型介质的区域。

2型介质的波动方程必需允许实值指数函数解答(上升指数函数或下降指数函数),而1型介质的波动方程那么必需允许行进波解答。

在光学里,1型介质可能是玻璃,而2型介质可能是真空。

在量子力学里,从粒子运动这方面来讲,1型介质区域是粒子总能量大于位能的区域,而2型介质是粒子总能量小于位能的区域(称为位势垒)。

假假设条件适当,从1型介质区域入射至2型介质区域,行进波的波幅会穿透过2型介质区域,再以进行波的形式,显现于第二个1型介质区域。

在量子力学里,穿透过的波幅能够合乎物理地说明为行进粒子。

遵守薛定谔波动方程,穿透波幅的绝对值平方和入射波幅的绝对值平方的比率给出了粒子隧穿的透射系数,也确实是其透射概率。

关于遵守其它种波动方程的光波、微波、绳波、声波等等,穿透波幅能够物理地说明为行进能量,而穿透波幅的绝对值平方和入射波幅的绝对值平方的比率那么给出了穿透能量和入射能量的比率。

隧穿粒子的能量不变,只有量子幅降低。

因此,找到粒子的概率也会降低。

这些"类似隧穿现象"发生的尺寸与行进波的波长有关。

关于电子来讲,2型介质区域的厚度通常只有几纳米。

相较之下,关于一个隧穿出原子核的α粒子来讲,厚度会是超小;关于光波来讲,尽管2型介质区域的厚度超大,类似现象仍旧会发生。

参考以下图:

图1隧穿效应展现图

认真观看薛定谔波动方程。

假假设粒子能够被视为一个局域化(localized)于一点的物体,那么粒子在介质区域内运动的行为是由粒子的动能设定的。

在1型介质区域内,动能是正值的;而在2型介质区域内,动能是负值的。

这现象并非会造成任何矛盾。

量子力学不许诺粒子局域化于一点。

粒子的波函数必是有些散开的("非局域的"),而非局域的物体,其动能的期望值必是正值的。

思路推导:

试探一个入射波,遇处处于x轴坐标0与a之间的位势Uo。

入射波的一部份会反射归去,成为反射波;另一部份那么会穿透过位势,成为透射波。

那么,在位势垒的左侧与右边,波函数是不同的。

在位势垒的内部,依照WKB近似,波函数大约为sin(8*pi*x4)。

一个粒子穿透过位势垒的概率等于透射率,当取所有物理参数都超大于普朗克常数的经典极限,那么,透射系数正确地变成零,也确实是说,粒子无法穿透过位势垒。

隧穿效应的例子

1)冷发射现象:

冷发射是指低温时在外静电场阻碍下电子从金属中飞出的现象。

依照金属势阱中的电子气,在温度T为0K时,电子从低能级开始,在不违背泡利原理的原那么下填充,一直填充到能级Ef。

Ef成为费米能级。

因为温度对电子气的能态散布阻碍是极小的,因此咱们就考虑T为0K的情形。

一样情形下,真空为无穷厚的势垒,要使电子从金属中溢出,必需传递给它能量,能量值不能小于脱出功:

W=Vo-Ef。

2)接触电势差:

两种不同性质的金属相接触会产生电势差的现象可用隧道效应来讲明。

依照前述的金属势阱模型,两种金属应有不同的脱出功和费米能级。

假设使它们相互接触,那么它们为有限宽的势垒分开。

因为金属Ⅰ中部份电子的能级高于金属Ⅱ的费米能级,这些电子有可能隧道跃迁到金属Ⅱ的自由能级上,明白金属的填充能级上限相等才停止。

3)α衰变:

α衰变确实是因为α粒子摆脱了本来不可能摆脱的强力的束缚而“逃出”原子核。

扫描隧道显微镜是量子隧穿效应的要紧应用之一。

扫描隧道显微镜能够克服一般光学显微镜像差的限制,通过穿遂电子扫描物体表面,从而分辨大大小于光波长的物体。

4)扫描隧道显微镜:

一般的显微镜无法观看到很多微小尺寸的物体;可是,扫描隧道显微镜能够清楚地观看到这些物体的细节。

扫描隧道显微镜克服了一般显微镜的极限问题(像差限制,波长限制等等)。

它能够用隧穿电子来扫描一个物体的表面。

3成立模型描述

微观粒子隧穿势垒是一种随机现象,可用蒙特卡罗方式对其进行模拟,运用蒙特卡罗方式处置此问题,第一依照要处置问题的规律,构建一个概率模型。

然后依据概率模型进行随机抽样,得出一组按已知散布的随机序列。

最后依据这一随机序列,借助运算机程序设计语言或图形软件,实现微观粒子势垒隧穿的动态随机进程的模拟。

概率模型的构建

设质量为m,能量为E的粒子,沿x轴正方向射向方势垒。

方势垒的势函数为U(x)=﹛Uo(0a)﹜;粒子的波函数H知足定态薛定谔方程

求解出透射率D和反射率R。

透射率D是透射波概率流密度与入射流概率密度之比,是粒子隧穿势垒的概率。

反射率R是反射波概率流密度与入射流概率密度之比,是粒子被势垒反射的概率。

粒子射向势垒,有两种可能,隧穿或反射,且两种可能的概率之和为1。

数据的生成

依照上述粒子入射方势垒的概率模型进行随机抽样,便取得随机变量的散布函数的随机数序列。

在蒙特卡罗方式中,有多种方式实现按已知散布的随机抽样。

关于此隧穿效应问题的特点,采纳离散型随机变量的直接抽样法。

据此,产生随机序列的步骤如下:

1)运算机产生一个0至1均匀散布的随机数r;2)假设知足1-D

那个随机数序列即是实现势垒隧穿动态随机进程的依据。

势垒隧穿动态随机进程的实现

利用上述生成的只有0和1两种取值的随机数序列,借助运算机程序设计语言或图形软件,就可实现势垒隧穿动态随机进程的模拟。

那个地址采纳编写Matlab程序的方式来实现。

在编写Matlab程序时,用Matlab的rand函数产生0至1之间均匀散布的随机数,利用EraseMode属性的Xor值,可实现对图形对象本身的擦除,从而实现动画模式,程序设计流程图如下

图2方势垒隧穿效应动态随机进程程序流程图

图3GUI界面的图

 

图4入射粒子能量大于势垒高度时的GUI界面显示图

图5入射粒子能量小于势垒高度时的GUI界面显示图

4模型组成模块功能描述

势垒隧穿的动态随机进程Matlab程序的功能描述如下:

functionvarargout=untitledhan(varargin)

%UNTITLEDHANM-filefor

%UNTITLEDHAN,byitself,createsanewUNTITLEDHANorraisestheexisting

%singleton*.

%H=UNTITLEDHANreturnsthehandletoanewUNTITLEDHANorthehandleto

%theexistingsingleton*.

%UNTITLEDHAN('CALLBACK',hObject,eventData,handles,...)callsthelocal

%functionnamedCALLBACKinwiththegiveninputarguments.

%UNTITLEDHAN('Property','Value',...)createsanewUNTITLEDHANorraisesthe

%existingsingleton*.Startingfromtheleft,propertyvaluepairsare

%appliedtotheGUIbeforeuntitledhan_OpeningFunctiongetscalled.An

%unrecognizedpropertynameorinvalidvaluemakespropertyapplication

%stop.Allinputsarepassedtountitledhan_OpeningFcnviavarargin.

%*SeeGUIOptionsonGUIDE'sToolsmenu.Choose"GUIallowsonlyone

%instancetorun(singleton)".

%Seealso:

GUIDE,GUIDATA,GUIHANDLES

%Copyright2002-2003TheMathWorks,Inc.

%Edittheabovetexttomodifytheresponsetohelpuntitledhan

%LastModifiedbyGUIDE21-Jan-202022:

42:

31

%Begininitializationcode-DONOTEDIT

gui_Singleton=1;

gui_State=struct('gui_Name',mfilename,...

gui_Singleton',gui_Singleton,...

gui_OpeningFcn',@untitledhan_OpeningFcn,...

gui_OutputFcn',@untitledhan_OutputFcn,...

gui_LayoutFcn',[],...

gui_Callback',[]);

ifnargin&&ischar(varargin{1})

=str2func(varargin{1});

end

ifnargout

[varargout{1:

nargout}]=gui_mainfcn(gui_State,varargin{:

});

else

gui_mainfcn(gui_State,varargin{:

});

end

%Endinitializationcode-DONOTEDIT

%---Executesjustbeforeuntitledhanismadevisible.

functionuntitledhan_OpeningFcn(hObject,eventdata,handles,varargin)

%Thisfunctionhasnooutputargs,seeOutputFcn.

%hObjecthandletofigure

%eventdatareserved-tobedefinedinafutureversionofMATLAB

%handlesstructurewithhandlesanduserdata(seeGUIDATA)

%varargincommandlineargumentstountitledhan(seeVARARGIN)

%Choosedefaultcommandlineoutputforuntitledhan

=hObject;

%Updatehandlesstructure

guidata(hObject,handles);

%UIWAITmakesuntitledhanwaitforuserresponse(seeUIRESUME)

%uiwait;

%Outputsfromthisfunctionarereturnedtothecommandline.

functionvarargout=untitledhan_OutputFcn(hObject,eventdata,handles)

%varargoutcellarrayforreturningoutputargs(seeVARARGOUT);

%hObjecthandletofigure

%eventdatareserved-tobedefinedinafutureversionofMATLAB

%handlesstructurewithhandlesanduserdata(seeGUIDATA)

%Getdefaultcommandlineoutputfromhandlesstructure

varargout{1}=;

functionedit1_Callback(hObject,eventdata,handles)

%hObjecthandletoedit1(seeGCBO)

%eventdatareserved-tobedefinedinafutureversionofMATLAB

%handlesstructurewithhandlesanduserdata(seeGUIDATA)

%Hints:

get(hObject,'String')returnscontentsofedit1astext

%str2double(get(hObject,'String'))returnscontentsofedit1asadouble

edit1=str2double(get(hObject,'String'));

ifisnan(edit1)

set(hObject,'String',0);

errordlg('Inputmustbeanumber','Error');

end

data=getappdata(gcbf,'metricdata');

=edit1;

setappdata(gcbf,'metricdata',data);%用以实现输入势垒高度的回调函数

%Executesduringobjectcreation,aftersettingallproperties.

functionedit1_CreateFcn(hObject,eventdata,handles)

%hObjecthandletoedit1(seeGCBO)

%eventdatareserved-tobedefinedinafutureversionofMATLAB

%handlesempty-handlesnotcreateduntilafterallCreateFcnscalled

%Hint:

editcontrolsusuallyhaveawhitebackgroundonWindows.

%SeeISPCandCOMPUTER.

Ifispc&&isequal(get(hObject,'BackgroundColor'),get(0,'defaultUicontrolBackgroundColor'))

set(hObject,'BackgroundColor','white');

end

functionedit2_Callback(hObject,eventdata,handles)

%hObjecthandletoedit2(seeGCBO)

%eventdatareserved-tobedefinedinafutureversionofMATLAB

%handlesstructurewithhandlesanduserdata(seeGUIDATA)

%Hints:

get(hObject,'String')returnscontentsofedit2astext

%str2double(get(hObject,'String'))returnscontentsofedit2asadouble

edit2=str2double(get(hObject,'String'));

ifisnan(edit2)

set(hObject,'String',0);

errordlg('Inputmustbeanumber','Error');

end

data=getappdata(gcbf,'metricdata');

=edit2;

setappdata(gcbf,'metricdata',data);%用以实现输入势垒宽度的回调函数

%---Executesduringobjectcreation,aftersettingallproperties.

functionedit2_CreateFcn(hObject,eventdata,handles)

%hObjecthandletoedit2(seeGCBO)

%eventdatareserved-tobedefinedinafutureversionofMATLAB

%handlesempty-handlesnotcreateduntilafterallCreateFcnscalled

%Hint:

editcontrolsusuallyhaveawhitebackgroundonWindows.

%SeeISPCandCOMPUTER.

ifispc&&isequal(get(hObject,'BackgroundColor'),get(0,'defaultUicontrolBackgroundColor'))

set(hObject,'BackgroundColor','white');

end

functionedit3_Callback(hObject,eventdata,handles)

%hObjecthandletoedit3(seeGCBO)

%eventdatareserved-tobedefinedinafutureversionofMATLAB

%handlesstructurewithhandlesanduserdata(seeGUIDATA)

%Hints:

get(hObject,'String')returnscontentsofedit3astext

%str2double(get(hObject,'String'))returnscontentsofedit3asadouble

edit3=str2double(get(hObject,'String'));

ifisnan(edit3)

set(hObject,'String',0);

errordlg('Inputmustbeanumber','Error');

end

data=getappdata(gcbf,'metricdata');

=edit3;

setappdata(gcbf,'metricdata',data);%用以实现输入粒子质量的回调函数

%Executesduringobjectcreation,aftersettingallproperties.

functionedit3_CreateFcn(hObject,eventdata,handles)

%hObjecthandletoedit3(seeGCBO)

%eventdatareserved-tobedefinedinafutureversionofMATLAB

%handlesempty-handlesnotcreateduntilafterallCreateFcnscalled

%Hint:

editcontrolsusuallyhaveawhitebackgroundonWindows.

%SeeISPCandCOMPUTER.

Ifispc&&isequal(get(hObject,'BackgroundColor'),get(0,'defaultUicontrolBackgroundColor'))

set(hObject,'BackgroundColor','white');

end

functionedit4_Callback(hObject,eventdata,handles)

%hObjecthandletoedit4(seeGCBO)

%eventdatareserved-tobedefinedinafutureversionofMATLAB

%handlesstructurewithhandlesanduserdata(seeGUIDATA)

%Hints:

get(hObject,'String')returnscontentsofedit4astext

%str2double(get(hObject,'String'))returnscontentsofedit4asadouble

edit4=str2double(get(hObject,'String'));

ifisnan(edit4)

set(hObject,'String',0);

errordlg('Inputmustbeanumber','Error');

end

data=getappdata(gcbf,'metricdata');

=edit4;

setappdata(gcbf,'metricdata',data);%用以实现输入粒子能量的回调函数

%Executesduringobjectcreation,aftersettingallproperties.

functionedit4_CreateFcn(hObject,eventdata,handles)

%hObjecthandletoedit4(seeGCBO)

%eventdatareserved-tobedefinedinafutureversionofMATLAB

%handlesempty-handlesnotcreateduntilafterallCreateFcnscalled

%Hint:

editcon

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 表格模板 > 表格类模板

copyright@ 2008-2022 冰豆网网站版权所有

经营许可证编号:鄂ICP备2022015515号-1