【摘要】：隨著現(xiàn)代電子科技的高速發(fā)展,電子電路的集成度越來越高,電子電路的尺寸也不斷的向著納米級別甚至更低的尺度發(fā)展。高集成度電路一直是現(xiàn)代工程研發(fā)和仿真運(yùn)算的熱門領(lǐng)域。在對于納米級別或者亞納米級別電路或者器件進(jìn)行研發(fā)設(shè)計時,由于量子效應(yīng)逐漸凸顯,僅使用傳統(tǒng)的麥克斯韋方程組的單物理建模手段已經(jīng)無法滿足實際需求。因此,針對量子效應(yīng)在納米尺寸仿真占比重的問題,能夠同時考量電磁效應(yīng)與量子效應(yīng)的建模方法的研究引起了電子科學(xué)學(xué)科學(xué)者們的廣泛關(guān)注。對于電磁效應(yīng)而言,麥克斯韋方程組能夠系統(tǒng)而準(zhǔn)確的描述電磁場的變化和影響。傳統(tǒng)的電磁效應(yīng)建模的出發(fā)點大多數(shù)是麥克斯韋方程組,由此衍生出諸多經(jīng)典理論。同時,由于單純理論解析對于實際問題來說,求解往往非常困難。隨著計算機(jī)技術(shù)的飛速發(fā)展,研究者們將計算機(jī)運(yùn)算技術(shù)與麥克斯韋方程組相結(jié)合,開創(chuàng)了經(jīng)典的時域有限差分、矩量法、有限元法等對實際問題而言有效且通用的計算電磁學(xué)算法理論,為工程設(shè)計研發(fā)做出了巨大的貢獻(xiàn)。對于量子效應(yīng)而言,薛定諤方程是系統(tǒng)解釋微觀粒子運(yùn)動、易于與計算機(jī)運(yùn)算結(jié)合的理論公式。通過求解薛定諤方程可以得到微觀粒子的本征能量、本征頻率和波函數(shù)分布等參數(shù)。而對于納米器件量子效應(yīng)的研究,其有效切入點之一便是求解器件的電子的本征能量、本征態(tài)等參數(shù),從而推導(dǎo)計算出器件的量子效應(yīng)對于電荷密度、電流大小的影響。如何準(zhǔn)確而高效的將電磁效應(yīng)中麥克斯韋方程組與量子效應(yīng)中薛定諤方程耦合,是研究電磁-量子多物理場問題的關(guān)鍵。電子受到電磁場影響產(chǎn)生運(yùn)動,而電子運(yùn)動本身代表了波函數(shù)的改變和新的電流產(chǎn)生。電子產(chǎn)生的電流同時產(chǎn)生新的電磁場,影響著原先電磁場。解決這個相互作用問題便是這個領(lǐng)域的核心思想�；谶@個思想,當(dāng)前研究者們對電磁-量子多物理問題取得了長足的進(jìn)展,主要分為兩種算法流派。對于第一種算法,因為受到電磁場影響的含時薛定諤方程的求解需求磁矢量、電標(biāo)量的參與,故其核心思想是:1.通過計算麥克斯韋方程組得到當(dāng)前時刻的電場和磁場值2.通過得到的電磁場值,運(yùn)用磁矢勢、電標(biāo)量的定義以及洛倫茲規(guī)范,計算出所需的磁矢勢、電標(biāo)量3.代入磁矢勢、電標(biāo)量,計算受到電磁場影響的含時薛定諤方程得到系統(tǒng)的波函數(shù)4.通過波函數(shù)計算出電子運(yùn)動產(chǎn)生的電流,作為額外項代入下一課時麥克斯韋方程組中計算,從而達(dá)到耦合電磁系統(tǒng)與量子系統(tǒng)的作用。這種算法直接使用傳統(tǒng)麥克斯韋方程組進(jìn)行迭代,在每一時間都需要進(jìn)行電磁場與磁矢勢、電標(biāo)量的置換運(yùn)算,然后代入薛定諤方程中,所以導(dǎo)致算法的耦合度不高、消耗的計算機(jī)內(nèi)存大、仿真所需時間偏長。同時,如果提高仿真空間的精度、大小,計算機(jī)內(nèi)存的消耗、計算性能的要求會急劇增加。如果與辛結(jié)構(gòu)、隱式差分等提高精度的方法相結(jié)合,更是對計算機(jī)性能和仿真時間需求的巨大挑戰(zhàn)。另一種流派的核心思想是通過一系列近似替換,將受電磁場影響的薛定諤方程近似表示為只包含電場的形式,同時忽略掉量子系統(tǒng)對于電磁系統(tǒng)的電流影響項。這樣算法的步驟相對前者會有很大的簡化:1.計算當(dāng)前時刻的電場磁場值2.計算只包含電場項、與磁矢勢電標(biāo)量無關(guān)的薛定諤方程近似形式3.代入當(dāng)前時刻的值計算下一時刻各種值。這種算法的仿真步驟簡單,對于計算機(jī)內(nèi)存和計算性能要求不高,計算速度也相對快。但是由于基于各種近似和忽略電磁效應(yīng)與量子效應(yīng)的相互影響,與實際情況有誤差,而且隨著仿真時間的增加,這個誤差會不斷累積放大,對于精確的求解納米級別的器件的長時間、復(fù)雜仿真情況,往往力不從心。針對上述的各種情況,本文主要通過理論研究和數(shù)值計算方法探索新的解決電磁效應(yīng)與量子效應(yīng)多物理場仿真的方法。本文的主要研究工作為以下幾個方面:(1)簡單介紹了電磁計算學(xué)中常用的數(shù)值方法的種類,包括矩量法、有限元法和時域有限差分法等。以麥克斯韋方程組為例,著重對本論文使用的時域有限差分的實行步驟進(jìn)行介紹,系統(tǒng)分析了該方法的數(shù)值穩(wěn)定性和數(shù)值色散特性。討論了時域有限差分方法的差分方法,分析了不同差分格式下離散格式的精度和計算性能等問題。(2)論文從麥克斯韋方程組出發(fā),通過對磁矢勢、電標(biāo)量的理論研究,推導(dǎo)出了能夠直接迭代的磁矢勢、電標(biāo)量方程組。該方程組能夠避開計算電場、磁場,直接迭代磁矢勢、電標(biāo)量的性質(zhì),更使用于各種需要磁矢勢和電標(biāo)量的建模場合。并通過時域有限差分方法,對提出的磁矢勢和電標(biāo)量方程組進(jìn)行離散差分,給出了其易于計算機(jī)實現(xiàn)的形式結(jié)構(gòu)。(3)分別以復(fù)數(shù)坐標(biāo)延伸和完美匹配吸收理論為基礎(chǔ),推導(dǎo)出適用于磁矢勢和電標(biāo)量迭代方程組的完美匹配吸收邊界。經(jīng)過仿真實驗與理論解對比,證明了所提出的吸收邊界能夠有效吸收邊界磁矢勢、電標(biāo)量,在長時間仿真中能夠保證不產(chǎn)生影響仿真結(jié)果的反射。(4)通過時域有限差分方法對薛定諤方程進(jìn)行定性研究,給出了薛定諤方程的離散差分形式。同時引入辛結(jié)構(gòu)進(jìn)一步耦合薛定諤方程,從理論層面證明了辛結(jié)構(gòu)對薛定諤方程的適用性,并給出了辛結(jié)構(gòu)中辛傳播子的求解方法。然后通過仿真實驗、理論分析值對分別使用普通二階時域有限差分方法、高階時域有限差分方法、引入辛結(jié)構(gòu)的時域有限差分方法的薛定諤離散數(shù)值方程進(jìn)行測試,證明了辛結(jié)構(gòu)對于量子系統(tǒng)仿真中長時間仿真、密集網(wǎng)格劃分的情況具有更穩(wěn)定、仿真結(jié)果更精確的優(yōu)勢。(5)從理論層面推導(dǎo)出了受到電磁效應(yīng)影響的薛定諤方程形式,針對其僅與磁矢勢和電標(biāo)量強(qiáng)相關(guān)的特點,提出了基于電磁勢能和薛定諤修正方程耦合的算法。該算法分為兩個部分:電磁部分使用提出的磁矢勢和電標(biāo)量方程取代傳統(tǒng)的麥克斯韋方程組;量子效應(yīng)方面使用受到電磁效應(yīng)影響的薛定諤方程。通過量子效應(yīng)計算出的電流項、磁矢勢和電標(biāo)量影響薛定諤方程這兩個方面成功耦合了電磁-量子多物理系統(tǒng)。并在此基礎(chǔ)上,對整個算法引入辛結(jié)構(gòu),提出了辛結(jié)構(gòu)下的電磁勢能方程-薛定諤方程耦合算法,使之具有了辛算法的長時間仿真穩(wěn)定性、量子系統(tǒng)能量保一性的優(yōu)點。(6)針對傳統(tǒng)量子態(tài)控制電流計算算法耦合度不高、計算誤差大的問題,以磁矢勢和電標(biāo)量方程組、薛定諤方程為基礎(chǔ),結(jié)合波函數(shù)變換和量子態(tài)控制電流計算理論,提出了基于電磁勢能的新型量子態(tài)控制電流計算算法。值得注意的是,在該算法的電磁勢能部分,我們通過引入洛倫茲規(guī)范,進(jìn)一步改進(jìn)了磁矢量和電標(biāo)量方程組,使之計算結(jié)構(gòu)更簡潔、計算速度更快,電標(biāo)量和磁矢勢處于同一取樣時刻,代入薛定諤方程中更為簡便。傳統(tǒng)的算法將受電磁場影響的薛定諤方程近似表示為只包含電場的形式,同時忽略掉量子系統(tǒng)對于電磁系統(tǒng)的電流影響項,帶來了理論近似誤差,無法通過改變計算精度消除。而提出的算法克服了這個缺點,從理論上來講,不存在近似誤差。最后通過仿真實驗對比提出的算法和傳統(tǒng)量子態(tài)控制電流計算算法,定性定量分析實驗結(jié)果,證明了提出的算法計算更為準(zhǔn)確、控制能力更強(qiáng)的優(yōu)點。綜上所述,本論文主要著重于研究有效解決電磁效應(yīng)-量子效應(yīng)多物理場仿真的算法,通過電磁勢能和辛算法的研究,提出了耦合度高、計算精確高的算法結(jié)構(gòu)。論文的主要創(chuàng)新點如下:(1)論文提出了磁矢勢和電標(biāo)量的迭代方程組,給出了其時域有限差分格式。同時系統(tǒng)的研究和給出了對于這個迭代方程組的完美匹配層算法。使得磁矢勢和電標(biāo)量方程組能夠如同傳統(tǒng)的麥克斯韋方程組仿真一樣,適用于各種需要模擬無限大空間和無反射空間的實際仿真情況。(2)通過耦合新型磁矢勢電標(biāo)量迭代方程組和受電磁效應(yīng)影響的薛定諤方程,構(gòu)建了解決電磁效應(yīng)-量子效應(yīng)多物理問題的計算結(jié)構(gòu)。同時理論證明了辛結(jié)構(gòu)對于提出的算法的可行性,將辛算法與之結(jié)合,創(chuàng)立了一種新型的基于電磁勢能的電磁-量子辛算法。(3)論文建立了不基于忽略量子影響假設(shè)的、使用電磁勢能方程組的量子態(tài)控制電流算法,并改進(jìn)提出了更為簡潔、計算效率更高的電磁勢能方程組形式。該算法能夠準(zhǔn)確的計算激發(fā)帶電粒子從基態(tài)至激發(fā)態(tài)所需要的電流的精確值,并從仿真角度證明所計算的控制電流的控制性能和穩(wěn)定性。
【學(xué)位授予單位】：山東大學(xué)
【學(xué)位級別】：博士
【學(xué)位授予年份】：2019
【分類號】：O413;TN701
【圖文】：

