孫蘇珍,王冬嶺

(西北大學(xué)數(shù)學(xué)學(xué)院,陜西 西安 710127)

1 引言

經(jīng)典Schr?dinger方程是量子力學(xué)中的一個(gè)基本方程,也是量子力學(xué)的一個(gè)基本假定,它描述了微觀粒子的狀態(tài)隨時(shí)間變化的規(guī)律.半經(jīng)典量子力學(xué)界于經(jīng)典物理和量子物理之間,運(yùn)用到了量子力學(xué)量子化的概念,但同時(shí)并沒(méi)有摒棄經(jīng)典物理中的一些理念,是經(jīng)典力學(xué)過(guò)渡到量子力學(xué)中的產(chǎn)物.在半經(jīng)典量子力學(xué)中,物質(zhì)運(yùn)動(dòng)狀態(tài)有一條近似的經(jīng)典軌跡,同時(shí)也會(huì)受到不可忽略的量子效應(yīng)擾動(dòng).從數(shù)學(xué)方程角度來(lái)看,由于單個(gè)粒子的狀態(tài)高度集中在某條近似經(jīng)典軌跡的附近,波函數(shù)會(huì)在局部形成強(qiáng)烈的高振蕩,即半經(jīng)典 Schr?dinger方程的解具有局部高振蕩特征.對(duì)于半經(jīng)典Schr?dinger方程的研究,文獻(xiàn)[1-7]在這方面做了很多工作.

整數(shù)階Laplace方程基于經(jīng)典的布朗運(yùn)動(dòng)假設(shè),分?jǐn)?shù)階Laplace方程基于粒子運(yùn)動(dòng)的 Lévy過(guò)程,Lévy過(guò)程有利于刻畫粒子的非局部特征,于是許多學(xué)者開(kāi)始重視研究分?jǐn)?shù)階 Laplace方程.文獻(xiàn) [8]利用 Rieze位勢(shì)給出了 (??)s的定義,至此人們開(kāi)始對(duì)分?jǐn)?shù)階Laplace方程展開(kāi)了豐富的研究并獲得了大量成果.文獻(xiàn)[12]給出了分?jǐn)?shù)階Laplace算子基本解,為之后研究分?jǐn)?shù)階Laplace方程解的性質(zhì)奠定了基礎(chǔ).文獻(xiàn)[9]給出了(??)s的最大值原理,文獻(xiàn)[10]研究了分?jǐn)?shù)階Laplace方程解的正則性.目前,分?jǐn)?shù)階Laplace算子在自然科學(xué)領(lǐng)域已經(jīng)有了重要的應(yīng)用.分?jǐn)?shù)階導(dǎo)數(shù)是一個(gè)非局部算子,與整數(shù)階微分方程相比,能更好地描述具有長(zhǎng)程效應(yīng)的物理現(xiàn)象和動(dòng)態(tài)過(guò)程.

首先簡(jiǎn)單介紹幾種分?jǐn)?shù)階Laplace的定義方法[11]:

(i)定義在Rn上的分?jǐn)?shù)Laplace算子是一個(gè)具有非局部性的擬微分算子,形式如下:

其中α∈(1,2)是任意實(shí)數(shù),Cn,a是依賴于空間維數(shù)n和α的常數(shù),該算子定義在S空間,即由Rn光滑速降函數(shù)構(gòu)成的Schwarz空間.

(ii)在S空間中,也可以用Fourier變換來(lái)等價(jià)地定義分?jǐn)?shù)階Laplace算子:

其中F(u)是u的Fourier變換.當(dāng)α=2時(shí),分?jǐn)?shù)Laplace算子退化成經(jīng)典的Laplace算子,當(dāng)1<α<2時(shí),該分?jǐn)?shù)階Laplace算子為非局部算子.

(iii)還可以將算子延拓到廣義函數(shù)空間

其中該空間內(nèi)積定義如下

(iv)還可以通過(guò)考慮在Rn×[0,+∞)的以下問(wèn)題來(lái)定義Laplace算子:

其中C=C(n,a)是一個(gè)合適的正常數(shù),詳見(jiàn)文獻(xiàn)[12].本文用Fourier變換來(lái)定義分?jǐn)?shù)階Laplace算子.

