專利名稱:一種基于Abbe矢量成像模型的光源-掩模同步優(yōu)化方法
技術(shù)領(lǐng)域:
本發(fā)明涉及一種基于Abbe (阿貝)矢量成像模型的光源-掩模同步優(yōu)化方法,屬于光刻分辨率增強(qiáng)技術(shù)領(lǐng)域。
背景技術(shù):
當(dāng)前的大規(guī)模集成電路普遍采用光刻系統(tǒng)進(jìn)行制造。光刻系統(tǒng)主要包括照明系統(tǒng)(包括光源和聚光鏡)、掩模、投影系統(tǒng)及晶片四部分。光源發(fā)出的光線經(jīng)過聚光鏡聚焦后入射至掩模,掩模的開口部分透光;經(jīng)過掩模后,光線經(jīng)由投影系統(tǒng)入射至涂有光刻膠的晶片上,這樣掩模圖形就復(fù)制在晶片上。
目前主流的光刻系統(tǒng)是193nm的ArF深度紫外光刻系統(tǒng),隨著光刻技術(shù)節(jié)點(diǎn)進(jìn)入45nm-22nm,電路的關(guān)鍵尺寸已經(jīng)遠(yuǎn)遠(yuǎn)小于光源的波長。因此光的干涉和衍射現(xiàn)象更加顯著,導(dǎo)致光刻成像產(chǎn)生扭曲和模糊。為此光刻系統(tǒng)必須采用分辨率增強(qiáng)技術(shù),用以提高成像質(zhì)量。光源-掩模協(xié)同優(yōu)化(source mask optimization,簡稱SM0)是一種重要的光刻分辨率增強(qiáng)技術(shù)。SMO利用光源及掩模之間的相互作用,通過改變光源明暗圖形、掩模圖形以及在掩模上添加細(xì)小的輔助圖形的方法,達(dá)到提高光刻成像分辨率的目的。較之傳統(tǒng)的分辨率增強(qiáng)技術(shù)(如光學(xué)鄰近效應(yīng)校正(optical proximity correction,簡稱0PC)等),SMO在掩模優(yōu)化過程中引入光源變量,增大了優(yōu)化自由度,從而能夠更為有效的提高光刻系統(tǒng)的分辨率。光源-掩模同步優(yōu)化(simultaneous source mask optimization,簡稱 SISM0)方法是實(shí)現(xiàn)SMO的重要方法之一。其特點(diǎn)是在每一次優(yōu)化迭代中,同時(shí)更新光源和掩模的像素值,從而同步改變光源圖形和掩模圖形。另一方面,為了進(jìn)一步提高光刻系統(tǒng)成像分辨率,目前業(yè)界普遍采用浸沒式光刻系統(tǒng)。浸沒式光刻系統(tǒng)為在投影物鏡最后一個(gè)透鏡的下表面與晶片之間添加了折射率大于I的液體,從而達(dá)到擴(kuò)大數(shù)值孔徑(numerical aperture,簡稱NA),提高成像分辨率的目的。由于浸沒式光刻系統(tǒng)具有高NA (NA > I)的特性,而當(dāng)NA > 0. 6時(shí),電磁場的矢量成像特性對光刻成像的影響已經(jīng)不能被忽視,因此對于浸沒式光刻系統(tǒng)其標(biāo)量成像模型已經(jīng)不再適用。為了獲取精確的浸沒式光刻系統(tǒng)的成像特性,必須采用基于矢量成像模型的SMO技術(shù),對浸沒式光刻系統(tǒng)中的光源和掩模進(jìn)行優(yōu)化。相關(guān)文獻(xiàn)(IEEE International Conference of Electron Devices andSolid-State Circuits EDSSC,2010,I 4)針對部分相干成像系統(tǒng),提出了一種較為高效的基于梯度的SISMO優(yōu)化方法。但是以上方法基于光刻系統(tǒng)的標(biāo)量成像模型,因此不適用于高NA的光刻系統(tǒng)。同時(shí),由于光源面上不同位置光線的入射角度不同,其對投影系統(tǒng)的作用存在差異,但是現(xiàn)有技術(shù)沒有考慮投影系統(tǒng)對光源面上不同點(diǎn)光源入射光線的響應(yīng)差異。因此采用現(xiàn)有方法獲取空間像與實(shí)際存在較大的偏差,進(jìn)而影響SISMO方法的優(yōu)化效果O
發(fā)明內(nèi)容
本發(fā)明的目的是提供一種基于Abbe矢量成像模型的SISMO方法。該方法采用基于矢量模型的SISMO技術(shù)對光源明暗圖形和掩模圖形進(jìn)行優(yōu)化,其可同時(shí)適用于具有高NA的浸沒式光刻系統(tǒng)以及具有低NA的干式光刻系統(tǒng)。實(shí)現(xiàn)本發(fā)明的技術(shù)方案如下一種基于Abbe矢量成像模型的SISMO方法,具體步驟為步驟101、將光源初始化為大小為NsXNs的光源圖形J,將掩模圖形M初始化為大小為NXN的目標(biāo)圖形Z,其中Ns和N為整數(shù);步驟102、設(shè)置初始光源圖形J上發(fā)光區(qū)域的像素值為1,不發(fā)光區(qū)域的像素值為O ;設(shè)定 NsXNs 的變量矩陣 Qs :當(dāng) J(xs, ys) = I 時(shí),Ω8( ) = |;τ ;當(dāng) J(xs, ys) = O 時(shí),
= 淇中J(xs,ys)表示光源圖形上各像素點(diǎn)(xs,ys)的像素值;設(shè)置初始掩模圖形M上開口部分的透射率為1,阻光區(qū)域的透射率為O ; 設(shè)定NXN的變量矩陣ΩΜ :當(dāng)M(x, y) = I時(shí),ΩΜ(x,少)=|疋;當(dāng)M(x, y) = O時(shí),
= ;其中M(x, y)表示掩模圖形上各像素點(diǎn)(X, y)的透射率;步驟103、將目標(biāo)函數(shù)D構(gòu)造為目標(biāo)圖形之與當(dāng)前光源圖形和掩模圖形對應(yīng)的光刻
膠中成像之間的歐拉距離的平方,即,其中之(x,_y)為目標(biāo)圖形
各像素點(diǎn)的像素值,Z(x, y)表示利用Abbe矢量成像模型計(jì)算當(dāng)前光源圖形和掩模圖形對應(yīng)的光刻膠中成像各像素點(diǎn)的像素值;步驟104、計(jì)算目標(biāo)函數(shù)D對于變量矩陣Qs的梯度矩陣VD(Qs),將光源圖形上各像素點(diǎn)的像素值之和Jsum近似為給定常數(shù),得到梯度矩陣ν·0(Ω8)的近似值▽次Qs);計(jì)算目標(biāo)函數(shù)D對于變量矩陣ΩΜ的梯度矩陣VO(Qm);步驟105、利用最陡速降法更新變量矩陣Qs,更新Ω^Ω8-\χ%(Ω8),其中%為預(yù)先設(shè)定的光源優(yōu)化步長,獲取對應(yīng)當(dāng)前Qs的光源圖形j,
J(Wj = |[1 + cosi2s(X,·^)];利用最陡速降法更新變量矩陣ΩΜ,更新,其中、為預(yù)先設(shè)定的掩模優(yōu)化步長,獲取對應(yīng)當(dāng)前Ωμ的掩模圖形M,M(x,>0 = |[1 + cosQm(x,>0]。獲取
對應(yīng)當(dāng)前M的二值掩模圖形,一般情況下tm = O. 5 ;
步驟206、基于光刻膠近似模型,根據(jù)空間像I計(jì)算光源圖形和掩模圖形對應(yīng)的光刻膠中的成像Z。本發(fā)明所述步驟203中針對單個(gè)點(diǎn)光源(xs,ys)獲取該點(diǎn)光源照明時(shí)對應(yīng)晶片位置上的空間像I (xs,ys)的具體過程為設(shè)定光軸的方向?yàn)閦軸,并依據(jù)左手坐標(biāo)系原則以z軸建立全局坐標(biāo)系(x,y,z);步驟301、針對單個(gè)點(diǎn)光源(xs,ys),計(jì)算點(diǎn)光源發(fā)出的光波在掩模上NXN個(gè)子區(qū)域的近場分布E ;其中,E為NXN的矢量矩陣,其每個(gè)元素均為一 3X1的矢量,表示全局坐標(biāo)系中掩模的衍射近場分布的3個(gè)分量;步驟302、根據(jù)近場分布E獲取光波在投影系統(tǒng)入瞳后方的電場分布Eebnt(a,/ ),其中,E=nt(a,#)為NXN的矢量矩陣,其每個(gè)兀素均為一 3X1的矢量,表不全局坐標(biāo)系中入瞳后方的電場分布的3個(gè)分量;步驟303、設(shè)光波在投影系統(tǒng)中傳播方向近似與光軸平行,進(jìn)一步根據(jù)入瞳后方的電場分布E=nt( ,/0)獲取投影系統(tǒng)出瞳前方的電場分布ΕΓ(α’,;0');其中,出瞳前方的電場分布ΕΤ(α',/ ’)為NXN的矢量矩陣,其每個(gè)元素均為一 3X1的矢量,表示全局坐標(biāo)系中出瞳前方的電場分布的3個(gè)分量;步驟304、根據(jù)投影系統(tǒng)出瞳前方的電場分布,獲取投影系統(tǒng)出瞳后方的電場分布;步驟305、利用沃爾夫Wolf光學(xué)成像理論,根據(jù)出瞳后方的電場分布獲取晶片上的電場分布Ewaf'并根據(jù)Ewafw獲取點(diǎn)光源對應(yīng)晶片位置上空間像I (xs,ys)。有益效果較之傳統(tǒng)的光源單獨(dú)優(yōu)化方法和掩模單獨(dú)優(yōu)化方法,本發(fā)明涉及的SISMO方法利用光源及掩模之間的相互作用,在掩模優(yōu)化過程中引入光源變量,增大了優(yōu)化自由度,從而能夠更為有效的提高光刻系統(tǒng)的分辨率。其次,本發(fā)明利用Abbe矢量成像模型描述光刻系統(tǒng)的成像過程,考慮了電磁場的矢量特性,優(yōu)化后的光源圖形和掩模圖形不但適用于小NA的情況,也適用于NA > O. 6的情況。再次,本發(fā)明利用優(yōu)化目標(biāo)函數(shù)的梯度信息,結(jié)合最陡速降法對光源圖形和掩模圖形進(jìn)行優(yōu)化,優(yōu)化效率高。最后,本發(fā)明將光源面柵格化成多個(gè)點(diǎn)光源,針對不同點(diǎn)光源分別計(jì)算其對應(yīng)的空間像,具有精確度高的優(yōu)點(diǎn),該方法可適用于不同形狀的光源,且滿足45nm及以下技術(shù)節(jié)點(diǎn)的光刻仿真需求。
圖I為本發(fā)明基于Abbe矢量成像模型的SISMO方法的流程圖。圖2為點(diǎn)光源發(fā)出光波經(jīng)掩模、投影系統(tǒng)后在晶片位置上形成空間像的示意圖。圖3為初始光源、初始掩模及其對應(yīng)的光刻膠中成像的示意圖。圖4為基于Abbe矢量成像模型的光源單獨(dú)優(yōu)化結(jié)果、初始掩模圖形及其對應(yīng)的光刻膠中成像的示意圖。圖5為初始光源圖形、基于Abbe矢量成像模型的掩模單獨(dú)優(yōu)化結(jié)果及其對應(yīng)的光 刻膠中成像的示意圖。圖6為采用基于Abbe矢量成像模型的SISMO方法優(yōu)化后的光源圖形、掩模圖形及其對應(yīng)的光刻膠中成像的示意圖。
具體實(shí)施例方式下面結(jié)合附圖進(jìn)一步對本發(fā)明進(jìn)行詳細(xì)說明。本發(fā)明的原理當(dāng)光線通過掩模在光刻膠中成像與目標(biāo)圖形相同或近似時(shí),則光刻系統(tǒng)中印制在晶片上的圖形具有很高的分辨率。因此本發(fā)明將SISMO的優(yōu)化目標(biāo)函數(shù)D構(gòu)造為目標(biāo)圖形與光源圖形和掩模圖形所對應(yīng)的光刻膠中成像之間的歐拉距離的平方;如
目標(biāo)圖形的大小為ΝΧΜΡ = Τ^Ρ(υ)-Ζ(υ))2丨為目標(biāo)圖形中各點(diǎn)的像素
值,Z(x,y)為光源圖形和掩模圖形所對應(yīng)的光刻膠中成像的像素值,Z (x,y)與之的值均為O或1,本發(fā)明用圖形或圖像上各像素點(diǎn)的中心坐標(biāo)來表示此像素點(diǎn)的位置。如圖I所示,本發(fā)明基于Abbe矢量成像模型的SISMO方法,具體步驟為步驟101、將光源初始化為大小為NsXNs的光源圖形J,將掩模圖形M初始化為大小為NXN的目標(biāo)圖形Z,其中Ns和N為整數(shù)。步驟102、設(shè)置初始光源圖形J上發(fā)光區(qū)域的像素值為1,不發(fā)光區(qū)域的像素值為
O;設(shè)定 NsXNs 的變量矩陣 Ω5 :當(dāng) J(xs, ys) = I 時(shí),fts(xs,凡)= |;τ ;當(dāng) J(xs, ys) = O 時(shí),
Qs(xs,凡)= 其中J(xs,ys)表示光源圖形上各像素點(diǎn)(xs,ys)的值。設(shè)置初始掩模圖形M上開口部分的透射率為I,阻光區(qū)域的透射率為O ;設(shè)定NXN的變量矩陣QM:iM(x,y)=I 時(shí),Ωμ(χ,>0 = |;γ ;當(dāng) M(x, y) = O 時(shí),ΩΜ0,>0 = ^·7Γ ;其中 M(x, y)表示掩模圖形上各
像素點(diǎn)(X,y)的透射率;步驟103、將目標(biāo)函數(shù)D構(gòu)造為目標(biāo)圖形與當(dāng)前光源圖形和掩模圖形對應(yīng)的光刻膠中成像之間的歐拉距離的平方,即,其中之(A7)為目標(biāo)圖形
各像素點(diǎn)的像素值,Z(x, y)表示利用Abbe矢量成像模型計(jì)算當(dāng)前光源圖形和掩模圖形對應(yīng)的光刻膠中成像各像素點(diǎn)的像素值。
本發(fā)明利用Abbe矢量成像模型計(jì)算當(dāng)前光源圖形和掩模圖形所對應(yīng)的光刻膠中成像的具體步驟為變量預(yù)定義如圖2所示,設(shè)定光軸的方向?yàn)閦軸,并依據(jù)左手坐標(biāo)系原則以z軸建立全局坐標(biāo)系(x,y,z);設(shè)部分相干光源面上任一點(diǎn)光源的全局坐標(biāo)為(xs,ys,zs),由該點(diǎn)光源發(fā)出并入射至掩模的平面波的方向余弦為Us,Y s),則全局坐標(biāo)與方向余弦之間的關(guān)系為as=xs, NAm, Ps=ys- NAm,八=Costsin'1 (NAm ·」x2s+y2s)]其中,NAni為投影系統(tǒng)物方數(shù)值孔徑。設(shè)掩模上任一點(diǎn)的全局坐標(biāo)為(X,y,z),基于衍射原理,從掩模入射至投影系統(tǒng)入瞳的平面波的方向余弦為(α,β,Y),其中(α,β,Y)是掩模(物面)上全局坐標(biāo)系(X,1,z)進(jìn)行傅里葉變換后的坐標(biāo)系。設(shè)晶片(像面)上任一點(diǎn)的全局坐標(biāo)為(xw,yw,zw),從投影系統(tǒng)出瞳入射至像面的平面波的方向余弦為(α ',β ',Υ'),其中(α ',β ',Y ')是晶片(像面)上全局坐標(biāo)系(xw,yw, zw)進(jìn)行傅里葉變換后的坐標(biāo)系。全局坐標(biāo)系與局部坐標(biāo)系之間的轉(zhuǎn)換關(guān)系建立局部坐標(biāo)系(ei, e| |), e I軸為光源發(fā)出光線中TE偏振光的振動(dòng)方向,eN軸為光源發(fā)出光線中TM偏振光的振動(dòng)方向。波矢量為^ =,由波矢量和光軸構(gòu)成的平面稱為入射面,TM偏振光的振動(dòng)方向在入射面內(nèi),TE偏振光的振動(dòng)方向垂直于入射面。則全局坐標(biāo)系與局部坐標(biāo)系的轉(zhuǎn)換關(guān)系為
權(quán)利要求
1.一種基于Abbe矢量成像模型的光源-掩模同步優(yōu)化方法,其特征在于,具體步驟為 步驟101、將光源初始化為大小為NsXNs的光源圖形J,將掩模圖形M初始化為大小為NXN的目標(biāo)圖形之,其中Ns和N為整數(shù); 步驟102、設(shè)置初始光源圖形J上發(fā)光區(qū)域的像素值為1,不發(fā)光區(qū)域的像素值為0 ;設(shè)定 NsXNs 的變量矩陣 Qs :當(dāng) J(xs, Ys) = I 時(shí),i2s(xs,凡)= |;r ;當(dāng) J(xs, Ys) = 0 時(shí),^sixs,ys) = ^淇中J(Xs,ys)表示光源圖形上各像素點(diǎn)(xs,ys)的像素值;設(shè)置初始掩模圖形M上開口部分的透射率為1,阻光區(qū)域的透射率為0 ;設(shè)定NXN的變量矩陣Qm:當(dāng)M(x,y) = I 時(shí),= ;當(dāng) M(x, y) = 0 時(shí),QmO,少) = |;r ;其中 M(x, y)表示掩模圖形上各像素點(diǎn)(x,y)的透射率; 步驟103、將目標(biāo)函數(shù)D構(gòu)造為目標(biāo)圖形與當(dāng)前光源圖形和掩模圖形對應(yīng)的光刻膠中成像之間的歐拉距離的平方,即,其中之(y)為目標(biāo)圖形各像素點(diǎn)的像素值,Z(x, y)表示利用Abbe矢量成像模型計(jì)算當(dāng)前光源圖形和掩模圖形對應(yīng)的光刻膠中成像各像素點(diǎn)的像素值; 步驟104、計(jì)算目標(biāo)函數(shù)D對于變量矩陣Qs的梯度矩陣VD(Qs),將光源圖形上各像素點(diǎn)的像素值之和Jsuffl近似為給定常數(shù),得到梯度矩陣VD(Qs)的近似值▽次Qs);計(jì)算目標(biāo)函數(shù)D對于變量矩陣Qm的梯度矩陣VG(Qm);步驟105、利用最陡速降法更新變量矩陣Qs,更新Qs為%-SnsXV^fls),其中%為預(yù)先設(shè)定的光源優(yōu)化步長,獲取對應(yīng)當(dāng)前Qs的光源圖形J,J(U) = |[l + C0Sfis(xs,凡)];利用最陡速降法更新變量矩陣W,更新 為^-\X VD(Qm),其中% 為預(yù)先設(shè)定的掩模優(yōu)化步長,獲取對應(yīng)當(dāng)前Qm的掩模圖形M,M(x,_y) = |[l + COSQM(x,>0];獲取對應(yīng)當(dāng)前M的二值掩模圖形=,tm為預(yù)定參量;…,當(dāng)M(x,少) 步驟106、計(jì)算當(dāng)前光源圖形J和二值掩模圖形Mb對應(yīng)的目標(biāo)函數(shù)D的值;當(dāng)該值小于預(yù)定閾值或者更新變量矩陣^^和的次數(shù)達(dá)到預(yù)定上限值時(shí),進(jìn)入步驟107,否則返回步驟104 ; 步驟107,終止優(yōu)化,并將當(dāng)前光源圖形J和掩模圖形Mb確定為經(jīng)過優(yōu)化后的光源圖形和掩模圖形。
2.根據(jù)權(quán)利要求I所述基于Abbe矢量成像模型的光源-掩模同步優(yōu)化方法,其特征在于,所述步驟103中利用Abbe矢量成像模型計(jì)算當(dāng)前光源圖形和掩模圖形對應(yīng)的光刻膠中成像的具體步驟為 步驟201、將掩模圖形M柵格化為NXN個(gè)子區(qū)域; 步驟202、將光源圖形J柵格化為NsXNs個(gè)子區(qū)域;步驟203、針對單個(gè)點(diǎn)光源(xs,ys),獲取該點(diǎn)光源照明時(shí)對應(yīng)晶片位置上的空間像I (xs,ys); 步驟204、判斷是否已經(jīng)計(jì)算出所有點(diǎn)光源對應(yīng)晶片位置上的空間像,若是,則進(jìn)入步驟205,否則返回步驟203 ; 步驟205、根據(jù)阿貝Abbe方法,對各點(diǎn)光源對應(yīng)晶片位置上的空間像I (xs,ys)進(jìn)行疊力口,獲取部分相干光源照明時(shí),晶片位置上的空間像I ; 步驟206、基于光刻膠近似模型,根據(jù)空間像I計(jì)算光源圖形和掩模圖形對應(yīng)的光刻膠中的成像。
3.根據(jù)權(quán)利要求2基于Abbe矢量成像模型的光源-掩模同步優(yōu)化方法,其特征在于,所述步驟203中針對單個(gè)點(diǎn)光源(xs,ys)獲取該點(diǎn)光源照明時(shí)對應(yīng)晶片位置上的空間像I(xs, ys)的具體過程為 設(shè)定光軸的方向?yàn)閦軸,并依據(jù)左手坐標(biāo)系原則建立全局坐標(biāo)系;(a,¢, y )是掩模上全局坐標(biāo)系(X,1,z)進(jìn)行傅里葉變換后的坐標(biāo)系,(a ' , ^ / , Y / )是晶片上全局坐標(biāo)系(xw, yw, zw)進(jìn)行傅里葉變換后的坐標(biāo)系; 步驟301、針對單個(gè)點(diǎn)光源(xs,ys),計(jì)算點(diǎn)光源發(fā)出的光波在掩模上NXN個(gè)子區(qū)域的近場分布E ;其中,E為NXN的矢量矩陣,其每個(gè)兀素均為一 3X1的矢量,表不全局坐標(biāo)系中掩模的衍射近場分布的3個(gè)分量; 步驟302、根據(jù)近場分布E獲取光波在投影系統(tǒng)入瞳后方的電場分布Ent(a,灼,其中,E= (a,#)為NXN的矢量矩陣,其每個(gè)元素均為一 3 X I的矢量,表示全局坐標(biāo)系中入瞳后方的電場分布的3個(gè)分量; 步驟303、設(shè)光波在投影系統(tǒng)中傳播方向近似與光軸平行,進(jìn)一步根據(jù)入瞳后方的電場分布奶獲取投影系統(tǒng)出瞳前方的電場分布W);其中,出瞳前方的電場分布Efrt( ',f)為NXN的矢量矩陣,其每個(gè)元素均為一 3X1的矢量,表示全局坐標(biāo)系中出瞳前方的電場分布的3個(gè)分量; 步驟304、根據(jù)投影系統(tǒng)出瞳前方的電場分布Errt(C^f),獲取投影系統(tǒng)出瞳后方的電場分布EebjVf); 步驟305、利用沃爾夫Wolf光學(xué)成像理論,根據(jù)出瞳后方的電場分布E〖xt( ’,f)獲取晶片上的電場分布Ewaf'并根據(jù)Ewafw獲取點(diǎn)光源對應(yīng)晶片位置上空間像I (xs,ys)。
全文摘要
本發(fā)明提供一種基于Abbe(阿貝)矢量成像模型的光源-掩模同步優(yōu)化方法,本方法設(shè)置光源圖形像素值和掩模中開口部分以及阻光部分的透射率,設(shè)置變量矩陣ΩS和ΩM,將目標(biāo)函數(shù)D構(gòu)造為目標(biāo)圖形與當(dāng)前光源和掩模對應(yīng)的光刻膠中成像之間的歐拉距離的平方;利用變量矩陣ΩS和ΩM以及目標(biāo)函數(shù)D引導(dǎo)光源和掩模圖形的同步優(yōu)化過程。較之傳統(tǒng)的光源單獨(dú)優(yōu)化或掩模單獨(dú)優(yōu)化,本發(fā)明所涉及的方法能夠更為有效的提高光刻系統(tǒng)的分辨率。同時(shí)采用本發(fā)明優(yōu)化后的光源和掩模不但適用于小NA的情況,也適用于NA>0.6的情況。另外本發(fā)明利用優(yōu)化目標(biāo)函數(shù)的梯度信息,結(jié)合最陡速降法對光源圖形和掩模圖形進(jìn)行優(yōu)化,優(yōu)化效率高。
文檔編號G03F1/76GK102707582SQ20121019978
公開日2012年10月3日 申請日期2012年6月18日 優(yōu)先權(quán)日2012年6月18日
發(fā)明者李艷秋, 董立松, 韓春營, 馬旭 申請人:北京理工大學(xué)