果和理論分析結(jié)果，觀察算法的精確性以及吸收邊界的有效性。那么對于電磁勢逡逑方程，我們同樣使用經(jīng)典的點源仿真從理論和仿真觀察這組電磁勢更新方程的性逡逑能。本小節(jié)會從仿真圖結(jié)果、仿真值與理論對比、計算累計誤差三個方面分析提逡逑出的電磁勢仿真算法以及電磁勢卷積完美匹配層理論的精確性和有效性。逡逑將一個偶極子源放在無限大空間的真空中，其源為Ｊ＝ｉｓｉｎ（２；ｒ／0的ｘ軸方向逡逑電流，頻率為０．１ＧＨｚ。仿真網(wǎng)格劃分為１６８，１６８，邋１６８的三維空間；外層的完美逡逑匹配層的厚度為２０個網(wǎng)格；每個網(wǎng)格劃分的單位是０．３米；仿真時間長為１０００００逡逑個時間步長。那么通過電磁場理論我們可以寫出其理論傳播解：逡逑Ａ邋＝邋ｘ＾－ｅ：ｋｒｅ＇￣－邐（３－６８）逡逑Ａｎｒ逡逑７７逡逑０邋＝邋ｃｏｓ６——ｅｉｋｒｅ＇２邐（３－６９）逡逑４＾ｒ逡逑仿真結(jié)果如下圖：逡逑ｐｈｉ邐ｐｈｉ逡逑