選擇譜方法(Galerkin譜方法[13]和擬譜方法[1-2])作為求解 Schr?dinger方程的首要數(shù)值方法,主要是因?yàn)樽V方法有高精度,并且Schr?dinger方程本身具有譜算子特征.分裂算法是微分方程計(jì)算中常用到的一種方法,它將原方程分裂成幾個(gè)更簡(jiǎn)單的方程,然后通過(guò)某種組合方式來(lái)構(gòu)造出原方程的數(shù)值算法.分裂算法和擬譜算法的結(jié)合,可以同時(shí)具備顯格式,高精度,保結(jié)構(gòu)的特點(diǎn).由于Schr?dinger方程在空間方向是線性的,可以直接用譜方法進(jìn)行空間離散變成線性常微分方程,并且勢(shì)函數(shù)可以精確求解,所以用分裂譜方法來(lái)求解Schr?dinger方程有很明顯的優(yōu)勢(shì).

本文主要通過(guò)分裂譜算法把文獻(xiàn)[14]中關(guān)于半經(jīng)典整數(shù)階Schr?dinger方程解的高振蕩行為推廣到空間分?jǐn)?shù)階Schr?dinger方程,并通過(guò)數(shù)值比較研究空間分?jǐn)?shù)階階數(shù)α對(duì)解的高振蕩行為的影響.

2 整數(shù)階半經(jīng)典 Schr?dinger方程

考慮d(d≥1)維線性 Schr?dinger（LS）方程:

其中u=u(x,t)是關(guān)于空間變量x和時(shí)間變量t的波函數(shù),i是虛數(shù)單位,~是普朗克常數(shù),V=V(x)是勢(shì)函數(shù),(x,t)∈R+×?,??Rd是波函數(shù)在空間中的定義域.在線性 Schr?dinger方程中,當(dāng)常數(shù)很小時(shí),即 0

半經(jīng)典的 Schr?dinger方程在許多應(yīng)用中起著核心作用,是量子力學(xué)的基本模型[15].具有時(shí)間依賴勢(shì)的Schr?dinger方程在量子物理學(xué)和理論化學(xué)中具有重要意義,考慮具有時(shí)間依賴勢(shì)的半經(jīng)典Schr?dinger方程的初值問(wèn)題[14]

初值條件u(x,0)=u0(x),其中V是一個(gè)實(shí)值的,具有時(shí)間依賴的勢(shì)函數(shù),0<ε?1是半經(jīng)典參數(shù).

(1)式描述了一個(gè)由電子和一個(gè)更大質(zhì)量的粒子或粒子系統(tǒng)組成的系統(tǒng),例如原子核或原子晶格,具有相互作用勢(shì)V,ε是電子和質(zhì)量更大的系統(tǒng)質(zhì)量之比的平方根.由于方程(1)是線性的,可直接對(duì)空間變量半離散化,將其轉(zhuǎn)化為線性O(shè)DE方程,然后通過(guò)計(jì)算矩陣指數(shù)來(lái)求解.但其中半經(jīng)典參數(shù)ε的微小變化會(huì)引起解的快速震蕩,因此給數(shù)值離散帶來(lái)了很大的困難,解決這些困難既需要算法的靈活性也需要對(duì)原方程解的結(jié)構(gòu)進(jìn)行一定的分析和理解.一個(gè)很好的方法是用分裂譜方法,因?yàn)镾chr?dinger方程在空間方向是線性的,可以直接用譜方法進(jìn)行離散變成線性方程,并且分裂譜方法對(duì)整數(shù)階和分?jǐn)?shù)階都有優(yōu)勢(shì),從而避免了分?jǐn)?shù)階算子的非局部特征差分離散導(dǎo)致的滿矩陣帶來(lái)的計(jì)算量增加.

3 空間分?jǐn)?shù)階半經(jīng)典 Schr?dinger方程

4 數(shù)值試驗(yàn)與分析

4.1 一維例子

根據(jù)上面的數(shù)值方法,首先考慮一維的數(shù)值例子,其中V(x)=2x2,初值條件

取x=[?π,π],N=1000,t=[0,0.5],M=500,通過(guò)分裂譜方法對(duì)空間分?jǐn)?shù)階半經(jīng)典 Schr?dinger方程進(jìn)行數(shù)值計(jì)算,下圖為數(shù)值解u(x,t)及Re(u(x,t)),Im(u(x,t))的圖像,其中取各種參數(shù)ε,α進(jìn)行比較分析.

