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

    Exploring the latest Pantheon SN Ia dataset by using three kinds of statistics techniques

    2021-04-28 02:27:04ShuangWangandXiaolinLuo
    Communications in Theoretical Physics 2021年4期

    Shuang Wang and Xiaolin Luo

    School of Physics and Astronomy,Sun Yat-Sen University,Guangzhou 510297,China

    Abstract In this work,we explore the cosmological consequences of the latest Type Ia supernova(SN Ia)dataset,Pantheon,by adopting the wCDM model.The Pantheon dataset currently contains the largest number of SN Ia samples,which contains 1048 supernovae on the redshift range 0 <z <2.3.Here we take into account three kinds of SN Ia statistics techniques,including:(1)magnitude statistics(MS),which is the traditional SN Ia statistics technique;(2)flux statistics(FS),which is based on the flux-averaging (FA) method; and (3) improved flux statistics (IFS),which combines the advantages of MS and FS.It should be mentioned that the IFS technique needs to scan the(zcut,Δz)parameters plane,where zcut and Δz are redshift cut-off and redshift interval of FA,respectively.The results are as follows.(1) Using the SN dataset only,the best FA recipe for IFS is(zcut,Δz)=(0.1,0.08);(2)comparing to the old SN dataset,JLA,adopting the Pantheon dataset can reduce the 2σ error bars of equation of state w by 38%,47%and 53%for MS,FS and IFS,respectively; (3) FS gives closer results to other observations,such as Baryon acoustic oscillations and cosmic microwave background;(4)compared with FS and IFS,MS more favors a Universe that will end in a ‘big rip’.

    Keywords: Dark Energy,Type Ia supernova,Cosmological Observations

    1.Introduction

    The Type Ia supernova (SN Ia),which is a sub-category of cataclysmic variable stars [1],has played a crucial role in discovering the cosmic acceleration [2,3].So far,SN Ia is still one of the most useful and powerful tools to investigate the nature of dark energy (DE) [4–12].

    In the recent 10 years,many high-quality supernova dataset have been released,such as ‘Union’ [13],‘Constitution’ [14],‘SDSS’[15],‘Union2’[16],‘SNLS3’[17],‘Union2.1’[18]and‘JLA’ [19].In 2018,the latest ‘Pantheon’ dataset [20],which contains 1048 SN Ia at the redshift range 0 <z <2.3,was released.Using the Pantheon sample,Scolnic et al had given the observational constraints on the wCDM and CPL models [20].

    On the other hand,along with the rapid growth of the number of SN Ia discovered,the studies on the systematic uncertainties of SN Ia have drawn more and more attention.It has been proved that [21],the classic SN Ia statistics method(hereafter we will call it ‘magnitude statistic’ (MS)) suffers from various systematic uncertainties,such as the calibration errors [17],the host–galaxy extinction [22–24],the gravitational lensing[25,26],different light-curve fitters[27]and the redshift evolution of SN color-luminosity parameter β[28–34].Therefore,the control of the systematic uncertainties of SN Ia have become one of the biggest challenges in SN cosmology.

    In order to reduce the systematic uncertainties of SN Ia,some interesting statistics techniques of SN Ia are proposed in the literature.For examples,in 2000,Wang proposed a new analysis technique,called flux-averaging (FA),to reduce the systematic errors caused by the weak lensing effect of SN Ia[35].The FA technique focus on the observed flux data of SN Ia,and then average these flux data at some redshift bins with same width.Hereafter,we will call this statistics method of SN Ia as ‘flux statistic’ (FS).The FS can reduce several systematic uncertainties of SN Ia [36–38],but it will lead to larger error bars of model parameters.In 2013,One of the present authors and Wang[28]proposed an improved version of flux-averaging.This new statistics method combines the advantages of MS and FS,and thus can reduce the systematic uncertainties and the error bars of model parameters at the same time[39–41].Hereafter,we will call this latest statistics method of SN Ia as ‘improved flux statistics’ (IFS).

    In this work,we will explore the Latest Pantheon SN Ia dataset by using all the three statistic techniques of SN Ia (i.e.MS,FS and IFS).It should be emphasized that,in the previous studies about the Pantheon samples[42–45],only MS was taken into account.On the other hand,in our previous studies[39–41],The cosmology fits are always performed by combining SN Ia samples with other observations,such as Cosmic microwave background(CMB) and Baryon acoustic oscillations (BAO).In this work,we will mainly focus on the cosmological constraints given by the Pantheon SN Ia Dataset alone.

    The paper is organized as follows.In section.2,we will introduce the methodology used in this work.In particular,we will show how to calculate the χ2function of SN Ia data,for the case of adopting MS,FS,IFS,respectively.In section.3,we will show the results of our studies.We will discuss the differences between cosmological consequences given by Pantheon dataset and by previous SN samples,the differences among the cosmological constraints given by MS,FS,IFS,as well as the ultimate fate of the Universe.Finally,the discussion and conclusion are shown in section.4

    2.Methodology

    In this section,we introduce how to calculate the χ2function of SN Ia data,for the case of adopting MS,FS,IFS,respectively.

    2.1.Magnitude statistics

    As shown in[20],adopting MS,the χ2function of SN Ia data can be expressed as

    Here the Δμ ≡μobs?μth,where μobsis the observational distance modulus of SN,given by [20]:

    where mBis the observed peak magnitude in the rest frame of the B band,M is the absolute B-band magnitude of a fiducial SN Ia,α is the coefficient of the relation between luminosity and stretch,X1describes the time stretching of the light curve,β is the coefficient of the relation between luminosity and color,and C describes the supernova color at maximum brightness.Furthermore,ΔMis a distance correction based on the host-galaxy mass of the SN Ia and ΔBis a distance correction based on predicted biases from simulation.

    The theoretically distance modulus of SN Ia μthcan be expressed as

    Here zcmbis the CMB rest frame redshift,zhelis the heliocentric redshift and dLis the luminosity distance of SN Ia,given by:

    r(z) is given by:

    where c is the speed of light,H0is the current value of the Hubble parameter H(z) and E(z)≡H(z)/H0.

    For simplicity,in this work we only consider the wCDM model (i.e.DE equation of state (EOF) w is a constant parameter) in a flat Universe.Based on the Friedmann equation,we can get:

    where Ωr,Ωmand Ωderepresent the current fractional densities of radiation,matter and dark energy,respectively.The radiation density parameter Ωris given by [46]:

    wherezeq=2.5 ×Tcmb=2.7255 K and h is the reduced Hubble constant.In the case of only adopting SN samples,we set the radiation density parameter Ωr=0.

    In addition,Cov is the total covariance matrix,which is given by:

    where the statistical matrix Dstatonly has the diagonal components,it includes the distance error of each SN Ia as follows:

    whereis the photometric error of the SN Ia distance,is the distance uncertainty from the mass step correction,is the uncertainty from the peculiar velocity uncertainty and redshift measurement uncertainty in quadrature,is the uncertainty from stochastic gravitational lensing,is the intrinsic scatter,andis the uncertainty from the distance bias correction.Furthermore,Csysis the systematic covariance for each SN Ia.One can find the more details about the uncertainty matrix Cov in [17].

    2.2.Flux statistics

    FA divides the whole redshift region into some bins with the same width.The segment points of various bins are zi=Δz·i,where Δz is the width of each bin and i=1,2,3,…,n.

    As shown in[35],adopting FS,the χ2function of SN Ia data can be expressed as:

    where

    The observational flux-averaged distance modulus is calculated by:

    where the ‘a(chǎn)bsolute luminosities’{ L(zcmb)}are

    with the ‘fluxs’ distance modulus

    On the other hand,the theoretical prediction is given by:

    For more details about the FA technique,see [35].

    2.3.Improved flux statistics

    IFS introduces a new parameter,i.e.redshift cut-off zcut.For the case of z <zcut,the χ2function is calculated by using MS;for the case of z ≥zcut,the χ2function is calculated by using FS.It means that:

    Comparing to MS,IFS introduces two new parameters,i.e.redshift cut-off zcutand the width of redshift bin Δz.Here we require that zcut=0.1·i,i=0,1,2,…,8;while Δz=0.01·j,j=4,5,6,…,11.

    Based on the JLA samples,[39]scanned the whole (zcut,Δz)plane and found that(zcut,Δz)=(0.6,0.06)will give the tightest DE constrains.But this result was obtained by using combined observational data,which included CMB,BAO and SN Ia.In this work,we mainly focus on the SN Ia data.Therefore,using the Pantheon sample alone,we will scan the parameter space of (zcut,Δz).

    It should be mentioned that,compared with the conventional approach,the IFS approach is very helpful to reduce the systematic uncertainties of SN Ia.For example,as shown in [28],adopting the IFS approach can significantly reduce the evolution of SN Ia color-luminosity parameter β.

    2.4.Other observational data

    In addition to the SN Ia samples,some other cosmological observations,such as CMB [47,48,46] and BAO [49–51],also play important roles in exploring the nature of DE.Therefore,for comparison,we also take CMB and BAO data into account.

    In this work,for CMB,we use the distance priors data extracted from Planck 2015[52].For BAO,we adopt the data from BOSS DR12 [53],which provides six data points of H(z) and DA(z) at z=0.38,0.51 and 0.61.For more details about calculating χ2of CMB and BAO,see [40].

    3.Cosmology fit results

    The cosmology fits of this work are performed by using COSMOMC package[54].Moreover,to access the ability of constraining DE for various SN Ia statistics techniques,we also take into account the quantity figure of merit (FoM)[55,56],which is the inverse of the area enclose by the 2σ confidence level (CL) contour of (w,Ωm),for the wCDM model.Therefore:

    where Cov(f1,f2,f3,…) is the covariance matrix of the chosen set of DE parameters.

    In this section,first of all,we discuss the best recipe for IFS.Then,we compare the cosmological consequence of the Pantheon dataset with the results of the JLA dataset.Finally,we compare the results of MS,FS and IFS,respectively.

    3.1.Searching the best FA recipe for IFS

    In this subsection,we scan the(zcut,Δz)plane to find the best FA recipe for IFS.As mentioned above,we require that zcut=0.1·i,i=0,1,2,…,8;while Δz=0.01·j,j=4,5,6,…,11.For each set of (zcut,Δz),we perform Markov chain Monte Carlo analysis by using the wCDM model.Then,we compute the corresponding values of the FoM,which are given by equation (19).

    A 3D graph of the values of FoM,given by different sets of (zcut,Δz),is shown in figure 1.It is clear that different values of (zcut,Δz) will give different FoMs.Based on this figure,we find that the best FA recipe of IFS is (zcut,Δz)=(0.1,0.08)with the FoM=278.24(denoted by a black dot).Hereafter,we use this recipe for all IFS techniques.

    Figure 1.3D graph of the values of FoM,which are given by different sets of(zcut,Δz),for the wCDM model.The black dot represents the best FA recipe (zcut,Δz)=(0.1,0.08) for IFS,which gives the FoM=278.34.

    Figure 2.FoM values given by different zcut,for wCDM.The solid red line represents the results constrained only by SN samples.The blue dashed line represents results constrained by the combined SN+BAO+CMB data.

    It should be mentioned that the best fitting results of(zcut,Δz) are quite different from the results of [39],which give(zcut,Δz)=(0.6,0.06).The reason is that the results of [39]are obtained by using the combined SN+BAO+CMB data,while the results of the current paper are obtained by using the SN data only.As a cross check,we also scan the parameter space of (zcut,Δz) by using the combined Pantheon+BAO+CMB data.For this case,the obtained results are very similar to the results of [39].This shows the stability of our method.

    From the above figure,one can see that varying zcutwill produce larger influence on the value of FoM than changing Δz.In figure 2,we give the results of FoM given by different zcut.The solid red line denotes the results given by using the SN data alone,while the dashed blue line represents the results given by using the combined SN+CMB+BAO data.For the case of using the SN data alone,the values of FoM rapidly decrease at the region zcut>0.1.For the case of using the combined SN+CMB+BAO data,zcut=0.2 will yield the maximal value of FoM.In other words,using the SN data alone will give a smaller zcut.

    3.2.The differences between the cosmological consequences of Pantheon and JLA

    In this subsection,we compare the differences between the cosmological consequences of Pantheon and JLA.For complete analysis,all the three statistics techniques,including MS,FS and IFS,are taken into account.

    In figure 3,we present the 1D marginalized probability distributions of w,which are produced by the Pantheon and JLA datasets,respectively.One can find that,for the case of using MS technique,the results of w given by the Pantheon and JLA data are quite different.for the case of using FS and IFS,the Pantheon and JLA data will give similar results for w.

    More details are shown in table 1.For 1σ CL,comparing to the case of using JLA data,using Pantheon data will decrease the the error bars of EoS w by 45%,43% and 56%,for MS,FS and IFS,respectively.For 2σ CL,compared with the case of using JLA data,using Pantheon data will decrease the the error bars of EoS w by 38%,47% and 53%,for MS,FS and IFS,respectively.In addition,using Pantheon data can also increase the values of FoM by 373%,127% and 153%,for MS,FS and IFS,respectively.These results show that compared with the JLA data,the Pantheon data can provide the much tighter DE constraints.

    3.3.The differences between the cosmological consequences of MS,FS and IFS

    In this part,we compare the differences among the cosmological consequences of MS,FS and IFS.

    Figure 3.1D marginalized probability distributions of w for three SN Ia statistics techniques.The solid black lines denote the result from Pantheon and the dashed red lines represent the result from JLA.For all the statistics techniques,the new dataset,Pantheon,gives tighter constraint on w than that constrained by the JLA dataset.

    In figure 4,we present the 1D marginalized probability distributions of w.The solid black,dashed red and dotted blue lines denote the results given by MS,FS and IFS,respectively.As a comparison,we also constrain the wCDM model by adopting the combined CMB+BAO data,which is represented by the dash-dotted cyan line.From this figure,one can see that the IFS can give the tightest constraint among the three SN Ia statistics techniques.In addition,comparing with other SN Ia statistics techniques,FS yields a more similar marginalized probability distribution of w to that given by the CMB+BAO dataset.

    In figure 5,we plot the 2σ error bars of w for three SN Ia statistics techniques.The solid black,dashed red and dotted blue lines represent the results of MS,FS and IFS,respectively.One can find that using IFS will yield the tightest constraint on w.Moreover,using FS will give a smallest lower limit of EoS w,which is less than ?1.As will be discussed in the next subsection,this will lead to a ‘cosmic doomsday’.

    4.Conclusions and discussions

    In this work,we investigate the cosmological consequences of the latest SN Ia dataset,Pantheon,by using the wCDM model.As a comparison,we also consider the JLA dataset.Moreover,three kinds of SN Ia statistics techniques,including MS,FS and IFS,are taken into account.In addition,we mainly focus on the cosmological consequences given by using SN data alone.

    First,based on the wCDM model,we scan the parameter space of(zcut,Δz)to determine the best recipe for IFS.Then,we compare the difference between the fitting results given by the Pantheon and JLA data.Finally,by using SN Ia samples alone,we compare the differences among the cosmological consequences given by MS,FS and IFS.

    Our main results are as follows:

    ? For IFS,we find the best FA recipe(zcut,Δz)=0.1,0.08),which can give the largest FoM=278.34 (see figure 1).Comparing to the case of adopting the combined SN+CMB+BAO dataset,using SN data alone will yield a smaller value of zcut(see figure 2).

    ? Comparing to the JLA dataset,the Pantheon dataset can give tighter DE constraints(see figure 3).Specifically,the Pantheon data can decrease the 2σ error bars of w by 38%47% and 53%,for MS,FS and IFS,respectively.In addition,using Pantheon data can also increase the values of FoM by 373%,127%and 153%,for MS,FS and IFS,respectively (see table 1).

    ? FS gives closer results to other observations,such as BAO and CMB (see figure 4).In addition,among the three SN Ia statistics techniques,IFS yields the tightest constraint on w (see figure 5).

    In this paper,we only discuss a specific DE model,i.e.the wCDM,which has a constant EoS w.It is interesting to consider the cases of adopting various dynamics DE models,such as quintessence [59],Chaplygin gas [60],holographic DE [61],agegraphic DE [62],Yang–Mills condensate [63],Chevalliear–Polarski–Linder parameterization [64,65] and binned parameterization [66–68].

    In addition,in this study we only compare the cosmological consequences of the SN observation with the CMBand BAO observations.It would be interesting to compare the results of the SN observation with some other cosmological observations,such as weak gravitational lensing [69],abundance of galaxy clusters [70],the Alcock–Paczynski effect[71],direct H0measurement [72] and the cosmic age test[73,74].This will be done in future works.

    Table 1.The FoM and CL of w for the JLA and Pantheon datasets.

    Figure 4.1D marginalized probability distributions of w produced by the Pantheon dataset.The solid black,dashed red and dotted blue lines denote the results given by MS,FS and IFS,respectively.As a comparison,we also constrain the wCDM model by adopting the combined CMB+BAO data,which is represented by the dash-dotted cyan line.

    Figure 5.2σ error bars of EoS w given by three SN Ia statistics techniques.The solid black,dashed red and dotted blue lines denote the results given by MS,FS and IFS,respectively.It’s clear that the IFS yields the tightest constraint on w.

    Acknowledgments

    SW is supported by the National Natural Science Foundation of China under Grant No.11405024 and the Fundamental Research Funds for the Central Universities under Grant No.16lgpy50.

    ORCID iDs

    av在线观看视频网站免费| 久久午夜亚洲精品久久| 婷婷色av中文字幕| 日韩在线高清观看一区二区三区| 91午夜精品亚洲一区二区三区| 91麻豆精品激情在线观看国产| 国产黄色小视频在线观看| 禁无遮挡网站| 最近最新中文字幕大全电影3| 天天一区二区日本电影三级| www日本黄色视频网| 99九九线精品视频在线观看视频| 色视频www国产| 国产在线男女| 亚洲精品自拍成人| 国产一区二区三区在线臀色熟女| 久99久视频精品免费| 国内久久婷婷六月综合欲色啪| 日本熟妇午夜| 26uuu在线亚洲综合色| 久久综合国产亚洲精品| 在线免费观看不下载黄p国产| 啦啦啦观看免费观看视频高清| 欧美不卡视频在线免费观看| 国产精品av视频在线免费观看| 日本爱情动作片www.在线观看| 国产日本99.免费观看| 国产av麻豆久久久久久久| 成年女人看的毛片在线观看| 国产精品一区二区三区四区久久| 亚洲国产精品成人久久小说 | 亚洲精品日韩av片在线观看| a级毛片a级免费在线| 日本一二三区视频观看| 日本三级黄在线观看| 观看免费一级毛片| 在线a可以看的网站| 亚洲人成网站在线播放欧美日韩| 在线免费观看的www视频| 美女内射精品一级片tv| 国产又黄又爽又无遮挡在线| 精品一区二区三区人妻视频| 精品久久国产蜜桃| 99热全是精品| 亚洲,欧美,日韩| 国产av在哪里看| 国产精品女同一区二区软件| 性色avwww在线观看| 成人二区视频| 麻豆成人av视频| 蜜臀久久99精品久久宅男| 久久精品影院6| 国内精品宾馆在线| 波多野结衣巨乳人妻| 久久久久久九九精品二区国产| 波多野结衣高清作品| 给我免费播放毛片高清在线观看| av专区在线播放| 99久久久亚洲精品蜜臀av| 亚洲av免费在线观看| 一本久久精品| 国产精品野战在线观看| 国产免费男女视频| 国产精品国产高清国产av| 国产69精品久久久久777片| 国产69精品久久久久777片| 国内久久婷婷六月综合欲色啪| 成人一区二区视频在线观看| 少妇熟女aⅴ在线视频| 成人高潮视频无遮挡免费网站| 欧美一区二区精品小视频在线| 精品久久久久久久久亚洲| 亚洲国产欧洲综合997久久,| 成人漫画全彩无遮挡| 一本一本综合久久| 在线播放国产精品三级| 国产激情偷乱视频一区二区| 伊人久久精品亚洲午夜| www.av在线官网国产| 精华霜和精华液先用哪个| 日本爱情动作片www.在线观看| 一级av片app| 亚洲,欧美,日韩| 成人三级黄色视频| 男女视频在线观看网站免费| 亚洲不卡免费看| 国产精品一区二区三区四区免费观看| avwww免费| 狠狠狠狠99中文字幕| 美女高潮的动态| 久久久久久久久久成人| 热99在线观看视频| 99久久精品国产国产毛片| 黄色欧美视频在线观看| 天天躁夜夜躁狠狠久久av| 国产精品三级大全| 少妇熟女欧美另类| 两性午夜刺激爽爽歪歪视频在线观看| 一区福利在线观看| 日本一二三区视频观看| 欧美最黄视频在线播放免费| 国产三级中文精品| 一个人免费在线观看电影| 啦啦啦啦在线视频资源| 伦理电影大哥的女人| 色噜噜av男人的天堂激情| 我要搜黄色片| 国产一级毛片在线| 波野结衣二区三区在线| www.色视频.com| 国产成人影院久久av| 午夜久久久久精精品| 小说图片视频综合网站| 又粗又硬又长又爽又黄的视频 | 国产成人精品一,二区 | 欧美极品一区二区三区四区| 免费av不卡在线播放| 九九久久精品国产亚洲av麻豆| 一区二区三区免费毛片| 小蜜桃在线观看免费完整版高清| 在线观看一区二区三区| 99riav亚洲国产免费| 国产高清三级在线| av天堂在线播放| 观看美女的网站| 成人亚洲精品av一区二区| 国语自产精品视频在线第100页| 色综合亚洲欧美另类图片| 亚洲欧美清纯卡通| 亚洲,欧美,日韩| 久久精品夜夜夜夜夜久久蜜豆| 少妇猛男粗大的猛烈进出视频 | 国产一区二区亚洲精品在线观看| 99久国产av精品| 国产白丝娇喘喷水9色精品| 免费av毛片视频| 国产亚洲av嫩草精品影院| 国产欧美日韩精品一区二区| 18禁黄网站禁片免费观看直播| 成人午夜高清在线视频| 精品人妻一区二区三区麻豆| 91午夜精品亚洲一区二区三区| 在线观看美女被高潮喷水网站| 99热这里只有精品一区| 日韩欧美 国产精品| 狂野欧美激情性xxxx在线观看| 啦啦啦啦在线视频资源| 亚洲中文字幕一区二区三区有码在线看| 国产精品日韩av在线免费观看| 桃色一区二区三区在线观看| 三级经典国产精品| 亚洲一级一片aⅴ在线观看| 国产淫片久久久久久久久| 国产午夜精品一二区理论片| 久久久久久伊人网av| 精品欧美国产一区二区三| 赤兔流量卡办理| 深夜a级毛片| 亚洲av熟女| 又爽又黄a免费视频| 久久热精品热| 亚洲aⅴ乱码一区二区在线播放| 国产久久久一区二区三区| 精品熟女少妇av免费看| 亚洲无线在线观看| av黄色大香蕉| 国产一区二区三区av在线 | 男人舔奶头视频| 国产精品女同一区二区软件| 真实男女啪啪啪动态图| 国产真实伦视频高清在线观看| 国产高清有码在线观看视频| 午夜精品一区二区三区免费看| ponron亚洲| 欧美日韩乱码在线| 国产色婷婷99| 哪个播放器可以免费观看大片| 久久午夜福利片| 国产精品精品国产色婷婷| 干丝袜人妻中文字幕| 小说图片视频综合网站| 青青草视频在线视频观看| 国产精品久久电影中文字幕| av.在线天堂| 一级毛片电影观看 | 99热这里只有是精品50| 日本色播在线视频| 久久午夜亚洲精品久久| 看十八女毛片水多多多| 亚洲熟妇中文字幕五十中出| 精品99又大又爽又粗少妇毛片| 波野结衣二区三区在线| 免费大片18禁| 国产在线精品亚洲第一网站| 黄片无遮挡物在线观看| 好男人在线观看高清免费视频| 国产亚洲精品av在线| 久久精品夜夜夜夜夜久久蜜豆| 99在线视频只有这里精品首页| 色5月婷婷丁香| 久久精品夜夜夜夜夜久久蜜豆| 国产淫片久久久久久久久| 亚洲,欧美,日韩| 亚洲av男天堂| 亚洲欧美精品综合久久99| 亚洲中文字幕一区二区三区有码在线看| 国产高清视频在线观看网站| 国产一区二区在线av高清观看| 日韩,欧美,国产一区二区三区 | 欧美高清成人免费视频www| 91av网一区二区| 日韩欧美精品免费久久| 国产午夜精品论理片| 国产高清激情床上av| 在线免费十八禁| 岛国毛片在线播放| 久久久久久久午夜电影| 亚洲第一电影网av| 亚洲自拍偷在线| 高清毛片免费观看视频网站| 精品少妇黑人巨大在线播放 | 春色校园在线视频观看| 成人永久免费在线观看视频| 亚洲欧美成人综合另类久久久 | 青春草视频在线免费观看| 久久99精品国语久久久| 国产亚洲精品久久久com| 亚洲成av人片在线播放无| 九九久久精品国产亚洲av麻豆| 亚洲成人久久爱视频| 一个人看视频在线观看www免费| 免费av毛片视频| 一级黄片播放器| 大香蕉久久网| 91久久精品国产一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 自拍偷自拍亚洲精品老妇| 99riav亚洲国产免费| 97超碰精品成人国产| www.色视频.com| 人妻系列 视频| 国产老妇女一区| 一级毛片电影观看 | 久久人妻av系列| 亚洲精品久久久久久婷婷小说 | 免费在线观看成人毛片| 国产成人精品久久久久久| 一夜夜www| av在线老鸭窝| 99久久精品热视频| 国产一级毛片在线| 男人舔奶头视频| 一本久久精品| 国产精品一区二区在线观看99 | 此物有八面人人有两片| 给我免费播放毛片高清在线观看| 欧美不卡视频在线免费观看| 欧美一区二区亚洲| 啦啦啦韩国在线观看视频| 成人三级黄色视频| 国产中年淑女户外野战色| 一级黄色大片毛片| 欧美日韩乱码在线| 日韩在线高清观看一区二区三区| 亚洲国产欧美人成| 又黄又爽又刺激的免费视频.| 3wmmmm亚洲av在线观看| 午夜激情福利司机影院| 九草在线视频观看| 色综合亚洲欧美另类图片| 国产白丝娇喘喷水9色精品| 成人毛片60女人毛片免费| 一边亲一边摸免费视频| 99久久成人亚洲精品观看| 国产三级在线视频| 麻豆成人午夜福利视频| 亚洲人成网站在线播放欧美日韩| 欧美在线一区亚洲| 久久亚洲国产成人精品v| 五月玫瑰六月丁香| 99久久成人亚洲精品观看| 亚洲乱码一区二区免费版| 一本一本综合久久| 欧美三级亚洲精品| 中国美女看黄片| 亚洲成av人片在线播放无| 99热只有精品国产| 在线播放国产精品三级| 亚洲精品日韩av片在线观看| 国模一区二区三区四区视频| 精品久久久久久久久久免费视频| 亚洲无线观看免费| 亚洲欧美成人综合另类久久久 | 中文字幕制服av| 亚洲精品乱码久久久久久按摩| 99riav亚洲国产免费| 精品少妇黑人巨大在线播放 | 国产精品一区www在线观看| 三级男女做爰猛烈吃奶摸视频| 青春草亚洲视频在线观看| 真实男女啪啪啪动态图| 国产蜜桃级精品一区二区三区| 国产精品一二三区在线看| 九草在线视频观看| 久久久久久伊人网av| 亚洲国产欧美人成| 久久午夜福利片| 好男人在线观看高清免费视频| 国产成人精品一,二区 | 国产成人aa在线观看| 欧美变态另类bdsm刘玥| 久久久午夜欧美精品| 一个人免费在线观看电影| 国产精品久久电影中文字幕| 成人特级黄色片久久久久久久| 亚洲18禁久久av| 久久久久网色| kizo精华| 26uuu在线亚洲综合色| 国产精华一区二区三区| 国产精品人妻久久久久久| 亚洲精华国产精华液的使用体验 | 日本一本二区三区精品| 亚洲综合色惰| 久久草成人影院| 日韩人妻高清精品专区| 日韩亚洲欧美综合| 永久网站在线| 禁无遮挡网站| 在线天堂最新版资源| 国产精品久久视频播放| 一级黄色大片毛片| 亚洲精品久久国产高清桃花| .国产精品久久| 日日摸夜夜添夜夜爱| 欧美激情国产日韩精品一区| 免费观看人在逋| 国产成人午夜福利电影在线观看| 欧美一区二区国产精品久久精品| 可以在线观看毛片的网站| av又黄又爽大尺度在线免费看 | 成人一区二区视频在线观看| 国产成人精品婷婷| 国内精品久久久久精免费| 久久久久久伊人网av| 蜜桃久久精品国产亚洲av| 深夜a级毛片| 国产黄色视频一区二区在线观看 | 变态另类丝袜制服| 联通29元200g的流量卡| 久久精品国产自在天天线| 岛国在线免费视频观看| 欧美日本亚洲视频在线播放| 中文精品一卡2卡3卡4更新| 岛国在线免费视频观看| 在线播放国产精品三级| 草草在线视频免费看| 校园春色视频在线观看| 色噜噜av男人的天堂激情| 如何舔出高潮| 国产一区二区三区av在线 | 永久网站在线| 国产成人影院久久av| 黄色一级大片看看| 两个人视频免费观看高清| 日韩成人伦理影院| 成人漫画全彩无遮挡| 天堂av国产一区二区熟女人妻| 深夜精品福利| 国产日韩欧美在线精品| 久久久成人免费电影| 一卡2卡三卡四卡精品乱码亚洲| 国产三级中文精品| 国产中年淑女户外野战色| 日韩人妻高清精品专区| 秋霞在线观看毛片| 久久久午夜欧美精品| 最近的中文字幕免费完整| 国产高潮美女av| 欧美色视频一区免费| 国产成人精品久久久久久| 美女cb高潮喷水在线观看| 变态另类成人亚洲欧美熟女| 狠狠狠狠99中文字幕| 成人美女网站在线观看视频| 12—13女人毛片做爰片一| 男女视频在线观看网站免费| 人人妻人人澡人人爽人人夜夜 | 老女人水多毛片| 国产 一区 欧美 日韩| 午夜福利在线观看免费完整高清在 | 热99在线观看视频| 午夜久久久久精精品| 精品国产三级普通话版| 久久国内精品自在自线图片| 搞女人的毛片| 久久久久久久久大av| 乱系列少妇在线播放| 伦理电影大哥的女人| 国产成人a∨麻豆精品| 波野结衣二区三区在线| 能在线免费看毛片的网站| 久久久久久久久中文| 国产高清激情床上av| 久99久视频精品免费| 色综合色国产| 国产精品久久久久久av不卡| 在线a可以看的网站| 亚洲va在线va天堂va国产| 欧美最新免费一区二区三区| 中国美女看黄片| 国产老妇女一区| 国产精华一区二区三区| 两个人的视频大全免费| 国产人妻一区二区三区在| 你懂的网址亚洲精品在线观看 | 亚洲图色成人| 菩萨蛮人人尽说江南好唐韦庄 | 波多野结衣巨乳人妻| 丰满乱子伦码专区| 91狼人影院| 日本免费a在线| 国产伦一二天堂av在线观看| 狠狠狠狠99中文字幕| 国产黄片美女视频| 成人无遮挡网站| av天堂中文字幕网| 久久国内精品自在自线图片| 黄色日韩在线| 99在线人妻在线中文字幕| 日韩制服骚丝袜av| 亚洲va在线va天堂va国产| 麻豆成人av视频| 久久人人爽人人爽人人片va| 午夜精品国产一区二区电影 | 青春草视频在线免费观看| www日本黄色视频网| 欧美色视频一区免费| 看十八女毛片水多多多| or卡值多少钱| 亚洲在线自拍视频| 欧美日韩精品成人综合77777| 国产黄色视频一区二区在线观看 | 人妻夜夜爽99麻豆av| 美女内射精品一级片tv| 久久久欧美国产精品| 99热只有精品国产| 女的被弄到高潮叫床怎么办| 日本黄大片高清| 亚洲七黄色美女视频| 看免费成人av毛片| videossex国产| 久久精品91蜜桃| 久久久久九九精品影院| 日韩欧美 国产精品| 三级毛片av免费| 日本与韩国留学比较| 如何舔出高潮| 99热网站在线观看| 国产精品野战在线观看| 噜噜噜噜噜久久久久久91| 国产精品久久久久久精品电影小说 | av卡一久久| 欧美最新免费一区二区三区| 日本欧美国产在线视频| 亚洲激情五月婷婷啪啪| 国产一级毛片在线| 99热这里只有是精品50| 亚洲18禁久久av| 99久久精品一区二区三区| a级毛片a级免费在线| 欧美激情国产日韩精品一区| 性插视频无遮挡在线免费观看| 九九在线视频观看精品| 日韩亚洲欧美综合| 国产极品天堂在线| 校园春色视频在线观看| 波多野结衣高清作品| 亚洲人成网站在线观看播放| 一级毛片久久久久久久久女| 亚洲经典国产精华液单| 久久久色成人| 国产亚洲5aaaaa淫片| 免费看日本二区| 一区福利在线观看| 国产日韩欧美在线精品| 中国国产av一级| 亚洲国产精品合色在线| 亚洲精品久久国产高清桃花| 欧美一区二区精品小视频在线| 久久99热6这里只有精品| 看非洲黑人一级黄片| 一本精品99久久精品77| 国产视频首页在线观看| 网址你懂的国产日韩在线| 好男人在线观看高清免费视频| 国产亚洲5aaaaa淫片| 亚洲国产精品久久男人天堂| 国产成人福利小说| 人妻制服诱惑在线中文字幕| 亚洲av免费高清在线观看| 只有这里有精品99| 尤物成人国产欧美一区二区三区| 国产一区二区亚洲精品在线观看| 国产探花极品一区二区| 99久久精品国产国产毛片| 2022亚洲国产成人精品| 亚洲第一电影网av| 九九热线精品视视频播放| 欧美激情久久久久久爽电影| 毛片女人毛片| 亚洲中文字幕一区二区三区有码在线看| 好男人视频免费观看在线| 老师上课跳d突然被开到最大视频| 一个人看视频在线观看www免费| 淫秽高清视频在线观看| 联通29元200g的流量卡| 亚州av有码| 亚洲av电影不卡..在线观看| 一区福利在线观看| 国产在线精品亚洲第一网站| 久久久精品欧美日韩精品| 婷婷精品国产亚洲av| 久久久国产成人免费| 国产精品无大码| 久久久a久久爽久久v久久| 91久久精品国产一区二区三区| 久久久久国产网址| 最近手机中文字幕大全| 中文字幕av成人在线电影| 中文字幕免费在线视频6| 国产精品,欧美在线| 国产精品电影一区二区三区| 国产黄片视频在线免费观看| 在线观看免费视频日本深夜| 在线免费十八禁| 欧美日本亚洲视频在线播放| 在线播放国产精品三级| 国产一区二区亚洲精品在线观看| 成人高潮视频无遮挡免费网站| 国产成人福利小说| 伊人久久精品亚洲午夜| 国产亚洲精品久久久com| 只有这里有精品99| 级片在线观看| 日韩强制内射视频| 在线观看免费视频日本深夜| 精品无人区乱码1区二区| 欧美日韩精品成人综合77777| 老女人水多毛片| 欧美一区二区亚洲| 永久网站在线| 特大巨黑吊av在线直播| 男女边吃奶边做爰视频| 国产精品.久久久| 亚洲欧美中文字幕日韩二区| 日韩精品有码人妻一区| 在线国产一区二区在线| 嫩草影院入口| 悠悠久久av| 国产91av在线免费观看| 久久久精品94久久精品| 人体艺术视频欧美日本| 啦啦啦啦在线视频资源| 欧美精品国产亚洲| 成年版毛片免费区| 狂野欧美白嫩少妇大欣赏| 色综合亚洲欧美另类图片| 欧美日本视频| 国产亚洲av嫩草精品影院| 欧美日韩精品成人综合77777| 亚洲丝袜综合中文字幕| 欧美成人a在线观看| 国产乱人视频| 亚洲国产精品国产精品| 国产一区二区三区在线臀色熟女| 国产熟女欧美一区二区| 男人狂女人下面高潮的视频| 哪里可以看免费的av片| 国产精品.久久久| 中文字幕精品亚洲无线码一区| 亚洲av.av天堂| 亚洲欧美成人精品一区二区| 亚洲四区av| 看非洲黑人一级黄片| 波多野结衣巨乳人妻| 国产成人影院久久av| 亚洲精品日韩av片在线观看| 免费一级毛片在线播放高清视频| 色播亚洲综合网| 日日啪夜夜撸| 亚洲在线自拍视频| 夜夜爽天天搞| 麻豆久久精品国产亚洲av| 久久久久性生活片| 26uuu在线亚洲综合色| 久久久久久久久中文| 最近手机中文字幕大全| 久久欧美精品欧美久久欧美| 亚洲av电影不卡..在线观看| 国产精品无大码| 在线播放国产精品三级| 久久鲁丝午夜福利片| 噜噜噜噜噜久久久久久91| 尾随美女入室| 综合色丁香网| 日本成人三级电影网站| 国模一区二区三区四区视频| 高清午夜精品一区二区三区 | 老熟妇乱子伦视频在线观看| 亚洲四区av| 亚洲精品久久久久久婷婷小说 | 麻豆成人午夜福利视频| 2021天堂中文幕一二区在线观| 亚洲精品乱码久久久v下载方式|