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

    用中子星限制暗物質(zhì)粒子散射截面?

    2021-10-09 06:38:50王弋塵繆雪麗邵立晶
    天文學報 2021年5期
    關(guān)鍵詞:散射截面熱化玻色子

    王弋塵 繆雪麗 邵立晶

    (1 北京大學物理學院北京 100871)

    (2 北京大學科維理天文與天體物理研究所北京 100871)

    (3 中國科學院國家天文臺北京 100012)

    1 引言

    暗物質(zhì)在現(xiàn)代宇宙學中扮演著不可或缺的重要角色.暗物質(zhì)作為一種宇宙學假說得到了多項獨立證據(jù)的支持,包括星系旋轉(zhuǎn)曲線[1]、星系團的觀測[2]、宇宙微波背景輻射[3]等,這些證據(jù)幾乎涵蓋了宇宙學的各個尺度.因此,在現(xiàn)階段的研究中,暗物質(zhì)的存在得到了科學家們的廣泛認可.

    隨著支持暗物質(zhì)存在的證據(jù)越來越豐富,更多的研究者把注意力放在了探究暗物質(zhì)的組成和存在形式上.上世紀80年代,人們對暗物質(zhì)的存在形式主要持兩種觀點[4],即弱相互作用大質(zhì)量粒子(Weakly Interacting Massive Particle,WIMP)和暈族大質(zhì)量致密天體(Massive Compact Halo Object,MACHO);后者包括各種低亮度的致密天體,有黑洞、中子星、白矮星、褐矮星和自由行星等.上世紀90年代,MACHO研究項目期望以MACHO來解釋星系和星系團中對比實際觀測而言理論預測的質(zhì)量缺失[5],但是后續(xù)的觀測結(jié)果把銀河系暗物質(zhì)暈質(zhì)量中MACHO的占質(zhì)比限制到了8%以下[6],進而否定了暗物質(zhì)主要由MACHO組成的假設(shè).所以,人們認為暗物質(zhì)的主要組成部分可能是WIMP或其他暗物質(zhì)候選者.WIMP作為暗物質(zhì)粒子的重要候選者,只參與引力和弱核力相互作用,其相互作用在電弱標度,屬于超出粒子物理標準模型的新粒子[7].WIMP速度較低,因此不太容易克服引力相互作用,更容易聚集成團.

    目前,探測暗物質(zhì)的方法主要分為直接探測和間接探測兩種.直接探測是指通過探測穿過地球的暗物質(zhì)粒子流與地球上實驗裝置中的靶物質(zhì)的相互作用來研究暗物質(zhì)粒子的物理性質(zhì).2010年,DAMA/LIBRA (Dark Matter/Large Sodium Iodide Bulk for Rare Processes)的結(jié)果顯示[8],實驗探測到了暗物質(zhì)粒子與放置在地表深處的碘化鈉晶體中靶原子的相互作用的年調(diào)制信號.DAMA/LIBRA的實驗認為該信號的年調(diào)制特性是由于地球在銀河系暗物質(zhì)背景下的軌道運動引起的.同年,CoGeNT合作組也聲稱獲得了年調(diào)制的暗物質(zhì)信號[9].文獻[10]對DAMA/LIBRA和CoGeNT的結(jié)果進行分析,認為兩個實驗的觀測結(jié)果能用質(zhì)量為7 GeV、與核子的散射截面為2×10?40cm2的暗物質(zhì)粒子來解釋(本文取光速c=1,下同).

    2018年,同樣使用碘化鈉作為靶材料但靈敏度比DAMA/LIBRA更高的COSINE-100沒有得到DAMA/LIBRA的結(jié)果[11].并且,目前靈敏度最高的XENON1T實驗也沒有觀測到任何顯著的暗物質(zhì)粒子事件[12].XENON1T實驗對暗物質(zhì)粒子與核子的散射截面σnX的探測靈敏度很高,對于非自旋依賴(Spin Independent,SI)的WIMP,在質(zhì)量mX≈30 GeV的情況下,σSI的探測靈敏度達到了4.7×10?47cm2[12].因此,XENON1T實驗沒有發(fā)現(xiàn)WIMP信號的結(jié)果對WIMP的mX和σnX給出了較強的限制.此外,并沒有發(fā)現(xiàn)WIMP的XENON1T實驗卻在2020年發(fā)現(xiàn)了中性輕子暗物質(zhì)粒子候選者軸子(Axion)的疑似信號[13].未來,更新升級的XENONnT實驗[14]和正在準備中的PandaX-4T (Particle and Astrophysical Xenon Detector-4T)實驗[15]的探測靈敏度將比XENON1T實驗更高,但是通過直接探測實驗的方法找到WIMP蹤跡的前景尚未可知.

    此外,對于自旋依賴和非自旋依賴的暗物質(zhì)粒子,直接探測實驗給出的限制有數(shù)量級上的區(qū)別.對于相同質(zhì)量的暗物質(zhì)粒子,直接探測實驗給出的自旋依賴比非自旋依賴散射截面的限制要弱3–4個數(shù)量級(與靶材料的種類有關(guān)).這是因為,對于非自旋依賴相互作用,原子核中的核子對總散射截面的貢獻可以疊加;而對于自旋依賴的相互作用,由于原子核內(nèi)不同核子自旋會相互抵消,使得總散射截面下降,進而導致對截面探測靈敏度的下降[16].因此,直接探測實驗對自旋依賴的暗物質(zhì)粒子散射截面的探測靈敏度上的不足,還需要其他探測方法補充.

    除直接探測外,另一種重要的探測暗物質(zhì)的方法是間接探測.間接探測是指以宇宙為實驗室,突破地球?qū)嶒灥木窒扌?結(jié)合天文觀測和理論模型預測,對暗物質(zhì)粒子的物理性質(zhì)給出限制[17].間接探測的一個重要思路就是考察暗物質(zhì)粒子與天體的相互作用.具體而言,就是尋找難以被直接探測的暗物質(zhì)粒子在特定天體環(huán)境下的可觀測效應(yīng),或?qū)ふ野滴镔|(zhì)粒子對天體演化的可觀測的影響.

    天體與暗物質(zhì)的相互作用可以捕獲暗物質(zhì)粒子,從而吸積暗物質(zhì)粒子[18].以往研究已考察過暗物質(zhì)粒子和地球表面物質(zhì)或和太陽物質(zhì)的相互作用[19–20].中子星是更加致密的天體,具有極端的引力環(huán)境,極高的粒子密度和極強的電磁場,暗物質(zhì)粒子在其中參與的物理過程可能會更加豐富和極端.

    文獻[17]首次考慮了中子星吸積暗物質(zhì)可能會導致暗物質(zhì)自引力塌縮形成黑洞并吞噬中子星的模型.基于該模型,后續(xù)研究對致密天體白矮星和中子星進行了更加細致的討論,給出了在暗物質(zhì)粒子不同質(zhì)量mX下對應(yīng)的暗物質(zhì)粒子-中子散射截面σnX的限制[16].在強簡并條件下,玻色子可以發(fā)生玻色-愛因斯坦凝聚(Bose-Einstein Condensate,BEC),而費米子具有較高費米簡并壓.暗物質(zhì)粒子為玻色子或費米子對應(yīng)了其本身在被中子星吸積后的不同物理過程,二者有顯著差異,所以區(qū)分暗物質(zhì)粒子是玻色子或費米子是必要的.文獻[17]的模型對非對稱玻色子暗物質(zhì)的參數(shù)空間給出了較強的限制[21–22].而對于費米子暗物質(zhì),如果假定費米子暗物質(zhì)存在吸引自相互作用,那么也能對暗物質(zhì)的參數(shù)空間給出較強限制[23–24].

    本文基于中子星吸積暗物質(zhì)發(fā)生引力塌縮形成黑洞吞噬中子星的物理模型[17,22],給出了在不同的暗物質(zhì)粒子質(zhì)量mX下對暗物質(zhì)粒子–中子散射截面σnX的限制.本文考慮的暗物質(zhì)模型是WIMP模型,且考慮暗物質(zhì)粒子是玻色子的情況.本文還討論了暗物質(zhì)粒子有無自相互作用以及有無湮滅條件對給出的暗物質(zhì)參數(shù)的影響.對暗物質(zhì)參數(shù)進行限制的系統(tǒng),本文選取了已發(fā)現(xiàn)的中子星—鄰近太陽系的PSR J0437–4715和銀心附近的PSR J1745–2900,也考慮了兩個在具有高暗物質(zhì)背景密度的球狀星團中心和銀心附近可能存在的年老中子星.我們將文中結(jié)果與XENON1T實驗探測的結(jié)果比較可知,本文采用的方法可以給出更強的限制,從而可以促進人們對暗物質(zhì)的理解.

    2 物理模型

    我們采用中子星吸積暗物質(zhì)在自引力塌縮下形成黑洞并吞噬中子星的物理模型[17,22],并考慮如下物理過程:(1)中子星吸積暗物質(zhì)粒子;(2)暗物質(zhì)粒子熱化并聚集在中子星中心附近一個小尺度范圍內(nèi);(3)當中子星內(nèi)部聚集的暗物質(zhì)粒子達到一定數(shù)目,暗物質(zhì)粒子發(fā)生自引力塌縮,形成小型黑洞;(4)黑洞經(jīng)過演化最終吞噬中子星.概括來說,我們把整個過程簡稱為吸積、熱化、塌縮和吞噬.下面對這幾個過程給出具體的說明.

    2.1 吸積

    一顆中子星周圍分布著暗物質(zhì)粒子,其中一部分暗物質(zhì)粒子的運動軌道將穿過星體.暗物質(zhì)粒子和中子星上以中子為主的重子物質(zhì)會發(fā)生散射,使暗物質(zhì)粒子損失能量,進而被中子星的引力場捕獲.當粒子的軌道束縛在中子星內(nèi)部時,可以認為暗物質(zhì)粒子被中子星吸積.

    對于主序星,也存在上述過程.文獻[18]討論了太陽對于暗物質(zhì)粒子的吸積,這是我們選取模型的吸積速率方程的雛形.文中考慮暗物質(zhì)粒子運動服從麥克斯韋-玻爾茲曼分布,對不同能量和角動量的軌道進行積分,得到天體在單位時間內(nèi)吸積的暗物質(zhì)粒子個數(shù)[18],

    其中,nX是天體所處環(huán)境的背景暗物質(zhì)粒子數(shù)密度,是暗物質(zhì)粒子的速度彌散度,G是萬有引力常數(shù),MS和RS分別是天體的質(zhì)量和半徑,vesc是逃逸速度,α是無量綱因子,其定義為α=2mXmnβ/(+),其中mn是中子質(zhì)量,β是徑向位置的函數(shù),fσ是與散射截面有關(guān)的因子,具體表達式為[18]:

    其中,σcrit=mnS是天體的特征截面參數(shù),ρB是星體本身物質(zhì)密度,σnX是暗物質(zhì)粒子–中子散射截面.(2)式是對所有可能的軌道取平均值,最后通過積分得到.通過計算fσ因子,我們能得到穿過天體的暗物質(zhì)粒子發(fā)生散射次數(shù)的期望值.

    對于中子星吸積暗物質(zhì)的情況,文獻[25]給出單位時間內(nèi)中子星吸積的暗物質(zhì)粒子數(shù),

    其中,

    式中NB為中子星的總重子數(shù),ρX是中子星所在位置的暗物質(zhì)密度.(3)式和(1)式相對應(yīng),整理可得相似形式.兩個表達式的區(qū)別在于:一、(3)式在討論暗物質(zhì)碰撞的能量損失時考慮了狹義相對論,因此得到了B2=的參數(shù)項,該項僅在mX>106GeV時會顯著影響結(jié)果;二、由于中子星表面的中子處于一定程度的簡并態(tài),存在泡利阻塞(Pauli blocking),能量過低的碰撞因無法使中子能量達到更高能級而被禁止,因此有限制因子ξ=min;三、(3)式中截面相關(guān)的因子fσ由(4)式給出,其中σsat=(0.45NBξ),其中RNS是中子星的半徑.(4)與(2)式的區(qū)別在于,(4)式考慮了廣義相對論效應(yīng),即在施瓦西度規(guī)下進行計算.

    實際的暗物質(zhì)粒子積累數(shù)并不是勻速增加的.這是因為被中子星吸積的暗物質(zhì)粒子可能會在中子星內(nèi)發(fā)生湮滅和衰變.我們定義暗物質(zhì)粒子的壽命為τ,湮滅截面為〈σav〉,則中子星吸積的暗物質(zhì)粒子數(shù)量NX的變化率可以寫成Facc與衰變速率Fdecay以及湮滅速率Fanni的和[26],

    其中,Facc由(3)式得到.在不考慮湮滅情況下中子星吸積的暗物質(zhì)粒子數(shù)為[27]

    2.2 熱化

    暗物質(zhì)粒子在被吸積進中子星后,并不會迅速形成穩(wěn)定的結(jié)構(gòu).暗物質(zhì)粒子將在中子星內(nèi)部繼續(xù)運動,和中子星內(nèi)部物質(zhì)通過碰撞繼續(xù)交換能量,直至暗物質(zhì)粒子的溫度與中子星內(nèi)部溫度相同.暗物質(zhì)粒子和中子通過不斷相互作用達到熱平衡的過程就是熱化過程.暗物質(zhì)粒子在非相對論的中子費米氣體中的熱化時間為[30]

    熱化時間tth主要和暗物質(zhì)粒子的質(zhì)量,暗物質(zhì)粒子–中子散射截面和中子星的溫度相關(guān).當mX?mn或mX?mn時,tth很小,熱化會在較短的時間內(nèi)完成;如果mX~mn,要在相同的時間完成熱化,需要更大的σnX.最后一項TNS是中子星的核心溫度,一般取為105–106K.

    當熱化時間tth與中子星壽命相比,滿足

    這樣的情況下,暗物質(zhì)會在較短的時間內(nèi)熱化并聚集在中子星內(nèi)部.如果不滿足(12)式,即便中子星吸積了足夠數(shù)量的暗物質(zhì)粒子,但暗物質(zhì)粒子的聚集程度很低,在一定的聚集區(qū)域內(nèi),暗物質(zhì)粒子的質(zhì)量之和遠小于暗物質(zhì)粒子聚集區(qū)域內(nèi)存在的非暗物質(zhì)的總質(zhì)量.這樣一來,如果暗物質(zhì)無法聚集到滿足塌縮所需要的條件,那么后續(xù)的塌縮和吞噬過程都不會發(fā)生.我們稱(12)式為熱化條件.

    是否滿足熱化條件,不僅取決于暗物質(zhì)粒子的物理性質(zhì),也取決于中子星的年齡.年齡tNS越大的中子星,留給暗物質(zhì)粒子完成熱化的時間越充足.考慮一個年齡為tNS~1 Gyr的年老中子星,若mX~1 GeV,達成熱化條件所需的散射截面要滿足σnX>10?53cm2,這已是在所考慮的質(zhì)量區(qū)間上熱化條件對散射截面最嚴格的要求,相比于目前直接觀測給出的對散射截面的限制而言是比較容易滿足的.但考慮一個tNS~105yr的中子星,若mX~1 GeV,所需的散射截面為σnX>10?48cm2,熱化條件將較難滿足.

    若暗物質(zhì)粒子在中子星內(nèi)部完成熱化過程,達到熱平衡.假定自相互作用并不顯著,我們可以通過維力定理計算暗物質(zhì)粒子在中子星內(nèi)部主要的位置分布,即其典型熱化半徑rth[30]:

    其中k為玻爾茲曼常數(shù),ρB=1.3×1018kg·cm?3為中子星物質(zhì)的密度.可以注意到,當mX越大,熱化半徑rth越小.計算表明,對于暗物質(zhì)mX~1 GeV的情況,熱化半徑并不小.暗物質(zhì)粒子在這樣的聚集程度下發(fā)生引力塌縮還很困難,即正常情況下,一個年老中子星累積的暗物質(zhì)粒子質(zhì)量在r MN,此處的MN是指熱化半徑?jīng)Q定的球體內(nèi)的中子星重子物質(zhì)的總質(zhì)量.

    2.3 塌縮

    對于有自相互作用的玻色子,我們考慮自相互作用為λ|?|4的情況[30–32],這種自相互作用表現(xiàn)為排斥力.達到錢德拉塞卡極限發(fā)生引力塌縮時暗物質(zhì)粒子數(shù)應(yīng)滿足[22],

    其中NChand是達到錢德拉塞卡極限發(fā)生引力塌縮時的臨界暗物質(zhì)粒子數(shù),普朗克質(zhì)量mPl=1.22×1019GeV.當λ=0時,NChand=,該式回到無自相互作用的錢德拉塞卡極限;當λ >時,(14)式括號中的第2項占主導,此時有NChand

    除此之外,玻色子暗物質(zhì)要發(fā)生引力塌縮還需要滿足吸積的暗物質(zhì)的總質(zhì)量相對中子星原本物質(zhì)占主導的條件.具體而言,在熱化半徑rth決定的球體以內(nèi),暗物質(zhì)的總質(zhì)量Macc=NXmX應(yīng)該大于中子星原有物質(zhì)的質(zhì)量MN,即[16],

    這個條件為暗物質(zhì)無BEC的質(zhì)量占主導時的情況.當然,如果聚集的玻色子暗物質(zhì)粒子發(fā)生BEC,其聚集的空間尺度會大大縮小,輕易滿足(15)式的條件,使得后續(xù)塌縮過程容易進行.發(fā)生BEC時暗物質(zhì)粒子數(shù)需要滿足的條件為[33]

    其中NBEC為發(fā)生BEC時的臨界暗物質(zhì)粒子數(shù).我們假設(shè)PB≈0.3ρB為中子星非暗物質(zhì)部分提供的壓強.下文的計算表明,考慮暗物質(zhì)發(fā)生BEC將比考慮暗物質(zhì)不發(fā)生BEC給出更強的限制.

    對于費米子,要發(fā)生自引力塌縮,需要引力超過費米簡并壓帶來的排斥作用,即[16],

    其中kF是費米波矢,?是約化普朗克常數(shù),r和V分別為中子星吸積的費米子暗物質(zhì)聚集的半徑和體積.可以從(17)式中得到自引力塌縮所需的暗物質(zhì)粒子數(shù)Nfermi=這樣得到的Nfermi比無自相互作用(λ=0)的玻色子暗物質(zhì)的錢德拉塞卡極限倍,而對于總暗物質(zhì)粒子數(shù)為(14)式的自相互作用的玻色子暗物質(zhì),雖然NChand在0下也對mPl/mX因子有3次方的依賴關(guān)系,但由于,造成玻色子暗物質(zhì)塌縮極限NChand降低,使得這種情況下的玻色子暗物質(zhì)的塌縮極限比費米子暗物質(zhì)的塌縮極限更加容易達到.除非假定費米子暗物質(zhì)存在較強的吸引自相互作用來對抗費米簡并壓,不然無法對費米子暗物質(zhì)給出較強的限制.關(guān)于假定費米子暗物質(zhì)粒子存在吸引的自相互作用的研究,請見文獻[23];基于星系暗物質(zhì)暈結(jié)構(gòu)相關(guān)的研究,請見文獻[34].本文主要考慮暗物質(zhì)粒子是玻色子的情形.

    2.4 吞噬

    中子星內(nèi)部吸積的暗物質(zhì)粒子自引力塌縮形成小型黑洞后,黑洞的后續(xù)演化還需要滿足一定的條件才能吞噬整個中子星.演化過程中主要存在3種效應(yīng):黑洞對中子星原有物質(zhì)的吸積、黑洞的霍金輻射、黑洞對暗物質(zhì)粒子的吸積.由此,黑洞質(zhì)量的變化率為[22]

    其中MBH是黑洞的質(zhì)量.等式右邊第1項是黑洞對中子星原本物質(zhì)的Bondi吸積速率,第2項是霍金輻射速率,第3項是黑洞對暗物質(zhì)粒子的吸積速率;vs是中子星內(nèi)部的聲速,一般取為vs=0.1.黑洞的質(zhì)量隨時間的變化率隨著MBH單調(diào)遞增.由此可知,只要黑洞形成的初始時刻的質(zhì)量變化滿足

    那么后續(xù)吞噬過程就可以發(fā)生,并且黑洞將在短時標內(nèi)吞噬中子星.否則,黑洞的質(zhì)量將不斷減小,無法吞噬中子星.

    3 參數(shù)選取與計算結(jié)果

    只要明確了暗物質(zhì)粒子形成黑洞吞噬中子星所需要滿足的條件,即可選擇中子星系統(tǒng),計算并得到在確定的中子星系統(tǒng)下滿足塌縮條件對應(yīng)的暗物質(zhì)粒子質(zhì)量mX和暗物質(zhì)粒子–中子散射截面σnX的范圍.我們實際選取的中子星系統(tǒng)是仍舊存在的,并沒有被暗物質(zhì)塌縮形成的黑洞吞噬,這意味著我們可以通過實際觀測得到的中子星系統(tǒng)來排除滿足塌縮條件的暗物質(zhì)的參數(shù)空間.由此,我們可以對暗物質(zhì)的參數(shù)空間進行限制.

    3.1 中子星的參數(shù)

    本文考慮了4個不同的中子星系統(tǒng),參數(shù)選取如表1所示.此外考慮到我們計算精度有限,因此不需要過于精確的參數(shù),所以我們采用的系統(tǒng)都采用了中子星的典型參數(shù),如中子星質(zhì)量MNS=1.5M⊙,半徑RNS=11 km.

    表1 所選中子星系統(tǒng)的參數(shù)Table 1 Parameters of the selected neutron star systems

    中子星吸積的暗物質(zhì)總數(shù)量NX與中子星所處環(huán)境的暗物質(zhì)密度ρX和中子星的年齡tNS正相關(guān),因此處于較高暗物質(zhì)密度環(huán)境和較大年齡的中子星更適合用來限制暗物質(zhì)粒子的參數(shù)空間.表1中前兩個為實際觀測到的中子星,其中PSR J0437–4715是一顆距離地球150 pc的鄰近中子星,具有較大的tNS和較低的ρX;PSR J1745–2900是一顆距銀心0.1 pc的磁星,具有較大的ρX和較低的tNS.除了真實觀測到的中子星之外,我們表1中還列出了兩個假想的系統(tǒng),分別是考慮未來可能會在球狀星團M4的中心或銀心(GC)附近觀測到的年齡為tNS~1 Gyr的中子星.對表格中的后兩種系統(tǒng)而言,它們兼具大年齡和高暗物質(zhì)密度的條件,因此可以給出對暗物質(zhì)參數(shù)空間更強的限制.當然,這兩類系統(tǒng)的發(fā)現(xiàn),依賴于未來射電望遠鏡的具體情況.

    3.2 BEC對限制的影響

    假定暗物質(zhì)粒子是玻色子,暗物質(zhì)粒子形成黑洞吞噬中子星的全部過程能發(fā)生需要滿足的條件包括:

    (i)不考慮BEC,此時需要滿足的條件為以下4個公式:(12)、(14)、(15)和(19)式;

    (ii)考慮BEC,此時需要滿足的條件為以下4個公式:(12)、(14)、(16)和(19)式.

    假定玻色子暗物質(zhì)粒子僅存在很小的λ|?|4型的自相互作用,我們?nèi)ˇ?10?30.我們分別研究PSR J0437–4715和PSR J1745–2900,計算考慮和不考慮BEC兩種情況.

    這兩個系統(tǒng)給出的限制如圖1所示.深色和淺色實線上方包絡(luò)的深色和淺色區(qū)域分別為考慮和不考慮BEC得到的暗物質(zhì)參數(shù)空間限制區(qū)域(即被排除的參數(shù)空間).點虛線下方包絡(luò)的區(qū)域為無熱化區(qū)域,如果無熱化區(qū)域和限制區(qū)域發(fā)生重疊,那么重疊的區(qū)域不會被排除,仍舊是暗物質(zhì)粒子可能存在的參數(shù)空間.可以發(fā)現(xiàn),考慮BEC,計算得到在σnX=10?45cm2時,兩個系統(tǒng)對應(yīng)排除的暗物質(zhì)的質(zhì)量范圍分別是10?1GeV

    圖1 考慮BEC與否的計算結(jié)果對比圖.深色實線上方包絡(luò)的深色區(qū)域為考慮BEC后排除的暗物質(zhì)粒子的參數(shù)空間;淺色實線上方包絡(luò)的淺色區(qū)域為不考慮BEC時排除的暗物質(zhì)粒子的參數(shù)空間.點虛線下方包絡(luò)的區(qū)域為無熱化區(qū)域.對PSR J0437–4715和PSR J1745–2900這兩個不同熱化時間的系統(tǒng),如果無熱化區(qū)域和限制排除的區(qū)域發(fā)生重疊,那么重疊的區(qū)域不會被排除,仍舊是暗物質(zhì)粒子可能存在的參數(shù)空間.Fig.1 Comparison of the excluded parameter space of dark matter particles with or without BEC.The deep regions above deep solid lines are excluded parameter space of dark matter particles with BEC,and the light regions above light solid lines are excluded without BEC.The regions below dot-dash lines are the non-thermalization regions.For PSRs J0437–4715 and J1745–2900 whose thermalization times are different,if a non-thermalization region overlaps with an excluded region,the overlapped region is still valid.

    3.3 結(jié)果和分析

    我們考慮玻色子暗物質(zhì)粒子的λ|?|4自相互作用耦合常數(shù)λ的不同取值(λ={0,10?30,10?25,10?15}),考慮有無明顯的湮滅截面(〈σav〉={0,10?42cm3·s?1}),基于表1中的4個系統(tǒng)進行計算,得到的結(jié)果如圖2所示.

    圖2 基于4個中子星系統(tǒng)在考慮BEC情況下對暗物質(zhì)粒子的限制.曲線的顏色從深到淺,分別對應(yīng)λ={0, 10?30,10?25, 10?15};實線和點線分別對應(yīng)有無湮滅截面〈σav〉={0, 10?42 cm3·s?1}.圖中每條顏色深淺不同的實線和點線向上包絡(luò)的區(qū)域?qū)?yīng)的是排除的暗物質(zhì)參數(shù)空間.點虛線向下包絡(luò)的區(qū)域為無熱化區(qū)域.如果無熱化區(qū)域和限制排除的區(qū)域發(fā)生重疊,那么該重疊的區(qū)域不會被排除,仍舊是暗物質(zhì)粒子可能存在的參數(shù)空間.圖中深色虛線是運行278.8 d的XENON1T直接探測實驗對WIMP的限制曲線(90%置信水平),深色虛線上方的參數(shù)空間被XENON1T實驗排除[12].圖中的淺色和中淺色虛線分別為正在準備中的XENONnT直接探測實驗[14]和PandaX-4T直接探測實驗[15]預測的對WIMP的探測靈敏曲線(90%置信水平).Fig.2 Constraints on dark matter particles calculated from four different neutron stars with BEC.Solid curves,from dark to light,correspond to different self-interaction parameters,λ={0,10?30,10?25,10?15}.Solid curves and dotted curves correspond to 〈σav〉={0,10?42 cm3·s?1}respectively.The region above each curve is the excluded parameter space of dark matter particles.The regions below dot-dash lines are non-thermalization regions.If a non-thermalization region overlaps with an excluded region,the overlapped region is still valid.The parameter space above the deep dashed curve is excluded by the 278.8-day data collected with XENON1T direct search experiment at 90% confidence level[12].The light dashed curve and the medium light dashed curve correspond to the projected sensitivity of XENONnT direct search experiment[14] and PandaX-4T direct search experiment[15] at 90%confidence level respectively,and the two experiments are in preparation.

    由圖2可知,在不考慮自相互作用和湮滅截面時,4個系統(tǒng)對相同質(zhì)量的暗物質(zhì)粒子的散射截面均有較強的限制作用.這是因為對玻色子暗物質(zhì)而言,無排斥的相互作用和無湮滅的情況,會使得塌縮需要的暗物質(zhì)粒子數(shù)更少,因此更容易滿足錢德拉塞卡極限.對單個具體的中子星系統(tǒng)PSR J0437–4715的計算結(jié)果對參數(shù)空間的限制不強,主要是因為暗物質(zhì)密度ρX太低.在無自相互作用和無湮滅截面情形時,在10?3GeV10?45cm2;在mX~101GeV時達到了對σnX最強的限制,排除了σnX>10?52cm2.PSR J1745–2900是一個年齡較小的磁星.雖然ρX相對PSR J0437–4715增大了5個數(shù)量級,但是限制區(qū)域沒有明顯擴大.一方面是因為PSR J1745–2900吸積的暗物質(zhì)粒子總數(shù)被其40 kyr的年齡壓低,另一方面是因為無熱化區(qū)域隨年齡降低而增大,導致有部分無熱化區(qū)域與限制區(qū)域交疊,使得交疊部分限制失效.由于更大的ρX和tNS的中子星可以對暗物質(zhì)的參數(shù)空間有更好的限制,我們考慮未來可能會在球狀星團M4的中心或者銀心附近發(fā)現(xiàn)的年齡在Gyr量級的中子星.我們預測給出了可能對玻色子暗物質(zhì)參數(shù)空間的限制.其中最好的限制結(jié)果是來自于在銀心附近可能發(fā)現(xiàn)的年老中子星,在10?4GeV10?47cm2;但其最強限制因為受無熱化的影響,部分區(qū)域限制無效.對銀心處的中子星系統(tǒng)而言最好的限制是來自于λ=10?30,〈σav〉=0的情況,在mX~102GeV時,達到σnX~10?54cm2.下一代大型射電天文望遠鏡—平方公里陣列(Square Kilometre Array,SKA)—由于建在南非和澳大利亞,能對南天區(qū)進行全面的巡天,有望在銀心附近發(fā)現(xiàn)中子星,對暗物質(zhì)研究起到促進作用.

    4 分析與討論

    4.1 與直接探測結(jié)果的比較

    XENON1T實驗的精度相比DAMA/LIBRA和CoGeNT更高,但沒有測到顯著的暗物質(zhì)粒子事件,相對其他直接探測實驗而言給出了對暗物質(zhì)參數(shù)空間更強的限制.文獻[12]給出了在90%的置信水平下,運行278.8 d的XENON1T實驗對WIMP的限制曲線,為圖2中的深色虛線.

    由圖2所示,XENON1T對mX=101–103GeV的暗物質(zhì)粒子–核子的散射截面σnX有好的探測靈敏度.在這個質(zhì)量范圍內(nèi)對不同的暗物質(zhì)粒子質(zhì)量mX,對應(yīng)的σnX的探測靈敏度在10?45–10?47cm2范圍內(nèi).由于XENON1T實驗實際未探測到暗物質(zhì)粒子的信號,所以靈敏度曲線以上的區(qū)域就是XENON1T排除的暗物質(zhì)參數(shù)空間.對比圖2內(nèi)中子星系統(tǒng)和XENON1T實驗給出的限制曲線,中子星系統(tǒng)給出的對暗物質(zhì)參數(shù)的限制區(qū)域隨自相互作用強度變化和有無湮滅截面而改變,但總體而言,中子星系統(tǒng)每個曲線排除的區(qū)域是明顯大于XENON1T實驗排除的暗物質(zhì)參數(shù)空間范圍.

    本文也對比了正在準備中XENONnT直接探測實驗[14]和PandaX-4T直接探測實驗[15],這兩個實驗預測給出的對WIMP探測靈敏度要比XENON1T更高,其靈敏度由圖2中淺色和中淺色虛線所示.這兩個實驗都對mX=101–103GeV的暗物質(zhì)粒子–核子的散射截面σnX有更好的探測靈敏度,且它們的最大探測靈敏度比XENON1T實驗有兩個數(shù)量級左右的提升,但是總體上的結(jié)果仍沒有明顯優(yōu)于本文采用的間接探測方法給出的限制.

    并且,XENON1T等直接探測實驗給出的靈敏度是針對非自旋依賴截面而言的,對于自旋依賴截面,XENON1T等直接探測實驗的靈敏度還要進一步降低至少兩三個數(shù)量級.因此可以認為本文采用的間接探測方法的限制強度對于玻色子暗物質(zhì)而言已經(jīng)可以與XENON1T等直接探測實驗的靈敏度媲美并有所超越.換而言之,在文章假設(shè)的模型框架下,按照本文給出的限制,基于本文的假設(shè),以現(xiàn)有的實驗條件,通過XENON1T等直接探測實驗找到為玻色子的WIMP暗物質(zhì)的難度不低.

    當然,文中的模型與分析都是初步的,并且也依賴于暗物質(zhì)是否為費米子或者玻色子等一些具體的物理假設(shè).另外,相關(guān)假設(shè)都不可避免地含有天體物理環(huán)境參量所帶來的不確定性.所以,直接探測實驗與間接探測實驗在這個意義上是互相補充,相互促進的.

    4.2 自相互作用和湮滅對限制的影響

    由圖2可知,自相互作用和湮滅截面在計算中產(chǎn)生了不同的影響.更大的λ使得限制區(qū)域向mX更大的方向移動,同時伴隨著對σnX限制的減弱.這是因為,一方面由于錢德拉塞卡極限NChand∝,在mX太小時NChand會變大,塌縮條件更難達成,從而決定了小質(zhì)量區(qū)域的限制;mX太大時,NChand會被壓低,從而使黑洞的初始質(zhì)量下降,更容易因霍金輻射蒸發(fā)從而無法吞噬中子星,這一效應(yīng)則決定了大質(zhì)量區(qū)域的限制.隨著λ的增大,塌縮條件更難達成,為抵消λ增大的影響,mX增大,限制區(qū)域的左邊界向大質(zhì)量方向移動;另一方面,隨著λ的增大,更大的NChand使得黑洞初始質(zhì)量MBH更大,削弱了霍金輻射引起的蒸發(fā),讓黑洞更容易吞噬中子星,從而使得限制區(qū)域的右邊界也向大質(zhì)量方向移動.

    對于〈σav〉=10?42cm3·s?1的情況,會使得參數(shù)空間的限制減弱.更進一步,如果〈σav〉?10?42cm3·s?1,限制區(qū)域會消失.這是因為,較大的〈σav〉會增大湮滅速率,使得吸積和湮滅的動態(tài)平衡在更小的NX處發(fā)生,參見(8)式.NX上限的降低將導致塌縮更難發(fā)生,由此限制區(qū)域變小.而對于〈σav〉≤10?42cm3·s?1或tNS~105yr而言,湮滅速率尚不足以在中子星年齡內(nèi)和吸積速率達到動態(tài)平衡,這種情況下湮滅對限制強度的影響會不明顯.

    考慮N ≈FacctNS,對于有高吸積速率且年輕的中子星系統(tǒng),如PSR J1745–2900,湮滅截面的影響會更弱.因為PSR J1745–2900有很大的Facc和很小的tNS,湮滅速率更難升高至吸積速率的水平,難以達到吸積和湮滅的動態(tài)平衡.相反,對于低吸積速率,年老的中子星系統(tǒng),如PSR J0437–4715,如果受到湮滅截面的影響,系統(tǒng)很可能在極短的時標內(nèi)達到吸積和湮滅的動態(tài)平衡,使得后續(xù)中子星的年齡增加對于改善限制來說沒有太大的幫助.

    未來若有高吸積速率年老系統(tǒng)的觀測會幫助進一步提升我們對暗物質(zhì)參數(shù)空間的限制.考慮到可能存在湮滅截面的作用,尋找適合檢驗模型的中子星系統(tǒng)的重心應(yīng)該放在高吸積速率上.當然,年齡對熱化區(qū)域的影響也必須考慮在內(nèi).

    5 結(jié)論

    本文基于前人研究的中子星吸積暗物質(zhì)導致自引力塌縮形成黑洞并吞噬中子星的物理模型[22],用已發(fā)現(xiàn)的中子星(包括鄰近太陽系的PSR J0437–4715和銀心附近的PSR J1745–2900)和未來可能發(fā)現(xiàn)的球狀星團中心和銀心附近的年老中子星,計算并給出了不同暗物質(zhì)粒子質(zhì)量mX對應(yīng)的暗物質(zhì)粒子–中子散射截面σnX的限制.

    對比中子星吸積的玻色子暗物質(zhì)粒子是否發(fā)生BEC的情況,我們看到,考慮BEC會使得mX和σnX的限制范圍各增寬3個數(shù)量級以上.所以,考慮BEC會增強本模型給出的限制.不考慮自相互作用和湮滅截面的條件下,文中考慮的4個系統(tǒng)對應(yīng)的限制均在10?3GeV10?45GeV的玻色子暗物質(zhì).隨著系統(tǒng)暗物質(zhì)密度ρX的增大和中子星年齡tNS的提高,系統(tǒng)給出的限制將進一步增強,對于銀心附近的年老中子星,在mX~102GeV時,最強限制可達到10?54cm2的量級.此外,由于更小的σnX落在了無熱化區(qū)域內(nèi),限制強度將無法進一步提升.與直接探測實驗XENON1T對比,我們發(fā)現(xiàn)該間接探測方法對暗物質(zhì)粒子mX-σnX參數(shù)空間的限制強度超過了XENON1T給出的限制.與正在準備中的XENONnT和PandaX-4T實驗相比,本文采用的間接探測方法也有一定優(yōu)勢.但我們使用的方法給出的這些限制是基于一定假設(shè)條件的,包括:(1)暗物質(zhì)粒子是玻色子;(2)湮滅截面小;(3)不存在太強的排斥性自相互作用;(4)可能發(fā)生BEC等.這些條件需要根據(jù)暗物質(zhì)模型與中子星物理進行進一步的考察.

    考慮暗物質(zhì)粒子的λ|?|4型排斥自相互作用和湮滅截面會對不同系統(tǒng)對應(yīng)的限制強度造成不同程度的影響.考慮自相互作用會使限制區(qū)域向mX-σnX平面的右上側(cè)移動.這樣一來,有限制的區(qū)域?qū)?yīng)的mX更大,但對不同mX下σnX的限制有所減弱.湮滅截面對限制強度的影響更為復雜.湮滅截面較大時,中子星能在較短時標下達到吸積和湮滅的動態(tài)平衡,暗物質(zhì)密度不再增長,從而使得限制變?nèi)?因此湮滅截面對低吸積速率的年老中子星影響較大,而對高吸積速率的年輕中子星影響較小.高吸積速率的年輕中子星受熱化條件的影響更大.以上這些討論為未來在這個研究方向如何選擇合適的中子星系統(tǒng)有一定的啟示.并且我們期待未來可能通過SKA等觀測在銀心附近能觀測到年老中子星,這樣的系統(tǒng)的發(fā)現(xiàn)有助于增強本文使用的模型對暗物質(zhì)參數(shù)空間的限制,從而幫助人們進一步理解暗物質(zhì).

    猜你喜歡
    散射截面熱化玻色子
    W玻色子超重了
    大自然探索(2023年7期)2023-11-14 13:08:00
    紀念希格斯玻色子發(fā)現(xiàn)10周年
    科學(2022年4期)2022-10-25 02:43:30
    捕捉希格斯玻色子之手
    科學(2022年4期)2022-10-25 02:42:42
    長波紅外光學系統(tǒng)混合被動無熱化設(shè)計
    LHCb =8 TeV的Drell-Yan-Z→e+e-數(shù)據(jù)對部分子分布函數(shù)的影響
    基于微波倍頻源太赫茲頻段雷達散射截面測量
    雷達學報(2018年1期)2018-04-04 01:56:56
    熱化學方程式的書寫和理解
    115In中子非彈性散射截面的實驗測量及蒙特卡羅修正
    核技術(shù)(2016年4期)2016-08-22 09:05:22
    基于MTF的紅外鏡頭無熱化性能評估方法研究
    平板形目標的量子雷達散射截面計算
    制導與引信(2014年4期)2014-05-25 00:34:37
    国产欧美日韩精品一区二区| 免费看a级黄色片| 日韩三级伦理在线观看| 欧美成人a在线观看| 免费在线观看成人毛片| 赤兔流量卡办理| 大话2 男鬼变身卡| 有码 亚洲区| 亚洲精品国产av成人精品| 女的被弄到高潮叫床怎么办| 蜜桃久久精品国产亚洲av| 国产精品人妻久久久久久| 人妻一区二区av| 亚洲欧美一区二区三区黑人 | 午夜激情欧美在线| 欧美成人精品欧美一级黄| 尾随美女入室| 国产黄频视频在线观看| 日本av手机在线免费观看| 国产精品美女特级片免费视频播放器| 久久久久久久久中文| 边亲边吃奶的免费视频| 久久精品综合一区二区三区| 国产精品美女特级片免费视频播放器| 亚洲精品一二三| 亚洲av日韩在线播放| 丰满人妻一区二区三区视频av| 男女视频在线观看网站免费| 欧美日韩亚洲高清精品| 欧美+日韩+精品| 2022亚洲国产成人精品| 日韩成人av中文字幕在线观看| 国产乱来视频区| 国产激情偷乱视频一区二区| 汤姆久久久久久久影院中文字幕 | 国产 一区 欧美 日韩| 一个人看的www免费观看视频| 亚洲精品成人久久久久久| 欧美一区二区亚洲| 亚洲一区高清亚洲精品| 欧美激情久久久久久爽电影| 亚洲国产精品专区欧美| 好男人视频免费观看在线| 熟妇人妻不卡中文字幕| 免费看光身美女| 91av网一区二区| 欧美性感艳星| 免费观看性生交大片5| 国产人妻一区二区三区在| 精品久久久噜噜| av.在线天堂| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲av免费高清在线观看| 免费黄频网站在线观看国产| ponron亚洲| 夫妻午夜视频| 超碰97精品在线观看| 欧美日本视频| 在线观看一区二区三区| 性色avwww在线观看| 久久久a久久爽久久v久久| 九九在线视频观看精品| 少妇的逼好多水| 简卡轻食公司| 国产激情偷乱视频一区二区| 91av网一区二区| 女人十人毛片免费观看3o分钟| 国产精品伦人一区二区| 白带黄色成豆腐渣| 亚洲熟妇中文字幕五十中出| 好男人在线观看高清免费视频| 成人美女网站在线观看视频| 中文字幕av成人在线电影| 少妇的逼好多水| 亚洲欧美一区二区三区国产| 熟妇人妻久久中文字幕3abv| 亚洲国产精品专区欧美| 在线 av 中文字幕| 中文资源天堂在线| 日本午夜av视频| 69av精品久久久久久| 熟妇人妻不卡中文字幕| 国产免费一级a男人的天堂| 七月丁香在线播放| 亚洲四区av| 成人一区二区视频在线观看| 国产精品无大码| 免费观看性生交大片5| 国产91av在线免费观看| 国产午夜福利久久久久久| 国产成人精品婷婷| 美女高潮的动态| 不卡视频在线观看欧美| 丰满人妻一区二区三区视频av| 日本欧美国产在线视频| 国产成人一区二区在线| 国产成人免费观看mmmm| 午夜激情福利司机影院| 岛国毛片在线播放| 草草在线视频免费看| 久久午夜福利片| videos熟女内射| 国产精品一区二区性色av| 麻豆国产97在线/欧美| 99热全是精品| 国产精品一区二区在线观看99 | 国产v大片淫在线免费观看| 人妻制服诱惑在线中文字幕| 日韩欧美三级三区| 成年免费大片在线观看| 大香蕉久久网| 又粗又硬又长又爽又黄的视频| 国产美女午夜福利| 尾随美女入室| 欧美bdsm另类| 日本一二三区视频观看| 欧美日本视频| 久久久久久久久大av| 99久久精品国产国产毛片| 成人一区二区视频在线观看| 只有这里有精品99| 久久6这里有精品| 亚洲国产精品国产精品| 麻豆久久精品国产亚洲av| 久久久久久国产a免费观看| 韩国av在线不卡| 日韩国内少妇激情av| 精品久久久久久电影网| 免费在线观看成人毛片| 青春草视频在线免费观看| 久久久久网色| 国产午夜精品一二区理论片| 免费播放大片免费观看视频在线观看| 卡戴珊不雅视频在线播放| 日日啪夜夜爽| 亚洲av日韩在线播放| 免费观看a级毛片全部| 视频中文字幕在线观看| 亚洲av电影不卡..在线观看| 亚洲精品一二三| 国产一区二区在线观看日韩| 性色avwww在线观看| 啦啦啦韩国在线观看视频| 日产精品乱码卡一卡2卡三| 久久6这里有精品| 最近最新中文字幕大全电影3| 日韩欧美 国产精品| 搡老乐熟女国产| 亚洲国产成人一精品久久久| 国产乱人视频| 亚洲婷婷狠狠爱综合网| 久久久欧美国产精品| 美女大奶头视频| 精品久久久久久久末码| 亚洲精品一二三| 国内揄拍国产精品人妻在线| 狂野欧美白嫩少妇大欣赏| 特级一级黄色大片| 国产亚洲91精品色在线| 韩国高清视频一区二区三区| 免费观看在线日韩| 久久99热这里只频精品6学生| 又爽又黄a免费视频| 永久免费av网站大全| 免费大片18禁| 亚洲在线观看片| 国产极品天堂在线| 最近手机中文字幕大全| 中文字幕av在线有码专区| 日本色播在线视频| 国产淫片久久久久久久久| 国产成人一区二区在线| 久久人人爽人人爽人人片va| 亚洲不卡免费看| 亚洲第一区二区三区不卡| 国产高清国产精品国产三级 | 国产淫片久久久久久久久| 日日摸夜夜添夜夜添av毛片| 日韩欧美三级三区| 女的被弄到高潮叫床怎么办| 久久精品人妻少妇| 亚洲欧美清纯卡通| 伊人久久精品亚洲午夜| 中文欧美无线码| 六月丁香七月| 欧美人与善性xxx| 精品久久久久久电影网| 免费看av在线观看网站| 国产爱豆传媒在线观看| 一区二区三区乱码不卡18| 国产成人午夜福利电影在线观看| 国产精品一及| 免费在线观看成人毛片| 亚洲av日韩在线播放| 观看美女的网站| 乱码一卡2卡4卡精品| 婷婷色麻豆天堂久久| 少妇猛男粗大的猛烈进出视频 | 天天躁夜夜躁狠狠久久av| 最近中文字幕高清免费大全6| 亚洲在线自拍视频| 中文欧美无线码| 毛片女人毛片| 女人久久www免费人成看片| 中文精品一卡2卡3卡4更新| 一级av片app| av.在线天堂| 国产综合精华液| 亚洲av成人精品一区久久| 日本与韩国留学比较| 亚洲va在线va天堂va国产| 国产伦理片在线播放av一区| 2018国产大陆天天弄谢| 国产在线一区二区三区精| 日韩精品青青久久久久久| 美女内射精品一级片tv| 91精品一卡2卡3卡4卡| 久久精品夜夜夜夜夜久久蜜豆| 国产淫片久久久久久久久| 国产精品一区二区三区四区免费观看| 最近视频中文字幕2019在线8| av线在线观看网站| freevideosex欧美| 99久久精品国产国产毛片| 国产成人freesex在线| 久久人人爽人人爽人人片va| 午夜久久久久精精品| 亚洲精品,欧美精品| 国产在视频线精品| 看十八女毛片水多多多| 嫩草影院精品99| 最近最新中文字幕大全电影3| 日韩欧美精品v在线| 久久国内精品自在自线图片| 久久久久久久久久久免费av| 少妇裸体淫交视频免费看高清| 99热这里只有是精品50| 亚洲熟女精品中文字幕| 亚洲最大成人手机在线| 日韩国内少妇激情av| 国产国拍精品亚洲av在线观看| 久久久久久久大尺度免费视频| 纵有疾风起免费观看全集完整版 | 中文资源天堂在线| 国产精品日韩av在线免费观看| 亚洲精品乱码久久久v下载方式| 亚洲精品亚洲一区二区| 国产精品久久久久久精品电影小说 | 午夜老司机福利剧场| 欧美人与善性xxx| 一级片'在线观看视频| 欧美成人午夜免费资源| 亚洲精品久久久久久婷婷小说| 禁无遮挡网站| 在线a可以看的网站| 视频中文字幕在线观看| 成年女人在线观看亚洲视频 | 精品久久久精品久久久| 成年版毛片免费区| 午夜亚洲福利在线播放| 国产亚洲精品av在线| 亚洲自偷自拍三级| 国产在视频线在精品| 97人妻精品一区二区三区麻豆| 亚洲成人久久爱视频| 99视频精品全部免费 在线| 少妇裸体淫交视频免费看高清| 欧美潮喷喷水| av又黄又爽大尺度在线免费看| 成年av动漫网址| 日韩精品青青久久久久久| 大香蕉97超碰在线| 亚洲成人av在线免费| 欧美一区二区亚洲| 欧美日本视频| or卡值多少钱| a级一级毛片免费在线观看| 日韩大片免费观看网站| 亚洲精品乱码久久久v下载方式| 人人妻人人澡欧美一区二区| 国产黄色视频一区二区在线观看| 国产免费福利视频在线观看| 春色校园在线视频观看| av天堂中文字幕网| 卡戴珊不雅视频在线播放| 亚洲性久久影院| 日本猛色少妇xxxxx猛交久久| 99久久人妻综合| 97人妻精品一区二区三区麻豆| 在线免费观看不下载黄p国产| 夜夜看夜夜爽夜夜摸| 久久6这里有精品| 尾随美女入室| 欧美性感艳星| 久久亚洲国产成人精品v| 国产国拍精品亚洲av在线观看| 精品人妻一区二区三区麻豆| 亚洲美女视频黄频| 久久草成人影院| 蜜臀久久99精品久久宅男| 在线免费十八禁| 国产av码专区亚洲av| 一区二区三区四区激情视频| 中文资源天堂在线| av在线天堂中文字幕| 国产欧美日韩精品一区二区| 亚洲内射少妇av| 国产精品麻豆人妻色哟哟久久 | av专区在线播放| 国产欧美另类精品又又久久亚洲欧美| 午夜久久久久精精品| 国内揄拍国产精品人妻在线| 成人亚洲精品av一区二区| 国产黄频视频在线观看| 久久精品久久久久久久性| 亚洲精品自拍成人| 97精品久久久久久久久久精品| 寂寞人妻少妇视频99o| 性色avwww在线观看| 网址你懂的国产日韩在线| 校园人妻丝袜中文字幕| 国产一区二区三区综合在线观看 | 有码 亚洲区| 国产成人精品一,二区| 亚洲国产精品专区欧美| 天堂影院成人在线观看| 日韩大片免费观看网站| 99热这里只有是精品50| 亚洲熟妇中文字幕五十中出| 青春草视频在线免费观看| 爱豆传媒免费全集在线观看| 又黄又爽又刺激的免费视频.| 亚洲美女视频黄频| 日日啪夜夜爽| 丰满人妻一区二区三区视频av| 精品亚洲乱码少妇综合久久| 高清日韩中文字幕在线| 久久久久久久亚洲中文字幕| 久久人人爽人人爽人人片va| 99久久精品一区二区三区| 天堂√8在线中文| 久久久久久久亚洲中文字幕| 又粗又硬又长又爽又黄的视频| 国产一区有黄有色的免费视频 | 女的被弄到高潮叫床怎么办| 在线观看一区二区三区| 直男gayav资源| 婷婷六月久久综合丁香| 色尼玛亚洲综合影院| 国产不卡一卡二| 国产精品蜜桃在线观看| 成人毛片a级毛片在线播放| 校园人妻丝袜中文字幕| 国产探花在线观看一区二区| 久久久久精品久久久久真实原创| 在现免费观看毛片| 久久久久久久久久成人| 亚洲熟妇中文字幕五十中出| 狠狠精品人妻久久久久久综合| av在线观看视频网站免费| 亚洲精品成人av观看孕妇| 91午夜精品亚洲一区二区三区| 久99久视频精品免费| 日韩欧美 国产精品| 麻豆成人午夜福利视频| 男女啪啪激烈高潮av片| 日韩一区二区三区影片| 哪个播放器可以免费观看大片| 国产一级毛片在线| 少妇被粗大猛烈的视频| 久久精品国产鲁丝片午夜精品| 草草在线视频免费看| 亚洲精品一区蜜桃| 免费不卡的大黄色大毛片视频在线观看 | 蜜臀久久99精品久久宅男| 非洲黑人性xxxx精品又粗又长| 国产乱人视频| 色视频www国产| av播播在线观看一区| 超碰97精品在线观看| 久久国产乱子免费精品| 高清毛片免费看| 婷婷六月久久综合丁香| 日本一本二区三区精品| 亚洲欧洲日产国产| av在线亚洲专区| 日韩成人av中文字幕在线观看| 国产探花极品一区二区| 男女视频在线观看网站免费| 在线天堂最新版资源| 亚洲最大成人手机在线| 国产乱人视频| 免费看av在线观看网站| 成人国产麻豆网| 亚洲美女搞黄在线观看| 国产 亚洲一区二区三区 | 又爽又黄a免费视频| 我的女老师完整版在线观看| 成人综合一区亚洲| 国产精品日韩av在线免费观看| 久久久久久久久中文| 青春草亚洲视频在线观看| 欧美日本视频| 亚洲aⅴ乱码一区二区在线播放| 一个人看的www免费观看视频| 成人无遮挡网站| 一级毛片黄色毛片免费观看视频| 好男人在线观看高清免费视频| 国产爱豆传媒在线观看| 麻豆av噜噜一区二区三区| 国产真实伦视频高清在线观看| 一本一本综合久久| 夫妻性生交免费视频一级片| 免费观看的影片在线观看| 国产毛片a区久久久久| 精品一区在线观看国产| 亚洲人成网站在线观看播放| 久久综合国产亚洲精品| 久久久久久久久久久免费av| 2021少妇久久久久久久久久久| 身体一侧抽搐| 少妇熟女欧美另类| 美女主播在线视频| videos熟女内射| 国产国拍精品亚洲av在线观看| 中文字幕av在线有码专区| 免费看美女性在线毛片视频| 国产精品一区二区三区四区免费观看| 亚洲精品国产av成人精品| 1000部很黄的大片| 韩国高清视频一区二区三区| 亚洲欧洲国产日韩| 夜夜爽夜夜爽视频| 一级二级三级毛片免费看| 亚洲最大成人中文| 丰满少妇做爰视频| 亚洲精品aⅴ在线观看| 午夜激情欧美在线| 国产一区二区在线观看日韩| 国产色爽女视频免费观看| 国产成人福利小说| 99re6热这里在线精品视频| 久久韩国三级中文字幕| 午夜福利在线观看免费完整高清在| 成年人午夜在线观看视频 | 亚洲国产精品成人久久小说| 人人妻人人澡欧美一区二区| 99久久人妻综合| 五月玫瑰六月丁香| 中文资源天堂在线| 亚洲高清免费不卡视频| 亚洲精品aⅴ在线观看| 99热这里只有是精品在线观看| 精品久久久久久久久亚洲| 精品亚洲乱码少妇综合久久| 一级a做视频免费观看| 亚洲av免费在线观看| 日韩电影二区| 一个人免费在线观看电影| 看十八女毛片水多多多| 精品一区二区三区视频在线| 国产精品无大码| 直男gayav资源| 国产av不卡久久| 18+在线观看网站| 成人午夜高清在线视频| 国产一级毛片七仙女欲春2| 在线免费观看不下载黄p国产| 97人妻精品一区二区三区麻豆| 精品人妻一区二区三区麻豆| 久久精品综合一区二区三区| 久久久久国产网址| 久久久色成人| 亚洲欧美日韩无卡精品| 精品人妻一区二区三区麻豆| 精品一区在线观看国产| 女人久久www免费人成看片| 男人和女人高潮做爰伦理| 免费高清在线观看视频在线观看| 色尼玛亚洲综合影院| 日韩av免费高清视频| 欧美激情国产日韩精品一区| 一区二区三区乱码不卡18| 中文字幕av成人在线电影| 晚上一个人看的免费电影| 97超视频在线观看视频| 日韩在线高清观看一区二区三区| 六月丁香七月| 男女啪啪激烈高潮av片| 熟妇人妻不卡中文字幕| 精品久久久久久久久av| 欧美激情久久久久久爽电影| 久久久久网色| 精华霜和精华液先用哪个| 国产精品麻豆人妻色哟哟久久 | 夫妻性生交免费视频一级片| 亚洲欧美一区二区三区国产| 久久精品夜夜夜夜夜久久蜜豆| 日韩视频在线欧美| 嫩草影院新地址| 精品国产一区二区三区久久久樱花 | 97人妻精品一区二区三区麻豆| 老女人水多毛片| 嫩草影院精品99| .国产精品久久| 亚洲欧美中文字幕日韩二区| 3wmmmm亚洲av在线观看| 大又大粗又爽又黄少妇毛片口| 在线观看美女被高潮喷水网站| 美女脱内裤让男人舔精品视频| 国产中年淑女户外野战色| 日韩伦理黄色片| 久久精品熟女亚洲av麻豆精品 | 免费av观看视频| 毛片一级片免费看久久久久| 91久久精品电影网| 如何舔出高潮| 少妇的逼好多水| 免费大片黄手机在线观看| 好男人在线观看高清免费视频| 久久国产乱子免费精品| 欧美97在线视频| 久久精品久久精品一区二区三区| 麻豆久久精品国产亚洲av| 国产精品伦人一区二区| 国产成人精品婷婷| www.av在线官网国产| 亚洲国产色片| 少妇被粗大猛烈的视频| 美女被艹到高潮喷水动态| 一二三四中文在线观看免费高清| 欧美精品一区二区大全| 久久久久久久大尺度免费视频| 午夜激情久久久久久久| 黄色配什么色好看| 久99久视频精品免费| 国内揄拍国产精品人妻在线| 小蜜桃在线观看免费完整版高清| 欧美xxxx黑人xx丫x性爽| 久久久午夜欧美精品| 极品少妇高潮喷水抽搐| 18禁在线播放成人免费| 欧美xxⅹ黑人| 亚洲国产成人一精品久久久| 国产成人福利小说| 黄色日韩在线| 韩国av在线不卡| 色吧在线观看| 十八禁网站网址无遮挡 | 亚洲精品色激情综合| 日日摸夜夜添夜夜添av毛片| 男插女下体视频免费在线播放| 黄片无遮挡物在线观看| 国产色婷婷99| 卡戴珊不雅视频在线播放| 亚洲精品亚洲一区二区| videos熟女内射| 综合色丁香网| 亚洲av日韩在线播放| 色视频www国产| av又黄又爽大尺度在线免费看| 天天一区二区日本电影三级| 亚洲国产色片| 国产伦精品一区二区三区视频9| 80岁老熟妇乱子伦牲交| 精品久久久精品久久久| 亚洲熟女精品中文字幕| 可以在线观看毛片的网站| 性插视频无遮挡在线免费观看| 日韩av在线免费看完整版不卡| 尤物成人国产欧美一区二区三区| 听说在线观看完整版免费高清| 天堂中文最新版在线下载 | 国产一区有黄有色的免费视频 | 国产成人午夜福利电影在线观看| 一本久久精品| 亚洲精品成人久久久久久| 天堂av国产一区二区熟女人妻| 欧美成人a在线观看| 亚洲精品,欧美精品| 国产免费又黄又爽又色| 卡戴珊不雅视频在线播放| 午夜福利在线在线| 国产成人a区在线观看| 69av精品久久久久久| freevideosex欧美| 国产成人精品一,二区| 卡戴珊不雅视频在线播放| 非洲黑人性xxxx精品又粗又长| 亚洲精品aⅴ在线观看| 久久精品国产亚洲av天美| 综合色av麻豆| 欧美日韩一区二区视频在线观看视频在线 | 丝袜美腿在线中文| 神马国产精品三级电影在线观看| 老师上课跳d突然被开到最大视频| 观看美女的网站| 22中文网久久字幕| 精品久久久久久久久久久久久| 亚洲av一区综合| 午夜福利高清视频| 国产综合懂色| 69av精品久久久久久| 国产精品国产三级专区第一集| 国产精品麻豆人妻色哟哟久久 | 欧美精品一区二区大全| 又爽又黄a免费视频| 久久99热这里只频精品6学生| 男女边摸边吃奶| 国产高清三级在线| 亚洲国产精品成人综合色| av在线观看视频网站免费| 免费人成在线观看视频色| 最后的刺客免费高清国语| 欧美日韩在线观看h|