圖 1 取 ε=10?2,α=2時(shí)對(duì)應(yīng)的 Re(u(x,t)),Im(u(x,t),u(x,t)圖像

圖 2 取 ε=10?2,α=1.8時(shí)對(duì)應(yīng)的 Re(u(x,t)),Im(u(x,t),u(x,t)圖像

圖3 取ε=10?2,α=1.6時(shí)對(duì)應(yīng)的Re(u(x,t)),Im(u(x,t),u(x,t)圖像

圖4 取ε=10?2,α=1.2時(shí)對(duì)應(yīng)的Re(u(x,t)),Im(u(x,t),u(x,t)圖像

圖5 取 ε=10?4,α=2時(shí)對(duì)應(yīng)的 Re(u(x,t)),Im(u(x,t),u(x,t)圖像

圖6 取ε=10?4,α=1.6時(shí)對(duì)應(yīng)的Re(u(x,t)),Im(u(x,t),u(x,t)圖像

從圖1-圖6可以看到,通過(guò)與整數(shù)階薛定諤方程的比較,當(dāng)α不斷變小時(shí),波函數(shù)在空間方向的振蕩頻率增高,振蕩寬度變小,峰值變大.說(shuō)明空間分?jǐn)?shù)階算子的階數(shù)α對(duì)于解的高振蕩特征具有明顯而直接的影響.空間分?jǐn)?shù)階算子的非局部特征具有使得解的振蕩聚集的效應(yīng).改變勢(shì)函數(shù)為V(x)=1+cos(x),其他條件和上面例子相同,分別取不同的α,對(duì)應(yīng)計(jì)算結(jié)果的圖像如下,從圖中不難發(fā)現(xiàn)和上面類似的結(jié)果,即空間分?jǐn)?shù)階算子對(duì)于解的高振蕩產(chǎn)生顯著的影響和聚集效應(yīng).

圖7 取ε=10?1,Re(u(x,t))在α=2,1.95,1.9的圖像

圖 8 取 ε=10?1,Re(u(x,t))在 α=1.85,1.8,1.75的圖像

圖 9 取 ε=10?1,Re(u(x,t))在 α=1.7,1.6,1.2的圖像

圖10 取 ε=10?3,α=2時(shí)對(duì)應(yīng)的 Re(u(x,t)),Im(u(x,t),u(x,t)圖像

圖 11 取 ε=10?3,α=1.6時(shí)對(duì)應(yīng)的 Re(u(x,t)),Im(u(x,t),u(x,t)圖像

圖 12 取 ε=10?3,α=1.2時(shí)對(duì)應(yīng)的 Re(u(x,t)),Im(u(x,t),u(x,t)圖像

4.2 二維例子

進(jìn)一步考慮二維空間分?jǐn)?shù)階半經(jīng)典Schr?dinger方程的例子

其中 ?=[?4π,4π]2×(0,T],取T=4,N=250,M=400,同樣用分裂譜方法對(duì)其進(jìn)行數(shù)值計(jì)算,分別取α=2,1.6,1.2,下圖為數(shù)值解及其等高線圖像.

圖13 取 ε=10?1,u(x,t)在α=2,1.6,1.2的數(shù)值解

圖14 取ε=10?1,u(x,t)在α=2,1.6,1.2數(shù)值解的等高線圖

圖15 取 ε=10?2,u(x,t)在α=2,1.6,1.2的數(shù)值解

圖16 取ε=10?2,u(x,t)在α=2,1.6,1.2數(shù)值解的等高線圖

圖17 取 ε=10?6,u(x,t)在α=2,1.6,1.2的數(shù)值解

圖18 取ε=10?6,u(x,t)在α=2,1.6,1.2數(shù)值解的等高線圖

對(duì)于二維情況,從計(jì)算結(jié)果可以看出,空間分?jǐn)?shù)階算子對(duì)于解的振蕩也產(chǎn)生了明顯而直接的影響.當(dāng)參數(shù)ε相對(duì)較大時(shí),空間分?jǐn)?shù)階算子對(duì)于解的振蕩的區(qū)域有著非常敏感的影響,隨著α的變小,振蕩區(qū)域開(kāi)始擴(kuò)散,振蕩的區(qū)域變大.但是當(dāng)參數(shù)ε相對(duì)較小時(shí),空間分?jǐn)?shù)階算子對(duì)于解的振蕩的區(qū)域的敏感性降低,隨著α的變小,振蕩區(qū)域變小,振蕩特征逐步減弱,此時(shí)的分?jǐn)?shù)階算子表現(xiàn)出光滑效應(yīng),即逐步消除振蕩的效果.當(dāng)參數(shù)繼續(xù)變小取ε=10?6時(shí),數(shù)值解基本上不再具有明顯的振蕩特征.