（ａ）邐（ｂ）逡逑圖３－２邋（ａ）為磁矢量ｘ分量在ｘｙ平面的強(qiáng)度圖；（ｂ）為磁矢ｘ分量在ｘｙ平面的強(qiáng)度圖。逡逑２ｘ邋１０邋９邐—一逡逑；—ａｎａｌｙｔｉｃａｌ邋ｉ逡逑：？邋ｓｉｍｕａｌｔｅｄ邋：逡逑１邐Ａ邋／？逡逑0｜＂廣．Ｗ／邐＼邐＼邐！邋Ｉ邋／邋＼邋／＂＇ｖ＇逡逑“Ｖｓ邋！邐Ｖ邋ｖ逡逑－１邐ｉ逡逑－２邐！逡逑■３邐丨ｉ逡逑．３邐Ｉ邋Ｉ逡逑ｉ邋Ｉ逡逑：邋１逡逑５邐１０邐１５邐２０邐２５邐３０邐３５邐４０邐４５邐５０逡逑Ａｘ｛．邋１６８／２，１６８＾）逡逑（ａ）逡逑０．８邐＇邐？邐ｉ邐：￣■邐ｈ逡逑！—ａｎａｌｙｔｉｃａｌ邋：逡逑ｆｖ邐丨、ｓｉｍｕａｌｔｅｄ；逡逑0６＿邐卜邐＾逡逑0．４邐ｉ＇邐！＼逡逑ｉｉ邋Ｍ逡逑°－２邐／，邋｜Ｊ邋ｉ邋Ａ逡逑ａ邋；ｒ－邋；邋？邐ｉ邋；邋＼逡逑ｖ邋＼；邋；邋ｉ邋；．邋■邋Ｊ邋Ｖ邋ｖ逡逑？０．２邐Ｖ邋Ｈ邋；邐Ｖ逡逑＼ｊ邋ｊ逡逑－０．４邐？逡逑４邋ｔ逡逑－０．６逡逑１逡逑？０８邐－邐１５邐２0＂￣邐２５￣邐３０邐３５邐４０邐４５邐５０逡逑ｐｈｉ（．邋１６８／２，１６８／２）逡逑（ｂ）逡逑圖３－３邋（ａ）磁矢量的ｘ分量仿真與理解計算值比較；（ｂ）電標(biāo)量的ｘ分量仿真與理解計算值逡逑比較。逡逑２８逡逑

