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

    Ultra-wide Bandwidth Observations of 19 Pulsars with Parkes Telescope

    2022-09-02 12:24:26ZuRongZhouJingBoWangNaWangGeorgeHobbsandShuangQiangWang

    Zu-Rong ZhouJing-Bo WangNa WangGeorge Hobbsand Shuang-Qiang Wang

    1 Xinjiang Astronomical Observatory,Chinese Academy of Sciences,Urumqi,Xinjiang 830011,China; wangjingbo@xao.ac.cn, na.wang@xao.ac.cn

    2 University of Chinese Academy of Sciences,19A Yuquan Road,Beijing 100049,China

    3 Xinjiang Key Laboratory of Radio Astrophysics,150 Science1-Street,Urumqi,Xinjiang 830011,China

    4 Key Lab of Radio Astronomy,Chinese Academy of Sciences,Beijing 100101,China

    5 CSIRO Astronomy and Space Science,PO Box 76,Epping,NSW 1710,Australia

    6 CAS Key Laboratory of FAST,NAOC,Chinese Academy of Sciences,Beijing 100101,China

    Abstract Flux densities are basic observation parameters to describe pulsars.In the most updated pulsar catalog,24%of the listed radio pulsars have no flux density measurement at any frequency.Here,we report the first flux density measurements,spectral indices,pulse profiles,and correlations of the spectral index with pulsar parameters for 19 pulsars employing the Ultra-Wideband Low receiver system installed on the Parkes radio telescope.The results for spectral indices of 17 pulsars are in the range between?0.6 and?3.10.The polarization profiles of thirteen pulsars are shown.There is a moderate correlation between the spectral index and spin frequency.For most pulsars detected,the signal-to-noise ratio of pulse profile is not high,so DM,Faraday rotation measure,and polarization cannot be determined precisely.Twenty-nine pulsars were not detected in our observations.We discuss the possible explanations for why these pulsars were not detected.

    Key words: (stars:) pulsars: general–stars: neutron–methods: data analysis

    1.Introduction

    Pulsars have been discovered for more than 53 yr(Hewish1968),though its radiation mechanism has not yet been fully understood(Jankowski et al.2018).One of the widely accepted models is that the radio emission regions are limited to the open polar cap inside the light-cylinder radius(Ruderman&Sutherland1975).However,the specific radiation area,physical process,and other details about pulsar radio emission are still unclear.Studying the radio spectra of pulsars is helpful to understand their emission mechanism.Unfortunately,only a few hundred pulsars have determined spectra.Only a tiny proportion of pulsars have been studied in a relatively wide frequency range(e.g.,Dai et al.2015).

    Measurements of flux density at multiple frequencies are needed when we determine the spectra of pulsars.The Australia Telescope National Facility (ATNF) Pulsar Catalog(Manchester et al.2005) provides a database for pulsar observational parameters.Pulsar flux densities in this catalog are relatively well-known near 1400 and 400 MHz,where most pulsars have been found,but few are known at other frequencies.We have referenced the most updated catalog,version 1.64,and out of the 2872 radio pulsars in the catalog,680 do not have flux density values in any radio waveband.Moreover,97.4%have no historical flux density measurements near 2 GHz,74.3% have no recorded flux density values near 400 MHz,and 32.5% have no recorded measurements of flux density close to 1400 MHz.The pulsar parameters in the catalog are collected by different telescopes and different generations of receivers and backends,and each has its own system error.Significant differences were found between multiple measurements(Levin et al.2013).Therefore,it is vital to obtain the absolute flux density calibrated measurements of pulsars(Jankowski et al.2018).Efforts to measure radio spectra began in earnest with Rankin et al.(1970) and Sieber (1973),followed by Malofeev &Malov (1980) and Izvekova et al.(1981).They measured the spectra at low frequencies near 100 MHz and below,finding that most pulsars have cliffy spectra that can be expressed as a simple power-law.Some pulsar spectra deviate from a simple power-law at low frequencies and show a turn-over(Rankin et al.1970),while some show a highfrequency cut-off in the form of a spectrum steepening or a break in the spectrum(Sieber1973).Lorimer et al.(1995) examined the spectra of 280 pulsars and obtained a mean spectral index of?1.6.A study was conducted of the spectral properties of 441 radio pulsars observed with the Parkes telescope (Jankowski et al.2018) and found about 79% of these pulsars could be classified as simple power-law spectra.

    People realized that pulsars were usually highly polarized long ago(Lyne &Smith1968).The long-term stability of the pulse profile and its complexity in terms of total intensity and polarization are also important characteristics of pulsars.The mean pulse profile and polarization properties help in understanding the geometry of the star,the pulse emission mechanism,and the beaming of pulsar radiation.Mean pulse profiles often have double or triple components,leading to different descriptions of the origin of the beamed emission(Backer1976),and the spectral index often varies from one component to the next (e.g.,Backer1972).Polarization properties of pulsars are normally described in terms of the four Stokes parameters.Many pulsars display a systematic variation across the pulse profile of the position angle (PA).The observed PA swings in many pulsars like the “S”-shaped curve are explained by the rotating vector model (Radhakrishnan&Cooke1969).Of course,the observed PA variations are not always continuous and smooth.Both normal pulsars and millisecond pulsars (MSPs) often can be seen with a high degree of linear polarization,and orthogonal-mode PA jumps(see e.g.,Yan et al.2011).Linear polarization is usually relatively stronger than circular polarization.Circular polarization most often has a sense reversal near the profile midpoint,usually associated with the core or central component of the profile(Rankin1983).

    In addition to the shape properties of the mean pulse profile,the width of the mean pulse profile is also used to describe the profile characteristics.The angular width of the mean pulse is also known as the width of the observed beam.W50andW10are pulse widths that are often mentioned and studied (see e.g.,Gould &Lyne1998).Pulsars exhibit a diverse frequency dependence of average pulse profile.Phillips &Wolszczan(1992) demonstrate the whole profile width and spacing between each component of three pulsars increase with the decrease of observation frequency from 4800 to 50 MHz using the Arecibo 305 m radio telescope.The increasing pulse width trend is very obvious in the frequency band below 1 GHz.Xilouris et al.(1996) extended the measurements of profile width to 32 GHz and combined their measurements with published values at lower frequencies,showing a better fit withW50of six pulsars.TheW10values versus the observing frequencies of 150 normal pulsars were well fitted with the Thorsett relationship(Thorsett1991;Chen &Wang2014).

    Lately,the Parkes 64 m radio telescope has been equipped with an ultra-wideband,low-frequency receiver.The frequency range of the receiver system is continuous from 704 to 4032 MHz.Moreover,the receiver has excellent sensitivity and polarization properties (Zhang et al.2019).The highsensitivity,well-calibrated and wide frequency observations with the Parkes telescope could efficiently provide us a systematic and uniform sample of pulsar flux densities(Hobbs et al.2020).In this work,we describe observations of 19 pulsars using the Ultra-Wideband Low(UWL)receiver system in Section2.Our results are shown in Section3.We discuss the results and give a concluding summary in Section4.

    2.Observations and Data Reduction

    We referenced the ATNF Pulsar Catalog version 1.59 to identify all radio pulsars without a flux density measurement near 1400 MHz that can be observed by the Parkes telescope.This list provides 505 pulsars.Then we removed all pulsars only detected in a survey that has a poorly defined sensitivity limit (for instance,some surveys are simply defined in the catalog as “miscellaneous” and very long observations with large telescopes may have detected the pulsars).After applying this selection,we have 327 pulsars remaining.Based on the known survey sensitivities,we know we can obtain a detectable signal close to 1400 MHz with observation durations of 30 minutes for all these pulsars.Still,some will be detectable at a high signal-to-noise ratio (S/N) with significantly shorter integration times.In particular,54 of these pulsars were discovered in the Molonglo surveys or with the Parkes 70 cm receiver(Manchester et al.1978).Such pulsars can be detected with only a few minute integration time.

    These 54 pulsars were observed with Parkes using the UWL receiver.All the details of the receiver and backends are described in Manchester et al.(2013) and Hobbs et al.(2020).Medusa and Parkes Digital Filter Banks 4 (PDFB4) are simultaneously used for the signal pre-processor systems in the observation.The frequencies of these two backends are centered at 2368 MHz and 1369 MHz,respectively.The total recording bandwidth of Medusa (3328 MHz) was subdivided into 3328 frequency channels,while PDFB4 of 256 MHz was subdivided into 1024 frequency channels.In the meantime,each pulsar period was divided into 1024 phase bins with each band.The data for all observations were recorded in a fold mode in 30 s and de-dispersed online coherently.Note that PDFB4 has no flux calibrator.

    The data were processed with the PSRCHIVE software package(Hotan et al.2014).We removed 5% of the band edges and manually excised data affected by narrow-band and impulsed radio frequency interference (RFI).We used the PSRCHIVE program PAZ to eliminate RFI automatically.Then we utilized PAZI to inspect the pulse profiles visually and remove subintegrations or frequency channels affected by RFI manually.We also relied on observations of the radio galaxy 3C 218(Hydra A)to transform the measured intensities to absolute flux densities,using on-and off-source pointing to measure the obvious brightness of the noise diode as a function of radio frequency.The pulsar observations were calibrated with their associated calibration files using PAC to transform the polarization products to Stokes parameters,to flatten the bandpass,and to calibrate the pulse profiles in flux density units.We formed a noise-free standard template from our observations using PAAS and then employed PSRFLUX to obtain the average flux density.All raw data of 54 pulsars can be downloaded from the CSIRO data archive(Hobbs et al.2011).

    3.Results

    We obtained data sets of 54 pulsars for which no flux densities were previously published near 1400 MHz,and havesuccessfully detected 25 pulsars.Only 14 of 25 pulsars were detected at 1369 and 2368 MHz at the same time.We gained pulse profiles and flux densities of 19 pulsars at 2368 MHz.No flux calibrator is available for the other six pulsars at 1369 MHz.

    3.1.Flux Density Measurements

    One of the fundamental properties of any astronomical source is its flux density,so we recommend that these flux densities should be included in the next version of the pulsar catalog.Our main results are listed in Table1,in which the pulsar name,pulse period,dispersion measure (DM),observation frequency,Modified Julian Date(MJD),observation times,observation length,S/N of the pulse profile,average flux density obtained at 2368 MHz,and width of pulse profiles at 50% (W50) are given in column order.Nineteen pulsars which are observed at 2368 MHz in Table1are mainly young pulsars with characteristic ages mostly between 106and 107yr,except for PSR J0348+0432 and PSR J2222?0137,which are in binary systems with a white dwarf (WD) companion with different masses.The observation time of most pulsars is about 10 minutes.(TheS2368represents flux density at 2368 MHz.)In Table1,flux densities range from 0.23 to 2.06 mJy.It is noted thatW50ranges from 3°to 21°.The error ofW50was estimated by determining how the width changes when the 50% flux density cuts across the profile move up or down by the baseline root mean square (rms) noise level (Dai et al.2015).W50andW10of 17 pulsars from literature at other frequencies are given in TableB2.TheW50values for PSR J1833?6023 at two frequencies are identical.TheW50of 10 pulsars decreases with frequency,but that of six pulsars increases with frequency.

    Table 1Pulsar Observational Parameters

    3.2.Spectral Indices

    The obtained flux densities can be used to measure the spectral indices.Many pulsars have been observed at one or more frequencies of 400,600,800,1400,2000,and 3000 MHz(e.g.,Dai et al.2015;Jankowski et al.2018).To estimate spectral indices,we divide the UWL data sets into three subbands and measure the flux density for each of them.Their center frequencies are close to 1400,2000,and 3000 MHz,with sub-bandwidths of 400 MHz,400 MHz,and 600 MHz,respectively.The results are listed in Table2asS1400,S2000,andS3000.Except for six pulsars (PSR J1057?4754,PSR J1805+0306,PSR J1810?5338,PSR J1854?1421,PSR J1903?0632,and PSR J2222?0137),the flux density of 13 other pulsars cannot be measured in all the three sub-bands since too many frequency channels are removed.The pulse profiles of only one or two sub-bands can be seen for some pulsars.Previously publishedS150,S200,S300,S400,S600,S700,S800,andS6000for these pulsars are used to estimate the spectral index,and the references are given in the last column of Table1.We assume a simple power-law of the formSν=bxα,where,ν is the center frequency and ν0=1.4 GHz a constant reference frequency.The fit parameters are the spectral index α and a constantb,and their values are provided in Table2.The spectral behavior of the 17 pulsars can be described well by a simple power-law over the frequency range considered.The spectral indices of all the 17 pulsars are between?0.6 and?3.10.For PSR J1810?5338,a broken power-law is better to describe the spectrum with a 600 MHz cut-off frequency,and α1is the spectral index before and α2the one after the break.The spectral plots of 18 pulsars and power-law fits to the data are depicted in Figure1.According to previous statistics,the mean spectral index of pulsars with a simple power-law spectrum is about?1.60.The majority of our spectral indices are flatter than the mean values from Jankowski et al.(2018),but our measurements are not particularly unusual as there are many other pulsars that also have similar spectral indices.As our sample is relatively small,we do not think our pulsar sample is particularly inclined to flat spectrum pulsars.

    Table 3Correlation of Spectral Index α with lg|x| for Different Pulsar Parameters x

    3.3.Polarization Profiles

    Pulse profiles at 2368 MHz for all the 19 pulsars are presented for the first time.The profiles of 20 among 25 pulsars were published in the European Pulsar Network (EPN)database near 1.4 GHz for the first time except for five pulsars(PSR J1057?4754,PSR J1604?7203,PSR J1625?4048,PSR J1728?0007,and PSR J1805+0306).7 http://www.epta.eu.org/epndbThe pulse profiles at 1369 MHz and 2368 MHz of most pulsars are similar in our observations.Most of the average pulse profiles in FiguresA1andA2have a single peak,and the pulse profiles are narrow.Figures2–4display the calibrated polarization profiles for 13 pulsars in our sample.A rotation measure (RM) fitting during the calibration is considered during the calibration of polarization.Only six pulsars show strong polarization characteristics at two frequencies simultaneously,while the other seven pulsars only have polarization profiles at 1369 or 2368 MHz.Polarization cannot be measured precisely for other pulsars due to the limited S/N.Results for individual pulsars are described in detail as follows.

    3.3.1.PSR J0348+0432

    PSR J0348+0432,a binary pulsar in a 2.46 hr orbit with a low-mass ELL1 WD companion,was first discovered in the Green Bank Telescope(GBT)350 MHz drift scan survey.The pulse profile of the pulsar at 350 MHz has three components(Lynch et al.2013),while Figure3shows a single peak at 2368 MHz and is linearly polarized.There is significant right-hand circular polarization and slight shallow rotation in the PA variation across the profile.

    3.3.2.PSR J1157?5112

    PSR J1157?5112 was discovered in the first high-frequency survey of intermediate Galactic latitudes.It is a 44 ms pulsar and the first recycled pulsar with an ultramassive (M>1.14Msolar) WD companion(Edwards et al.2001).The polarization profile for PSR J1157?5112 shows a single peak at 1369 MHz in Figure4.The profile is linearly polarized,and there is little or no variation in the PA.No significant circular polarization is observed across the whole profile.

    3.3.3.PSR J1420?5416(PSR B1417?54)

    PSR J1420?5416 was first discovered at 408 MHz in the Molonglo survey(Manchester et al.1978).Figure3affirms the pulse profile of this pulsar consists of two closely spaced components at 2368 MHz.The leading component is less linearly polarized than the trailing component,and the PA exhibits a continuous decrease across the trailing component.There is significant right-hand circular polarization.

    3.3.4.PSR J1615?2940(PSR B1612?29)

    Similar to PSR J1420?5416,PSR J1615?2940 was also discovered in the second Molonglo pulsar survey (Manchester et al.1978).Figure3depicts a pulse profile with two components at 2368 MHz.The pulse profile is very low in linear polarization,and thus there is very little PA measured.There is a sense reversal of the circular polarization from righthand to left-hand under the trailing part of the profiles.

    3.3.5.PSR J1728?0007(PSR B1726?00)

    PSR J1728?0007 was found in a survey for short period pulsars(Stokes et al.1985).Weisberg et al.(1999) published its profile at 1418 MHz,but the linear polarization is too weak to measure.Figure3demonstrates the pulse profile consists of two components at 2368 MHz.The pulsar is linearly polarized with a continuously increasing PA across the profile.

    3.3.6.PSR J1746?2856

    PSR J1746?2856,a highly dispersed pulsar,was discovered in the direction of the Galactic Center at 3.1 GHz with the Parkes radio telescope (Johnston et al.2006).The mean pulse profiles and polarization parameters for this pulsar at 2368 MHz are shown in Figure3.The degree of linear polarization is low,with very little PA measured.There are probably four or five distinct components,but it is impossible to classify this pulsar without additional information.

    Figure 1.Spectra for 18 pulsars and power-law fits to the data (blue lines).The red bars give the uncertainties of the flux densities at the different frequencies.

    Figure 2.Average polarization profiles for PSR J1805+0306,PSR J1810?5338,PSR J1833?6023,PSR J1848?1952,PSR J1903?0632 and PSR J1854?1421 at 1369 and 2368 MHz.The total intensity is shown in black,while the linear polarization and circular-polarization are shown in red and blue respectively.

    Figure 2.(Continued.)

    Figure 3.Average polarization profiles for PSR J0348+0432,PSR J1420?5416,PSR J1615?2940,PSR J1728?0007,PSR J1746?2856,and PSR J2222?0137 at 2368 MHz.The total intensity is shown in black,while linear polarization and circular polarization are shown in red and blue respectively.

    Figure 4.Average polarization profile for PSR J1157?5112 at 1369 MHz.

    3.3.7.PSR J1805+0306(PSR B1802+03)

    Polarization profiles at 2368 and 1369 MHz for PSR J1805+0306 are presented in Figure2.The profiles display two components at both frequencies.The leading part of the profile has relatively high fractional linear polarization,whereas the trailing part of the profile is essentially unpolarized at both frequencies.There is a shallow rotation of the PA through the leading component at 2368 MHz.The polarization profiles at 1369 MHz are in good agreement with those presented by Weisberg et al.(1999) at Arecibo 1418 MHz.

    3.3.8.PSR J1810?5338(PSR B1806?53)

    As is visible in Figure2,the mean pulse profiles of PSR J1810?5338 are almost the same at the two frequencies.There are three pulse components of the pulsar,which are similar to that at 660 MHz(Jiang &You2014).The leading and trailing components are weaker than the central component of the pulse profile.The central component is less linearly polarized than the leading and trailing components at both frequencies.Flat increasing PA swings through the overall profiles can be seen both at 2368 and 1369 MHz.There is a sense reversal of the circular polarization from right-hand to left-hand under the whole profiles at both frequencies.

    3.3.9.PSR J1833?6023(PSR B1828?60)

    The mean pulse profiles and polarization parameters for PSR J1833?6023 at 2368 and 1369 MHz are similar,which are shown in Figure2.The pulsar was first discovered at a frequency of 408 MHz by the second Molonglo pulsar survey(Manchester et al.1978).PSR J1833?6023 has a profile with only one narrow component,and the profile is linearly polarized.There is a strong sense reversal of the circular polarization from right-hand to left-hand under the overall pulse,and a steep decreasing PA swings across it.

    3.3.10.PSR J1848?1952(PSR B1845?19)

    PSR J1848?1952 was first discovered by an extensive survey of pulsars which had been undertaken using observations at the Molonglo Radio Observatory and the Australian National Radio Astronomy Observatory,Parkes.The observing frequency for both the Molonglo and Parkes observations was 408 MHz(Manchester et al.1978).As displayed in Figure2,the mean pulse profile of PSR J1848?1952 at both 2368 and 1369 MHz consists of two close sharp components.The PA variations are different,but the linear and circular polarizations are almost the same at two frequencies.The overall pulse has obvious linear polarization,and the trailing component is less linearly polarized than the leading component.A high degree of right-handed circular polarization was observed across the whole profile.There is a continuous decreasing PA across the profile at 2368 MHz,while the PA variation is complex with a shallow rotation swing through it at 1369 MHz.

    3.3.11.PSR J1854?1421(PSR B1851?14)

    Similar to PSR J1848?1952,PSR J1854?1421 was first also discovered by the second Molonglo pulsar survey(Manchester et al.1978).The pulse profile of the pulsar in Figure2exhibits a single peak at both 2368 and 1369 MHz and is linearly polarized.There is a continuous decreasing PA across the profile at 1369 MHz,while the PA variation is discontinuously decreasing with a PA jump at 2368 MHz.No significant circular polarization is observed across the whole profile.

    3.3.12.PSR J1903?0632(PSR B1900?06)

    Mean pulse profiles and polarization parameters for PSR J1903?0632 at 2368 and 1369 MHz are very similar,as displayed in Figure3.The pulsar was discovered during a systematic search for pulsars at a frequency of 408 MHz,carried out with the Jodrell Bank Mark IA radio telescope in 1972 (Davies et al.1972).Similar to 408 MHz,the pulse profile of the pulsar also has two pulse components at two frequencies.The trailing component is less linearly polarized than the leading component.PSR J1903?0632 shows a small amount of circular polarization from left-hand to right-hand under the overall profile.The PA variation across two components appears continuous and has a negative slope.

    3.3.13.PSR J2222?0137

    PSR J2222?0137,a 2.44 day binary pulsar with a massive CO WD companion,was first discovered in the GBT 350 MHz drift scan survey.The mean pulse profile and polarization parameters for PSR J2222?0137 at 2368 MHz given inFigure3have far more details than previously published results at 820 MHz (Boyles et al.2013).The pulse profile is very similar at both frequencies,with only one pulse component.The linear polarization is low under the overall profile.There is a hint of circular polarization from left-hand to right-hand against the profile.The PA variation across the profile is complex with regions of increasing and decreasing PAs,and a PA jump can be seen at phase 0.504.

    3.4.Correlations of Spectral Index with Pulsar Parameters

    We test for a correlation between spectral index α and lg|x|for all the 18 pulsars,wherexis one of the pulsar parameters below.The parameters are spin frequency ν,spin-down rate ˙ν,the magnetic field at the light cylinder radiusBLC,the characteristic age тc,the surface magnetic fieldBsurf,the spin-down luminosity ˙E,pulse periodP,and period derivative ˙Pof the pulsar.We took all values of pulsar parameters from the ATNF Pulsar Catalog.Most are covariant because these quantities depend on basic pulsar parameters,such as pulse period and its derivative.

    We first measured the correlation by visual inspection,then computed the Spearman rank correlation coefficient to characterize its strength.We test all pulsars in our single power-law data set first,then the isolated pulsars.There are only two pulsars in the binary system.The correlation coefficients,correspondingp-values,and the number of pulsarsN,for which correlations are computed,are listed in Table3.

    We find a moderate correlation between the spectral index and the pulse period for normal pulsars.All the other combinations show no correlation.For isolated pulsars,we find a mildly negative correlation of spectral index with the surface magnetic field.

    4.Discussion and Conclusions

    The minimum detectable flux density (Smin) of the UWL system can be evaluated using the radiometer equation(Manchester et al.1996)

    where α=S/N,β is the factor of digitization and other processing losses,Tsysis the sum of receiver temperature and sky temperature (K),Gis the gain of the telescope,Npis the number of polarizations(two in this case),Δν is the observing bandwidth (MHz),Tis the observation length (s),Wis the effective pulse width in time units,andPis the pulse period(Dewey et al.1985).According to Equation (1),the minimum detectable flux density of the UWL system varies from pulsar to pulsar with differentW,P,andT.Here,we adopt α=10,β=1.5 (Manchester et al.1996),G=0.64 K Jy?1(Edwards et al.2001),Tsys=22 K,and Δν=3328 MHz(Hobbs et al.2020).The minimum detectable flux density(Smin)of the UWL system ranges from 0.02 to 0.76 mJy for pulsars with different duty cycles.

    However,54 pulsars have been observed and only 25 have been detected.Among the 29 pulsars undetected,six pulsars are MSPs with periods from 1.7 to 3.7 ms,and the periods of the other 23 pulsars range from 44 ms to 4.3 s.We estimated the flux densities at 2368 MHz assuming the spectral index α=?1.6 for undetected pulsars with only one published flux density (the flux densities of PSR J1231?1411 and J2256?1024 were estimated separately).As shown in TableB1,all the estimatedS2368of 22 pulsars are potentially detectable by the UWL system.

    1.Equation (1) is not suitable for estimating the sensitivity of the UWL system.Because the bandwidth of the UWL system is very wide,the flux densities of pulsars in the bandwidth change significantly and decrease in the highfrequency part.

    2.The flux density of a radio pulsar decreases rapidly with increasing observing frequency.According to previous statistics,the spectral index of pulsars varies significantly.As expressed in TableB1,the estimated flux density at 2368 MHz for some pulsars is significantly higher than the detection threshold of the UWL system.Those undetected pulsars may have a steeper radio spectrum.

    3.The bandwidth of the UWL system cannot be fully utilized.Part of the bandwidth is contaminated by RFI and has to be removed.Therefore,the sensitivity of the UWL system is not as high as we calculated using Equation (1).In fact,most of the frequency channels from 700 to 1000 MHz have been removed,and the radio emission of pulsars is bright at the low end of the UWL system.

    4.The radio emission of pulsars is strongly affected by propagation through the interstellar medium (Kumamoto et al.2021).The typical observing length of our observations is only 10 minutes,but the timescale of diffractive scintillation could be hours,which is much longer than our observing length.

    The correlation coefficients for pulsars in our sample,along with previous results from Zhao et al.(2019),Jankowski et al.(2018),Han et al.(2016),and Lorimer et al.(1995),are listed in Table3.Jankowski et al.(2018) obtained Spearman rank correlations of the spectral index with various pulsar parameters for 276 pulsars with simple power-law spectra.They found moderate correlations with ˙ν,BLC,and ˙Efor young pulsars.Han et al.(2016)and Lorimer et al.(1995)ascertained a very weak correlation for the spectral index with ˙E.For pulsars in our sample,we find moderate correlations between spectral index with ν,which are in the sense of steeper spectra for fast-rotating pulsars.No consistent correlation is found across the different samples,and results from disparate samples are distinct.

    Table B1Estimated Flux Densities at 2368 MHz (S2368) for 22 Pulsars

    Table B2Pulse Widths from Literature for 17 Pulsars

    In this work,we have carried out ultra-wide bandwidth observations of 19 pulsars with the Parkes telescope.Flux density measurements,spectral properties,polarization profiles,and pulse widths at 2368 MHz or 1369 MHz have been presented.The non-detection of polarization in 12 out of 25 pulsars is perhaps not be surprising.As depicted in FigureA1,the S/N of some pulse profiles (PSRs J0057?7201,J1721?1939,J1816?5548,etc.) is low.It is not very likely to measure polarization with a low S/N.Other pulse profiles have a medium S/N.The degree of linear polarization for these pulsars may be low.To measure the pulse profile,polarization and RM with high precision and study their evolution with frequency,longer observations for these pulsars are needed.We also note that there is a large number of pulsars in the catalog that still do not have such basic information as the flux density near 1.4 GHz.Observing each pulsar for a longer time or using telescopes with larger observation apertures (such as FAST)at a specific frequency is needed when we continue this work in the future.

    Acknowledgments

    The Parkes radio telescope is part of the ATNF which is funded by the Commonwealth of Australia for operation as a National Facility managed by CSIRO.This paper includes archived data obtained through the CSIRO Data Access Portal(http://data.csiro.au).This work is supported by the National Natural Science Foundation of China (Grant No.NSFC12041304),National SKA Program of China (No.2020SKA0120100),Youth Innovation Promotion Association of Chinese Academy of Sciences,National Key Research and Development Program of China (No.2017YFA0402602),the CAS Jianzhihua project,and Heaven Lake Hundred-Talent Program of Xinjiang Uygur Autonomous Region of China.

    Appendix A

    Figures of Averaged Pulse Profiles

    Figure A1.Averaged pulse profiles at 2368 MHz were obtained with the Parkes telescope for seven pulsars.

    Figure A2.Averaged pulse profiles of 13 pulsars at a center frequency of 1369 MHz were obtained with the Parkes telescope.

    Figure A2.(Continued.)

    Averaged pulse profiles at 2368 MHz or 1369 MHz were obtained with the Parkes telescope for 16 pulsars,five of which were obtained at both frequencies.

    Appendix B

    Estimated Flux Densities at 2368 MHz (S2368) of 22 Pulsars and Pulse Widths at other Frequencies for 17 Pulsars

    Estimated flux densities of 22 Pulsars were estimated at 2368 MHz (S2368).Pulse widths for 17 pulsars were measured previously at other frequencies.

    亚洲男人天堂网一区| 法律面前人人平等表现在哪些方面| 99国产精品99久久久久| 国产日韩欧美视频二区| 乱人伦中国视频| 他把我摸到了高潮在线观看 | 美女福利国产在线| 在线观看www视频免费| 国产一区有黄有色的免费视频| 91精品三级在线观看| 午夜福利免费观看在线| 汤姆久久久久久久影院中文字幕| 一区二区av电影网| 老司机午夜十八禁免费视频| 精品一区二区三区av网在线观看 | 高清视频免费观看一区二区| 欧美黑人精品巨大| 人妻 亚洲 视频| 日韩欧美三级三区| 亚洲第一av免费看| 亚洲精品国产一区二区精华液| 日本av免费视频播放| 国产欧美日韩一区二区精品| 精品一区二区三区四区五区乱码| 三级毛片av免费| 日韩大片免费观看网站| 法律面前人人平等表现在哪些方面| 亚洲中文日韩欧美视频| 一区二区三区乱码不卡18| 欧美国产精品va在线观看不卡| 亚洲国产欧美网| 国产一区二区 视频在线| 淫妇啪啪啪对白视频| 制服人妻中文乱码| 免费久久久久久久精品成人欧美视频| 亚洲欧洲精品一区二区精品久久久| 啦啦啦视频在线资源免费观看| 老司机在亚洲福利影院| 亚洲第一av免费看| 亚洲av第一区精品v没综合| 在线av久久热| 色播在线永久视频| 好男人电影高清在线观看| 岛国毛片在线播放| 后天国语完整版免费观看| 久久人妻av系列| 国产视频一区二区在线看| 国产免费av片在线观看野外av| av网站免费在线观看视频| 人妻久久中文字幕网| 免费日韩欧美在线观看| 少妇精品久久久久久久| 精品免费久久久久久久清纯 | 啦啦啦视频在线资源免费观看| 午夜视频精品福利| 99在线人妻在线中文字幕 | 国产1区2区3区精品| 每晚都被弄得嗷嗷叫到高潮| 国产在线视频一区二区| 波多野结衣一区麻豆| 精品国产乱子伦一区二区三区| 在线观看www视频免费| 国产精品二区激情视频| 欧美乱妇无乱码| 国产在线精品亚洲第一网站| 蜜桃国产av成人99| 国产欧美日韩精品亚洲av| 亚洲天堂av无毛| 精品国内亚洲2022精品成人 | 亚洲国产毛片av蜜桃av| 久久人妻福利社区极品人妻图片| 亚洲免费av在线视频| 大片免费播放器 马上看| 老熟女久久久| 亚洲av电影在线进入| 亚洲伊人久久精品综合| av电影中文网址| 色在线成人网| 免费观看人在逋| 91麻豆av在线| 国产精品国产高清国产av | 久久天躁狠狠躁夜夜2o2o| www.精华液| 一级a爱视频在线免费观看| 日本欧美视频一区| 老汉色av国产亚洲站长工具| 国产精品av久久久久免费| 中文字幕人妻丝袜制服| 狠狠狠狠99中文字幕| 欧美日韩亚洲综合一区二区三区_| 国产有黄有色有爽视频| 亚洲精华国产精华精| 国产精品美女特级片免费视频播放器 | 亚洲熟女毛片儿| 久久中文看片网| 久久精品亚洲av国产电影网| 亚洲欧美日韩高清在线视频 | 亚洲第一欧美日韩一区二区三区 | 成在线人永久免费视频| 亚洲成a人片在线一区二区| 国产在视频线精品| 欧美精品高潮呻吟av久久| 久久精品国产亚洲av高清一级| 老司机福利观看| 成年版毛片免费区| 日韩欧美一区视频在线观看| 国产精品久久久久久人妻精品电影 | 久久久久久人人人人人| 国产成人精品在线电影| 性少妇av在线| 午夜91福利影院| 国产一区有黄有色的免费视频| 欧美精品一区二区大全| 国产片内射在线| 亚洲精品中文字幕在线视频| 大香蕉久久网| aaaaa片日本免费| 一本大道久久a久久精品| 亚洲国产看品久久| 久久免费观看电影| 国产日韩一区二区三区精品不卡| 悠悠久久av| 国产精品一区二区在线观看99| 久久久久久亚洲精品国产蜜桃av| 久久精品亚洲熟妇少妇任你| 国产淫语在线视频| 99riav亚洲国产免费| 99精品在免费线老司机午夜| 亚洲av成人不卡在线观看播放网| 少妇裸体淫交视频免费看高清 | 久久精品91无色码中文字幕| 久久影院123| 婷婷成人精品国产| 成人亚洲精品一区在线观看| 欧美黄色淫秽网站| 五月天丁香电影| 欧美黑人精品巨大| 制服诱惑二区| 18禁裸乳无遮挡动漫免费视频| 91麻豆av在线| 丝袜美足系列| 黄色毛片三级朝国网站| 搡老熟女国产l中国老女人| 曰老女人黄片| 精品国产亚洲在线| 精品一品国产午夜福利视频| 中文字幕另类日韩欧美亚洲嫩草| 国产激情久久老熟女| 日韩一区二区三区影片| 日韩三级视频一区二区三区| 最黄视频免费看| 久久人人爽av亚洲精品天堂| 亚洲午夜精品一区,二区,三区| 欧美日韩黄片免| 久久精品国产a三级三级三级| 欧美+亚洲+日韩+国产| 久久亚洲精品不卡| 中文字幕另类日韩欧美亚洲嫩草| 少妇的丰满在线观看| 在线观看66精品国产| 精品少妇黑人巨大在线播放| 国产麻豆69| 久久久久网色| 亚洲成人国产一区在线观看| 国产极品粉嫩免费观看在线| 一区二区三区乱码不卡18| 欧美黑人精品巨大| 久9热在线精品视频| 欧美乱码精品一区二区三区| 精品国产一区二区三区四区第35| 最新的欧美精品一区二区| 50天的宝宝边吃奶边哭怎么回事| 十八禁高潮呻吟视频| 国产成人精品无人区| 精品第一国产精品| 亚洲专区中文字幕在线| 国内毛片毛片毛片毛片毛片| 男女高潮啪啪啪动态图| 久久婷婷成人综合色麻豆| 露出奶头的视频| 亚洲精品乱久久久久久| 热99国产精品久久久久久7| 91麻豆精品激情在线观看国产 | 久久精品91无色码中文字幕| 一级毛片精品| 亚洲国产看品久久| 欧美日韩视频精品一区| 咕卡用的链子| 久久精品国产亚洲av香蕉五月 | 一级毛片精品| 欧美在线黄色| 在线观看免费日韩欧美大片| 丝袜美足系列| 亚洲国产看品久久| 国产欧美日韩一区二区精品| 一本综合久久免费| 考比视频在线观看| 国产熟女午夜一区二区三区| 久久精品国产综合久久久| av视频免费观看在线观看| 精品熟女少妇八av免费久了| 亚洲一区中文字幕在线| 老汉色av国产亚洲站长工具| 国产单亲对白刺激| 亚洲欧美日韩另类电影网站| 亚洲av日韩精品久久久久久密| av又黄又爽大尺度在线免费看| 亚洲三区欧美一区| 国产男女超爽视频在线观看| 男女高潮啪啪啪动态图| 亚洲 国产 在线| 中文字幕人妻丝袜制服| 搡老乐熟女国产| 丝瓜视频免费看黄片| 91成年电影在线观看| 国产av精品麻豆| 亚洲人成电影免费在线| 蜜桃在线观看..| 高清黄色对白视频在线免费看| 国产精品欧美亚洲77777| 婷婷成人精品国产| 在线天堂中文资源库| 桃花免费在线播放| 亚洲 欧美一区二区三区| 一边摸一边做爽爽视频免费| 少妇被粗大的猛进出69影院| 1024视频免费在线观看| 首页视频小说图片口味搜索| 成人18禁高潮啪啪吃奶动态图| 久久久久久久国产电影| 亚洲av美国av| 中亚洲国语对白在线视频| 精品第一国产精品| 国产一区二区三区视频了| 亚洲视频免费观看视频| 亚洲人成电影免费在线| 超碰成人久久| a级片在线免费高清观看视频| 18禁黄网站禁片午夜丰满| 亚洲欧美一区二区三区久久| 欧美av亚洲av综合av国产av| av有码第一页| 成人国产av品久久久| 国产成人啪精品午夜网站| 9热在线视频观看99| 久久中文看片网| 亚洲伊人久久精品综合| 久热爱精品视频在线9| 两个人免费观看高清视频| 亚洲色图 男人天堂 中文字幕| 亚洲 欧美一区二区三区| 波多野结衣av一区二区av| 久久天躁狠狠躁夜夜2o2o| 国产99久久九九免费精品| 精品亚洲乱码少妇综合久久| 亚洲欧美一区二区三区黑人| 日韩中文字幕视频在线看片| 两性午夜刺激爽爽歪歪视频在线观看 | 日韩 欧美 亚洲 中文字幕| 母亲3免费完整高清在线观看| 久久久久视频综合| 王馨瑶露胸无遮挡在线观看| 日韩视频在线欧美| 久久久久久久精品吃奶| 91大片在线观看| 欧美国产精品va在线观看不卡| 丝袜喷水一区| 中亚洲国语对白在线视频| 欧美日韩亚洲国产一区二区在线观看 | 久久中文字幕人妻熟女| 国产精品亚洲av一区麻豆| 亚洲国产欧美网| 国产高清激情床上av| 99riav亚洲国产免费| 亚洲欧美日韩高清在线视频 | 午夜福利在线观看吧| www.熟女人妻精品国产| 少妇 在线观看| 他把我摸到了高潮在线观看 | 亚洲欧美激情在线| 久久中文看片网| 别揉我奶头~嗯~啊~动态视频| 免费一级毛片在线播放高清视频 | 成人亚洲精品一区在线观看| 国产aⅴ精品一区二区三区波| 新久久久久国产一级毛片| 久久性视频一级片| www日本在线高清视频| 成人免费观看视频高清| 亚洲综合色网址| 精品亚洲乱码少妇综合久久| 后天国语完整版免费观看| 成年版毛片免费区| 国产av精品麻豆| 一级毛片女人18水好多| 一区二区日韩欧美中文字幕| 夜夜夜夜夜久久久久| 亚洲国产成人一精品久久久| 桃红色精品国产亚洲av| 欧美精品啪啪一区二区三区| 天天添夜夜摸| 窝窝影院91人妻| 考比视频在线观看| 精品久久蜜臀av无| 久久精品国产a三级三级三级| 香蕉国产在线看| 国产精品98久久久久久宅男小说| 丝袜美腿诱惑在线| 免费一级毛片在线播放高清视频 | 91九色精品人成在线观看| 日本黄色日本黄色录像| 国产人伦9x9x在线观看| 精品视频人人做人人爽| 首页视频小说图片口味搜索| 成人国产一区最新在线观看| 成年人黄色毛片网站| 一级黄色大片毛片| 老鸭窝网址在线观看| 中文字幕精品免费在线观看视频| 女性生殖器流出的白浆| 50天的宝宝边吃奶边哭怎么回事| 亚洲视频免费观看视频| 国产精品久久电影中文字幕 | 欧美激情 高清一区二区三区| 熟女少妇亚洲综合色aaa.| 无遮挡黄片免费观看| 精品熟女少妇八av免费久了| 国产无遮挡羞羞视频在线观看| 交换朋友夫妻互换小说| 一区二区三区乱码不卡18| av有码第一页| 久久青草综合色| 亚洲第一青青草原| 美女扒开内裤让男人捅视频| 波多野结衣av一区二区av| 国产精品秋霞免费鲁丝片| 飞空精品影院首页| 日本av免费视频播放| 亚洲三区欧美一区| 亚洲精品自拍成人| 手机成人av网站| 肉色欧美久久久久久久蜜桃| 久久 成人 亚洲| 久久久精品免费免费高清| 青草久久国产| 欧美一级毛片孕妇| 国产精品一区二区精品视频观看| 捣出白浆h1v1| av福利片在线| 日韩免费高清中文字幕av| 精品亚洲成国产av| 亚洲七黄色美女视频| 大片电影免费在线观看免费| 下体分泌物呈黄色| 亚洲久久久国产精品| av网站免费在线观看视频| 国产一区二区在线观看av| 国产精品国产av在线观看| 亚洲精品国产精品久久久不卡| 一级片免费观看大全| 日韩中文字幕视频在线看片| 少妇猛男粗大的猛烈进出视频| 老熟女久久久| 久久av网站| 亚洲精品自拍成人| 日韩中文字幕视频在线看片| 日日爽夜夜爽网站| 久久久久久人人人人人| 露出奶头的视频| 操美女的视频在线观看| 一本久久精品| 97人妻天天添夜夜摸| 国产欧美日韩精品亚洲av| 两人在一起打扑克的视频| 侵犯人妻中文字幕一二三四区| 美女高潮到喷水免费观看| 欧美日韩亚洲国产一区二区在线观看 | 精品熟女少妇八av免费久了| 欧美亚洲日本最大视频资源| 免费人妻精品一区二区三区视频| 午夜两性在线视频| 高清毛片免费观看视频网站 | 国产日韩一区二区三区精品不卡| 久久毛片免费看一区二区三区| www.熟女人妻精品国产| 亚洲成人国产一区在线观看| 久久影院123| 18禁国产床啪视频网站| 制服人妻中文乱码| 男女床上黄色一级片免费看| 亚洲 欧美一区二区三区| 19禁男女啪啪无遮挡网站| 99精品在免费线老司机午夜| 岛国毛片在线播放| 国产免费视频播放在线视频| 亚洲av电影在线进入| 在线观看舔阴道视频| 成在线人永久免费视频| 黄色毛片三级朝国网站| 久久久欧美国产精品| 精品国产乱码久久久久久小说| 色综合婷婷激情| 免费少妇av软件| 国产有黄有色有爽视频| 精品熟女少妇八av免费久了| 成人18禁高潮啪啪吃奶动态图| www日本在线高清视频| av天堂在线播放| 亚洲成av片中文字幕在线观看| 国产精品久久久久久人妻精品电影 | 国产成人av激情在线播放| videos熟女内射| 国产精品偷伦视频观看了| 亚洲成人手机| 久久天堂一区二区三区四区| 极品人妻少妇av视频| 成人特级黄色片久久久久久久 | 亚洲综合色网址| 人人妻,人人澡人人爽秒播| 日本撒尿小便嘘嘘汇集6| 波多野结衣av一区二区av| 欧美成人免费av一区二区三区 | 久久毛片免费看一区二区三区| 国产片内射在线| 午夜福利,免费看| 精品国内亚洲2022精品成人 | 国产成人欧美| 国产精品亚洲一级av第二区| 亚洲九九香蕉| av天堂久久9| 亚洲国产毛片av蜜桃av| bbb黄色大片| 亚洲一区中文字幕在线| 国产三级黄色录像| 亚洲avbb在线观看| 久久国产亚洲av麻豆专区| 一边摸一边抽搐一进一出视频| 中文字幕人妻丝袜一区二区| 欧美日韩福利视频一区二区| 日本av手机在线免费观看| av欧美777| 欧美精品一区二区免费开放| av免费在线观看网站| 母亲3免费完整高清在线观看| 侵犯人妻中文字幕一二三四区| 一区二区三区精品91| 999久久久精品免费观看国产| 欧美日韩一级在线毛片| 精品亚洲成国产av| 一级毛片电影观看| 91精品三级在线观看| 丰满迷人的少妇在线观看| 成年动漫av网址| 在线av久久热| 人人妻,人人澡人人爽秒播| 99精品久久久久人妻精品| 变态另类成人亚洲欧美熟女 | 性色av乱码一区二区三区2| 9色porny在线观看| 国产av一区二区精品久久| 99热国产这里只有精品6| 中文亚洲av片在线观看爽 | 91国产中文字幕| 色老头精品视频在线观看| 90打野战视频偷拍视频| 人人澡人人妻人| 十八禁高潮呻吟视频| 中亚洲国语对白在线视频| 在线观看舔阴道视频| 99在线人妻在线中文字幕 | 亚洲中文日韩欧美视频| 国产高清激情床上av| 日本vs欧美在线观看视频| 女同久久另类99精品国产91| 欧美激情极品国产一区二区三区| 亚洲精品乱久久久久久| 好男人电影高清在线观看| 最新的欧美精品一区二区| 99国产精品99久久久久| 高清视频免费观看一区二区| 日韩中文字幕视频在线看片| 婷婷成人精品国产| 亚洲五月婷婷丁香| 国产精品1区2区在线观看. | 久久久久久久久久久久大奶| 国产片内射在线| 狠狠精品人妻久久久久久综合| 国产高清国产精品国产三级| 蜜桃国产av成人99| 免费看十八禁软件| 色综合婷婷激情| 日韩熟女老妇一区二区性免费视频| 青青草视频在线视频观看| 国产真人三级小视频在线观看| 色精品久久人妻99蜜桃| 亚洲男人天堂网一区| 国产亚洲午夜精品一区二区久久| 在线观看一区二区三区激情| 亚洲五月婷婷丁香| 亚洲色图av天堂| 日韩欧美免费精品| 啦啦啦在线免费观看视频4| 美女午夜性视频免费| 女人久久www免费人成看片| 天天添夜夜摸| 一区二区av电影网| 国产高清视频在线播放一区| 欧美成狂野欧美在线观看| 欧美亚洲日本最大视频资源| 超色免费av| 天堂俺去俺来也www色官网| 亚洲五月婷婷丁香| 建设人人有责人人尽责人人享有的| 欧美精品人与动牲交sv欧美| 下体分泌物呈黄色| 黄频高清免费视频| 精品高清国产在线一区| 啦啦啦 在线观看视频| 正在播放国产对白刺激| 80岁老熟妇乱子伦牲交| 久久热在线av| 欧美激情极品国产一区二区三区| 精品高清国产在线一区| 午夜福利免费观看在线| 热re99久久国产66热| 成人免费观看视频高清| 色综合婷婷激情| 久久久国产精品麻豆| 夜夜骑夜夜射夜夜干| 日本a在线网址| 曰老女人黄片| 男女边摸边吃奶| 亚洲熟女毛片儿| 国产精品成人在线| 丰满饥渴人妻一区二区三| 欧美在线一区亚洲| 黄色视频在线播放观看不卡| 日韩欧美免费精品| 在线 av 中文字幕| 一级片免费观看大全| 成年版毛片免费区| 十八禁人妻一区二区| 中文亚洲av片在线观看爽 | 日本精品一区二区三区蜜桃| 人人澡人人妻人| 岛国毛片在线播放| 亚洲av片天天在线观看| 性色av乱码一区二区三区2| 美女国产高潮福利片在线看| 中文字幕人妻丝袜制服| 欧美 亚洲 国产 日韩一| 久久99一区二区三区| 老司机亚洲免费影院| 亚洲一码二码三码区别大吗| www日本在线高清视频| 满18在线观看网站| 亚洲欧美一区二区三区久久| 国产男女超爽视频在线观看| 黑人操中国人逼视频| 国产成人精品在线电影| 精品卡一卡二卡四卡免费| 19禁男女啪啪无遮挡网站| 久久99热这里只频精品6学生| 老汉色av国产亚洲站长工具| 日日爽夜夜爽网站| 国产一区二区激情短视频| 美女国产高潮福利片在线看| 精品国产乱码久久久久久小说| 日韩精品免费视频一区二区三区| 免费一级毛片在线播放高清视频 | 亚洲专区国产一区二区| 国产成人一区二区三区免费视频网站| 国产av一区二区精品久久| 丝袜美足系列| 免费看十八禁软件| 精品人妻1区二区| 亚洲情色 制服丝袜| 亚洲欧美一区二区三区黑人| 国产高清视频在线播放一区| 亚洲国产中文字幕在线视频| 日韩熟女老妇一区二区性免费视频| 成年版毛片免费区| 国产精品秋霞免费鲁丝片| 欧美老熟妇乱子伦牲交| 欧美 亚洲 国产 日韩一| 麻豆国产av国片精品| 日本欧美视频一区| 久久99热这里只频精品6学生| 69av精品久久久久久 | av网站在线播放免费| 成人免费观看视频高清| 国产视频一区二区在线看| 在线av久久热| 国产一区二区三区综合在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲欧美色中文字幕在线| xxxhd国产人妻xxx| 久久免费观看电影| 免费黄频网站在线观看国产| 啦啦啦中文免费视频观看日本| 久久久久久久久久久久大奶| 99久久人妻综合| 天堂8中文在线网| 电影成人av| 国产男靠女视频免费网站| 黄片播放在线免费| 人成视频在线观看免费观看| 韩国精品一区二区三区| 久久久久久久久免费视频了| 人成视频在线观看免费观看| 99国产极品粉嫩在线观看| 国产一区二区 视频在线| 狠狠狠狠99中文字幕| 日韩制服丝袜自拍偷拍| 精品一区二区三区av网在线观看 |