• <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

    亚洲男人的天堂狠狠| 国产精品一区二区三区四区久久| 无遮挡黄片免费观看| 在线观看一区二区三区| 又粗又爽又猛毛片免费看| 精品国产亚洲在线| 日韩中文字幕欧美一区二区| 白带黄色成豆腐渣| 三级国产精品欧美在线观看| 国产真人三级小视频在线观看| www日本黄色视频网| 99精品欧美一区二区三区四区| 国产高清激情床上av| 九九热线精品视视频播放| 欧美一区二区亚洲| 一个人观看的视频www高清免费观看| a在线观看视频网站| 在线观看免费午夜福利视频| 99热这里只有精品一区| 一级毛片高清免费大全| h日本视频在线播放| 亚洲无线观看免费| 亚洲中文字幕一区二区三区有码在线看| 中文字幕av成人在线电影| 精品人妻1区二区| 婷婷六月久久综合丁香| 久久精品人妻少妇| 婷婷六月久久综合丁香| 嫩草影视91久久| 91av网一区二区| 19禁男女啪啪无遮挡网站| 亚洲欧美日韩卡通动漫| 亚洲五月婷婷丁香| 亚洲欧美精品综合久久99| 久久久久精品国产欧美久久久| 成人一区二区视频在线观看| 欧美乱码精品一区二区三区| 国内精品一区二区在线观看| 又黄又粗又硬又大视频| h日本视频在线播放| 国产亚洲精品久久久久久毛片| 身体一侧抽搐| 18禁黄网站禁片免费观看直播| 18禁国产床啪视频网站| 99在线人妻在线中文字幕| 日韩欧美三级三区| 在线免费观看的www视频| 久久久精品欧美日韩精品| 亚洲五月天丁香| 国产中年淑女户外野战色| 久久精品国产99精品国产亚洲性色| 热99re8久久精品国产| 18禁黄网站禁片午夜丰满| 欧美中文日本在线观看视频| 国产精品久久电影中文字幕| 久久国产精品人妻蜜桃| 人人妻人人看人人澡| 亚洲成av人片在线播放无| 欧美不卡视频在线免费观看| 亚洲中文字幕一区二区三区有码在线看| 国产爱豆传媒在线观看| 国产成人系列免费观看| 欧美午夜高清在线| 亚洲内射少妇av| 俄罗斯特黄特色一大片| 免费看十八禁软件| 91久久精品国产一区二区成人 | 在线播放国产精品三级| 亚洲人成网站高清观看| 99久久成人亚洲精品观看| 啦啦啦韩国在线观看视频| 欧美乱码精品一区二区三区| 在线a可以看的网站| 波多野结衣高清无吗| 色综合亚洲欧美另类图片| 91麻豆精品激情在线观看国产| 变态另类丝袜制服| 制服丝袜大香蕉在线| 嫁个100分男人电影在线观看| 国产单亲对白刺激| 免费在线观看日本一区| 国产亚洲精品一区二区www| 伊人久久精品亚洲午夜| 波野结衣二区三区在线 | 十八禁网站免费在线| 香蕉丝袜av| 一边摸一边抽搐一进一小说| 最后的刺客免费高清国语| 麻豆国产av国片精品| 老司机福利观看| 国内精品久久久久精免费| 欧美日韩精品网址| 床上黄色一级片| 黄色片一级片一级黄色片| 亚洲国产欧美网| 国产野战对白在线观看| 免费搜索国产男女视频| 免费观看精品视频网站| av中文乱码字幕在线| 一个人免费在线观看电影| 国产日本99.免费观看| 午夜福利成人在线免费观看| 亚洲国产欧美网| 999久久久精品免费观看国产| 亚洲精品在线观看二区| 亚洲18禁久久av| 一个人观看的视频www高清免费观看| 很黄的视频免费| 午夜激情福利司机影院| 内地一区二区视频在线| 精品久久久久久久人妻蜜臀av| 欧美精品啪啪一区二区三区| 一进一出好大好爽视频| 色精品久久人妻99蜜桃| 国内精品久久久久久久电影| 欧美中文日本在线观看视频| 欧美一级毛片孕妇| 国产在线精品亚洲第一网站| 小蜜桃在线观看免费完整版高清| 搞女人的毛片| 一级毛片高清免费大全| 精品国内亚洲2022精品成人| 两性午夜刺激爽爽歪歪视频在线观看| 一级黄色大片毛片| 欧美av亚洲av综合av国产av| 日本精品一区二区三区蜜桃| 一个人免费在线观看电影| www.www免费av| 首页视频小说图片口味搜索| 日本一本二区三区精品| 亚洲色图av天堂| 午夜福利在线观看免费完整高清在 | 欧美大码av| 999久久久精品免费观看国产| 亚洲精品国产精品久久久不卡| 日本撒尿小便嘘嘘汇集6| 又黄又爽又免费观看的视频| 色综合亚洲欧美另类图片| 国产欧美日韩一区二区精品| 国产精品一区二区三区四区免费观看 | 日韩有码中文字幕| 免费无遮挡裸体视频| 88av欧美| 中出人妻视频一区二区| 国语自产精品视频在线第100页| 村上凉子中文字幕在线| 精品久久久久久久毛片微露脸| 99久久无色码亚洲精品果冻| 亚洲欧美日韩高清专用| 熟女少妇亚洲综合色aaa.| 看黄色毛片网站| 亚洲欧美日韩东京热| 亚洲 欧美 日韩 在线 免费| 欧美av亚洲av综合av国产av| 国产成人福利小说| 嫩草影院入口| 亚洲欧美激情综合另类| 国内久久婷婷六月综合欲色啪| 香蕉久久夜色| 日韩欧美 国产精品| 国产精品av视频在线免费观看| 偷拍熟女少妇极品色| 狂野欧美白嫩少妇大欣赏| 久久国产精品人妻蜜桃| 免费无遮挡裸体视频| 午夜福利在线在线| 老汉色av国产亚洲站长工具| 麻豆国产av国片精品| 最近最新中文字幕大全电影3| 国产精品综合久久久久久久免费| 老鸭窝网址在线观看| 亚洲在线观看片| 精品人妻1区二区| 免费搜索国产男女视频| 免费av观看视频| 国产精品香港三级国产av潘金莲| 51国产日韩欧美| 国产精品国产高清国产av| 日本精品一区二区三区蜜桃| 精品久久久久久成人av| 波野结衣二区三区在线 | 熟女人妻精品中文字幕| 一个人观看的视频www高清免费观看| 三级国产精品欧美在线观看| 午夜福利高清视频| 欧美bdsm另类| 香蕉丝袜av| 久99久视频精品免费| 丰满人妻熟妇乱又伦精品不卡| 一个人看视频在线观看www免费 | 亚洲激情在线av| 91久久精品国产一区二区成人 | 亚洲国产精品合色在线| 国产野战对白在线观看| 嫩草影院精品99| 女警被强在线播放| 亚洲欧美日韩高清在线视频| 中文亚洲av片在线观看爽| 午夜a级毛片| 色播亚洲综合网| 99久国产av精品| 又黄又爽又免费观看的视频| 中文字幕熟女人妻在线| 国产色婷婷99| 亚洲欧美日韩无卡精品| bbb黄色大片| 国产黄片美女视频| 免费大片18禁| 国内精品久久久久精免费| 97人妻精品一区二区三区麻豆| 午夜福利视频1000在线观看| 成人国产综合亚洲| 亚洲内射少妇av| 熟女少妇亚洲综合色aaa.| 成人av一区二区三区在线看| 一边摸一边抽搐一进一小说| 欧美激情久久久久久爽电影| 最近视频中文字幕2019在线8| 色尼玛亚洲综合影院| 亚洲成人久久性| 一二三四社区在线视频社区8| 波多野结衣高清无吗| 亚洲电影在线观看av| 亚洲熟妇熟女久久| 1024手机看黄色片| 欧美乱妇无乱码| 精品人妻偷拍中文字幕| 亚洲熟妇熟女久久| 国产亚洲av嫩草精品影院| 亚洲第一欧美日韩一区二区三区| 91在线观看av| 日本免费一区二区三区高清不卡| 亚洲精品美女久久久久99蜜臀| 午夜免费激情av| 免费人成在线观看视频色| 最好的美女福利视频网| 亚洲av熟女| 国产精品久久电影中文字幕| 天堂网av新在线| 日本免费a在线| 蜜桃亚洲精品一区二区三区| 黑人欧美特级aaaaaa片| av天堂中文字幕网| 国产成年人精品一区二区| 久久午夜亚洲精品久久| 亚洲在线自拍视频| 九九在线视频观看精品| 国产精品久久久人人做人人爽| 欧美性感艳星| 精品午夜福利视频在线观看一区| 免费人成在线观看视频色| 熟女人妻精品中文字幕| 男女那种视频在线观看| 全区人妻精品视频| 色综合站精品国产| 99热这里只有精品一区| 99久久久亚洲精品蜜臀av| a级毛片a级免费在线| 无人区码免费观看不卡| 日日夜夜操网爽| 亚洲真实伦在线观看| 91在线精品国自产拍蜜月 | 国产免费一级a男人的天堂| 欧美性猛交黑人性爽| 亚洲av不卡在线观看| 网址你懂的国产日韩在线| 国产欧美日韩精品亚洲av| av片东京热男人的天堂| www日本在线高清视频| 亚洲精品色激情综合| 亚洲欧美激情综合另类| 欧美日韩精品网址| 99在线人妻在线中文字幕| 亚洲第一欧美日韩一区二区三区| 国产午夜福利久久久久久| 久久久久久久久久黄片| 俄罗斯特黄特色一大片| 最近在线观看免费完整版| 精品国内亚洲2022精品成人| 久久欧美精品欧美久久欧美| 日本a在线网址| 亚洲国产精品久久男人天堂| bbb黄色大片| 久久久久久大精品| 可以在线观看的亚洲视频| 国产在线精品亚洲第一网站| 午夜福利高清视频| 亚洲狠狠婷婷综合久久图片| 国产97色在线日韩免费| 国产亚洲精品综合一区在线观看| 热99在线观看视频| 麻豆成人午夜福利视频| 啦啦啦韩国在线观看视频| 国产一级毛片七仙女欲春2| 亚洲激情在线av| 国产亚洲av嫩草精品影院| 老司机深夜福利视频在线观看| www日本在线高清视频| 18禁黄网站禁片免费观看直播| 女人十人毛片免费观看3o分钟| 女人被狂操c到高潮| 在线国产一区二区在线| 少妇丰满av| 亚洲欧美日韩高清专用| 成人一区二区视频在线观看| 国产欧美日韩一区二区三| 成年人黄色毛片网站| 好看av亚洲va欧美ⅴa在| 国产精品久久久久久亚洲av鲁大| 国产精品自产拍在线观看55亚洲| 亚洲精品456在线播放app | av在线天堂中文字幕| 夜夜夜夜夜久久久久| 99精品在免费线老司机午夜| 国产乱人视频| 亚洲性夜色夜夜综合| 亚洲av免费高清在线观看| 99国产精品一区二区三区| 亚洲成人中文字幕在线播放| 免费搜索国产男女视频| 国产视频内射| 亚洲欧美日韩东京热| 欧美一区二区国产精品久久精品| 一级黄片播放器| 久久久久国产精品人妻aⅴ院| 国产真实伦视频高清在线观看 | 老司机在亚洲福利影院| 久久精品影院6| 99riav亚洲国产免费| 国产av不卡久久| 国产麻豆成人av免费视频| 日日摸夜夜添夜夜添小说| 国产一级毛片七仙女欲春2| 99在线人妻在线中文字幕| 亚洲欧美激情综合另类| 国产亚洲精品av在线| 宅男免费午夜| 国产成人欧美在线观看| 亚洲 国产 在线| 亚洲性夜色夜夜综合| 国产色爽女视频免费观看| 亚洲精品一区av在线观看| 中亚洲国语对白在线视频| 免费观看精品视频网站| 搡老岳熟女国产| 毛片女人毛片| 国产欧美日韩精品一区二区| 精品国产美女av久久久久小说| 18禁黄网站禁片午夜丰满| 波多野结衣巨乳人妻| 欧美一区二区精品小视频在线| 母亲3免费完整高清在线观看| 日韩大尺度精品在线看网址| 九九在线视频观看精品| 日韩欧美一区二区三区在线观看| 久久精品91无色码中文字幕| 亚洲在线自拍视频| 在线播放国产精品三级| 亚洲人成网站高清观看| 一区二区三区高清视频在线| av视频在线观看入口| 免费在线观看日本一区| 国产主播在线观看一区二区| 亚洲激情在线av| 大型黄色视频在线免费观看| 3wmmmm亚洲av在线观看| 久久6这里有精品| 狂野欧美白嫩少妇大欣赏| 欧美乱码精品一区二区三区| 免费av毛片视频| 好男人在线观看高清免费视频| 欧美乱码精品一区二区三区| 99久久精品热视频| 三级毛片av免费| 亚洲国产精品合色在线| 长腿黑丝高跟| 久久精品国产99精品国产亚洲性色| 成年女人看的毛片在线观看| 国产精品女同一区二区软件 | 在线观看一区二区三区| 最近在线观看免费完整版| 国产毛片a区久久久久| 国产国拍精品亚洲av在线观看 | 欧美日韩精品网址| 国产一区二区激情短视频| 99久久精品一区二区三区| 国产一区二区在线av高清观看| a级一级毛片免费在线观看| 亚洲av日韩精品久久久久久密| 色综合婷婷激情| 国产精品自产拍在线观看55亚洲| 久久精品国产亚洲av涩爱 | 亚洲av不卡在线观看| 男女下面进入的视频免费午夜| 国产成人啪精品午夜网站| 男女做爰动态图高潮gif福利片| 国产成人啪精品午夜网站| 日韩欧美一区二区三区在线观看| 国产成人av激情在线播放| 亚洲精品美女久久久久99蜜臀| 国产激情欧美一区二区| 乱人视频在线观看| 日本熟妇午夜| 国模一区二区三区四区视频| 熟女电影av网| 欧美日韩一级在线毛片| 国产三级在线视频| 亚洲男人的天堂狠狠| 免费av毛片视频| 国产成人a区在线观看| 床上黄色一级片| 午夜日韩欧美国产| 狠狠狠狠99中文字幕| 婷婷亚洲欧美| 18禁裸乳无遮挡免费网站照片| 免费av不卡在线播放| 两个人视频免费观看高清| 午夜日韩欧美国产| 欧美中文日本在线观看视频| 韩国av一区二区三区四区| 日韩欧美一区二区三区在线观看| 欧美色视频一区免费| 午夜两性在线视频| 两个人视频免费观看高清| 午夜精品一区二区三区免费看| 免费看a级黄色片| 99热6这里只有精品| 夜夜看夜夜爽夜夜摸| 18美女黄网站色大片免费观看| 在线十欧美十亚洲十日本专区| 免费av毛片视频| 蜜桃亚洲精品一区二区三区| 午夜福利高清视频| 日韩欧美精品免费久久 | 看黄色毛片网站| 老汉色av国产亚洲站长工具| www.熟女人妻精品国产| 村上凉子中文字幕在线| 免费观看的影片在线观看| 国产精品国产高清国产av| 有码 亚洲区| 制服丝袜大香蕉在线| 亚洲人成网站高清观看| 久久精品国产亚洲av香蕉五月| 免费看a级黄色片| 免费搜索国产男女视频| 无限看片的www在线观看| 亚洲无线在线观看| 一级毛片女人18水好多| 国产精品,欧美在线| 午夜福利成人在线免费观看| 99热只有精品国产| 亚洲国产欧美网| 国产高清有码在线观看视频| 国产精品99久久99久久久不卡| 给我免费播放毛片高清在线观看| ponron亚洲| 欧美中文日本在线观看视频| 久久精品国产清高在天天线| av专区在线播放| 国内毛片毛片毛片毛片毛片| 在线观看66精品国产| 国产欧美日韩一区二区精品| 无遮挡黄片免费观看| svipshipincom国产片| 亚洲精品美女久久久久99蜜臀| 午夜日韩欧美国产| 夜夜爽天天搞| 欧美一区二区精品小视频在线| 最新在线观看一区二区三区| 免费av毛片视频| av专区在线播放| 性色av乱码一区二区三区2| 亚洲av二区三区四区| 成年人黄色毛片网站| 人人妻,人人澡人人爽秒播| 成人亚洲精品av一区二区| 国产精品免费一区二区三区在线| 中文字幕高清在线视频| 亚洲av熟女| 国产午夜福利久久久久久| 国内精品久久久久久久电影| 久久这里只有精品中国| 99热精品在线国产| 日本a在线网址| 日本黄色视频三级网站网址| xxxwww97欧美| 国产精品久久久人人做人人爽| 啦啦啦韩国在线观看视频| 少妇裸体淫交视频免费看高清| 变态另类成人亚洲欧美熟女| 国产成人av教育| 国产一区在线观看成人免费| 两个人看的免费小视频| 深爱激情五月婷婷| 麻豆成人av在线观看| 欧美成狂野欧美在线观看| 国产三级在线视频| 中文在线观看免费www的网站| 成人特级av手机在线观看| 757午夜福利合集在线观看| 国产亚洲av嫩草精品影院| 日本五十路高清| 亚洲国产精品久久男人天堂| 久久久久久国产a免费观看| 色老头精品视频在线观看| 久久九九热精品免费| 人人妻,人人澡人人爽秒播| 国产久久久一区二区三区| 久久精品影院6| 亚洲va日本ⅴa欧美va伊人久久| 午夜福利视频1000在线观看| 久久人人精品亚洲av| 在线播放国产精品三级| 一个人免费在线观看的高清视频| 黄色成人免费大全| 一级黄片播放器| 黄色女人牲交| 免费av观看视频| 最近最新中文字幕大全免费视频| 亚洲国产精品合色在线| a级毛片a级免费在线| 亚洲五月天丁香| 少妇高潮的动态图| 日日干狠狠操夜夜爽| 欧美日韩国产亚洲二区| 欧美一区二区亚洲| 国产精品亚洲美女久久久| 观看美女的网站| 日韩精品中文字幕看吧| 欧美成人a在线观看| 欧美区成人在线视频| 黄片大片在线免费观看| 999久久久精品免费观看国产| 久久久色成人| 欧美一区二区精品小视频在线| 嫩草影院入口| 搡老熟女国产l中国老女人| 观看美女的网站| 国产伦一二天堂av在线观看| 天堂动漫精品| 亚洲第一电影网av| 88av欧美| 精品久久久久久久人妻蜜臀av| 嫩草影院精品99| 久久精品国产亚洲av涩爱 | 淫秽高清视频在线观看| 女同久久另类99精品国产91| 国产97色在线日韩免费| 亚洲18禁久久av| 亚洲国产精品合色在线| 国产成人aa在线观看| 国产伦精品一区二区三区视频9 | 午夜免费观看网址| 99精品在免费线老司机午夜| 一级毛片高清免费大全| 搡老妇女老女人老熟妇| av天堂中文字幕网| 国产精华一区二区三区| 国内精品美女久久久久久| 在线观看免费午夜福利视频| 欧美又色又爽又黄视频| 国产成人欧美在线观看| 在线观看一区二区三区| 天堂影院成人在线观看| 亚洲国产欧洲综合997久久,| 亚洲国产欧美网| а√天堂www在线а√下载| 婷婷精品国产亚洲av| 黄片大片在线免费观看| 12—13女人毛片做爰片一| 中文字幕精品亚洲无线码一区| 最好的美女福利视频网| 熟女少妇亚洲综合色aaa.| 宅男免费午夜| 免费一级毛片在线播放高清视频| 极品教师在线免费播放| 免费无遮挡裸体视频| 亚洲最大成人手机在线| 日本免费一区二区三区高清不卡| 国产69精品久久久久777片| 啦啦啦免费观看视频1| 99热这里只有是精品50| 久久久久久久精品吃奶| 欧美又色又爽又黄视频| 男女那种视频在线观看| 女人十人毛片免费观看3o分钟| 欧美+亚洲+日韩+国产| 国产 一区 欧美 日韩| 狠狠狠狠99中文字幕| 伊人久久大香线蕉亚洲五| 国产 一区 欧美 日韩| 国产高清有码在线观看视频| 亚洲精品亚洲一区二区| 男女那种视频在线观看| 日韩av在线大香蕉| 男人舔奶头视频| 亚洲av成人精品一区久久| 欧美zozozo另类| 国产精品影院久久| 亚洲在线自拍视频| 母亲3免费完整高清在线观看| 亚洲精品一区av在线观看| 怎么达到女性高潮| 久久精品91无色码中文字幕| 日韩欧美精品v在线| 窝窝影院91人妻| 成人无遮挡网站| 少妇人妻一区二区三区视频| 亚洲av中文字字幕乱码综合| 亚洲五月婷婷丁香| 男女下面进入的视频免费午夜| 欧美日本视频| 国产真实乱freesex| 久久久成人免费电影|