ｐｈｉ（．邋１６８／２，１６８／２）逡逑（ｂ）逡逑圖３－３邋（ａ）磁矢量的ｘ分量仿真與理解計算值比較；（ｂ）電標(biāo)量的ｘ分量仿真與理解計算值逡逑比較。逡逑２８逡逑

【相似文獻(xiàn)】

相關(guān)期刊論文 前10條

1 陳宵瑋;孫建強(qiáng);王一帆;;三耦合薛定諤方程組的高階保能量方法[J];安徽師范大學(xué)學(xué)報(自然科學(xué)版);2017年03期

2 孫浩翔;;Numerov算法求解一維薛定諤方程研究[J];科學(xué)家;2017年17期

3 孫玉真;姚炳學(xué);;薛定諤方程的群不變解和守恒律[J];青島科技大學(xué)學(xué)報(自然科學(xué)版);2013年05期

4 陳皓,高明,汪青杰;用有限差分法解薛定諤方程[J];沈陽航空工業(yè)學(xué)院學(xué)報;2005年01期

5 Leop.Kadanoff;陳健寧;;論復(fù)雜性[J];世界科學(xué);1987年11期

6 薛立作;王興駿;;“薛定諤方程”教學(xué)探討[J];甘肅農(nóng)大學(xué)報;1988年03期

