• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    磁斯格明子拓撲特性及其動力學微磁學模擬研究進展?

    2018-08-02 05:48:12孔令堯
    物理學報 2018年13期
    關鍵詞:明子鐵磁霍爾

    孔令堯

    (中國工程物理研究院流體物理研究所,綿陽 621900)

    具有非平庸拓撲性的新型磁結構斯格明子,由于其拓撲穩(wěn)定性、尺寸小、低電流驅動等方面的顯著優(yōu)勢,有望應用于自旋電子學儲存器件.拓撲和凝聚態(tài)物理學的結合,使得斯格明子展現(xiàn)出很多有趣的拓撲物理現(xiàn)象,吸引了眾多的研究興趣,同時這些性質也是其電流驅動下動力學特點的重要影響因素.本文從斯格明子的拓撲物理學基礎及其自旋電子學器件應用相關動力學兩個方面介紹了相關研究進展.在拓撲物理基礎方面,介紹了斯格明子的拓撲霍爾效應、斯格明子霍爾效應以及自旋軌道轉矩等拓撲性質,由此討論了斯格明子的動力學性質及其計算方法;在動力學方面,從非均勻電流驅動生成斯格明子、電流驅動下的穩(wěn)定輸運、產生湮滅過程的人工控制幾個賽道存儲應用關心的問題簡要地介紹了相關微磁學模擬研究最新進展.

    1 引 言

    1962年,Skyrme[1]提出一種受拓撲保護的“粒子”,被稱為斯格明子 (skyrmion).1975年,Polyakov[2]在理論上預測了一種二維Heisenberg自旋系統(tǒng)中的亞穩(wěn)態(tài),其自旋結構與Skyrme提出的斯格明子一致,首次在磁性系統(tǒng)中提出磁斯格明子概念.直到2009年,Mühlbauer等[3]在體手性B20材料MnSi單晶樣品中成功地觀測到磁性斯格明子在倒空間中形成的中子衍射點,證實了磁性斯格明子的存在,隨后Yu等[4]利用洛侖茲透射電子顯微鏡在Fe0.5Co0.5Si單晶樣品中實現(xiàn)實空間中的斯格明子磁成像.隨著研究的深入,存在斯格明子的材料從體手性B20結構材料拓展至重金屬/鐵磁/絕緣體多層材料的界面手性材料中[5?8],界面手性材料中的斯格明子被稱為奈爾(Néel)型斯格明子,相對的體材料中的斯格明子稱為布洛赫(Bloch)型.

    作為實空間存在的拓撲性的準粒子,斯格明子展現(xiàn)了眾多的拓撲物理性質,如emergent場[9]、拓撲霍爾效應[10?13]、斯格明子霍爾效應[14,15]、自旋軌道轉矩[9,16]等,激起了眾多科技工作者的研究興趣.有著拓撲保護的斯格明子顯示出特別的穩(wěn)定性,并且可以被小到102A/cm2的電流驅動[17],比磁疇的驅動電流小了4個數(shù)量級[9].正是由于斯格明子的穩(wěn)定性、小尺寸、低電流驅動等特性,使之有潛力作為斯格明子賽道存儲器件[18,19].賽道存儲器是一種新型非易失性存儲器,把斯格明子和鐵磁態(tài)分別作為信息的“1”和“0”約束在賽道內,可以實現(xiàn)三維信息存儲,并通過電流驅動來進行讀取寫入[20,21].

    本文分別介紹了斯格明子的拓撲物理基礎及其動力學微磁學模擬相關研究進展.在拓撲物理基礎方面,介紹了斯格明子的結構能量、拓撲荷、拓撲霍爾效應、斯格明子霍爾效應以及自旋軌道轉矩等性質,由此討論了斯格明子的動力學計算方法;在動力學方面,從非均勻電流驅動生成斯格明子、電流驅動下的穩(wěn)定輸運、產生湮滅過程的人工控制等幾個賽道存儲應用關心的問題,簡要介紹了相關微磁學模擬研究的進展.

    2 斯格明子的拓撲物理基礎

    2.1 手性磁體中的斯格明子

    一般的傳統(tǒng)磁性材料都具有中心反演對稱性,因而呈現(xiàn)出鐵磁性、反鐵磁性等共線的磁結構.但在一類如MnSi[3]、Fe1?xCoxSi[3,22]、FeGe[23,24]、M nGe[25?27]、Cu2OSeO3[28,29]等中心反演對稱破缺的立方結構B20型化合物中,由于自旋、軌道、晶格等多種自由度的關聯(lián)與耦合同時存在,導致了非線性的DM(Dzyaloshinskii-Moriya)相互作用[30?32],

    Si,Sj是相鄰的自旋;Dij是DM作用矢量.在反演操作下Si和Sj交換后,能量的正負號發(fā)生變化,體現(xiàn)了空間反演對稱性的打破.(1)式對應的能量最小值的構型是Si與Sj在垂直于Dij的面內,根據(jù)Dij的符號構成手性相互垂直,這種構型就是典型的非共線磁構型中的最近鄰構型.包含DM相互作用的體系在外磁場以及溫度的作用下出現(xiàn)多種磁有序結構,例如螺旋磁(helical)、鐵磁、圓錐螺磁(conical)、斯格明子相[4,10,28,33?35]等. 特別地,斯格明子是一種存在于一定溫度、磁場窗口下的拓撲有序磁渦旋狀結構,圖1(a)和圖1(c)分別描述了布洛赫型和奈爾型斯格明子的二維結構.包含最近鄰海森伯相互作用、DM相互作用以及外磁場作用的哈密頓量可以寫為

    其中J是相鄰自旋海森伯交換相互作用系數(shù),g是Lande系數(shù),μB是玻爾磁矩.

    在B20化合物中由于空間反演對稱性破缺而產生的DM相互作用矢量是從一個自旋指向另外一個自旋的,這樣的DM相互作用的結果是相鄰自旋傾向于相互垂直排列.而海森伯交換相互作用(鐵磁)下相鄰自旋是平行的.在這兩者的共同影響下相鄰的自旋呈現(xiàn)出一定夾角的排列方式(圖1(b)所示),呈螺旋式(helical)旋轉.從一維來看由于相鄰自旋的旋轉角度是固定的,在一維上自旋構成了波狀結構自旋鏈(圖1(a)).二維中一定磁場存在時,自旋的取向傾向與外磁場平行,而相鄰自旋又傾向一定夾角的排列方式,此時在DM相互作用、鐵磁相互作用、外磁場相互作用的共同影響下,自旋形成了圖1(a)中的斯格明子結構,這種結構被稱為布洛赫型斯格明子,在一個斯格明子內部沿著直徑各個方向都呈現(xiàn)出波狀的結構.根據(jù)(2)式,不考慮外磁場時圖1(b)中的相鄰自旋的夾角的大小應當為θ≈D0/2J,圖1(a)中的波長也是斯格明子的直徑,為λ=4πaJ/D0,D0是DM作用矢量的模.

    圖1 典型的斯格明子構型[7] (a)布洛赫型斯格明子示意圖;(b)海森伯相互作用與體手性DM相互作用共同作用下的最近鄰自旋構型;(c)奈爾類斯格明子示意圖;(d)自旋軌道耦合導致的界面DM相互作用Fig.1.Skyrmion configuration and DM interaction[7]:(a)B loch type skyrmion;(b)neighbor spin configuration under Heisenberg interaction and bulk DM interaction;(c)Néel type skrymion;(d)interfacial DM interaction induced by spin orbital coup ling.

    多數(shù)體手性磁體中的斯格明子存在的溫度窗口非常有限,限制了其自旋電子學應用.隨著相關研究的深入,另一類手性磁體中的斯格明子被發(fā)現(xiàn),在二維重金屬-輕鐵磁金屬-絕緣體的異質結中,重金屬層中的自旋軌道耦合作用使得界面反演對稱被打破,斯格明子可以在這種界面手性磁體中存在,并且其溫度窗口可能拓寬至室溫[5?8].這種界面手性磁體中DM作用矢量Dij垂直于連接相鄰自旋的向量rij(圖1(e)),該作用矢量的符號決定了自旋結構是右手手性還是左手手性.界面DM相互作用使得相鄰自旋相互垂直,類似的在海森伯鐵磁作用的共同作用下,相鄰自旋形成如圖1(d)中紅色箭頭顯示的鐘擺式旋轉,最終形成的一維自旋波的結構如圖1(c)所示,圖1(c)中還顯示了這些作用下的斯格明子的結構,也就是是奈爾型斯格明子.(2)式決定了奈爾型斯格明子的直徑、相鄰自旋的夾角等參數(shù)都與布洛赫型的相同,但是兩者斯格明子中的相鄰自旋的旋轉方向是不同的.

    2.2 斯格明子的拓撲荷

    奈爾型的斯格明子與布洛赫型的斯格明子構型差異很大,然而奈爾構型的斯格明子中每個自旋旋轉90?就得到了布洛赫構型,兩者可以通過連續(xù)變形相互轉化,因此這兩種構型是拓撲相同的.斯格明子的磁結構投影到一個球面上,構成了一個由手性自旋構成的閉合球面,球面上的自旋都指向球面外,如圖2(a)所示,球面內的自旋構型類似于一個“磁單極”的效果[36],這樣的一個球面體現(xiàn)了斯格明子的拓撲性質,可以用拓撲荷(topological charge)

    來表征其拓撲性質[37,38],其中m是描述自旋取向的單位矢量,其意義是從某一點出發(fā)遍歷整個球面,這些點上的自旋連續(xù)變化的周期數(shù),這里的周期對應的是球面立體角4π.可以看出,對于每一個斯格明子,其拓撲荷為±1,拓撲荷不隨著其構型的擾動而變化.從斯格明子的構型來看,其內部的三角形網格上的三個自旋構成了一個立體張角,如圖2(b)所示.對應到到球面投影上可以發(fā)現(xiàn),所有的這些立體角的疊加就是這個球面的立體角4π,因此可以使用計算立體角的方式來計算斯格明子的拓撲荷[39]:

    其中m1,m2,m3為圖中相鄰的三個自旋方向的單位矢量;ρ是歸一化系數(shù),

    圖2 斯格明子的拓撲荷[40,41] (a)斯格明子的投影變換;(b)拓撲荷的計算Fig.2.The topological charge of skyrmion[40,41]:(a)Stereographic projection of a skyrmion;(b)the caculation of topological charge.

    由于斯格明子拓撲荷的符號是由DM作用矢量的符號決定的,而DM作用矢量由材料、結構決定,因此在同一個材料中一般不會同時存在拓撲荷分別為1和?1的斯格明子.考慮到對于平庸的自旋結構,如鐵磁、vortex等其拓撲荷為0,數(shù)值模擬中可以用拓撲荷來表征體系的斯格明子數(shù)目∫,并且在單個斯格明子的情況下,可以通過積分來計算該斯格明子的中心位置[40?42].另外還要注意的是,在一定的溫度范圍內,由于DM相互作用和溫度擾動的共同作用,溫度擾動將引入局域非結構的拓撲荷[40],這些局域拓撲荷的總和不為零,但不代表任何結構化的構型.

    單純地從實空間結構來理解斯格明子的拓撲性,哈密頓量給出了最近鄰相互作用,確定了相鄰的兩個自旋構型,是典型的短程有序,在一維自旋鏈中形成了自旋波的形式(圖1(a),圖1(c));然而到了二維時,當自旋的構型以斯格明子為基態(tài)時,整個體系通過斯格明子這種尺度上比最近鄰自旋大得多的有序結構來實現(xiàn)能量最低,可以說是一種中程有序.在不同的態(tài)與這種中程有序的態(tài)相互轉變中,總要克服一定的勢壘,因此這種額外的中程有序為體系提供了額外的穩(wěn)定性,也是斯格明子穩(wěn)定性的一種表現(xiàn).

    從拓撲的角度來分析,在絕熱近似和自旋連續(xù)轉動的條件下,非平庸拓撲自旋結構和平庸拓撲自旋結構無法相互轉換[43],意味著無法在絕熱近似下通過連續(xù)的自旋轉動變化來完成磁性斯格明子的產生和湮滅過程,因此斯格明子的拓撲自旋結構是受到拓撲保護的,增加了其穩(wěn)定性.在MnSi等的B20結構磁體的實驗中就發(fā)現(xiàn)可以長期存在的亞穩(wěn)態(tài)斯格明子[44,45],證明了拓撲穩(wěn)定性的存在.然而斯格明子的拓撲穩(wěn)定性也帶來了一些問題,比如控制斯格明子產生湮滅的人工調控手段成為一個問題,在數(shù)值模擬中,亞穩(wěn)態(tài)的存在也同樣使得計算過程變得更為復雜.

    2.3 拓撲霍爾效應與斯格明子霍爾效應

    根據(jù)洪德定則,載流子的自旋總是試圖與所在位置的局域電子的自旋保持平行.當一個自由電子穿過斯格明子時,雖然根據(jù)斯格明子的構型,自由電子離開時和入射時的自旋方向是平行的,然而在通過的過程中其自旋的方向是不斷翻轉的.將這樣的翻轉過程等效為一個磁場,就是自由電子通過斯格明子時感受到的虛擬emergent磁場bz[9],

    對比(3)式中計算的拓撲荷可以發(fā)現(xiàn),emergent磁場強度是正比于拓撲荷的bz∝Q,是與斯格明子的拓撲性相關的.當斯格明子以一定速度運動時,運動的emergent磁場還產生一個電場,圖3顯示了載流子在這些場的共同作用下將產生一個偏轉,這就是拓撲霍爾效應(topological Hall effect)[10?13,15].

    當自由電子流經斯格明子時受到拓撲霍爾效應產生偏轉,而斯格明子也在自由電子的反作用下,受到一個反向的作用力,這個作用力是馬格努斯力[38],磁結構向著另一側偏轉,這個效應是斯格明子霍爾效應(skyrmion Hall effect).霍爾效應是帶電粒子在磁場中運動引發(fā)的偏轉現(xiàn)象,電荷的正負性決定了該粒子向著相反的方向偏離.相比之下,拓撲荷為±1,可以被視為準粒子的斯格明子,在電流的自旋轉移矩驅動下運動中產生的斯格明子霍爾效應也有類似的現(xiàn)象,拓撲荷的符號決定了斯格明子拓撲霍爾效應的偏轉方向[14,15].

    圖3 拓撲霍爾效應斯格明子霍爾效應[38]Fig.3.Topological Hall effect and skyrmion Hall effect[38].

    M agnon描述了自旋在有限溫的激發(fā)下偏離其平衡位置的程度,可以視為攜帶了與局域自旋反平行的自旋.當自由電子流經斯格明子時產生了拓撲霍爾效應,同樣攜帶自旋方向信息的magnon流也就是熱流流經絕緣體斯格明子時也會產生類似的現(xiàn)象,magnon與emergent磁場相互作用使得magnon流偏移,就是拓撲magnon霍爾效應.當偏移的magnon在材料的邊界處積累時,就會產生一個可測量的溫度差,實驗[46]和數(shù)值模擬[47,48]都證實了這一現(xiàn)象.

    2.4 手性磁體中的拓撲磁結構與相圖

    前面討論了手性磁體中相鄰自旋構型的取向,在DM作用和海森伯交換作用的共同競爭下,斯格明子存在于一定的磁場窗口內,注意到斯格明子之間區(qū)域的自旋構型是鐵磁態(tài),也說明了相鄰自旋鐵磁排列與螺旋排列之間的能量差不是很大,在外磁場的調節(jié)下以及有限溫的情況下將呈現(xiàn)出豐富的磁結構構型.

    在二維體手性材料里主要有鐵磁態(tài)、斯格明子態(tài)、螺旋態(tài).圖4(a)顯示了在二維體手性磁體中隨著磁場溫度變化體系磁結構的相圖.低磁場下是螺旋態(tài),結構如圖4(b)所示,隨磁場的上升螺旋態(tài)逐漸變?yōu)樗垢衩髯討B(tài).在適當?shù)拇艌鰷囟认滤垢衩髯訒帕谐擅芏逊e的結構,如圖4(c),類似于晶體中的原子排列,所以稱為斯格明子晶體態(tài)[4].在準一維納米條帶中,類似的堆積結構也被稱為packed斯格明子態(tài)[49].隨著納米條帶寬度的降低,同樣的堆積結構逐漸變成一維斯格明子鏈.相圖中介于鐵磁態(tài)和斯格明子晶體態(tài)之間的斯格明子與鐵磁的混合態(tài),也可以稱之為斯格明子氣體態(tài);當體系中只有一個斯格明子時可以稱為單獨斯格明子態(tài).在界面手性材料中,螺旋態(tài)變成了一種類似的螺線(spiral)態(tài),如圖4(d)所示,螺線態(tài)與螺旋態(tài)結構的主要差異是相鄰自旋旋轉方向的不同.

    圖4 二維手性磁體中的磁結構 (a)二維體手性材料中的相圖[42];(b)螺旋態(tài)結構:(c)斯格明子晶體態(tài);(d)螺線態(tài)Fig.4.Magnetic structure in 2D chiral magnets:(a)Phase diagram in 2D bulk chiral magnetism[42];(b)structure of helical state;(c)structure of skyrmioncrystall;(d)structure of spiral state in interfacial chiral magnetism.

    數(shù)值模擬中可以用蒙特卡羅方法模擬退火過程[50]來計算不同磁場、溫度下的基態(tài),一般經典的Metropolis[51]算法就可滿足要求,overrelaxation算法、非線性共軛梯度法等可以用來加速計算[40,52,53].計算基態(tài)時可以在高溫(幾倍J的溫度)下使用隨機態(tài)作為初始態(tài),逐步降低溫度至零溫,一般每個溫度計算5×104個循環(huán)就能很好地收斂[54].此外在計算基態(tài)時要特別注意體系的格點數(shù)目,前文中已經討論了在零磁場時,螺旋態(tài)、螺線態(tài)的周期都是確定的,網格數(shù)量應當是周期的整數(shù)倍,因此在蒙特卡羅計算中網格數(shù)量應當與J,D等參數(shù)相匹配.

    2.5 斯格明子的動力學

    斯格明子的自旋拓撲性質使之成為一個獨特的研究實空間拓撲輸運現(xiàn)象的平臺[3,4,9,11,38,43,55].特別地,在金屬斯格明子材料中外加電流導致自旋轉移矩效應[56,57],由此驅動斯格明子的運動,并且電流對斯格明子的驅動非常有效率,臨界電流密度在102A/cm2左右,比一般的鐵磁疇壁的驅動電流小了4個數(shù)量級[20],同時對材料中的缺陷有強的抗干擾能力[17].正是由于斯格明子較小的尺寸、較低的電流操作性能,使之有潛力成為斯格明子賽道存儲器件.

    在界面斯格明子材料中,當電流流經重金屬層時,由于強的自旋軌道耦合產生自旋霍爾效應,如圖5所示.向上和向下的自旋分離,從而向上層注入自旋極化電流[6,58],被稱為自旋軌道轉矩效應[9,16].自旋軌道轉矩效應的作用下,金屬層中的自旋角動量注入上層輕鐵磁層中,最終推動了斯格明子的運動[16,59?62],也是一種有效地控制斯格明子動力學的手段.

    圖5 界面手性磁材料中的自旋軌道轉矩[63,64]Fig.5.Spin Orbit Torque in interfacial chiral magnets[63,64].

    自旋轉移矩效應和自旋軌道轉矩效應都有利于斯格明子被電流驅動,根據(jù)Landau-Lifshitz-Gilbert(LLG)方程可以給出斯格明子在電流輸運下的動力學方程的形式[64]:

    其中γ是磁回轉率(gyromagnetic ratio),α是磁阻尼系數(shù),Ms是飽和磁化率.方括號中的第一項是電流j1對應的自旋轉移矩項;第二項是自旋軌道轉矩項,θSH是霍爾效應系數(shù),tf是薄膜厚度,je是流經重金屬層的電流,這一項是針對多層膜結構中的奈爾型斯格明子計算其在電流驅動下的動力學過程.(7)式中,等效場Heff=??F/?(m)包含了自旋間的相互作用以及外磁場的作用,F是磁自由能.(2)式中給出了格點模型中的哈密頓量,在微磁學模擬中一般采用連續(xù)的模型,相應的自由能的形式可以寫為

    方括號中的三項分別為海森伯交換相互作用、DM相互作用和外磁場的作用所對應的能量密度.其中DM相互作用密度,在體相互作用的情況下可以寫為EDMbulk=D0m·(?×m);而界面相互作用時,其能量密度為EDMinter=D0(mz(?·m)?(m·?)mz),這里的D0為相應的DM作用矢量的模.求導后可以給出有效場的形式:

    右端的三項分別為海森伯交換相互作用、DM相互作用以及外磁場作用對有效場的貢獻,分別為

    有效場的形式中采用了歸一化的磁矩m以及相應的參數(shù),實際計算中需要注意參數(shù)的量綱以及正負號的選取等問題.

    值得注意的是自旋轉移矩項中的j1并非只針對電流,也可以是熱流或者magnon流[42,65].前文中已經討論了magnon可以被視為與局域自旋反平行的自旋,熱的傳播可以被視為magnon流的輸運.因此當magnon流經斯格明子時,類似于自由電子為載流子時的自旋轉移矩效應,也將產生移矩效應,可以在絕緣體材料中驅動斯格明子的移動[66].理論上magnon很好地解釋了熱流下的斯格明子動力學,而數(shù)值模擬上magnon并不是一個容易描述的物理量,可以采用隨機場L來描述不同溫度下自旋的熱擾動,含隨機力的LLG方程[67,68]為

    對于隨機場L有

    這里的a是計算網格的邊長.這里要注意含隨機力的LLG方程(14)需要采用Heun算法來計算[67].

    微磁學模擬中對動力學方程(7)的計算,在時間離散上一般采用四階龍格-庫塔算法,空間離散多數(shù)采用有限差分處理,也可以使用有限元方法,相比之下有限元方法在空間離散上更為精確,但是計算開銷非常大.常用的計算軟件有MuMax3[69],OOMMF[70]等,MuMax提供了基于CUDA的GPU并行計算能力,比OOMMF計算更快但是更難增加自定義模塊.

    3 斯格明子自旋電子器件應用相關動力學研究進展

    3.1 作為賽道存儲應用的斯格明子

    賽道存儲器是一種新型非易失性存儲器[20],傳統(tǒng)的磁存儲如硬盤等是考慮移動存儲介質從而在固定位置讀取存儲介質中的磁信息,而賽道存儲則考慮在一個準一維條帶上用磁疇來存儲信息,使用電流驅動磁疇運動到相應的寫入讀取位置來實現(xiàn).經過幾代的發(fā)展,目前已經使用了自旋軌道轉矩來輔助電流驅動磁疇運動、相互之間反鐵磁的兩條賽道來加速磁疇運動[21].賽道存儲器的一個賽道單元上可以存儲很多個數(shù)據(jù),并且可以實現(xiàn)三維記錄,所以賽道存儲器有望實現(xiàn)超高密度的存儲.

    影響賽道存儲器性能最重要的三個問題包括:一是提高疇壁的移動速度,可以加快數(shù)據(jù)的讀取速度;二是過高的臨界電流密度會產生大的焦耳熱以及功耗,需要進一步降低電流密度;三是提高賽道存儲器的存儲密度,也就是單根納米帶單元中能有更多的磁疇壁、兩根納米帶之間能有更小的間距.

    具有穩(wěn)定拓撲結構的斯格明子的尺寸非常小,界面手性磁體中斯格明子的大小是5 nm[71],滿足了賽道存儲器件的密度要求;驅動斯格明子運動的臨界電流只有102A/cm2[17],比驅動磁疇所需的電流密度小了4個數(shù)量級[9],解決了焦耳熱問題.此外斯格明子在納米條帶中的輸運過程有跨越障礙[7]等特點.這些特性使得斯格明子是制作超高密度賽道存儲器非常好的選擇[18,19].

    作為賽道存儲器件應用就要求斯格明子在室溫具有穩(wěn)定的電流操作性、可控的產生湮滅過程以及方便的可探測性.目前室溫下存在的斯格明子已經被觀測到[14],產生湮滅過程的人工控制也有很多相應的手段[41,72?77].然而綜合來看,室溫下產生湮滅控制仍然是一個挑戰(zhàn)[37],此外,方便可集成的探測手段相關研究仍在起步階段[78].

    3.2 非均勻電流下斯格明子的生成

    賽道存儲器性能的優(yōu)化要求盡可能快地驅動斯格明子以加快數(shù)據(jù)的讀取速度,意味著要使用較強的電流密度.斯格明子雖然也能通過幾何限制區(qū)域,然而界面手性材料中,幾何受限導致的電流密度不均勻時,微磁學模擬發(fā)現(xiàn)電流驅動奈爾磁疇條帶通過幾何限制區(qū)域時會動態(tài)地生成斯格明子[79,80],實驗上也成功地觀察到這一現(xiàn)象[63],在體手性材料中微磁學模擬也預測了類似現(xiàn)象[81].這里主要討論界面手性材料中電流驅動下動態(tài)生成斯格明子的原因[64].

    前文已經討論了在奈爾型斯格明子的界面手性材料的重金屬-輕鐵磁金屬-絕緣體的異質結的多層膜中,電流通過重金屬層時,由于強的自旋軌道耦合而產生自旋霍爾效應,使得重金屬層中的自旋角動量注入上層鐵磁金屬中,從而產生自旋軌道轉矩.自旋軌道轉矩對上層金屬中磁疇壁的作用如圖6(a)所示,其對應的有效霍爾場Bsh

    其中是自旋矢量;薄膜的模量單位矢量;是驅動電流方向;是鐵磁層的厚度,Ms是飽和磁化強度,自旋霍爾系數(shù)θsh=Js/Jc是重金屬層中自旋流密度Js與電流密度Jc的比值.圖6(a)中顯示了手性磁疇中對應Bsh的強度與自旋矢量和驅動電流的方向有關.因此在包含條帶狀磁疇的納米線中,圖6(b)區(qū)域②,③中,自旋矢量的取向平行于所以自旋軌道轉矩對這兩個區(qū)域的磁疇壁沒有驅動作用.而區(qū)域①中的Bsh最大,所以這里自旋轉矩效應對磁結構的驅動效果最大.

    如圖6(c)所示,在引入幾何受限結構后,電流分布不再均勻,在左右兩邊分別是發(fā)散、匯聚的構型.在左側入口處電流的y方向分量如圖6(d)所示,其結果是產生了圖6(e)中的等效力,在這個等效力的作用下,整個磁疇條帶向右移動的同時,在右側出口處逐漸收緊,最終生成了圖6(e)中的斯格明子結構.

    圖6 電流驅動下幾何限制結構動態(tài)生成斯格明子[63,64] (a)自旋軌道轉矩對左半邊中的有效霍爾場;(b)條帶中奈爾磁疇對應的不同區(qū)域;(c)通過幾何限制位置時的電流分布;(d),(e)電流y分量的分布以及對應的等效作用力;(f)斯格明子的形成Fig.6.Skyrmion creation in geometrical constriction[63,64]:(a)The effective spin Hall field act on left side domain wall by spin orbital torque;(b)the three region of Néel domain wall in the stripe;(c)the distribution of electrical current in the geometrical constriction region;(d),(e)the in homogeneoous current and effective field around the neck region;(f)the genaration of skyrmion.

    3.3 電流驅動下斯格明子的穩(wěn)定輸運

    前面討論了拓撲霍爾效應和斯格明子霍爾效應作為斯格明子在電流驅動下的拓撲性的體現(xiàn),帶來了豐富的斯格明子動力學性質.然而,作為信息載體的斯格明子在賽道中電流驅動之下會出現(xiàn)一個垂直于電流方向的速度,最終結果是偏離到賽道邊緣被釘扎或者在大的電流驅動下在邊緣處湮滅[82?86].這兩種情況都是在賽道存儲應用中不希望遇到的,相關的理論以及數(shù)值模擬工作研究了怎樣減小或者避免斯格明子拓撲效應.例如,通過增強垂直于賽道邊緣的磁各向異性來更好地限制斯格明子在賽道內運動[84];類斯格明子構型skyrmionium因為拓撲荷為0,電流驅動時沒有斯格明子霍爾效應[84,87,88];反鐵磁斯格明子[84,86]、反鐵磁作用的斯格明子雙層膜結構[82,86]在電流驅動下因為兩個正好相反的斯格明子霍爾效應相互抵消而呈現(xiàn)出垂直于電流方向的速度為零的運動特性.

    圖7 斯格明子霍爾效應的消除 (a)反鐵磁斯格明子[83];(b)反鐵磁耦合的雙層斯格明子[82]Fig.7.Elimination of skyrmion Hall effect:(a)Antiferromagnetic skyrmion[83];(b)antiferromagnetic coupled Skyrmion bilayer[82].

    反鐵磁斯格明子中可以分解出拓撲荷分別為1,?1的兩個子斯格明子,如圖7(a)所示,這兩個子斯格明子受到的馬格努斯力正好相反,最終反鐵磁斯格明子在電流驅動下不顯現(xiàn)出霍爾效應.反鐵磁的斯格明子在實驗上難以實現(xiàn),然而反鐵磁耦合的斯格明子雙層結構在Co/Ru/Co異質結上實現(xiàn)了[89].類似于反鐵磁斯格明子,這種結構也由于馬格努斯合力為零而不顯現(xiàn)斯格明子霍爾效應,不同點在于將拓撲荷為1,?1的兩個斯格明子拆成兩層以反鐵磁作用耦合在一起,如圖7(b)所示.反鐵磁帶來了凈磁矩為零從而使得探測這樣的斯格明子結構成為難題,理論上預測了隨著自旋結構重排,反鐵磁斯格明子也可以顯示出拓撲霍爾效應[90],從而可能被探測到.

    3.4 斯格明子產生湮滅過程中的人工控制

    賽道存儲應用要求作為信息載體的斯格明子的產生湮滅過程完全可控.受到拓撲保護的斯格明子在較低的溫度下的產生湮滅過程有一個很大的勢壘[91],而較高的溫度下斯格明子本身就有著很強的隨機性[78].因此,斯格明子的產生湮滅過程的人工控制手段就成為了一個熱門話題.

    從圖4(a)的相圖來看,產生湮滅斯格明子的過程可以分為兩類:一類是從螺旋態(tài)或者螺線態(tài)到斯格明子之間的相變,另一類是從鐵磁態(tài)到斯格明子態(tài)之間的相變.螺旋態(tài)或者螺線態(tài)的構型比較接近斯格明子,所以它們之間的相變相對容易,靜態(tài)下合適的磁場溫度中,helical態(tài)的每個條帶在有邊界的參與時成為兩端封閉、拓撲荷為1的條帶狀磁疇構型,在磁場提高時會逐漸收縮或者從中間打斷生成斯格明子[92],在沒有邊界參與時磁疇將逐漸演化直至打斷成斯格明子;動態(tài)下前文已經討論的非均勻電流生成斯格明子就是典型的生成方法.

    然而賽道存儲器材中用有無斯格明子來代表儲存的數(shù)據(jù)是“1”和“0”,因此從鐵磁態(tài)到斯格明子之間的產生湮滅過程是最為重要的.不幸的是,受到拓撲保護的斯格明子與鐵磁態(tài)之間的相互轉變過程是非連續(xù)的[43],除了(2)式中適當?shù)拇艌龃翱谕膺€需要有外部的擾動來輔助拓撲成核,這些擾動有垂直方向的電流/極化電流脈沖[41,72,73]、激光脈沖[74,75]、磁場脈沖[76]、局部加熱[77]等.

    對于布洛赫型斯格明子,垂直電流輔助產生湮滅的作用分為三類[41]:電流的感生磁場、自旋極化電流的自旋轉矩效應、電流帶來的焦耳熱.其中熱效應對應著隨機擾動,在不同溫度擾動下有相應的概率完成拓撲成核過程從而輔助斯格明子的轉變.而感生磁場是漩渦狀的,對比布洛赫型的斯格明子可以發(fā)現(xiàn)漩渦狀的磁場會使得自旋的構型從鐵磁態(tài)向著接近斯格明子構型的方向轉變,因此有助于斯格明子的形成.而極化電流輔助斯格明子成核時,極化方向如圖8(a)所示,和垂直方向有一定夾角,極化電流的作用是使得一部分區(qū)域的自旋沿著面內排布從而使整個自旋構型更接近斯格明子.圖8(b)中顯示了奈爾型斯格明子在STM針尖極化電流輔助下的產生湮滅過程示意圖[71],注意到圖中顯示的極化電流的極化方向是垂直于平面的,與布洛赫型的有著顯著差異.

    斯格明子的產生湮滅過程在拓撲上對應了拓撲成核過程,既然體系無法自發(fā)的完成轉變過程,外界擾動提供了拓撲成核點來輔助轉變過程.布洛赫型的斯格明子的結構特征顯示出中心點周圍自旋的面內分量呈現(xiàn)漩渦狀排列,而拓撲成核點也應當有類似的構型,因此其對應的外加擾動更多的呈現(xiàn)出漩渦狀的特點;同時奈爾型的斯格明子的構型特征是中心點附近區(qū)域的自旋在z方向的分量反平行于周圍自旋,因此對應的外界擾動更多地呈現(xiàn)出垂直于其所在平面的特點.

    圖8 自旋極化電流產生斯格明子示意圖 (a)布洛赫型[41];(b)奈爾型[71]Fig.8.Generation of skyrmion by spin polarized current:(a)B loch type[41];(b)Néel type[71].

    4 結論與展望

    斯格明子相關研究是一個迅速發(fā)展的領域,本文從賽道存儲方面簡要介紹了一些電流驅動斯格明子、產生湮滅人工控制手段及其相關的拓撲物理性質.在這一領域仍然有很多應用與基礎問題需要繼續(xù)探索,例如可集成的便捷探測手段、電流驅動下霍爾效應的消除、低磁場室溫下斯格明子的控制、拓撲成核過程的拓撲保護性等,這些問題的解決可以使斯格明子的賽道存儲應用推向實際.

    猜你喜歡
    明子鐵磁霍爾
    關于兩類多分量海森堡鐵磁鏈模型的研究
    這個大童話講貓(五)名偵探西尼·霍爾
    機器人兒子
    你好,鐵磁
    哎喲,不怕!
    好孩子畫報(2016年5期)2016-05-14 09:54:26
    你好,鐵磁
    學生天地(2016年27期)2016-04-16 05:15:41
    機器人兒子
    機器兒子
    離子推力器和霍爾推力器的異同
    太空探索(2015年6期)2015-07-12 12:48:42
    一維交替鐵磁-反鐵磁耦合的海森堡鏈[Mn(N3)2(pybox)]n
    亚洲美女搞黄在线观看| 免费高清在线观看日韩| 又大又黄又爽视频免费| 天美传媒精品一区二区| 18+在线观看网站| 久久 成人 亚洲| 国产有黄有色有爽视频| 丰满饥渴人妻一区二区三| 高清av免费在线| 不卡视频在线观看欧美| 亚洲综合精品二区| 日韩精品有码人妻一区| 日韩一区二区三区影片| 少妇人妻 视频| 伊人久久国产一区二区| 久久久久人妻精品一区果冻| 日本黄色日本黄色录像| 一二三四中文在线观看免费高清| 在线观看一区二区三区激情| av在线观看视频网站免费| 亚洲精品成人av观看孕妇| 亚洲av男天堂| 夫妻性生交免费视频一级片| 男人舔女人的私密视频| 综合色丁香网| 寂寞人妻少妇视频99o| 999久久久国产精品视频| 中国国产av一级| 午夜av观看不卡| 国产精品不卡视频一区二区| 久久久久网色| 国产一区亚洲一区在线观看| 久久女婷五月综合色啪小说| 成年av动漫网址| 热re99久久国产66热| 超色免费av| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 校园人妻丝袜中文字幕| 最黄视频免费看| 不卡视频在线观看欧美| 香蕉精品网在线| 久久精品国产综合久久久| 国产精品人妻久久久影院| 欧美日韩av久久| 欧美xxⅹ黑人| 国产成人a∨麻豆精品| 亚洲国产欧美在线一区| 多毛熟女@视频| 蜜桃国产av成人99| 欧美日韩国产mv在线观看视频| 99久久人妻综合| 欧美日韩国产mv在线观看视频| 在线观看三级黄色| 在线观看免费日韩欧美大片| 久久久久国产网址| 在线观看一区二区三区激情| 建设人人有责人人尽责人人享有的| 久久久久久久精品精品| 一级毛片 在线播放| 母亲3免费完整高清在线观看 | www.av在线官网国产| 巨乳人妻的诱惑在线观看| 日韩中文字幕欧美一区二区 | 韩国精品一区二区三区| 午夜老司机福利剧场| 精品午夜福利在线看| 免费高清在线观看视频在线观看| 在线精品无人区一区二区三| 精品亚洲乱码少妇综合久久| 国产精品99久久99久久久不卡 | 中文欧美无线码| 人妻一区二区av| 久久久久久人妻| 2021少妇久久久久久久久久久| 美女xxoo啪啪120秒动态图| 国产精品亚洲av一区麻豆 | 夫妻性生交免费视频一级片| 蜜桃在线观看..| 高清在线视频一区二区三区| 婷婷成人精品国产| 看免费av毛片| 午夜激情久久久久久久| 亚洲天堂av无毛| 亚洲经典国产精华液单| 欧美精品高潮呻吟av久久| 99久国产av精品国产电影| 国产综合精华液| 男女国产视频网站| 精品国产一区二区久久| 国产熟女欧美一区二区| 久久人人爽人人片av| 青青草视频在线视频观看| 女人精品久久久久毛片| 成人国产av品久久久| 欧美日韩av久久| 国产亚洲欧美精品永久| 精品福利永久在线观看| 人妻 亚洲 视频| 18禁动态无遮挡网站| 日韩 亚洲 欧美在线| 在线亚洲精品国产二区图片欧美| 国产精品熟女久久久久浪| 成人漫画全彩无遮挡| 18禁裸乳无遮挡动漫免费视频| 街头女战士在线观看网站| 80岁老熟妇乱子伦牲交| 曰老女人黄片| 在线观看国产h片| 精品久久久久久电影网| 国产在线视频一区二区| 国产日韩欧美亚洲二区| 久久精品夜色国产| 亚洲精品日本国产第一区| 精品少妇久久久久久888优播| 成人毛片60女人毛片免费| 91精品三级在线观看| 丁香六月天网| av国产久精品久网站免费入址| 欧美日韩精品成人综合77777| 永久免费av网站大全| 国产av码专区亚洲av| 国产成人一区二区在线| 国产日韩欧美在线精品| 飞空精品影院首页| av女优亚洲男人天堂| 2021少妇久久久久久久久久久| 亚洲精品av麻豆狂野| 男女国产视频网站| 国产精品免费大片| 国产精品亚洲av一区麻豆 | 午夜免费男女啪啪视频观看| 午夜影院在线不卡| 久久精品国产鲁丝片午夜精品| 亚洲精品第二区| 亚洲四区av| 曰老女人黄片| 国产精品久久久久成人av| 99九九在线精品视频| 久久韩国三级中文字幕| 久久国产亚洲av麻豆专区| 成人漫画全彩无遮挡| 男女高潮啪啪啪动态图| 男女高潮啪啪啪动态图| 极品人妻少妇av视频| 亚洲国产av新网站| 日韩精品有码人妻一区| 免费观看a级毛片全部| 欧美亚洲 丝袜 人妻 在线| 男女啪啪激烈高潮av片| 制服丝袜香蕉在线| 满18在线观看网站| 免费观看av网站的网址| 狠狠婷婷综合久久久久久88av| 国产成人精品在线电影| 欧美黄色片欧美黄色片| 国产av一区二区精品久久| 国产精品麻豆人妻色哟哟久久| videos熟女内射| 十八禁高潮呻吟视频| www.自偷自拍.com| 久久久久久人妻| 久久久久久久国产电影| 自线自在国产av| 国产精品蜜桃在线观看| 久久久久久久国产电影| 亚洲美女搞黄在线观看| av网站在线播放免费| 亚洲人成电影观看| a级毛片在线看网站| 久久免费观看电影| 人成视频在线观看免费观看| 美女高潮到喷水免费观看| 国产一区二区三区综合在线观看| 国产欧美亚洲国产| 成年人免费黄色播放视频| 在线天堂中文资源库| 老汉色av国产亚洲站长工具| www.av在线官网国产| 欧美黄色片欧美黄色片| a级片在线免费高清观看视频| 欧美激情 高清一区二区三区| 91精品三级在线观看| 在线观看免费视频网站a站| 午夜免费观看性视频| 欧美97在线视频| 在线观看免费视频网站a站| 久久精品熟女亚洲av麻豆精品| 少妇的逼水好多| 熟女av电影| 捣出白浆h1v1| 久久精品国产亚洲av涩爱| 亚洲久久久国产精品| 90打野战视频偷拍视频| 亚洲男人天堂网一区| 免费看不卡的av| 亚洲国产精品999| 久久鲁丝午夜福利片| 国产成人精品无人区| 国产麻豆69| 亚洲av男天堂| 亚洲伊人久久精品综合| 18禁裸乳无遮挡动漫免费视频| 99久久中文字幕三级久久日本| 少妇的丰满在线观看| 又粗又硬又长又爽又黄的视频| 精品一区二区三区四区五区乱码 | 中文字幕人妻熟女乱码| 又大又黄又爽视频免费| 亚洲综合色网址| 两个人看的免费小视频| av国产久精品久网站免费入址| 欧美日韩av久久| 三上悠亚av全集在线观看| 国产一区有黄有色的免费视频| 男人添女人高潮全过程视频| 蜜桃在线观看..| 成人免费观看视频高清| 天堂中文最新版在线下载| 久久精品亚洲av国产电影网| 日本wwww免费看| 欧美少妇被猛烈插入视频| 国产成人免费无遮挡视频| 色婷婷久久久亚洲欧美| 美女视频免费永久观看网站| 久久精品国产a三级三级三级| 美女中出高潮动态图| 91精品三级在线观看| 亚洲欧美精品自产自拍| 只有这里有精品99| 美女国产视频在线观看| 欧美另类一区| 咕卡用的链子| 一区二区三区激情视频| 伊人久久国产一区二区| 亚洲精华国产精华液的使用体验| 免费大片黄手机在线观看| av.在线天堂| 青草久久国产| 日韩一卡2卡3卡4卡2021年| 成人国产av品久久久| 亚洲美女视频黄频| 黄色毛片三级朝国网站| 七月丁香在线播放| 国产精品人妻久久久影院| 成人影院久久| 黄色怎么调成土黄色| 国产免费一区二区三区四区乱码| 春色校园在线视频观看| 国产片特级美女逼逼视频| 婷婷色综合www| 免费黄色在线免费观看| 少妇熟女欧美另类| 在线亚洲精品国产二区图片欧美| av在线app专区| 考比视频在线观看| 国产不卡av网站在线观看| 亚洲男人天堂网一区| 在线天堂最新版资源| 18在线观看网站| 精品一区二区三区四区五区乱码 | 国产成人免费观看mmmm| 亚洲久久久国产精品| 免费高清在线观看日韩| 亚洲欧美精品自产自拍| 狠狠婷婷综合久久久久久88av| 国产一区二区三区av在线| 精品一区二区免费观看| 高清av免费在线| 少妇猛男粗大的猛烈进出视频| 免费日韩欧美在线观看| 免费黄色在线免费观看| 男女午夜视频在线观看| 少妇精品久久久久久久| 日韩在线高清观看一区二区三区| 精品国产一区二区三区久久久樱花| 国产免费福利视频在线观看| 亚洲天堂av无毛| 有码 亚洲区| 久久韩国三级中文字幕| 最近2019中文字幕mv第一页| 久久久久久人人人人人| 免费黄网站久久成人精品| 日韩精品免费视频一区二区三区| 女性生殖器流出的白浆| 美女大奶头黄色视频| 精品一区二区三区四区五区乱码 | 青春草亚洲视频在线观看| 9191精品国产免费久久| 热re99久久精品国产66热6| 欧美成人午夜免费资源| 另类亚洲欧美激情| 免费观看av网站的网址| 美女主播在线视频| 啦啦啦在线观看免费高清www| 人妻少妇偷人精品九色| 久久精品熟女亚洲av麻豆精品| 岛国毛片在线播放| 两性夫妻黄色片| 91久久精品国产一区二区三区| 国产成人免费无遮挡视频| 性色av一级| 赤兔流量卡办理| 成年美女黄网站色视频大全免费| 日韩不卡一区二区三区视频在线| av一本久久久久| 2021少妇久久久久久久久久久| av在线app专区| 99国产精品免费福利视频| 天天操日日干夜夜撸| 十分钟在线观看高清视频www| av国产精品久久久久影院| 国产精品免费视频内射| 国产亚洲av片在线观看秒播厂| 久久久久精品人妻al黑| 最新的欧美精品一区二区| 久久久国产一区二区| 一本久久精品| 亚洲美女搞黄在线观看| 亚洲精华国产精华液的使用体验| 考比视频在线观看| 黄频高清免费视频| 搡女人真爽免费视频火全软件| 亚洲av.av天堂| 自拍欧美九色日韩亚洲蝌蚪91| 国产成人精品婷婷| 国产精品一区二区在线不卡| 中文乱码字字幕精品一区二区三区| 男人舔女人的私密视频| 91精品三级在线观看| 2022亚洲国产成人精品| 色哟哟·www| 一级毛片我不卡| 日韩制服丝袜自拍偷拍| 欧美+日韩+精品| 一级,二级,三级黄色视频| 一区二区av电影网| 亚洲,一卡二卡三卡| 麻豆精品久久久久久蜜桃| 亚洲三级黄色毛片| 你懂的网址亚洲精品在线观看| 亚洲av成人精品一二三区| 一级毛片电影观看| 国产片特级美女逼逼视频| 自线自在国产av| 亚洲精品美女久久av网站| freevideosex欧美| 激情五月婷婷亚洲| 九色亚洲精品在线播放| 久久久a久久爽久久v久久| 亚洲成色77777| 亚洲欧美成人综合另类久久久| 亚洲色图综合在线观看| 999精品在线视频| 国产伦理片在线播放av一区| 一级片免费观看大全| 国产亚洲欧美精品永久| av卡一久久| 中文字幕制服av| 亚洲精品一二三| 国产精品 国内视频| av免费在线看不卡| videossex国产| 新久久久久国产一级毛片| 亚洲欧美成人精品一区二区| 国产极品天堂在线| 2018国产大陆天天弄谢| 不卡av一区二区三区| 在线观看一区二区三区激情| 人体艺术视频欧美日本| 日韩一本色道免费dvd| 一边亲一边摸免费视频| 天天躁日日躁夜夜躁夜夜| 999精品在线视频| 亚洲人成网站在线观看播放| 国产精品一区二区在线观看99| 69精品国产乱码久久久| 最近中文字幕高清免费大全6| 亚洲av综合色区一区| 午夜福利一区二区在线看| 精品人妻偷拍中文字幕| 免费久久久久久久精品成人欧美视频| 亚洲欧美一区二区三区久久| 精品99又大又爽又粗少妇毛片| 在线天堂中文资源库| 国产精品一国产av| 日韩成人av中文字幕在线观看| 妹子高潮喷水视频| 精品人妻熟女毛片av久久网站| av免费观看日本| 日日摸夜夜添夜夜爱| 自线自在国产av| 亚洲经典国产精华液单| 极品人妻少妇av视频| 午夜老司机福利剧场| 日本91视频免费播放| 人人妻人人澡人人爽人人夜夜| 国产熟女午夜一区二区三区| av电影中文网址| 91成人精品电影| 国产熟女欧美一区二区| 日韩三级伦理在线观看| 女性被躁到高潮视频| 日韩一区二区视频免费看| 电影成人av| 中文字幕最新亚洲高清| 婷婷成人精品国产| 深夜精品福利| 哪个播放器可以免费观看大片| 久久综合国产亚洲精品| 久久久久久人人人人人| 波野结衣二区三区在线| 久久久久久伊人网av| 亚洲综合精品二区| 美女中出高潮动态图| 最新的欧美精品一区二区| 久久ye,这里只有精品| 久久久久国产精品人妻一区二区| 1024香蕉在线观看| 国产精品偷伦视频观看了| 少妇猛男粗大的猛烈进出视频| 久久国产精品男人的天堂亚洲| 国产精品无大码| 9色porny在线观看| 国产乱来视频区| kizo精华| 欧美bdsm另类| 午夜福利影视在线免费观看| h视频一区二区三区| 午夜av观看不卡| 亚洲国产欧美网| 亚洲第一av免费看| 日韩一本色道免费dvd| 9色porny在线观看| 熟妇人妻不卡中文字幕| 日本-黄色视频高清免费观看| 欧美中文综合在线视频| 国产无遮挡羞羞视频在线观看| 亚洲av电影在线进入| 一本色道久久久久久精品综合| 少妇 在线观看| 国产一区二区激情短视频 | 青草久久国产| 精品少妇久久久久久888优播| 女人高潮潮喷娇喘18禁视频| 最近最新中文字幕大全免费视频 | 天天躁狠狠躁夜夜躁狠狠躁| 国产精品 国内视频| 国产av国产精品国产| 伊人久久国产一区二区| 美女大奶头黄色视频| 久久97久久精品| 国产日韩欧美在线精品| 国产精品国产三级专区第一集| 久久韩国三级中文字幕| 国产色婷婷99| 大香蕉久久成人网| 99久久综合免费| av视频免费观看在线观看| 免费在线观看黄色视频的| www日本在线高清视频| 午夜福利,免费看| 久久这里只有精品19| 97人妻天天添夜夜摸| 五月开心婷婷网| 亚洲综合色惰| 国产精品免费视频内射| 老汉色av国产亚洲站长工具| 可以免费在线观看a视频的电影网站 | 午夜福利,免费看| 精品少妇一区二区三区视频日本电影 | 1024视频免费在线观看| 午夜免费鲁丝| 男女午夜视频在线观看| 久久久久久久久久久久大奶| 久久久久久久精品精品| 亚洲国产精品国产精品| 国产精品 国内视频| 久久精品久久精品一区二区三区| 免费黄频网站在线观看国产| 日产精品乱码卡一卡2卡三| 国产免费视频播放在线视频| 日韩视频在线欧美| 国产精品一国产av| 春色校园在线视频观看| 午夜福利,免费看| 七月丁香在线播放| 久久人人爽人人片av| 国产高清不卡午夜福利| 97在线视频观看| 韩国精品一区二区三区| 日韩av不卡免费在线播放| 少妇猛男粗大的猛烈进出视频| 日韩电影二区| 26uuu在线亚洲综合色| 欧美日韩av久久| 久久久久久久久久久免费av| √禁漫天堂资源中文www| 自线自在国产av| 日日摸夜夜添夜夜爱| 日韩一区二区视频免费看| 亚洲美女视频黄频| 一二三四在线观看免费中文在| 女人久久www免费人成看片| 国产激情久久老熟女| 免费看av在线观看网站| 亚洲欧美一区二区三区黑人 | 国产成人精品一,二区| 中文天堂在线官网| videossex国产| 久久久久久久久久人人人人人人| 黄片播放在线免费| 熟女av电影| 国产日韩一区二区三区精品不卡| 熟妇人妻不卡中文字幕| 亚洲欧美一区二区三区黑人 | 啦啦啦在线免费观看视频4| 亚洲五月色婷婷综合| 丝袜在线中文字幕| 男人添女人高潮全过程视频| 大陆偷拍与自拍| 亚洲欧洲精品一区二区精品久久久 | 成人毛片a级毛片在线播放| 91精品国产国语对白视频| 国产不卡av网站在线观看| 国产男人的电影天堂91| 18禁观看日本| 一区二区三区精品91| 欧美日韩av久久| 午夜福利在线免费观看网站| 人人妻人人添人人爽欧美一区卜| 秋霞伦理黄片| 成人国产av品久久久| 国产欧美亚洲国产| www.精华液| 免费黄频网站在线观看国产| 国产精品一国产av| 欧美精品亚洲一区二区| 你懂的网址亚洲精品在线观看| 高清不卡的av网站| 国产又色又爽无遮挡免| av天堂久久9| 久久精品国产鲁丝片午夜精品| 欧美精品一区二区免费开放| 免费日韩欧美在线观看| 伦理电影大哥的女人| 国产综合精华液| 赤兔流量卡办理| 国产一区亚洲一区在线观看| 飞空精品影院首页| 国产精品国产三级专区第一集| 五月天丁香电影| 国产高清不卡午夜福利| 啦啦啦在线免费观看视频4| 国产伦理片在线播放av一区| 国产亚洲午夜精品一区二区久久| 天天躁狠狠躁夜夜躁狠狠躁| 日韩不卡一区二区三区视频在线| 超碰97精品在线观看| 视频在线观看一区二区三区| 老司机亚洲免费影院| 男女下面插进去视频免费观看| 国语对白做爰xxxⅹ性视频网站| 国产白丝娇喘喷水9色精品| 高清不卡的av网站| 日韩,欧美,国产一区二区三区| 我的亚洲天堂| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲国产欧美日韩在线播放| 久久99一区二区三区| 成人漫画全彩无遮挡| 另类亚洲欧美激情| 各种免费的搞黄视频| 97在线人人人人妻| 久久久精品国产亚洲av高清涩受| 18禁动态无遮挡网站| 欧美人与性动交α欧美精品济南到 | 激情五月婷婷亚洲| 80岁老熟妇乱子伦牲交| 亚洲欧洲日产国产| 国产xxxxx性猛交| 七月丁香在线播放| 涩涩av久久男人的天堂| 色吧在线观看| 国产精品麻豆人妻色哟哟久久| 中文字幕色久视频| 日本欧美视频一区| 又粗又硬又长又爽又黄的视频| 国产又爽黄色视频| 中文精品一卡2卡3卡4更新| 精品国产一区二区三区四区第35| xxx大片免费视频| 精品视频人人做人人爽| 午夜激情久久久久久久| 日韩欧美一区视频在线观看| 国产亚洲最大av| 黄片小视频在线播放| 青青草视频在线视频观看| 日韩在线高清观看一区二区三区| 亚洲欧洲精品一区二区精品久久久 | 菩萨蛮人人尽说江南好唐韦庄| 久久精品国产鲁丝片午夜精品| 七月丁香在线播放| 欧美av亚洲av综合av国产av | 免费不卡的大黄色大毛片视频在线观看| 免费黄频网站在线观看国产| 一级毛片 在线播放| av在线老鸭窝| 2018国产大陆天天弄谢| 午夜福利一区二区在线看| 国产高清国产精品国产三级| 中文字幕人妻丝袜一区二区 | 成人漫画全彩无遮挡| 人人妻人人爽人人添夜夜欢视频| 一区二区三区四区激情视频| 性高湖久久久久久久久免费观看|