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

    活動(dòng)星系核動(dòng)量驅(qū)動(dòng)反饋機(jī)制的研究

    2016-07-11 01:35:38張小霞
    天文研究與技術(shù) 2016年1期
    關(guān)鍵詞:外流

    付 艷,張小霞

    (中國(guó)科學(xué)院國(guó)家天文臺(tái),北京 100012)

    ?

    活動(dòng)星系核動(dòng)量驅(qū)動(dòng)反饋機(jī)制的研究

    付艷,張小霞

    (中國(guó)科學(xué)院國(guó)家天文臺(tái),北京100012)

    摘要:觀測(cè)中發(fā)現(xiàn)星系中心超大質(zhì)量黑洞的質(zhì)量與星系核球速度彌散或星系核球質(zhì)量之間存在緊致關(guān)系。一般認(rèn)為這類關(guān)系是由活動(dòng)星系核反饋造成的,但是活動(dòng)星系核反饋過程的微觀物理實(shí)現(xiàn)并不清楚。甚至活動(dòng)星系核反饋的提出者(Silk)也對(duì)其有效性提出質(zhì)疑。利用星系的觀測(cè)數(shù)據(jù)重新考慮了活動(dòng)星系核的反饋過程,檢驗(yàn)了活動(dòng)星系核反饋的有效性。收集了一個(gè)包括29個(gè)星系的樣本,通過樣本星系的觀測(cè)數(shù)據(jù)對(duì)它們的勢(shì)場(chǎng)和由動(dòng)量反饋而致的氣體殼層在勢(shì)場(chǎng)中的運(yùn)動(dòng)做了更準(zhǔn)確的研究,發(fā)現(xiàn)絕大多數(shù)星系的動(dòng)量反饋是有效的。

    關(guān)鍵詞:活動(dòng)星系核反饋;動(dòng)量驅(qū)動(dòng)反饋;外流;康普頓冷卻

    星系中心普遍存在超大質(zhì)量的黑洞[1-3]。觀測(cè)發(fā)現(xiàn)星系中超大質(zhì)量黑洞的質(zhì)量MBH與星系核球速度彌散σ或者星系核球質(zhì)量Mb之間存在緊致關(guān)系,即MBH-σ關(guān)系或MBH-Mb關(guān)系[2,4-12]。這種緊致關(guān)系預(yù)示著超大質(zhì)量黑洞與星系相互影響共同演化。人們一般采用活動(dòng)星系核反饋機(jī)制解釋MBH-σ(或MBH-Mb)關(guān)系。活動(dòng)星系核反饋機(jī)制包括能量驅(qū)動(dòng)反饋機(jī)制[13-14]和動(dòng)量驅(qū)動(dòng)反饋機(jī)制[15-17]。能量驅(qū)動(dòng)反饋是指黑洞吸積釋放的總能量能夠有效地轉(zhuǎn)移為星系中氣體整體的熱能或動(dòng)能,并足以使這些氣體擺脫星系束縛而逃出星系,這一機(jī)制可導(dǎo)致MBH∝σ5。動(dòng)量反饋則是指由黑洞吸積導(dǎo)致的風(fēng)或外流將動(dòng)量傳遞給周圍氣體,并掃除這些氣體從而終止了黑洞吸積和星系成長(zhǎng),并預(yù)言MBH∝σ4。目前尚不清楚活動(dòng)星系核反饋究竟是能量還是動(dòng)量反饋機(jī)制主導(dǎo)以及反饋機(jī)制的具體微觀物理過程。

    根據(jù)星系的黑洞質(zhì)量和速度彌散測(cè)量,文[1]采用等溫球模型研究活動(dòng)星系核反饋的有效性。他們估計(jì)了氣體從星系中逃出時(shí)所需的能量約為Mgcσ(c和Mg分別為光速和星系中氣體的總質(zhì)量),并將其與黑洞吸積釋放的總能量ηMBHc2(η=0.1為質(zhì)能轉(zhuǎn)換或輻射效率)進(jìn)行對(duì)比,發(fā)現(xiàn)單純由黑洞吸積形成的風(fēng)可能不足以將星系中氣體全部吹出星系。因此活動(dòng)星系核的反饋可能無(wú)效。

    文[1]采用等溫球模型描述星系的質(zhì)量密度分布和勢(shì)場(chǎng)分布。等溫球分布在靠近星系核球中心位置密度與實(shí)際星系的密度分布相比過高,而且在星系外邊界密度并不收斂。因而文[1]可能過高地估計(jì)了氣體的逃逸速度ve(r)。另外文[1]采用(res-r)/ve(res)(res為氣體逃逸半徑)估計(jì)氣體逃逸時(shí)標(biāo)τ時(shí),過高估計(jì)了氣體的逃逸時(shí)標(biāo)。由于星系中氣體的逃逸速度隨星系中心的距離降低,因而文[1]有可能過高估計(jì)了氣體逃出活動(dòng)星系核所需的時(shí)間以及所需的能量。總之,對(duì)于文[1]計(jì)算的ηMBHc2/Mgcσ<1,其核活動(dòng)的反饋相應(yīng)未必真的無(wú)效。為解釋此問題,收集了一個(gè)星系樣本,每一個(gè)樣本星系有觀測(cè)給出的面亮度分布數(shù)據(jù)。根據(jù)這些觀測(cè)可以得到星系的密度分布輪廓,并估計(jì)星系的寄主暗物質(zhì)暈質(zhì)量及暗物質(zhì)密度分布和星系勢(shì)場(chǎng)。利用這些密度分布,重新考慮了在這些樣本星系中的動(dòng)量反饋過程,并檢驗(yàn)了動(dòng)量反饋過程的有效性。

    論文第1節(jié)引入觀測(cè)樣本以及有關(guān)星系面亮度分布的觀測(cè)數(shù)據(jù)[8,10-12,18-22];第2節(jié)根據(jù)面亮度分布得到樣本星系的物質(zhì)密度分布,并估計(jì)暗物質(zhì)暈的大小和密度分布;第3節(jié)重新考慮活動(dòng)星系核動(dòng)量反饋過程及其有效性;第4節(jié)探討氣體向外擴(kuò)散時(shí)康普頓冷卻對(duì)活動(dòng)星系核反饋的影響;第5節(jié)討論影響動(dòng)量反饋有效性的幾個(gè)因素;第6節(jié)得出結(jié)論。

    1星系樣本

    為檢驗(yàn)活動(dòng)星系核反饋的有效性,收集了29個(gè)星系的相關(guān)觀測(cè)數(shù)據(jù)(表1)。這29個(gè)星系核區(qū)部分(0.1″-10″)均由哈勃望遠(yuǎn)鏡觀測(cè)。哈勃觀測(cè)的星系面亮度分布I(R)可以由Nuker law擬合,即

    (1)

    式中,α、β和γ為Nuker law分布的擬合參數(shù);R為到星系中心的二維投影距離;rb和Ib分別為V波段系面亮度的截?cái)喟霃胶徒財(cái)喟霃教幍拿媪炼?。?dāng)R?rb時(shí),I(R)∝R-β;當(dāng)R?rb時(shí),I(R)∝R-γ。各樣本星系V波段的相關(guān)擬合參數(shù)[8,10-12,18-22]在表1中列出,文中選取的數(shù)據(jù)采用宇宙學(xué)一致模型,即哈勃常數(shù)H0為70 km s-1/Mpc,物質(zhì)相對(duì)密度為Ωm=0.27,暗能量相對(duì)密度為ΩΛ=0.73。

    對(duì)于核區(qū)之外部分,星系的面亮度分布可以由de Vaucouleurs給出,即

    (2)

    式中,Re為V波段星系的半光度半徑;Ie為Re處星系的面亮度(Re的數(shù)據(jù)來(lái)自文[19])。

    為研究活動(dòng)星系核反饋機(jī)制的有效性,需要知道星系中心黑洞的質(zhì)量。表1中一部分樣本星系中心黑洞的質(zhì)量直接來(lái)自測(cè)量值,而對(duì)于其它沒有黑洞直接測(cè)量的樣本星系則通過MBH-σ關(guān)系,即

    (3)

    給出黑洞質(zhì)量估計(jì)。參數(shù)(α,β)在Mbulge≤1011.5M⊙時(shí),取(8.40 ± 0.09, 5.08 ± 0.7);在Mbulge>1011.5M⊙時(shí),取(8.52 ± 0.47, 4.69 ± 2.69 )[10]。

    2星系的物質(zhì)分布和勢(shì)場(chǎng)分布

    星系內(nèi)的物質(zhì)包括星系中心的超大質(zhì)量黑洞、恒星、氣體和暗物質(zhì)暈。若假設(shè)星系為球?qū)ΨQ分布,除中心超大質(zhì)量黑洞外,星系內(nèi)的物質(zhì)密度分布可表述為

    (4)

    式中,ρ*(r)為恒星質(zhì)量密度分布;ρg(r)為氣體的密度分布;ρh(r)為暗物質(zhì)質(zhì)量密度分布。由于考慮橢圓星系或核球,相對(duì)于恒星成分和暗物質(zhì)成分,氣體成分可能并不重要(例如10%ρ*(r)),因此在勢(shì)場(chǎng)的計(jì)算中,忽略氣體的貢獻(xiàn)。

    2.1恒星質(zhì)量密度輪廓

    星系的質(zhì)量密度分布ρ*(r)可以通過對(duì)星系面亮度分布I作積分變換直接得到[23],即

    (5)

    式中,Υ指星系中恒星的質(zhì)量與星系光度比;R為到星系中心的投影距離;dI/dR是星系面亮度隨R變化的梯度。

    (6)

    根據(jù)(6)式可以求得星系的外半徑R0。通過對(duì)(5)式進(jìn)行積分可以得到半徑r以內(nèi)的恒星總質(zhì)量為

    (7)

    2.2暗物質(zhì)質(zhì)量密度輪廓

    星系的暗物質(zhì)質(zhì)量密度分布一般可以由Nav-

    arro-Frenk-White(NFW)密度輪廓給出[24],即

    (8)

    (9)

    式中,M*,0是恒星特征質(zhì)量;M1是暗物質(zhì)暈的特征質(zhì)量;γ、δ和β均為擬合參數(shù)。根據(jù)文[26],M*,0、M1、β、δ和γ分別為10.72、12.35、0.43、0.56和1.54。對(duì)于表1中M*,total>1011.5M⊙的5個(gè)星系,它們都存在于室女座星系團(tuán)中,其寄主暗物質(zhì)暈的質(zhì)量約為Mh~1.2×1015M⊙[27]。

    根據(jù)星系及其寄主暗物質(zhì)暈的質(zhì)量密度分布可以得到它們的勢(shì)場(chǎng)分布[23]:

    (10)

    圖2展示了只考慮恒星物質(zhì)或考慮暗物質(zhì)暈的情況下,星系中不同位置r處的氣體或粒子的逃逸速度ve(r),以及這兩種情況下星系外半徑R0處的逃逸速度ve(R0)。由圖2可知,當(dāng)考慮星系周圍暗物質(zhì)暈時(shí),星系氣體的逃逸速度更大,即氣體要擺脫星系和暗物質(zhì)暈的束縛更加困難。

    圖1星系質(zhì)量密度分布。紅色實(shí)線代表星系的恒星質(zhì)量密度分布,藍(lán)色點(diǎn)線代表包括恒星和暗物質(zhì)質(zhì)量密度分布,綠色虛線代表等溫球模型下星系質(zhì)量密度分布

    Fig.1Density profiles for those galaxies in Table 1. The red solid lines represent the model only considering stellar mass, the blue dotted lines represent those considering both stellar mass and dark matter mass, and the green dashed lines represent the profiles given by the isothermal model圖2星系中氣體或粒子的逃逸速度隨半徑的分布。紅色實(shí)線和藍(lán)色點(diǎn)線分別代表只考慮恒星物質(zhì)和考慮恒星物質(zhì)和暗物質(zhì)的情況;黑色實(shí)心點(diǎn)和黑色空心點(diǎn)分別代表只考慮恒星物質(zhì)和考慮恒星物質(zhì)及其寄主暗物質(zhì)暈的情況,星系外半徑R0處的逃逸速度

    Fig.2Escape velocitiesat different radii. The blue dotted lines show the results obtained by considering both the contribution from the stars and the dark matter halo to the potential, and the red solid lines show the results obtained by only considering the contribution of the stars to the potential. The open (Full) circles represent escape velocities at the galaxy outer boundaryR0by considering both stars and dark matter halo (or only stars)

    3動(dòng)量驅(qū)動(dòng)反饋

    動(dòng)量反饋可能是導(dǎo)致MBH-σe關(guān)系的主導(dǎo)機(jī)制[16-17]。由于核活動(dòng)性輸出的動(dòng)量~fpL(t)/c(或動(dòng)能)會(huì)將星系核心處和星系中的氣體擠壓至氣體殼層(r→r+dr)并向外驅(qū)離此氣體殼層,這里L(fēng)(t)是活動(dòng)星系核的光度,觀測(cè)表明fp為量級(jí)在1左右的常數(shù)因子[28],標(biāo)志外流的動(dòng)量輸出與輻射光子總動(dòng)量間的比例。若只考慮恒星物質(zhì)的影響,根據(jù)動(dòng)量守恒定律可有

    (11)

    M*(

    (12)

    式中,ξ=fp/fg。

    星系一般存在于暗物質(zhì)暈之中,若考慮暗物質(zhì)暈對(duì)氣體的引力束縛作用,上述動(dòng)量方程則應(yīng)替換為

    (13)

    式中,Mh(

    (14)

    在核活動(dòng)觸發(fā)后,給定L(t)的演化形式,星系核心的氣體由于核活動(dòng)的動(dòng)量輸出推至氣體殼層r(t)→r(t)+dr(t)中。此殼層隨時(shí)間演化可以由上述(12)式和(14)式給出。在核活動(dòng)結(jié)束時(shí),氣體殼層的半徑達(dá)到rf,速度為vf,若vf大于rf處的逃逸速度ve(rf),則認(rèn)為核活動(dòng)性的動(dòng)量反饋?zhàn)銐驈?qiáng)以至于可以將星系中的氣體驅(qū)離星系。若vf

    求解(12)式和(14)式得出部分星系中氣體能夠克服恒星(或包括暗物質(zhì)暈)的引力束縛所能達(dá)到的速度vf和能夠到達(dá)的位置rf。核活動(dòng)觸發(fā)后星系中氣體外流,樣本星系中部分星系的活動(dòng)星系核單位時(shí)間內(nèi)輸出的動(dòng)量大于氣體受到來(lái)自恒星(或包括暗物質(zhì))的引力束縛,氣體動(dòng)量會(huì)一直增加,當(dāng)核活動(dòng)終止時(shí),有些星系氣體已成功地逃出星系束縛(rf>R0)(或暗物質(zhì)暈的維里半徑rf>rvir);部分星系的氣體動(dòng)量先增加,活動(dòng)星系核演化到一定時(shí)間,單位時(shí)間內(nèi)來(lái)自中心引擎的動(dòng)量小于來(lái)自恒星(或包含暗物質(zhì))的引力束縛時(shí),氣體動(dòng)量減少。這些星系中有些在核活動(dòng)終止時(shí),氣體外流的速度vf依然大于氣體的逃逸速度ve(rf),氣體就能夠成功地?cái)[脫星系的束縛,但是這些星系中有些能夠逃出恒星的束縛,卻逃不出來(lái)自暗物質(zhì)暈的引力束縛。還有部分星系在核活動(dòng)未終止前,氣體殼層已經(jīng)減速到0,停止向外擴(kuò)散。

    圖3展示了核活動(dòng)終止時(shí),星系的尺度R0或暗物質(zhì)暈的維里半徑rvir與氣體殼層到達(dá)的位置rf之間的關(guān)系,圖中只展示了核活動(dòng)終止后氣體外流速度vf不為0的星系。

    圖4核活動(dòng)終止時(shí),氣體殼層外流速度與其所在處逃逸速度的比值。假定L(t)=LEdd(MBH)=const,(a)為只考慮

    4氣體冷卻

    在氣體殼層向外擴(kuò)散的過程中,文[16]指出康普頓冷卻對(duì)氣體殼層的影響可能很重要。冷卻非常迅速就有可能導(dǎo)致氣體停止向外擴(kuò)散,當(dāng)氣體殼層向外擴(kuò)散時(shí)會(huì)產(chǎn)生激波而被加熱,熱電子與光子發(fā)生康普頓散射,光子吸收電子的能量,使得氣體冷卻??灯疹D散射導(dǎo)致氣體冷卻的時(shí)標(biāo)為

    (15)

    其中E是電子的動(dòng)能量;dE/dt是電子能量耗散速率。若僅考慮簡(jiǎn)單的單電子逆康普頓散射,其輻射功率[32]為

    (16)

    (17)

    式中,v(t)和r(t)是活動(dòng)星系核演化到t時(shí)刻氣體殼層的速度和到星系中心的距離;mp是質(zhì)子質(zhì)量;me是電子的質(zhì)量;rkpc是r(t)以千秒差距為單位的距離值;M8=MBH/103M⊙。

    核活動(dòng)演化到t時(shí)刻,若對(duì)應(yīng)的康普頓冷卻時(shí)標(biāo)tc(t)遠(yuǎn)小于核活動(dòng)時(shí)標(biāo),則氣體會(huì)被有效地冷卻,氣體的動(dòng)量和能量將最終以光子的形式釋放,使氣體無(wú)法逃出星系;反之,若tc(t)遠(yuǎn)大于核活動(dòng)時(shí)標(biāo)甚至Hubble年齡,則氣體不能被有效冷卻,將繼續(xù)向外擴(kuò)散,最終當(dāng)核活動(dòng)終止時(shí),若氣體殼層的速度vf大于逃逸速度ve(rf)時(shí),氣體就能夠逃出星系。

    圖5展示了核活動(dòng)終止時(shí),被激波加速氣體外流的冷卻時(shí)標(biāo)tc。由圖5可知,核活動(dòng)終止后部分星系中的氣體能夠逃出星系的束縛,其康普頓冷卻時(shí)標(biāo)tc遠(yuǎn)大于活動(dòng)星系核演化時(shí)標(biāo)τs,康普頓冷卻對(duì)這部分星系的動(dòng)量反饋過程無(wú)影響。

    圖5核活動(dòng)終止時(shí),氣體殼層激波加熱后康普頓冷卻時(shí)標(biāo)。假定L(t)=LEdd(MBH)=const,(a)為只考慮恒星物質(zhì)

    5討論

    核活動(dòng)終止后,樣本星系中部分星系的氣體能夠逃出星系(或暗物質(zhì)暈),部分星系中的氣體在核活動(dòng)尚未終止已經(jīng)停止向外擴(kuò)散,與前者對(duì)比,后者星系恒星總質(zhì)量M*,total(或暗物質(zhì)暈質(zhì)量Mh)與中心黑洞的質(zhì)量MBH比值偏大。因此有可能是黑洞質(zhì)量估計(jì)偏小??紤]了黑洞的真實(shí)質(zhì)量稍大于測(cè)量給出的中心值(例如為測(cè)量值再加上一個(gè)標(biāo)準(zhǔn)偏差)這一情況下的活動(dòng)星系核的動(dòng)量反饋。

    6結(jié)果

    本文引入星系觀測(cè)數(shù)據(jù)重構(gòu)了星系的恒星質(zhì)量密度分布和其寄主暗物質(zhì)暈的密度輪廓,并據(jù)此重新考慮了動(dòng)量反饋過程。研究并討論了影響動(dòng)量反饋有效性的幾個(gè)因素,諸如星系中氣體質(zhì)量與恒星質(zhì)量之比f(wàn)g,外流動(dòng)量輸出與中心引擎輻射光子總動(dòng)量之比f(wàn)p以及星系周圍的暗物質(zhì)暈引力束縛的影響。發(fā)現(xiàn)動(dòng)量驅(qū)動(dòng)反饋機(jī)制在絕大多數(shù)星系演化過程中是有效的,從而修正了文[1]得到的動(dòng)量反饋無(wú)效這一結(jié)果。

    致謝:感謝陸由俊研究員在理論知識(shí)和數(shù)值計(jì)算方面的指導(dǎo),以及在論文寫作中的幫助。

    參考文獻(xiàn):

    [1]Silk J, Nusser A. The massive-black-hole-velocity-dispersion relation and the halo baryon fraction: a case for positive active galactic nucleus feedback[J]. The Astrophysical Journal, 2010, 725: 556-560.

    [2]Magorrian J, Tremaine S, Richstone D, et al. The demography of massive dark objects in galaxy centers[J]. The Astronomical Journal, 1998, 115(6): 2285-2305.

    [3]Richstone D, Ajhar E A, Bender R, et al. Supermassive black holes and the evolution of galaxies[J]. Nature, 1998, 395(6701): A14-A19.

    [4]Ferrarese L, Merritt D. A fundamental relation between supermassive black holes and their host galaxies[J]. The Astrophysical Journal, 2000, 539(1): L9-L12.

    [5]Gebhardt K, Bender R, Bower G, et al. A relationship between nuclear black hole mass and galaxy velocity dispersion[J]. The Astrophysical Journal, 2000, 539(1): L13-L16.

    [6]Tremaine S, Gebhardt K, Bender R, et al. The slope of the black hole mass versus velocity dispersion correlation[J]. The Astrophysical Journal, 2002, 574(2): 740-753.

    [7]H?ring N, Rix H W. On the black hole mass-bulge mass relation[J]. The Astrophysical Journal, 2004, 604(2): L89-L92.

    [8]Gültekin K, Richstone D O, Gebhardt K, el al. The M-σ and M-L relations in galactic bulges and determinations of their intrinsic scatter[J]. The Astrophysical Journal, 2009, 698(1): 198-221.

    [9]Mcconnell N J, Ma C P, Graham J R, et al. The black hole mass in the brightest cluster galaxy NGC 6086[J]. The Astrophysical Journal, 2011, 728(2): 100-120.

    [10]Mcconnell N J, Ma C P. Revisiting the scaling relations of black hole masses and host galaxy properties[J]. The Astrophysical Journal, 2013, 764(2): 184-197.

    [11]Graham A W, Scott N. The M BH-L spheroid relation at high and low masses, the quadratic growth of black holes, and intermediate-mass black hole candidates[J]. The Astrophysical Journal, 2013, 764(2): 151-169.

    [12]Kormendy J, Ho L C. Coevolution (or not) of supermassive black holes and host galaxies[J]. Annual Review of Astronomy and Astrophysics, 2013, 51(2): 511-653.

    [13]Cattaneo A. Quasars and galaxy formation[J]. Monthly Notices of the Royal Astronomical Society, 2001, 324(1): 128-140.

    [14]Tiziana D M, Volker S, Lars H. Energy input from quasars regulates the growth and activity of black holes and their host galaxies[J]. Nature, 2005, 433(7026): 604-607.

    [15]Fabian A C. The obscured growth of massive black holes[J]. Monthly Notices of the Royal Astronomical Society, 1999, 308: L39-L43.

    [16]King A. Black holes, galaxy formation, and the M B H-σ relation[J]. The Astrophysical Journal, 2003, 596(1): L27-L29.

    [17]Murray N, Quataert E, Thompson T A. On the maximum luminosity of galaxies and their central black holes: feedback from momentum-driven wind[J]. The Astrophysical Journal, 2005, 618: 569-585.

    [18]Faber S M, Tremaine S, Ajhar E, et al. The centers of early-type galaxies with HST. IV. central parameter relations[J]. The Astronomical Journal, 1997, 114: 1771-1796.

    [19]Lauer T R, Gebhardt K, Faber S M, et al. The centers of early-type galaxies with Hubble Space Telescope. VI. bimodal central surface brightness profiles[J]. The Astrophysical Journal, 2007, 664(1): 226-256.

    [20]Cappellari M, Scott N, Alatalo K, et al. The Atlas3D project-XV. benchmark for early-type galaxies scaling relations from 260 dynamical models: mass-to-light ratio, dark matter, Fundamental Plane and Mass Plane[J]. Monthly Notices of the Royal Astronomical Society, 2013, 432(3): 1709-1741.

    [21]Cappellari M, Remid R M, Alatalo K, et al. The Atlas3D project-XX. mass-size and mass-σ distributions of early-type galaxies: bulge fraction drives kinematics, mass-to-light ratio, molecular gas fraction and stellar initial mass function[J]. Monthly Notices of the Royal Astronomical Society, 2013, 432(3): 1862-1893.

    [22]Cappellari M, Emsellem E, Krajnovic'D, et al. The Atlas3D project-I. a volume-limited sample of 260 nearby early-type galaxies: science goals and selection criteria[J]. Monthly Notices of the Royal Astronomical Society, 2011, 413(2): 813-836.

    [23]Binney J, Tremaine S. Galactic dynamics[M]. 2nd ed. Princeton: Princeton University Press, 1987.

    [24]Navarro J F, Frenk C S, White S D M. The structure of cold dark matter halos[J]. The Astrophysical Journal, 1996, 462: 563-575.

    [25]Duffy A R, Schaye J, Kay S T, et al. Dark matter halo concentrations in the Wilkinson Microwave Anisotropy Probe year 5 cosmology[J]. Monthly Notices of the Royal Astronomical Society, 2008, 309: L64-L68.

    [26]Behroozi P S, Conroy C, Wechsler R H. A comprehensive analysis of uncertainties affecting the stellar mass-halo mass relation for[J]. The Astrophysical Journal, 2010, 717: 379-403.

    [27]Fouque P, Solanes J M, Sanchis T, et al. Structure, mass and distance of the Virgo cluster from a Tolman-Bondi model[J]. Astronomy & Astrophysics, 2001, 375(3): 770-780.

    [28]Pounds K A, King A R, Page K L, et al. Evidence of a high-velocity ionized outflow in a second narrow-line quasar PG 0844+349[J]. Monthly Notices of the Royal Astronomical Society, 2003, 346(4): 1025-1030.

    [29]Madau P, Haardt F, Dotti M. Super-critical growth of massive black holes from strellar-mass seeds[J]. The Astrophysical Journal, 2014, 784: L38-L41.

    [30]Page D N, Thorne K S. Disk-accretion onto a black hole. I. time-averaged structure of accretion disk[J]. The Astrophysical Journal, 1974, 191: 499-506.

    [31]Yu Q, Tremaine S. Observational constraints on growth of massive black holes[J]. Monthly Notices of the Royal Astronomical Society, 2002, 335(4): 965-976.

    [32]Ishibashi W, Courvoisier J L. The physical origin of the X-ray power spectral density break timescale in accreting black holes[J]. Astronomy & Astrophysics, 2012, 540(4): L2-L5.

    [33]Oser L, Naab T, Ostriker J P, et al. The cosmological size and velocity dispersion evolution of massive early-type galaxies[J]. The Astrophysical Journal, 2012, 744: 63-71.

    CN 53-1189/PISSN 1672-7673

    Research on the Momentum Driven Feedback Mechanism

    Fu Yan, Zhang Xiaoxia

    (National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100012 China, Email: fuyan@bao.ac.cn)

    Key words:AGN feedback; Momentum driven feedback; Outflow; Compton cooling

    Abstract:The observations have shown a tight correlation between the mass of supermassive black holes (MBH) at galactic centers and the stellar velocity dispersion (σ), or bulge stellar mass (Mb). It is believed that such relationships are caused by the feedback from the nucleus activity, but in the microscopic physics view the AGN feedback process proceeding is not clear. Recently, Silk & Nusser (2010) suspect the validity of the AGN feedback process by investigating a number of galaxies with observationally determined black hole mass and the stellar velocity dispersion. In the paper, we revisit the study of Silk & Nusser (2010) by considering more details about the galaxy properties. We combine a sample of 29 galaxies. Adopting the observed data of those galaxies, we calculate and reconsider the process of AGN feedback to test the effectiveness of AGN feedback in the potential of those galaxies, their density profile and the potentials, and dark matter. We find that the momentum driven feedback can be effective in most galaxies.

    基金項(xiàng)目:國(guó)家自然科學(xué)基金 (11373031) 資助.

    收稿日期:2015-04-27;

    修訂日期:2015-05-07

    作者簡(jiǎn)介:付艷,女,碩士. 研究方向:活動(dòng)星系核. Email: fuyan@bao.ac.cn

    中圖分類號(hào):P142.9

    文獻(xiàn)標(biāo)識(shí)碼:A

    文章編號(hào):1672-7673(2016)01-0031-13

    猜你喜歡
    外流
    基于Fluent的賽車翼板外流場(chǎng)設(shè)計(jì)與仿真
    人口外流成因及對(duì)策
    活力(2019年17期)2019-11-26 00:42:20
    人力資源流動(dòng)對(duì)農(nóng)村經(jīng)濟(jì)發(fā)展的影響
    從“外流”到“回流”“鄉(xiāng)村旅游”下的村落需要返璞歸真
    不外流
    對(duì)煤炭企業(yè)專業(yè)技術(shù)人才外流的思考
    商情(2017年44期)2017-12-27 02:07:12
    外流販毒高危預(yù)警模型初探
    基于非結(jié)構(gòu)混合網(wǎng)格人椅系統(tǒng)外流場(chǎng)數(shù)值模擬
    消費(fèi)“外流”現(xiàn)象何時(shí)改變?
    農(nóng)村金融資金流向內(nèi)在誘因及其影響研究
    亚洲精品亚洲一区二区| 欧美日韩福利视频一区二区| 亚洲国产中文字幕在线视频| 国产精品乱码一区二三区的特点| 国产三级中文精品| 淫秽高清视频在线观看| 三级毛片av免费| 国产精品一区二区免费欧美| 一区二区三区激情视频| 精品免费久久久久久久清纯| 嫩草影院入口| 日韩av在线大香蕉| 亚洲中文日韩欧美视频| 脱女人内裤的视频| 色综合婷婷激情| 日韩欧美在线乱码| 看黄色毛片网站| 搡老熟女国产l中国老女人| 中文资源天堂在线| 国产黄片美女视频| 国内精品美女久久久久久| 日本五十路高清| 搡女人真爽免费视频火全软件 | 精品久久久久久,| 久久久成人免费电影| 久久精品影院6| 51午夜福利影视在线观看| 99久久无色码亚洲精品果冻| 99久久成人亚洲精品观看| 老熟妇乱子伦视频在线观看| 一进一出好大好爽视频| tocl精华| 国产成人av激情在线播放| 在线天堂最新版资源| 国产精品久久久人人做人人爽| 久久久久久久精品吃奶| 级片在线观看| 欧美三级亚洲精品| 哪里可以看免费的av片| 在线免费观看不下载黄p国产 | 国内久久婷婷六月综合欲色啪| 757午夜福利合集在线观看| 日本与韩国留学比较| 中文字幕熟女人妻在线| 国产69精品久久久久777片| 午夜福利成人在线免费观看| 丰满的人妻完整版| 日韩免费av在线播放| ponron亚洲| 五月玫瑰六月丁香| 国产精品嫩草影院av在线观看 | 亚洲欧美日韩卡通动漫| 久久久久精品国产欧美久久久| 久久国产精品影院| 免费高清视频大片| av女优亚洲男人天堂| 亚洲在线观看片| 999久久久精品免费观看国产| 久久中文看片网| 亚洲五月婷婷丁香| 亚洲av中文字字幕乱码综合| 亚洲美女黄片视频| 国产精品久久久久久人妻精品电影| 国产精品久久久久久久久免 | 久久精品91蜜桃| 丰满人妻一区二区三区视频av | 一个人看的www免费观看视频| 国产不卡一卡二| 色播亚洲综合网| 成人18禁在线播放| 精品一区二区三区人妻视频| 90打野战视频偷拍视频| 国产亚洲欧美98| 亚洲一区二区三区不卡视频| av欧美777| 嫩草影院入口| 国产成人a区在线观看| 999久久久精品免费观看国产| 波多野结衣高清作品| 亚洲欧美日韩高清专用| 天堂av国产一区二区熟女人妻| 我的老师免费观看完整版| 国产精品三级大全| 黄色片一级片一级黄色片| 中文字幕av在线有码专区| 国产精品久久久久久亚洲av鲁大| 欧美成人性av电影在线观看| 制服人妻中文乱码| 久久久久久人人人人人| 搡老妇女老女人老熟妇| 级片在线观看| 色哟哟哟哟哟哟| 国产一区二区亚洲精品在线观看| 免费在线观看亚洲国产| 手机成人av网站| av欧美777| 亚洲精品久久国产高清桃花| 黑人欧美特级aaaaaa片| 久久欧美精品欧美久久欧美| 色播亚洲综合网| 中文字幕人妻熟人妻熟丝袜美 | 夜夜夜夜夜久久久久| 国产免费男女视频| 香蕉av资源在线| 成年人黄色毛片网站| 窝窝影院91人妻| 波多野结衣高清无吗| 国内毛片毛片毛片毛片毛片| svipshipincom国产片| 亚洲 欧美 日韩 在线 免费| 国产一区二区三区在线臀色熟女| 女人被狂操c到高潮| 亚洲国产精品合色在线| 国产成年人精品一区二区| 在线观看日韩欧美| 色老头精品视频在线观看| 国产免费av片在线观看野外av| 亚洲最大成人手机在线| 一级毛片高清免费大全| 男人舔女人下体高潮全视频| 欧美成人一区二区免费高清观看| 国产av麻豆久久久久久久| 午夜福利在线在线| 亚洲精华国产精华精| 国产色爽女视频免费观看| 成人av一区二区三区在线看| 69人妻影院| 高清在线国产一区| 日韩欧美 国产精品| 欧美性猛交黑人性爽| 国产成人欧美在线观看| 黄色成人免费大全| 麻豆一二三区av精品| av在线蜜桃| 国产精品日韩av在线免费观看| 国产中年淑女户外野战色| 91九色精品人成在线观看| eeuss影院久久| 在线观看午夜福利视频| 亚洲国产欧美网| 亚洲精品成人久久久久久| 最近最新免费中文字幕在线| 在线国产一区二区在线| 精品国产美女av久久久久小说| 久久久色成人| 日韩国内少妇激情av| 免费av毛片视频| 精品一区二区三区视频在线 | 久久久久久久精品吃奶| 国产探花在线观看一区二区| 午夜a级毛片| 美女cb高潮喷水在线观看| 亚洲av中文字字幕乱码综合| 欧美bdsm另类| 免费观看的影片在线观看| 真人做人爱边吃奶动态| 国产高清激情床上av| www.www免费av| 国产av麻豆久久久久久久| 在线视频色国产色| 亚洲国产色片| 岛国在线观看网站| 国产一区二区亚洲精品在线观看| 国产精品 欧美亚洲| 身体一侧抽搐| 一个人观看的视频www高清免费观看| 国产在线精品亚洲第一网站| www日本黄色视频网| www.www免费av| 麻豆一二三区av精品| av国产免费在线观看| 午夜福利视频1000在线观看| 熟妇人妻久久中文字幕3abv| 国产欧美日韩一区二区三| 操出白浆在线播放| 激情在线观看视频在线高清| 日本精品一区二区三区蜜桃| 欧美日韩亚洲国产一区二区在线观看| av在线蜜桃| 亚洲av一区综合| 精品久久久久久久久久久久久| 一区二区三区高清视频在线| 少妇裸体淫交视频免费看高清| 国产在视频线在精品| 999久久久精品免费观看国产| 亚洲精品成人久久久久久| 欧美午夜高清在线| 我要搜黄色片| 给我免费播放毛片高清在线观看| 亚洲在线自拍视频| 欧美日韩黄片免| 成人三级黄色视频| 国内久久婷婷六月综合欲色啪| 白带黄色成豆腐渣| 午夜福利欧美成人| 91在线精品国自产拍蜜月 | 午夜福利在线在线| 欧美日韩乱码在线| 青草久久国产| 国产日本99.免费观看| 亚洲av中文字字幕乱码综合| av在线蜜桃| 亚洲av成人av| 亚洲欧美日韩高清专用| av黄色大香蕉| 亚洲无线观看免费| 久久婷婷人人爽人人干人人爱| 欧美一区二区国产精品久久精品| 老司机午夜十八禁免费视频| 99精品在免费线老司机午夜| 亚洲欧美日韩东京热| 亚洲av成人不卡在线观看播放网| 国产精品一及| 午夜福利在线观看吧| 蜜桃久久精品国产亚洲av| 日韩精品中文字幕看吧| 3wmmmm亚洲av在线观看| 日日干狠狠操夜夜爽| 亚洲国产中文字幕在线视频| 欧美日韩精品网址| 婷婷六月久久综合丁香| 午夜免费激情av| 男女那种视频在线观看| av在线蜜桃| 亚洲内射少妇av| 欧美日韩一级在线毛片| 国产精品电影一区二区三区| 一进一出抽搐gif免费好疼| 国产成人福利小说| 无人区码免费观看不卡| 在线播放无遮挡| 亚洲美女黄片视频| tocl精华| 欧美区成人在线视频| 听说在线观看完整版免费高清| 国产 一区 欧美 日韩| 日韩欧美精品v在线| 久久午夜亚洲精品久久| 成人国产综合亚洲| 欧美一区二区亚洲| 少妇的逼水好多| 色综合亚洲欧美另类图片| 国产亚洲欧美98| 欧美日韩黄片免| 日韩av在线大香蕉| 内射极品少妇av片p| 国产精品一区二区三区四区久久| 99久久综合精品五月天人人| 久久久久久九九精品二区国产| а√天堂www在线а√下载| 窝窝影院91人妻| 久久99热这里只有精品18| 夜夜夜夜夜久久久久| 不卡一级毛片| 桃红色精品国产亚洲av| h日本视频在线播放| 亚洲狠狠婷婷综合久久图片| av专区在线播放| 香蕉av资源在线| 国产精品久久久久久久久免 | 亚洲欧美日韩无卡精品| 老司机午夜十八禁免费视频| 12—13女人毛片做爰片一| 亚洲精品影视一区二区三区av| 欧美日韩国产亚洲二区| 嫩草影院精品99| 国产成人av激情在线播放| 热99在线观看视频| 亚洲av电影不卡..在线观看| 欧美av亚洲av综合av国产av| 97超级碰碰碰精品色视频在线观看| 精品国产超薄肉色丝袜足j| 69av精品久久久久久| 18+在线观看网站| 美女高潮喷水抽搐中文字幕| 国产在视频线在精品| 97人妻精品一区二区三区麻豆| 日韩有码中文字幕| 1000部很黄的大片| 真人一进一出gif抽搐免费| 99久久99久久久精品蜜桃| 亚洲国产欧洲综合997久久,| 制服人妻中文乱码| 高清日韩中文字幕在线| 国产久久久一区二区三区| 人人妻人人澡欧美一区二区| 成人性生交大片免费视频hd| 精品国产超薄肉色丝袜足j| 天天躁日日操中文字幕| 老司机在亚洲福利影院| 欧美激情在线99| 在线播放无遮挡| 亚洲内射少妇av| 国产色爽女视频免费观看| 久久亚洲精品不卡| 精品久久久久久成人av| 国产不卡一卡二| 最新在线观看一区二区三区| 久久久久亚洲av毛片大全| 精品一区二区三区视频在线 | 精品无人区乱码1区二区| 国产一区二区三区视频了| 中国美女看黄片| 日韩精品中文字幕看吧| 久久久久久久精品吃奶| 色老头精品视频在线观看| 国内毛片毛片毛片毛片毛片| 一本综合久久免费| 亚洲无线在线观看| 中文字幕高清在线视频| 波多野结衣高清作品| 国产精品 欧美亚洲| 国产成人系列免费观看| 欧美丝袜亚洲另类 | 欧美最新免费一区二区三区 | 久久久国产成人免费| 国产精品久久久久久久电影 | 伊人久久大香线蕉亚洲五| 亚洲欧美日韩高清专用| av在线蜜桃| 亚洲五月天丁香| 波多野结衣巨乳人妻| 久久99热这里只有精品18| 成人性生交大片免费视频hd| 国产探花极品一区二区| 久久久久久大精品| 成人一区二区视频在线观看| 久久中文看片网| 国产高清videossex| avwww免费| 国产成人欧美在线观看| 最好的美女福利视频网| 又爽又黄无遮挡网站| 国产av在哪里看| 黄色片一级片一级黄色片| 成人av一区二区三区在线看| 可以在线观看的亚洲视频| 精品不卡国产一区二区三区| 色综合亚洲欧美另类图片| 久久欧美精品欧美久久欧美| 黄色视频,在线免费观看| 欧美极品一区二区三区四区| 午夜福利在线观看免费完整高清在 | 最好的美女福利视频网| 禁无遮挡网站| 黄色成人免费大全| 成人无遮挡网站| 国产精品野战在线观看| 国产精品久久久人人做人人爽| 内地一区二区视频在线| 他把我摸到了高潮在线观看| 国产色爽女视频免费观看| 宅男免费午夜| 怎么达到女性高潮| 国产一区二区激情短视频| av专区在线播放| 亚洲精品粉嫩美女一区| 国产激情欧美一区二区| 90打野战视频偷拍视频| 久久久久免费精品人妻一区二区| 午夜亚洲福利在线播放| 少妇高潮的动态图| 亚洲精品影视一区二区三区av| 精品福利观看| 欧美成人a在线观看| 亚洲av二区三区四区| 国产精品 国内视频| 免费无遮挡裸体视频| 男女视频在线观看网站免费| 久久香蕉精品热| 桃红色精品国产亚洲av| 精品一区二区三区人妻视频| 18禁黄网站禁片午夜丰满| 中文字幕av成人在线电影| 亚洲人与动物交配视频| 久久精品人妻少妇| 国产视频内射| 3wmmmm亚洲av在线观看| 黄色丝袜av网址大全| 国产欧美日韩一区二区三| 欧美大码av| 久久久久久久久久黄片| 欧美最黄视频在线播放免费| or卡值多少钱| 国产探花极品一区二区| 嫩草影院入口| 亚洲 欧美 日韩 在线 免费| 成人性生交大片免费视频hd| 日本五十路高清| 成人鲁丝片一二三区免费| 亚洲国产精品sss在线观看| 99久久综合精品五月天人人| 国产蜜桃级精品一区二区三区| 99久久成人亚洲精品观看| 国产91精品成人一区二区三区| 一本精品99久久精品77| 波多野结衣高清作品| 窝窝影院91人妻| 久久久国产精品麻豆| 亚洲午夜理论影院| 欧美成人性av电影在线观看| 三级男女做爰猛烈吃奶摸视频| 高清毛片免费观看视频网站| 国产99白浆流出| 在线观看日韩欧美| 亚洲乱码一区二区免费版| 性欧美人与动物交配| 五月伊人婷婷丁香| 欧美国产日韩亚洲一区| 国产高清视频在线观看网站| 精品久久久久久久人妻蜜臀av| a级毛片a级免费在线| 久久久久免费精品人妻一区二区| 69av精品久久久久久| av福利片在线观看| 欧美成狂野欧美在线观看| 国产高潮美女av| x7x7x7水蜜桃| 五月伊人婷婷丁香| 制服人妻中文乱码| 成年女人永久免费观看视频| 搡老熟女国产l中国老女人| 精品免费久久久久久久清纯| 三级男女做爰猛烈吃奶摸视频| 美女cb高潮喷水在线观看| 亚洲成人久久性| 国内精品久久久久久久电影| 亚洲精品日韩av片在线观看 | 欧美绝顶高潮抽搐喷水| 天天一区二区日本电影三级| 高清日韩中文字幕在线| 欧美一级毛片孕妇| 亚洲 国产 在线| 久久性视频一级片| 国产精品一区二区三区四区免费观看 | av福利片在线观看| 不卡一级毛片| 最近最新免费中文字幕在线| 又黄又粗又硬又大视频| 亚洲无线在线观看| 怎么达到女性高潮| 免费电影在线观看免费观看| 久久草成人影院| 97碰自拍视频| 中国美女看黄片| 色老头精品视频在线观看| 精品久久久久久久末码| 人妻夜夜爽99麻豆av| 非洲黑人性xxxx精品又粗又长| 午夜福利18| 国产久久久一区二区三区| 99精品欧美一区二区三区四区| 亚洲人成网站高清观看| 舔av片在线| www日本在线高清视频| 日韩欧美精品免费久久 | 国产私拍福利视频在线观看| 99热这里只有是精品50| 精品久久久久久久毛片微露脸| 99国产精品一区二区三区| 久久精品91无色码中文字幕| 性色av乱码一区二区三区2| a级毛片a级免费在线| 欧美高清成人免费视频www| 狠狠狠狠99中文字幕| 亚洲专区国产一区二区| 午夜精品久久久久久毛片777| 午夜福利在线在线| 法律面前人人平等表现在哪些方面| 97超视频在线观看视频| 亚洲成av人片在线播放无| 久久久久久久精品吃奶| 麻豆久久精品国产亚洲av| 亚洲成a人片在线一区二区| 在线观看一区二区三区| 欧美一级a爱片免费观看看| 嫩草影院精品99| 成年版毛片免费区| 麻豆成人午夜福利视频| 亚洲精品日韩av片在线观看 | 精品一区二区三区人妻视频| 最新在线观看一区二区三区| 亚洲激情在线av| av专区在线播放| 搡老妇女老女人老熟妇| 精品国产亚洲在线| 老司机午夜十八禁免费视频| 99久久精品热视频| 国产激情偷乱视频一区二区| 欧美一区二区亚洲| 熟女电影av网| 亚洲成人久久爱视频| 亚洲精品成人久久久久久| 黑人欧美特级aaaaaa片| 叶爱在线成人免费视频播放| 国产精品1区2区在线观看.| 一级作爱视频免费观看| 国产一区在线观看成人免费| 国产成人av激情在线播放| av在线天堂中文字幕| 亚洲精品一卡2卡三卡4卡5卡| 国产淫片久久久久久久久 | 亚洲,欧美精品.| 色综合婷婷激情| 天天躁日日操中文字幕| 亚洲专区中文字幕在线| 内地一区二区视频在线| 国产精品三级大全| 亚洲av五月六月丁香网| 操出白浆在线播放| 精品国产三级普通话版| a级一级毛片免费在线观看| 亚洲欧美一区二区三区黑人| 国产午夜福利久久久久久| 最新在线观看一区二区三区| 97超视频在线观看视频| 看片在线看免费视频| 有码 亚洲区| 亚洲欧美精品综合久久99| 午夜久久久久精精品| 久久这里只有精品中国| 亚洲精品456在线播放app | 美女被艹到高潮喷水动态| 亚洲欧美一区二区三区黑人| 国产免费一级a男人的天堂| 少妇的丰满在线观看| 亚洲精品色激情综合| 又黄又粗又硬又大视频| 国产精品嫩草影院av在线观看 | 亚洲av五月六月丁香网| 国产精品美女特级片免费视频播放器| 在线观看一区二区三区| 天堂av国产一区二区熟女人妻| 99久久无色码亚洲精品果冻| 成人国产一区最新在线观看| 日本在线视频免费播放| 精品久久久久久久毛片微露脸| 久久国产精品影院| av片东京热男人的天堂| 久久久久久久久久黄片| 伊人久久精品亚洲午夜| 最近在线观看免费完整版| 欧美在线黄色| 99久久精品一区二区三区| 夜夜爽天天搞| 日韩欧美国产一区二区入口| 国产高清videossex| 极品教师在线免费播放| 在线天堂最新版资源| 成人鲁丝片一二三区免费| 男女床上黄色一级片免费看| 99国产精品一区二区三区| 久9热在线精品视频| 午夜福利欧美成人| 看片在线看免费视频| 国产精品久久视频播放| 99国产极品粉嫩在线观看| 最新美女视频免费是黄的| 少妇裸体淫交视频免费看高清| 亚洲一区二区三区色噜噜| 亚洲av不卡在线观看| 久久久久久久午夜电影| 女同久久另类99精品国产91| 欧美3d第一页| АⅤ资源中文在线天堂| h日本视频在线播放| x7x7x7水蜜桃| 美女高潮喷水抽搐中文字幕| 又黄又爽又免费观看的视频| 狂野欧美白嫩少妇大欣赏| 国产欧美日韩一区二区三| 国产伦在线观看视频一区| 亚洲av免费高清在线观看| 中文字幕人成人乱码亚洲影| 看片在线看免费视频| 午夜老司机福利剧场| 亚洲自拍偷在线| 中文字幕久久专区| 天天躁日日操中文字幕| 免费看a级黄色片| 夜夜躁狠狠躁天天躁| 亚洲黑人精品在线| 日本五十路高清| 欧美一级a爱片免费观看看| 亚洲aⅴ乱码一区二区在线播放| 嫩草影视91久久| 国产又黄又爽又无遮挡在线| 亚洲七黄色美女视频| 久久精品人妻少妇| 国产 一区 欧美 日韩| 亚洲欧美一区二区三区黑人| 伊人久久大香线蕉亚洲五| 观看免费一级毛片| 亚洲av熟女| 日韩欧美国产一区二区入口| 久久久久久久亚洲中文字幕 | www.熟女人妻精品国产| 欧美日韩乱码在线| 午夜福利高清视频| 国产aⅴ精品一区二区三区波| 制服丝袜大香蕉在线| 成人av一区二区三区在线看| 久久久久精品国产欧美久久久| 一个人免费在线观看电影| 日本免费一区二区三区高清不卡| 亚洲天堂国产精品一区在线| 岛国在线免费视频观看| 久久久久亚洲av毛片大全| 午夜精品一区二区三区免费看| 国产精品 国内视频| 日韩成人在线观看一区二区三区| 免费一级毛片在线播放高清视频| 午夜免费激情av| 一区福利在线观看| 国产熟女xx|