7 成t ,萬唯實,曾謹(jǐn)言;薛定諤方程的解在r=0鄰域的行為[J];大學(xué)物理;1989年02期

8 竹葉青青;;薛定諤方程[J];第二課堂(C);2019年06期

9 李玉祥;溫學(xué)飛;;薛定諤方程組多峰解的存在性(英文)[J];Journal of Southeast University(English Edition);2012年04期

10 魏益煥;;關(guān)于推廣的薛定諤方程的一個注釋[J];渤海大學(xué)學(xué)報(自然科學(xué)版);2006年01期

相關(guān)會議論文 前10條

1 洪永興;陳文;傅卓佳;;局部近似特解方法在薛定諤方程中的應(yīng)用[A];中國力學(xué)大會-2015論文摘要集[C];2015年

2 王雪瑩;;淺談薛定諤方程引入的一種方法[A];2005年全國高校非物理類專業(yè)物理教育學(xué)術(shù)研討會論文集[C];2005年

3 朱加民;;高階非線性薛定諤方程的精確解研究[A];浙江省光學(xué)學(xué)會第九屆學(xué)術(shù)年會暨新型光電技術(shù)青年論壇論文集[C];2005年

4 胡先權(quán);羅光;馬燕;崔立鵬;;高次正冪與逆冪勢函數(shù)的疊加的徑向薛定諤方程的解析解研究[A];中國數(shù)學(xué)力學(xué)物理學(xué)高新技術(shù)交叉研究學(xué)會第十二屆學(xué)術(shù)年會論文集[C];2008年

5 吳曉飛;;修正高階非線性薛定諤方程的顯式行波解[A];浙江省光學(xué)學(xué)會第九屆學(xué)術(shù)年會暨新型光電技術(shù)青年論壇論文集[C];2005年

6 李學(xué)明;孫志剛;;一種新型精確求解三原子薛定諤方程的傳播子[A];中國化學(xué)會第30屆學(xué)術(shù)年會摘要集-第十八分會：電子結(jié)構(gòu)理論方法的發(fā)展與應(yīng)用[C];2016年

7 黃浪揚(yáng);;廣義高階非線性薛定諤方程的多辛算法[A];全國計算物理學(xué)會第六屆年會和學(xué)術(shù)交流會論文摘要集[C];2007年

8 馮赫;張逸竹;閻天民;江玉海;;網(wǎng)絡(luò)上量子演化的雙節(jié)點分解方法[A];第六屆全國計算原子與分子物理學(xué)術(shù)會議論文集[C];2016年

9 邢藝凡;肖欽文;儲迪陽;孫喜策;吳俊;;基于狀態(tài)轉(zhuǎn)移矩陣的薛定諤方程的解[A];第二十六屆中國控制會議論文集[C];2007年

10 郭永新;吳亞波;;分析力學(xué)與物理學(xué)的現(xiàn)代發(fā)展[A];慶祝中國力學(xué)學(xué)會成立50周年暨中國力學(xué)學(xué)會學(xué)術(shù)大會’2007論文摘要集（下）[C];2007年

相關(guān)重要報紙文章 前2條

1 國家社科基金項目“哲學(xué)視域下的‘多世界解釋’問題研究”負(fù)責(zé)人、山西大學(xué)教授 賀天平;多世界解釋：量子力學(xué)詮釋中的一支奇葩[N];中國社會科學(xué)報;2016年

2 本報記者 劉霞;“潛伏”在日常生活中的7個方程式[N];科技日報;2012年

相關(guān)博士學(xué)位論文 前10條

1 向琛;時域量子—電磁多物理場耦合算法的研究[D];山東大學(xué);2019年

2 田智鯤;依賴于時間的薛定諤方程的兩網(wǎng)格方法[D];湘潭大學(xué);2018年

3 朱露;耗散型薛定諤方程解的大時間性態(tài)[D];南京大學(xué);2014年

4 房祥東;變分法在擬線性薛定諤方程中的應(yīng)用[D];大連理工大學(xué);2014年

5 吳科;廣義擬線性薛定諤方程的駐波解[D];云南大學(xué);2015年

6 白潔靜;含Delta勢薛定諤方程的數(shù)值方法研究[D];南京師范大學(xué);2016年

7 李全清;一般擬線性薛定諤方程解的存在性，多重性及集中現(xiàn)象[D];云南大學(xué);2016年

8 王建云;薛定諤方程有限元超收斂研究[D];湘潭大學(xué);2016年

9 張慧星;全空間上幾類非線性橢圓動力系統(tǒng)解的存在性及相關(guān)動力學(xué)行為[D];中國礦業(yè)大學(xué);2012年

10 李雪陽;半經(jīng)典量子動力系統(tǒng)的數(shù)值方法[D];湘潭大學(xué);2014年

相關(guān)碩士學(xué)位論文 前10條

1 朱思琦;非線性薛定諤方程有限時間爆破解的構(gòu)造[D];中國科學(xué)技術(shù)大學(xué);2019年

2 黎楚欣;分?jǐn)?shù)薛定諤方程下光束傳輸及其相互作用研究[D];深圳大學(xué);2018年

3 王萍;多耦合非線性薛定諤方程的能量守恒格式[D];江西師范大學(xué);2018年

4 應(yīng)東健;三組分耦合非線性薛定諤方程的組合解[D];河北工業(yè)大學(xué);2016年

5 高志云;耦合非線性薛定諤方程在光纖領(lǐng)域中的應(yīng)用[D];東北大學(xué);2017年

6 張康;直波導(dǎo)中的高效率單向聲傳播的研究[D];南京大學(xué);2017年

7 曹萌;二維光滑有界區(qū)域上非線性磁薛定諤方程的集中解[D];華中師范大學(xué);2018年

8 唐志君;薛定諤方程支配系統(tǒng)的能控性[D];東北師范大學(xué);2007年

9 石廣斌;變分法在幾類薛定諤方程中的應(yīng)用[D];曲阜師范大學(xué);2011年

10 劉慧;變系數(shù)非線性耦合薛定諤方程解的相關(guān)研究[D];華北電力大學(xué);2015年

本文編號：2791715