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

    ARTICLE High-Resolution Threshold Photoelectron Spectroscopy by Vacuum Ultraviolet Laser Velocity-Map-Imaging Method?

    2016-04-08 06:35:44ZhouLuHongGaoYuntaoXuLeiYangChowShingLamYaniceBenitezNgDepartmentofChemistryUniversityofCaliforniaDavisCA95616USADatedReceivedonDecember2015AcceptedonDecember152015
    CHINESE JOURNAL OF CHEMICAL PHYSICS 2016年1期

    Zhou Lu,Hong Gao,Yun-tao Xu,Lei Yang,Chow-Shing Lam,Yanice Benitez,C.Y.Ng?Department of Chemistry,University of California,Davis CA 95616,USA(Dated:Received on December 1,2015;Accepted on December 15,2015)

    ?

    ARTICLE High-Resolution Threshold Photoelectron Spectroscopy by Vacuum Ultraviolet Laser Velocity-Map-Imaging Method?

    Zhou Lu,Hong Gao,Yun-tao Xu,Lei Yang,Chow-Shing Lam,Yanice Benitez,C.Y.Ng?Department of Chemistry,University of California,Davis CA 95616,USA
    (Dated:Received on December 1,2015;Accepted on December 15,2015)

    Key words:Photoionization,Threshold photoelectron,Velocity-map imaging,Radical

    ?Part of the special issue for“the Chinese Chemical Society’s 14th National Chemical Dynamics Symposium”.

    ?Author to whom correspondence should be addressed.E-mail: cyng@ucdavis.edu

    I.INTRODUCTION

    Due to the space charge e ff ect,ions cannot be prepared in high concentrations.As a result,most of the conventional methods for spectroscopic measurements of neutral species cannot be directly applied to the study of ions.Spectroscopic properties of cations, such as rotational and vibrational constants and ionization energies(IEs),can be directly measured by highresolution photoelectron spectroscopy.Thus,progress and development made in photoionization and photoelectron spectroscopy can be considered as advancement in ion spectroscopy.

    Modern high-resolution photoelectron spectroscopic techniques that are based on the detection of zero kinetic energy(ZEKE)photoelectrons have their origin in the development of threshold photoelectron(TPE) measurements.The TPE spectroscopy,also previously referred to as ZEKE spectroscopy,was initiated about a half of century ago when tunable laboratory vacuum ultraviolet(VUV)photoionization sources based on monochromatized discharge radiation became available.By tuning the VUV photon energy across a photoionization transition band,the most popular approach for TPE detection involves the measurement of ZEKE photoelectrons using a dc extraction electric fi eld and a steradiancy analyzer,which essentially consists of a metal capillary to de fi ne a small solid angle for discriminating kinetic photoelectrons.Since ZEKE photoelectrons can be collected e ffi ciently with a small electric fi eld,this approach has been proven to give excellent TPE detection sensitivity.Nevertheless,the inability to discriminate against kinetic electrons traveling in the same solid angle de fi ned by the steradiancy analyzer[1], along with the relatively low optical resolutions o ff ered by tunable laboratory VUV discharge sources,has generally limited the TPE resolution to about 20?30 meV (full-width at half-maximum,FWHM).

    When a high-resolution pulsed laser is used as the photoionization source,the employment of a delayed electric fi eld pulse(with respect to the application of the photoionization laser)for TPE extraction is expected to be e ff ective for dispersing kinetic photoelectrons,and thus improve the TPE resolution.This idea of using a pulsed laser photoionization source together with the delayed pulsed electric fi eld photoelectron detection has led to the development of the pulsed fi eld ionization (PFI)-ZEKE or PFI-photoelectron(PFI-PE)spectros-copic technique[2?4].The VUV-PFI-PE scheme was later understood to originate from the detection of ZEKE photoelectrons formed by Stark electric fi eld ionization of high-n(n>100)Rydberg states populated by pulsed VUV laser excitation.Thus,the development of the PFI-PE detection method can be considered to derive from the TPE measurement scheme.Subsequent to the discovery of the laser PFI-PE technique,Zhu and Johnson introduced a variation of the PFI detection method,namely,the mass analyzed threshold ionization (MATI)detection scheme by detecting PFI-photoions (PFI-PIs)instead of PFI-PEs formed by PFI of high-n (n>100)Rydberg states[5].Due to the need to separate prompt background ions from MATIs or PFI-PIs,the PFI-PI measurements are more di ffi cult to perform and the achievable energy resolution is generally lower compared to that for PFI-PE measurements.Nevertheless, the PFI-PI detection provides information about the mass identity as well as the internal energy state of the photoion formed from the neutral precursor molecule, and has shown to allow the preparation of rovibronically selected molecular ions for unimolecular and bimolecular dynamics studies[6?8].This aspect of PFI studies is not in the scope of the present experiment.

    At the present time,the laser PFI-PE and PFI-PI detection techniques[9]are recognized as state-of-theart photoionization methods,which generally o ff er the highest energy resolution of 1?4 cm?1(FWHM).Although the VUV-PFI-PE and VUV-PFI-PI schemes are in principle applicable to all molecular systems,the application of these VUV laser PFI methods for the study of important molecular species,such as transient reaction intermediates and radicals,remains nontrivial.Since the PFI signal is derived from PFI of high-n (n>100)Rydberg species populated by photoexcitation (PEX)and the application of the electric fi eld pulse for PFI detection is usually delayed up to a fewμs with respect to the laser excitation,these high-n(n>100) Rydberg species can undergo decay induced by stray electric fi elds existed at the PEX/PI region.Without a careful e ff ort to minimize stray electric fi elds,the PFIPE and PFI-PI signals can be reduced,particularly in high-resolution PFI studies.This is the major disadvantage of PFI measurements.

    A highly e ff ective method for the separation and thus suppression of kinetic electrons in TPE detection is the velocity-map imaging(VMI)photoelectron(VMI-PE) scheme,which has been demonstrated previously by using a tunable laboratory VUV discharge photoionization source[10],achieving an energy resolution limited mostly by the optical bandwidth of 3?30 meV (FWHM)of the VUV discharge sources used.In addition to the merit of velocity-focusing of photoelectrons, the VMI-TPE and VMI-PE detection schemes also allow the collection of all low energy photoelectrons because a relatively high dc extraction electric fi eld and a large imaging microchannel plate(MCP)detector can be used in VMI-PE measurements.That is,the VMIPE and VMI-TPE detection schemes do not su ff er from the signal-decay problem as in PFI measurements,and thus result in enhanced sensitivity for VMI-PE measurements.

    The nature of VMI measurements is that the velocity resolution(?v)depends only on?r,and remains constant across the image(i.e.,for all the velocities),where r is the radial distance measured in term of the center of the VMI-PE image on the imaging detector.Thus, the achievable bandwidth(?E)in term of photoelectron√kinetic energy(E)is expected to be proportional toE and can be signi fi cantly narrowed by lowering E[11].This relation suggests that the highest photoelectron resolution or the narrowest?E value can be achieved by detecting TPEs as the VUV laser energy is scanned across the photoelectron transition bands of interest.This detection scheme is referred to the VUVVMI-TPE method.

    Both the VMI-PE measurement of neutral photoionization[12]and the slow electron velocity-map imaging (SEVI)measurement[13,14]of anion phot √odetachment take advantage of the character of?E∝E for VMIPE detection to achieve higher electron energy resolution.However,in many aspects,the VMI-PE and VMITPE measurements are di ff erent from the SEVI measurements.Photoionization of neutral species slightly below the ionization threshold is mediated by excitation to high-n Rydberg states.Due to the continuity of oscillator strength,the ionization threshold for neutral photoionization exhibits the step-function behavior[15].The measured photoionization thresholds for neutral species can be lowered by the applied electric fi eld in the photoionization/photoexcitation(PI/PEX) region,making it necessary to make Stark shift corrections,which is particularly important in high resolution photoelectron measurements.Anionic species have no Rydberg states;and the threshold law for anion photodetachment is proportional to E(l+1/2),where l is the angular momentum of the outgoing photoelectron.That is,the yield of photodetachment can be zero at true thresholds.The threshold law can also limit the photodetachment cross section to very low values at photon energies near the threshold.

    Preliminary reports on the successful implementation and application of the VUV-VMI-TPE and VUV-VMIPE detection methods for the high resolution photoelectron study of propargyl radical(C3H3)have been communicated[12,16].In this work,we present more detailed considerations and procedures on the application of these methods using C3H3and allyl(C3H5)radicals and their cationsas examples. The comparison between the VUV-VMI-TPE spectra observed here and the VUV-PFI-PE spectra[17]reported previously has illustrated the excellent performance in sensitivity and resolution of the VUV-VMIPE and VUV-VMI-TPE methods.Because the IE of chlorobenzene(C6H5Cl)(IE(C6H5Cl)=9.072 eV)has been well-established by VUV-PFI-PE measurements[18,19]and is close to the IE(C3H3)=8.698 eV[12,17], we have also chosen C6H5Cl to demonstrate the VUVVMI-PE and VUV-VMI-TPE methods.

    II.EXPERIMENTAL CONSIDERATION

    In the present experimental study,two apparatuses were used,namely,the VUV laser PFI-PE apparatus and the VUV laser VMI-PE and VMI-TPE apparatus.The experimental arrangements and procedures employed in using these apparatuses for VUV-PFI-PE detections have been described in detail;and those for VUV-VMI-PE and VUV-VMI-TPE measurements have been brie fl y described in Refs.[12,16].Thus,a more detailed description for the VUV-VMI-PE and VUVVMI-TPE measurements is given below in the present study.

    Both the PFI and VMI apparatuses employed in this work were equipped with a tunable VUV laser system,which was generated by resonance-enhanced fourwave di ff erence-frequency mixing(2ω1?ω2)schemes in a pulsed Kr or Xe gas jet as the nonlinear medium, where ω1and ω2represent the respective ultraviolet (UV)and visible laser frequencies generated by a UV and a visible dye laser,respectively.The UV and visible dye lasers were pumped by a common injection seeded Nd:YAG laser(Spectra-physics,GCR-290)operated at 30 Hz.The UV ω1frequency was fi xed at 212.556 nm (222.560 or 249.629 nm)to match the two-photon(2ω1) resonances of the Kr(Xe)transitions.The visible ω2output was scanned in the range of 485?760 nm to generate the desired VUV(2ω1?ω2)output range of 8.1?9.4 eV as required for the present experiment[16, 20,21].Here,the tunable VUV(2ω1?ω2)output was selected by an o ff-axis plano-convex MgF2lens before entering the PEX/PI region;and its intensity was monitored by a Cu photoelectric detector.The optical bandwidth of the VUV laser output was measured to be 0.12 cm?1(FWHM)for the VUV-PFI-PE apparatus [20,21]and 0.45 cm?1(FWHM)for the VUV-VMI-PE and VUV-VMI-TPE apparatus[22?26].

    A.The VUV-PFI-PE apparatus

    As described above,the VUV-PFI-PE apparatus consists of a pulsed(30 Hz)tunable VUV laser system as the photoionization source,a pulsed supersonic beam production system to introduce the gaseous sample,a time-of- fl ight(TOF)mass spectrometer for photoion detection,and a photoelectron spectrometer for PFIPE detection.The TOF ion mass spectrometer and the TOF electron spectrometer are situated above and below the PI/PEX center,such that the central axes of the TOF ion mass spectrometer and TOF electron spectrometer are perpendicular to the traveling direction of the molecular beam.

    Using the VUV-PFI-PE apparatus along with procedures described previously[20,21,27?29],we have measured the VUV-PFI-PE spectra of C6H5Cl near its ionization threshold.In this experiment,the C6H5Cl vapor pressure at room temperature(11.8 Torr)was seeded in He to a total stagnation pressure of 30 psi prior to supersonic expansion through a pulse valve (General valve,repetition rate=30 Hz).The resulting C6H5Cl pulsed beam was skimmed by a conical skimmer(diameter=2 mm)before intersecting perpendicularly with the VUV(2ω1?ω2)laser beam in the PI/PEX center.The timings of operating the pulsed valve and the VUV photoionization laser were synchronized by delayed pulse generators(Stanford Research System,DG535).

    B.The VUV-VMI-PE and VUV-VMI-TPE apparatus

    The VUV-VMI-PE and VUV-VMI-TPE apparatus was modi fi ed from the VUV laser VMI-photoion(VMIPI)apparatus,which has been described in detail[22, 23,30].Since electrons are susceptible to perturbation by the Earth and stray magnetic fi elds,we have shielded the electron-imaging lenses and the electron TOF tube of the apparatus by two layers ofμmetals.The VUVVMI-PE and VUV-VMI-TPE apparatus,in sequential order,consists of a pulsed molecular beam production system,a set of VMI aperture lenses,a 75-cm electron TOF tube,and a VMI-PE detector.The VMI-PE detector comprises of a dual set of electron-imaging MCPs (diameter=75 mm),a P47 phosphor screen,and a CCD camera.

    In the present experiment,the gas sample was seeded in He or Ar to a total stagnation pressure of 30 psi prior to supersonic expansion into the beam source chamber through a pulsed valve(Evan-Lavie Model: EL-5-2004,nozzle diameter=0.2 mm,30 Hz).The gas sample beam thus produced traveled along the central axes of the VMI lenses and the VMI-PE detector,and was skimmed by two conical skimmers(diameter=1 mm)before entering the PI/PEX region to intersect perpendicularly the VUV(2ω1?ω2)laser beam. The photoelectrons thus formed were velocity-focused and mapped onto the electron-imaging MCP detector [31].When collecting photoelectron images,the front MCP plate was grounded at all time,and a dc voltage of around+1.5 kV was applied to the second MCP plate.The amplitude of this dc voltage applied can alter the ampli fi cation of the photoelectron signal.Thus, the actual amplitude was adjusted depending on the signal level.A high voltage was applied to the P47 phosphor and the resulting fl uorescence was captured by a CCD camera and transferred to a computer using the DAVIS7 software provided by LaVision.The images recorded represent 3-dimensionsal(3D)Newtonian spheres mapped onto the 2D detector plane.The 3D distribution of photoelectrons can be reconstructedfrom the recorded 2D VMI-PE image by the inverse Abel transformation[32];and thus the kinetic energy distribution of photoelectrons can be obtained from the reconstructed 3D VMI-PE images.

    We have measured the VUV-VMI-PE and VUVVMI-TPE spectra of C6H5Cl,C3H3and C3H5near their ionization thresholds.A supersonically cooled radical beam source based on 193.3 nm laser photodissociation of propargyl chloride(C3H3Cl)and allyl bromide(C3H5Br)was used to prepare the C3H3and C3H5radicals,respectively[12,16,20,33].The C3H3Cl (C3H5Br)precursor sample was introduced into the beam source chamber as a gas jet of~10%C3H3Cl (C3H5Br)seeded in He.The 193 nm ArF excimer laser pulse(GAM laser,pulse energy≈10 mJ,repetition rate=30 Hz)intersects the precursor sample jet at the tip of the nozzle,inducing the formation of C3H3(C3H5)radicals from photodissocation of C3H3Cl (C3H5Br).It has been demonstrated previously that radicals,such as C3H3and C3H5formed in the high pressure region of the pulsed jet can undergo supersonic cooling by collisions with He or Ar atoms in the jet. The supersonically cooled C3H3(C3H5)radical beam passed through two conical skimmers(diameter=1 mm) before intersecting the tunable VUV laser radiation at the PI/PEX region of the VUV-VMI apparatus.

    The imaging MCP detector used has an active area with a diameter of 75 mm,and the distance between the PI/PEX center and the center of the MCP detector is 747 mm.The photoelectron signal arriving at the MCP is monitored by the 960 pixel×960 pixel area.For VMI-TPE measurements of C3H3,C3H5,and C6H5Cl, the smallest central image area of 0.5 mm×0.5 mm (3 pixels×3 pixels)was used to gate the VUV-VMITPE counts.Thus,the critical solid angle that characterizes the steradiancy TPE detector is very small, corresponding to a very high TPE energy resolution. The optimal resolution of the VMI-TPE technique can be achieved by lowering the electric fi eld applied at the PI/PEX region and reducing the pixel area at the center of the VMI-PE image during the scan of VUV laser energy.The VUV-VMI-TPE spectra obtained by scanning the VUV energy are normalized by the VUV intensities recorded by a Cu photoelectric detector.

    III.RESULTS AND DISCUSSIONS

    A main motivation of the present study of C6H5Cl, C3H3,and C3H5and their cations is to establish proper procedures for VUV-VMI-PE and VUV-VMITPE measurements as a general method for highresolution photoelectron spectroscopic studies.

    A.VUV-PFI-PE,VUV-VMI-PE,and VUV-VMI-TPE measurements for

    FIG.1(a)VUV-VMI-PEimagerecordedathν (VUV)=9.240 eV and F=33.2 V/cm.(b)VUV-VMI-PE spectrum for the formation of C6H5Cl+(X2B1)converted from the VMI-PE image shown in(a).

    The study of a stable molecule such as C6H5Cl is aimed to illustrate the important consideration and procedures involved in the application of the VUV-VMI-PE and VUV-VMI-TPE methods.As an example,we show in Fig.1(a)the VUV-VMI-PE image for the formation ofobtained at hν(VUV)=9.240 eV and F=33.2 V/cm.On the basis of inverse Abel transformation and the consideration of VUV photon energy,this photoelectron image was converted to the VUV-VMIPE spectrum foras shown in Fig.1(b). The photoelectron bands ofresolved in Fig.1(b)correspond to the VMI-PE rings observed in Fig.1(a).

    TheVUV-VMI-PEspectrumobservedat hν(VUV)=9.240eViscomparedtoVUV-VMIPE spectra formeasured at other hν(VUV)values of 9.110,9.150,9.185,9.200,9.220, 9.290,9.340,9.377,9.425,and 9.450 eV in Fig.2. The downward pointing red arrows marked the actual hν(VUV)energies used for individual VUV-VMI-PE measurements.The VMI-PE spectra reveal the excitation of two major vibrational progressions,nv6aand v7a+nv6a,n=0?2,as marked on top of the spectra in Fig.2.The present measurement yields 420±5 and 1120±5 cm?1for the v6aand v7avibrational frequencies,respectively,which are in good agreement with the previous(1+1′)PFI-PE measurements[19].

    FIG.2The VUV-VMI-PE spectra for C6H5Cl+(X2B1) converted from VUV-VMI-PE images measured at the hν(VUV)=9.110,9.150,9.185,9.200,9.220,9.240,9.290, 9.340,9.377,9.425,and 9.450 eV and F=33.2 V/cm. The red downward pointing arrows represent the actual hν(VUV)values used for individual VUV-VMI-PE measurements.The top energy scale is the VUV energy in eV;and the bottom energy scale is in cm?1measured with respect to the peak of C6H5Cl+origin band.The vibrational assignments are marked by red droplines on top of the fi gure. The excitation of two major vibrational progressions,nv6aand v7a+nv6a,n=0?2,are observed as marked on top of the fi gure.

    The comparison of Fig.2 shows that as hν(VUV)is decreased,the number of vibrational bands observable in the VUV-VMI-PE spectrum is reduced because of the reduced photoelectron kinetic energy span,which is determined by the di ff erence between the hν(VUV)employed and the origin band position of C6H5Cl+(X2B1). As pointed out above concerning the application of the VMI-PE method,the kinetic energy of photoelectrons ejected from a vibrational band is proportional to the square of photoelectron velocity,and thus to r2,where r is the radius of the VMI-PE ring observed on the imaging MCP detector.For a given hν(VUV),the outermost VMI-PE ring(radius=rmax)observed in this experiment correspond to the formation of the origin band of C6H5Cl+(X2B1).Thus,we expect that hν(VUV)used is proportional to rm2ax.This expectation is con fi rmed by the excellent linear plot of hν(VUV)versus rm2axplot observed in Fig.3.This linear plot yields an intercept of 9.070 eV at rmax=0.This value can be taken as the ionization threshold of C6H5Cl at a dc electric fi eld F=33 V/cm.To determine the IE(C6H5Cl)value from this intercept based on the VUV-VMI-PE measurements of Fig.2 requires the Stark shift correction. S tein rmce in t e hde Ib Ey ( Ct h 6eH 5VCUl )V=-7P 3F1I7-2P. 1E± m1.e6a csmur ?em 1h ena ts b [1e 8 e,n 1 d9 e]-, a simple energy calibration of the VUV-VMI-PE spectra would involve normalizing the peak position of the VUV-VMI-PE origin band to the latter IE(C6H5Cl) value.

    The comparison of the VUV-VMI-PE spectra in Fig.2 also reveals that as the photoelectron kinetic energy E for a selected vibrational band is decreased due to the lowering of the hν(VUV)value,the FWHMs

    0 40000120000 80000

    rmax/ Arb. unit

    9.10

    9.20

    9.30

    9.40

    9.50

    Photonenergy / eV

    2 FIG.3 The plot of hν(VUV)versus r2max,where rmaxin arbitrary units represents the radius of the outmost ring of the VMI-PE image for the C6H5Cl+(X2B1)origin band measured at selected hν(VUV)values of 9.110,9.150,9.185, 9.200,9.220,9.240,9.290,9.340,9.377,9.425,and 9.450 eV and F=33.2 V/cm.A linear plot is observed,yielding an intercept of 9.070 eV at rmax=0. for the vibrational bands appearing in the spectra become narrower,indicative of higher energy resolutions at lower E values.This observation is as pointed out above.The FWHM for the origin band obtained hν(VUV)=9.110 eV is 25 cm?1,and is the lowest value observed in Fig.2.Spectral simulation indicated that this FWHM value represents a limit set by overlapping rotational transitions and the rotational temperature achieved in supersonic expansion for the C6H5Cl molecular beam.

    The relation of?E∝√E observed in the VUV-VMIPE measurement of C6H5Cl indicates that the highest resolution can be achieved by using the VUV-VMI-TPE method.This method involves scanning the hν(VUV) energy across the vibrational band of interest and gating the TPEs arriving at the center of the VMI-PE detector.We have recorded the VUV-VMI-TPE spectrum for C6H5Cl+(X2B1)in the hν(VUV)range of 9.05?9.22 eV,and this spectrum is compared in Fig.4 with the VUV-PFI-PE spectrum for C6H5Cl+(X2B1) measured in the range of 9.05?9.42 eV.The VUV-PFIPE spectrum has been corrected for the Stark shift effect,and thus the IE(C6H5Cl)is marked by the peak of the VUV-PFI-PE origin band.Due to the Stark shift e ff ect,the vibrational peaks observed in the VUVVMI-TPE spectrum are red shifted with respect to the corresponding vibrational peak resolved in the VUVPFI-PE spectrum.The assignment of the vibrational bands is marked on top of the spectra in the fi gure.In addition to the strong v6aand v7avibrational bands, weak vibrational bands attributable to v1and v12vibrational excitations are also observed in these spectra.Previous simulation of the VUV-PFI-PE origin band indicates that the PFI-PE resolution achieved was 1.5?2.0 cm?1(FWHM).The essentially identical structures and FWHMs(≈25 cm?1)for the vibrational

    9.069.129.189.249.309.369.42

    Relativeintensity/A

    rb. unitυ1

    υ6a

    υ12υ7a+n*υ6a

    VUV energy / eV

    FIG.4(a)Comparison of the VUV-VMI-TPE spectrum (upper curve)with the VUV-PFI-PE spectrum(lower curve)for C6H5Cl+measured in the hν(VUV)range of 9.05?9.42 eV and F=33.2 V/cm.The VUV-VMI-PE spectrum has not been corrected for the Stark shift e ff ect, whereas the Stark shift for the VUV-PFI-PE spectrum has been corrected by normalizing the VUV-PFI-PE origin band of C6H5Cl+to the IE(C6H5Cl).The vibrational assignments are marked by red droplines on top of the spectra. bands observed in the VUV-VMI-TPE and the VUVPFI-PE spectra of Fig.4 support the conclusion that the photoelectron energy resolution attained in the VUVVMI-TPE measurement is comparable to that observed in the VUV-PFI-PE detection.

    Although the FWHMs of vibrational peaks resolved in the VUV-VMI-TPE and VUV-PFI-PE spectra of Fig.4 are about the same,the energy pro fi les observed for the VUV-PFI-PE peaks appear to be different from those of the VUV-VMI-TPE peaks.In order to examine the di ff erent energy pro fi le of the VUV-PFI-PE and VUV-VMI-TPE vibrational band for,we show in Fig.5 the magni fi ed VUV-PFI-PE and the VUV-VMI-TPE spectra for the origin band of.The IE(C6H5Cl)is marked by the downward pointing blue arrow or the VUV-PFIPE peak as the VUV-PFI-PE spectrum has been corrected for the Stark shift due to the pulsed electric f i eld applied in the PFI-PE measurement.The di ff erent spectral characteristics observed in the VUV-VMITPE and VUV-PFI-PE measurements have been brie fl y discussed previously[16].As shown in Fig.5,both of these spectra are asymmetric.The VUV-PFI-PE origin band forshows higher intensities on the low energy side of the origin band,which have been attributed to forced autoionization in VUV-PFIPE measurements[4].The VUV-VMI-TPE spectrum reveals higher intensities on the high energy side due to contribution from hot photoelectrons,resulting in a tailing structure toward higher energies.The discrimination of hot photoelectrons in the present detection of TPEs or ZEKE photoelectrons is achieved by the VMIPE arrangement,such that the TPEs are expected to

    FIG.5 A magni fi ed plot of the VUV-VMI-TPE spectrum (red curve)and the VUV-PFI-PE spectru+m(blue curve) for the origin vibrational band of C6H5Clin the energy rhaansg bee oefn 7 c3o0r5re0c?te7d32 f5o0r ctmhe? 1S.taTrkhe s hVifUtV e -ffP eFctI.-P TEh esp veecrttriucmaldownward pointing arrow,which passes through the peak of the VUV-PFI? -PE spectrum marks the IE(C6H5Cl).The peak position of the VUV-VMI-PE spectrum is 17 cm?1 lower than that of the VUV-PFI-PE spectrum,attributing to the Stark shift induced by the dc fi eld of F=33.2 V/cm used in the VUV-VMI-PE measurement.The dashed line gives the estimated pro fi le of the VMI-TPE spectrum.The intensities above the dashed line can be attributed as the hot photoelectron background[16].

    7305073100731507320073250 VUV / cm-1

    Electronintensity/A

    rb. unitarrive at the center of the imaging MCP detector.

    B.DC Stark shift in VMI-TPE measurements for C6H5Cl+(X2B1)

    A dc electric fi eld F=33.5 V/cm was applied at the PI/PEX region in the VMI-TPE measurement of theorigin band shown in Fig.4 and Fig.5. In an ideal situation of achieving completed discrimination of hot photoelectrons,the transmission for true TPEs is zero at the VUV energies above IE(C6H5Cl). Furthermore,we expect the energy pro fi le for VMI-PE detection is a symmetric function.We have estimated the energy pro fi le for the true VUV-VMI-TPE detection by drawing the dashed line in Fig.5 to separate the hot photoelectron contribution.That is,the signal above the dashed line are attributed to hot electrons travelling perpendicular to the MCP detector and arriving at the center of the VMI-PE image.The peak position of theorigin band observed in the VUV-VMITPE spectrum is found to be 17 cm?1lower than that found in the VUV-PFI-PE spectrum.This fi nding provides an estimate that the S√tark shift correction from the true IE(C6H5Cl)is?3.1F in cm?1,which is h√alf of the value predicated by classical formula of?6.1F in cm?1.

    FIG.6The VUV-VMI-TPE spectra of the origin vibrationalbandofC6H5Cl+intheenergyregionof 73075?73225 cm?1measuredatthedcelectric fi eld F=33.2,66.2,98.4,130.5,and 162.0 V/cm.The vertical blue arrow marks the IE(C6H5Cl).The dashed curves are the Gaussian function fi ts to individual VUV-VMI-TPE bands.The intensities above the dashed curves at the high energy side of the VUV-VMI-TPE peak can be attributed to contribution of hot photoelectron background.

    The PFI-PE method is based on the detection of ZEKE photoelectron originating from delayed PFI of high-n(n>100)Rydberg states by a delayed,PFI fi eld after laser excitation.It is well-established that the lowering of the ionization threshold in delayed P√FI detection is governed by the diabatic formula?4F in cm?1[34].To our knowledge,a detailed examination of the electric fi eld e ff ects on the VMI-PE and VMITPE methods has not been reported.We have performed a series of VUV-VMI-TPE measurements for theorigin band at the dc electric fi elds of F=33.2,66.2,98.4,130.5,and 162.0 V/cm.These spectra obtained without the Stark shift correction are depicted in Fig.6.The vertical line shown in the fi gure marks the IE(C6H5Cl)value determined in previous VUV-PFI-PE measurements.As shown in the fi gure, the increase of the dc electric fi eld F results in shifting of the VMI-TPE origin band to a lower hν(VUV) value together with broadening of the VMI-TPE origin band.All the VMI-TPE peaks observed in Fig.6 have a similar intensity pro fi le,i.e.,it drops rapidly on the low energy side and exhibits a tailing structure on the high energy side.We fi nd that the VMI-TPE peaks observed can be simulated well by a Gaussian function except the tailing structure toward higher energies as shown by the dashed curves of Fig.6.The intensities above the dashed curves can be attributed to hot photoelectrons.While the broader VMI-PE peaks observed at a higher F fi eld correspond to higher photoelectron intensities and thus better experimental sensitivity,the energy resolution observed is lower at a higher F value.

    FIG.7 Plot of the peak position of the VUV-VMI-TPE spectrum for the C6H5Cl+origin band versus the square root of the dc electric fi eld F used for photoelectron extraction.The solid line represents a linear least squares fi t to the plot,giving a slope of?3.10±0.01 cm?1/(V/cm?1)1/2and an intercept of 73171.83±0.13 cm?1at F=0.

    According to the Stark shift observed in the VUVVMI-TPE spectrum of Fig.5,the Stark shift induced by a F fi eld can be estimated by the formula, in cm?1.For the dc electric fi elds of F=33.2,66.2, 98.4,130.5,and 162.0 V/cm are predicted to be 35, 50,61,70,and 78 cm?1,respectively.These Stark shift values were found to be consistent with the VUVVMI-TPE peak positions observed in Fig.6.The Gaussian fi ts to individual VMI-PE peaks for F=33.2,66.2, 98.4,130.5,and 162.0 V/cm are shown as the dashed red curves in Fig.6.The center positions of the best fi t Gaussian curves are also in agreement with the peak positions of the V√MI-TPE spectra.The center peak positions of theF in Fig.7.The solid line is the linear least squares fi t to the data with a slope of?3.10±0.01 cm?1(V/cm)?1/2and an intercept of 7 iz3 a1t 7 io1 n . 8t3 h±r0es.1h3olcdm m?e1a.s u Sri n edce a tt h F e = l a0t,t e itr c vaanl u b ee i si d tehnet i ifio endas the IE(C6H5Cl).Thus,this least square analysis provided a more precise IE value compared to the values determined in previous VUV-PFI-PE measurements.Furthermore,the least squares fi t of F√ig.7 also provides a more precise formula,(?3.10±0.01)F in cm?1,for the Stark shift correction for the present VUV-VMI-TPE measurements.

    C.VUV-VMI-TPE measurements of C3H3and C3H5

    The high resolution photoelectron studies of C3H3and C3H5radicals have been made in previous VUVPFI-PE measurements[17,35,36].Recently,we have shown that high resolution studies of transient radicals, such as C3H3,can also be performed using the VUVVMI-PE and VUV-VMI-TPE methods[12,16].The present study provides further examinations of the dc electric fi eld e ff ect on the photoionization of C3H3and C3H5radicals.

    FIG.8 The VUV-VMI-TPE spectra for the origin vibrational band of C3H3+measured at the dc electric fi elds of F=3.2,16.8,26.8,and 33.5 V/cm.The C3H3sample was prepared by a pulsed supersonically cooled radical beam source based on 193 nm photodissociation of C3H3Cl.All VUV-VMI-PE spectra have been co√rrected for the dc Stark shift according to the formula?3.1F cm?1.

    Figure 8 depicts the VUV-VMI-TPE spectra for the origin band ofobtained at the applied dc electric fi elds of F=3.2,8.4,16.8,26.8,and 33.5 V/cm.The comparison of these spectra reveals that the FWHM for the VUV-VMI-TPE origin band increases as F is increased.All these spectra have been corrected for the dc Star√k shifts according to the formula derived above(?3.1F cm?1).At F=33.5 V/cm, the VMI-TPE spectrum reveals a left shoulder structure on the lower energy of the main peak,while a right shoulder structure is not noticeable.As the dc electric fi eld is decreased to 3.2 V/cm(the bottom trace in Fig.8),the right shoulder structure becomes discernible.The VUV-VMI-TPE spectrum obtained at F=3.2 V/cm is in excellent agreement with that reported recently by Gao et al.after taking into account the experimental uncertainties.The simulation (to be presented below)of the VUV-VMI-TPE spectrum obtained at F=3.2 V/cm yielded an IE(C3H3) (fr7o0m17 4th±e2 r cemce?n1t) V,wUhVic-hV Mis Iid-TenPtEic aml etoas tuhree mvaelnute bdyer iGveado et al.[16]and that from the VUV-PFI-PE measurements by Jacovella et al.[17]. Figur

    e 9 compares the VUV-PFI-PE spectra for the C3H5+(X1A1)origin band measured by Xi et al.[36] (top curve)and Gasser et al.[35](middle curve)with that recorded by the present VUV-VMI-TPE measurement(bottom curve).All three experiments used a pulsed supersonic molecular beam source to introduce the radical sample into the PI/PEX center.The VUVPFI-PE experiment of Xi et al.used a radical beam source based on pyrolysis,and thus the rotational temperature achieved was higher and thus the FWHM of the origin band is larger compared to that attained in the experiment of Grasser et al.and the present VUVVMI-TPE experiment,which used a supersonic radical beam source based on 193 nm photodissociation.The dc electric fi eld used was F=3.2 V/cm;and the VUVVMI-TPE spectrum of Fig.9 has been corrected for the Stark shift.The double peak structure resolved in the VUV-PFI-PE spectra has been shown to arise from partially resolved rotational contours based on spectral simulations to be presented below.

    FIG.9 Comparison of the VUV-MI-PE spectrum(bottom) for the C3H5+origin band with the VUV-PFI-PE spectra for the C3H5+origin band obtained by Xi et al.(top)and Gasser et al.[35](middle).The VUV-VMI-TPE spectrum was obtained by setting the dc electric fi eld F=3.2 V/cm. The VUV-VMI-TPE and the VUV-PFI-PE spectra have been corrected for the Stark shift e ff ect.

    D.Simulation of the origin bands for C3H3+(X1A1)and C3H5+(X1A1)

    The observation of partially resolved rotational contours in the VUV-VMI-TPE spectra allows rotational simulations of the VUV-VMI-TPE origin bands ofandas shown in Fig.10 and Fig.11,respectively.Although the semi-empirical simulation presented below cannot be considered as de fi nitive,it is useful for providing a more precise IE value and a more reliable spectral assignments.The semi-empirical simulation procedures have been discussed in detail previously,and thus only a brief description of the selection rules for the molecular symmetry group C2v(M)is given below[16,17,35?39].

    Since both the neutral and cationicandsystems belong to the molecular symmetry group C2v(M),the asymmetric top rotational energy levels can be labeled as,whereis the total angular momentum quantum number andare the projection quantum numbers of the neutral(cation).By using the previous reported rotational constants A′′,B′′,and C′′and[16,35,36],the Ray’s asymmetry parameters,de fi ned as(2B′′-A′′-C′′)/(A′′-C′′)for the neutral and

    FIG.10ComparisonoftheVUV-VMI-TPEspectrum for the origin band of C3H3+in the region of 70050?70300 cm?1with the best simulated spectrum, which is the sum of all rotational branch contributions (shown below the best simulated spectrum).The simulation is based on the selection rules that?Ka=odd and ?Kc=even and odd.The simulation uses a Gaussian instrumental pro fi le(FWHM=3.5 cm?1),assumes a rotational temperature of 35 K,and the rotational branch ratios of N:O:P:Q:R:S:T=2:3:5:3:5:3:2.

    C+)/(A+-C+)for the cation,yield the value of?0.998(?0.998)and?0.926(?0.905)for ofand).Since the asymmetric parameters are very close to?1,bothand C3H5can be considered as prolate symmetric top molecules.The rotational energy levels(Erot)of C3H3andare calculated by diagonzalizing the Hamiltonian matrix using the symmetric top rotational basis set.The intensityfor a rotational transition from a neutral state to a cationic level is proportional to the rotational population of the neutral molecule according to the Boltzmann distribution, where Trotrepresents the rotational temperature for the molecular sample,?Erotis the rotational energy measured with respect to the ground rovibronic state of the neutral molecule,and rBis a scaling parameter for a given rotational branch.In the electric dipole approximation for the selection rule,the conservation of the total angular momentum ignoring the electron spin gives:|?J|=|J+?J′′≤l+1,where l represents the orbital angular momentum quantum number of the outgoing photoelectron.The simulation takes into consideration of the rotational branches,N,O,P,Q,R,S,and T for ?J=0,±1,±2,and±3.

    Thegeneralsymmetryselectionrulesforhasbeendiscussedindetailpreviously[16,17].The symmetry analysis of the rovibronic wavefunctions for the neutral and cation shows that the rotational selection rules for the photoionization

    was recorded at F=3.2 V/cm, achieving a FWHM of 18 cm?1.The simulation is based on the selection rules that?Ka=odd,?Kc=odd for l=odd, and?Ka=even,?Kc=even for l=even,where l represents the orbital angular momentum quantum number of the outgoing photoelectron.The best simulated spectrum was obtained b assumes a Gaussian instrumental pro fi le (FWHM=1.0 cm?1),a rotational temperature of 20 K,and rotational branch ratios of N:O:P:Q:R:S:T=3:4:7:6:7:4:3.

    Figure 10 reveals the best overall simulated spectrum (second spectrum from the top)in comparison with the experimental VMI-TPE spectrum(top spectrum)obtained at F=3.2 V/cm.Since two H atoms are located symmetrically about the main molecular axis of C3H3, it is necessary for the neutral rotational population to take into account the nuclear spin statistics,yielding a ratio of 3:1 for even.The best simulated spectrum for the formation ofis the sum of the rotational contributions of the Q;P and R;O and S;and N and T branches shown below the best overall simulated spectrum in Fig.10.The corresponding rotational branching ratios that give the best simulation are N:O:P:Q:R:S:T=2:3:5:3:5:3:2.The simulation also provide a rotational temperature of Trot=35 K for the C3H3radical sample.The experimental VUVVMI-TPE spectrum was simulated by assuming a Gaussian instrumental energy function(FWHM=3.5 cm?1). The experimental spectrum reveals a main peak with the distinctive left and right shoulders.The main peak of theorigin band is contributed by the photoionization transitionwhich lies 9 cm?1higher than the IE(C3H3)de fi ned the photoionization transitionThe simulation yields the IE(C3H3)value of 70169 cm?1.Since F=3.2 V/cmis applied in the PI/PEX region during the TPE detection,the Stark shift correction is determined to be 5.4 cm?1.Thus,the IE(C3H3)is determined to be 70174.4±3.5 cm?1,which is in excellent agreement with the published IE values of 70175±2 and 70174.5±2.0 cm?1determined by Hong et al.[16]and Jacovella et al.[17]respectively.

    Figure 11 depicts the best overall simulated spectrum (second spectra from the top)for the comparison with the experimental VMI-TPE spectrum(top spectrum) obtained at F=3.2 V/cm.The best simulated spectrum represents the sum of the contributions of the N,O,P,Q,R,S,and T-branches with the branching ratios of N:O:P:Q:R:S:T=3:4:7:6:7:4:3.The simulation yielded an estimate of 20 K for the rotational temperature Trotof C3H5radicals.The experimental VUV-VMI-TPE spectrum was simulated by a Gaussian instrumental energy function(FWHM=1 cm?1).The experimental VUV-VMI-TPE spectrum for the origin band oreveals a discernible double peak structure,which is also observed in the previous PFIPE measurements[35,36].The double peak structure (separation≈5 cm?1)observed in the VUV-VMI-TPE curve is less pronounced compared to that found in the previous VUV-PFI-PE,and is mainly contributed by the P and R rotational branches.The simulation yields an ionization threshold of 65580 cm?1for C3H5.After taking into account the Stark shift of 5.4 cm?1due to the applied of F fi eld of 3.2 V/cm at the PI/PEX region,the IE(C3H3)is determined to be 65585±2 cm?1. This latter value is also in good accord with the previous PFI-PE measurements of 65584.6±2 cm?1determined by Xing et al.[36]and 65580.1±2.0 cm?1obtained by Gasser et al.[35].

    E.VUV-VMI-PE spectra for

    The VUV-VMI-PE and VUV-VMI-TPE detection schemes are complimentary methods,which provide high detection sensitivity as well as high energy resolution.Although the VUV-VMI-TPE method can achieve higher energy resolution,it requires the continuous scanning over a wide VUV energy range.The VUV-VMI-PE technique has the capability of detecting photoelectrons over a broad energy range with good energy resolutions,and can avoid the time consuming and challenging task of tuning the VUV laser over a wide energy range.

    FIG.12 The VMI-PE spectrum for C3H5+converted from VUV-VMI-PE images recorded at the respective VUV energies of 8.192,8.284,and 8.323 eV(marked by downward pointing arrows).The vibrational assignments are marked by drop lines.The top scale is VUV energy in eV and the bottom energy scale is measured with respect to the origin band of C3H5+.

    The application of the VUV-VMI-PE method on radical species has been demonstrated in the previous study of;and the VUV-VMI-PE study ofhas not been made.Figure 12 depicts the VUV-VMI-PE spectra forrecorded at F=33.5 V/cm and hν(VUV)=8.192,8.284,and 8.323 eV.The photoelectron kinetic energiesof the VMI-PE spectrum were calibrated by setting the peak position of the VMI-PE origin band ofto equal the IE(C3H5) value of 65585±2 cm?1.The VUV-VMI-PE spectra reveal discernible excitation of the v7+(C?C?C bending) vibrational progression at 436,741,876,and 1309 cm?1 above theorigin band[35,36].Similar to the comparisonof the VUV-VMI-PE spectra forshown in Fig.2,the VUV-VMI-PE spectra forof Fig.12 show that the VMI-PE resolution improves as the hν(VUV)is decreased closer to the ionization threshold of C3H5,resulting in narrower FWHMs of the vibrational bands for[11].As an example,the FWHM of thevibrational band was found to decrease from 88 cm?1to 66 cm?1 to 24 cm?1as the hν(VUV)energy is decreased from 8.323 eV to 8.284 eV to 8.192 eV.

    IV.CONCLUSION

    We have recorded the VUV-VMI-TPE and VUVVMI-PE spectra of,,andnear their ionization thresholds,aiming to establish the experimental procedures for general high-resolution photoelectron measurements of gaseous molecules.The comparison of these measurements with the corresponding VUV-PFI-PE measurementsshows that the energy resolutions achieved in the VUVVMI-TPE and VUV-PFI-PE detection are comparable. Since a relatively large dc fi eld can be used to extract all photoelectrons produced,we expect the VMI-PE method to have higher detection sensitivity than that observed in the PFI-PE method.The VUV-VMI-TPE spectra of C6H5Cl,C3H3,and C3H5have also been examined as a function of the applied dc electric fi eld F employed for photoelectron extraction.The analysis of the observed Stark shift for the VU√V-VMI-TPE measurements yields a Stark shift of?3.1F in cm?1, whichis half the value predicted by classical formula in cm?1.

    V.ACKNOWLEDGMENTS

    This work was supported by the National Science Foundation under CHE-0910488 and CHE-1462172.C. Y.Ng also acknowledges the support by the Chemical Sciences,Geosciences and Biosciences Division,Offi ce of Basic Energy Sciences,O ffi ce of Science,(US) Department of Energy(DOE)under Contract No.DEFG02-02ER15306.

    [1]T.Baer,W.B.Peatman,and E.Schlag,Chem.Phys. Lett.4,243(1969).

    [2]K.Mller-Dethlefs,M.Sander,and E.W.Schlag,Chem. Phys.Lett.112,291(1984).

    [3]K.Muller-Dethlefs and E.W.Schlag,Annual Rev. Phys.Chem.42,109(1991).

    [4]E.W.Schlag,ZEKE Spectroscopy,Edition Anglaise, Cambridge:Cambridge University Press,(1998).

    [5]L.Zhu and P.Johnson,J.Chem.Phys.94,5769(1991).

    [6]Y.C.Chang,H.Xu,Y.Xu,Z.Lu,Y.H.Chiu,D.J. Levandier,and C.Y.Ng,J.Chem.Phys.134,201105 (2011).

    [7]Y.Xu,B.Xiong,Y.C.Chang,and C.Y.Ng,J.Chem. Phys.137,241101(2012).

    [8]Y.Xu,B.Xiong,Y.C.Chang,and C.Y.Ng,J.Chem. Phys.139,024203(2013).

    [9]C.Y.Ng,Annu.Rev.Phys.chem.53,101(2002).

    [10]T.Baer and Y.Li,Int.J.Mass Spectrom.219,381 (2002).

    [11]C.J.Hammond and K.L.Reid,Phys.Chem.Chem. Phys.10,6762(2008).

    [12]H.Gao,Z.Lu,L.Yang,J.Zhou,and C.Y.Ng,J. Chem.Phys.137,161101(2012).

    [13]C.Nicole,I.Sluimer,F.Rosca-Pruna,M.Warntjes, M.Vrakking,C.Bordas,F.Texier,and F.Robicheaux, Phys.Rev.Lett.85,4024(2000).

    [14]A.Osterwalder,M.J.Nee,J.Zhou,and D.M.Neumark,J.Chem.Phys.121,6317(2004).

    [15]D.M.Neumark,J.Phys.Chem.A 112,13287(2008).

    [16]H.Gao,Y.Xu,L.Yang,C.S.Lam,H.Wang,J.Zhou, and C.Y.Ng,J.Chem.Phys.135,224304(2011).

    [17]U.Jacovella,B.Gans,and F.Merkt,J.Chem.Phys. 139,084308(2013).

    [18]X.Ripoche,P.Asselin,F.Piuzzi,and I.Dimicoli, Chem.Phys.175,379(1993).

    [19]T.G.Wright,S.I.Panov,and T.A.Miller,J.Chem. Phys.102,4793(1995).

    [20]C.S.Lam,H.Wang,Y.Xu,K.C.Lau,and C.Y.Ng, J.Chem.Phys.134,144304(2011).

    [21]H.Woo,K.C.Lau,J.Zhan,C.Y.Ng,C.L.Li,W. K.Li,and P.M.Johnson,J.Chem.Phys.119,7789 (2003).

    [22]H.Gao,Y.Song,W.M.Jackson,and C.Y.Ng,J. Chem.Phys.138,191102(2013).

    [23]H.Gao,Y.Song,L.Yang,X.Shi,Q.Yin,C.Y.Ng,and W.M.Jackson,J.Chem.Phys.135,221101(2011).

    [24]Z.Jingang,B.Jones,Y.Xueliang,W.Jackson,and C. Y.Ng,J.Chem.Phys.128,014305(2008).

    [25]X.Yang,J.Zhou,B.Jones,C.Y.Ng,and W.M.Jackson,J.Chem.Phys.128,084303(2008).

    [26]J.Zhou,K.C.Lau,E.Hassanein,H.Xu,S.X.Tian, B.Jones,and C.Y.Ng,J.Chem.Phys.124,034309 (2006).

    [27]Y.Hou,H.K.Woo,P.Wang,X.Xing,C.Y.Ng,and K.C.Lau,J.Chem.Phys.129,114305(2008).

    [28]K.C.Lau,H.Woo,P.Wang,X.Xing,and C.Y.Ng, J.Chem.Phys.124,224311(2006).

    [29]B.Reed,C.S.Lam,Y.C.Chang,X.Xing,D.S.Yang, and C.Y.Ng,Astrophys.J.693,940(2009).

    [30]Z.Lu,Y.C.Chang,Q.Z.Yin,C.Y.Ng,and W.M. Jackson,Science 346,61(2014).

    [31]A.T.J.B.Eppink and D.H.Parker,Rev.Sci.Instrum. 68,3477(1997).

    [32]C.J.Dasch,Appl.Opt.31,1146(1992).

    [33]S.Willitsch,J.M.Dyke,and F.Merkt,Helvetica Chim. Acta 86,1152(2003).

    [34]S.Pratt,J.Chem.Phys.98,9241(1993).

    [35]M.Gasser,A.M.Schulenburg,P.M.Dietiker,A.Bach, F.Merkt,and P.Chen,J.Chem.Phys.131,014304 (2009).

    [36]X.Xing,B.Reed,K.C.Lau,C.Y.Ng,X.Zhang,and G.B.Ellison,J.Chem.Phys.126,171101(2007).

    [37]R.Signorell and F.Merkt,Mol.Phys.92,793(1997). [38]S.Willitsch,A.Haldi,and F.Merkt,Chem.Phys.Lett. 353,167(2002).

    [39]S.Willitsch and F.Merkt,Int.J.Mass Spectrom.245, 14(2005).

    老司机在亚洲福利影院| 狂野欧美激情性xxxx| 色精品久久人妻99蜜桃| 欧美乱妇无乱码| 99精品欧美一区二区三区四区| 国产三级在线视频| 午夜福利18| 又粗又爽又猛毛片免费看| 国产精品,欧美在线| 欧美黑人欧美精品刺激| 欧美乱妇无乱码| 国产高清videossex| 岛国在线观看网站| av在线蜜桃| 宅男免费午夜| 亚洲精品久久国产高清桃花| 欧美日韩福利视频一区二区| 天堂av国产一区二区熟女人妻| 全区人妻精品视频| 黄片小视频在线播放| 国产精品久久久久久人妻精品电影| 夜夜夜夜夜久久久久| 天堂av国产一区二区熟女人妻| 欧美日韩黄片免| 18禁国产床啪视频网站| 日本一二三区视频观看| 国产色爽女视频免费观看| 免费av不卡在线播放| eeuss影院久久| 免费大片18禁| 色综合欧美亚洲国产小说| 欧美日本视频| 女人高潮潮喷娇喘18禁视频| 午夜亚洲福利在线播放| 老司机福利观看| 好看av亚洲va欧美ⅴa在| 又粗又爽又猛毛片免费看| 亚洲国产欧美人成| 老司机午夜十八禁免费视频| 老司机福利观看| 天天一区二区日本电影三级| 热99在线观看视频| 男插女下体视频免费在线播放| 国产av一区在线观看免费| 色老头精品视频在线观看| 亚洲国产精品sss在线观看| 噜噜噜噜噜久久久久久91| 狂野欧美白嫩少妇大欣赏| 久久久久久久亚洲中文字幕 | 色综合站精品国产| 每晚都被弄得嗷嗷叫到高潮| 欧美又色又爽又黄视频| 两个人的视频大全免费| 亚洲真实伦在线观看| 精品国产美女av久久久久小说| 高清日韩中文字幕在线| 91麻豆精品激情在线观看国产| 美女黄网站色视频| 国产精品三级大全| 一个人免费在线观看的高清视频| 久久6这里有精品| 老鸭窝网址在线观看| 最近在线观看免费完整版| x7x7x7水蜜桃| 99久国产av精品| 国产精品久久久久久久久免 | 欧美bdsm另类| 亚洲自拍偷在线| 日本五十路高清| 99久久精品热视频| 亚洲狠狠婷婷综合久久图片| 99久久无色码亚洲精品果冻| 麻豆国产av国片精品| 少妇的丰满在线观看| 精品久久久久久久久久免费视频| 天天一区二区日本电影三级| 国产三级中文精品| 一进一出好大好爽视频| 国产精品自产拍在线观看55亚洲| 中国美女看黄片| 首页视频小说图片口味搜索| 啦啦啦韩国在线观看视频| 国内精品一区二区在线观看| 欧美一区二区精品小视频在线| 国产高清激情床上av| 亚洲性夜色夜夜综合| 一本久久中文字幕| av专区在线播放| 伊人久久大香线蕉亚洲五| 亚洲av日韩精品久久久久久密| 亚洲人成伊人成综合网2020| 手机成人av网站| 一区福利在线观看| 日韩欧美精品免费久久 | 国产精品亚洲一级av第二区| 欧美+亚洲+日韩+国产| 久久婷婷人人爽人人干人人爱| 中国美女看黄片| 久久国产精品人妻蜜桃| 97碰自拍视频| 成人午夜高清在线视频| 欧美中文综合在线视频| 一区二区三区激情视频| 亚洲av日韩精品久久久久久密| 非洲黑人性xxxx精品又粗又长| 亚洲成av人片在线播放无| 最好的美女福利视频网| 午夜福利在线观看吧| 熟女少妇亚洲综合色aaa.| 精品福利观看| 日韩欧美三级三区| 亚洲成人久久爱视频| 在线观看舔阴道视频| 色在线成人网| 中亚洲国语对白在线视频| av黄色大香蕉| 午夜福利在线在线| 国产精品 国内视频| 午夜久久久久精精品| 国产99白浆流出| 亚洲精品国产精品久久久不卡| 亚洲欧美激情综合另类| 欧美绝顶高潮抽搐喷水| 18禁国产床啪视频网站| 日本撒尿小便嘘嘘汇集6| 免费看光身美女| 男女做爰动态图高潮gif福利片| 18美女黄网站色大片免费观看| 亚洲天堂国产精品一区在线| 亚洲精品一区av在线观看| 亚洲av第一区精品v没综合| 欧美日韩亚洲国产一区二区在线观看| 精品一区二区三区人妻视频| 天美传媒精品一区二区| 精品久久久久久久久久久久久| 亚洲中文字幕一区二区三区有码在线看| 亚洲国产欧美人成| 蜜桃久久精品国产亚洲av| 欧美区成人在线视频| 长腿黑丝高跟| 校园春色视频在线观看| 亚洲电影在线观看av| 国产精品三级大全| 亚洲精品乱码久久久v下载方式 | www国产在线视频色| 成年版毛片免费区| 女生性感内裤真人,穿戴方法视频| 免费在线观看亚洲国产| 日本成人三级电影网站| av在线天堂中文字幕| 精品99又大又爽又粗少妇毛片 | 又黄又爽又免费观看的视频| 免费高清视频大片| 亚洲av不卡在线观看| 国产精品久久久久久久电影 | 天天一区二区日本电影三级| 九色成人免费人妻av| 两人在一起打扑克的视频| 欧美日韩瑟瑟在线播放| 99在线人妻在线中文字幕| 国产视频内射| 欧美三级亚洲精品| 日本免费a在线| www.熟女人妻精品国产| 久久欧美精品欧美久久欧美| 免费看十八禁软件| 国产亚洲精品av在线| 国内揄拍国产精品人妻在线| 非洲黑人性xxxx精品又粗又长| 又黄又爽又免费观看的视频| 国产亚洲欧美在线一区二区| 精品久久久久久久久久久久久| 午夜福利高清视频| 久久6这里有精品| 国产成人a区在线观看| 久久久久亚洲av毛片大全| 91久久精品国产一区二区成人 | 国产真人三级小视频在线观看| 久久久久九九精品影院| 特级一级黄色大片| 国产高清视频在线观看网站| 亚洲久久久久久中文字幕| 身体一侧抽搐| 一个人观看的视频www高清免费观看| 精品久久久久久久久久免费视频| 色噜噜av男人的天堂激情| 美女cb高潮喷水在线观看| 色综合婷婷激情| av欧美777| 欧美色视频一区免费| 又黄又粗又硬又大视频| 丰满人妻一区二区三区视频av | 国产99白浆流出| 超碰av人人做人人爽久久 | 精华霜和精华液先用哪个| 日韩高清综合在线| 亚洲国产日韩欧美精品在线观看 | 99热精品在线国产| 日韩 欧美 亚洲 中文字幕| 成年人黄色毛片网站| 搡老熟女国产l中国老女人| 久久久久久久久大av| 少妇人妻精品综合一区二区 | 我的老师免费观看完整版| 香蕉丝袜av| 国产色爽女视频免费观看| 男女午夜视频在线观看| 欧美日韩一级在线毛片| 免费看光身美女| 亚洲精品美女久久久久99蜜臀| 97超级碰碰碰精品色视频在线观看| 在线观看66精品国产| 久久久久久久午夜电影| 亚洲天堂国产精品一区在线| 国产真人三级小视频在线观看| 国产欧美日韩一区二区精品| 亚洲va日本ⅴa欧美va伊人久久| 国产成年人精品一区二区| 最新在线观看一区二区三区| 日本黄大片高清| 我要搜黄色片| 国产精品野战在线观看| 欧美黄色片欧美黄色片| 欧洲精品卡2卡3卡4卡5卡区| 欧美又色又爽又黄视频| 在线观看午夜福利视频| 欧美激情久久久久久爽电影| 90打野战视频偷拍视频| 精品人妻一区二区三区麻豆 | 亚洲欧美日韩高清专用| 国产黄片美女视频| 99国产综合亚洲精品| 九色国产91popny在线| 亚洲国产精品999在线| 一本久久中文字幕| 成年版毛片免费区| 一夜夜www| 免费人成在线观看视频色| 少妇丰满av| 婷婷六月久久综合丁香| 中出人妻视频一区二区| 每晚都被弄得嗷嗷叫到高潮| 高清毛片免费观看视频网站| 欧美不卡视频在线免费观看| 欧美一区二区亚洲| 99久久精品一区二区三区| 久久这里只有精品中国| 国产亚洲欧美98| 国产极品精品免费视频能看的| 日本免费a在线| 99精品在免费线老司机午夜| 日韩欧美三级三区| 欧美日韩亚洲国产一区二区在线观看| 国产野战对白在线观看| 51国产日韩欧美| 全区人妻精品视频| 国产精品久久视频播放| 婷婷亚洲欧美| 亚洲av一区综合| 日韩大尺度精品在线看网址| 丰满人妻一区二区三区视频av | 国产精品影院久久| 最近在线观看免费完整版| 国产精品99久久99久久久不卡| 中文字幕熟女人妻在线| 最新美女视频免费是黄的| 91av网一区二区| 黄色片一级片一级黄色片| 精品不卡国产一区二区三区| 国产精品久久久久久人妻精品电影| 国产亚洲精品综合一区在线观看| 99在线人妻在线中文字幕| 黄色丝袜av网址大全| 亚洲自拍偷在线| 成人性生交大片免费视频hd| 麻豆久久精品国产亚洲av| 欧美黑人巨大hd| 在线观看舔阴道视频| 天堂√8在线中文| 日韩欧美在线二视频| 婷婷精品国产亚洲av在线| 免费观看的影片在线观看| 色老头精品视频在线观看| 天天添夜夜摸| 精品久久久久久久人妻蜜臀av| 国产精品,欧美在线| 午夜福利成人在线免费观看| 少妇人妻一区二区三区视频| 国产日本99.免费观看| 精品国产超薄肉色丝袜足j| 日韩av在线大香蕉| 两性午夜刺激爽爽歪歪视频在线观看| 天天添夜夜摸| 91在线观看av| 舔av片在线| 国产激情欧美一区二区| 亚洲国产欧洲综合997久久,| 两个人看的免费小视频| 免费av毛片视频| 老司机福利观看| 亚洲真实伦在线观看| 在线观看日韩欧美| 99久久久亚洲精品蜜臀av| 欧美成人一区二区免费高清观看| 激情在线观看视频在线高清| а√天堂www在线а√下载| 一区二区三区国产精品乱码| 成年女人永久免费观看视频| 九色国产91popny在线| 男女做爰动态图高潮gif福利片| 久久精品国产99精品国产亚洲性色| ponron亚洲| 亚洲av免费高清在线观看| av天堂中文字幕网| 亚洲精品影视一区二区三区av| 国产免费男女视频| 在线观看美女被高潮喷水网站 | 成人高潮视频无遮挡免费网站| 90打野战视频偷拍视频| 午夜日韩欧美国产| 午夜福利免费观看在线| 一级毛片高清免费大全| 久久久久久大精品| 中文字幕av在线有码专区| 国产精品一区二区三区四区久久| 国产精品一区二区免费欧美| 日本精品一区二区三区蜜桃| 精品久久久久久久久久免费视频| 五月玫瑰六月丁香| 免费观看精品视频网站| 十八禁网站免费在线| 尤物成人国产欧美一区二区三区| 精品国产美女av久久久久小说| 亚洲男人的天堂狠狠| 啦啦啦免费观看视频1| av片东京热男人的天堂| 最好的美女福利视频网| 久久久国产成人精品二区| 午夜福利在线在线| 夜夜夜夜夜久久久久| 欧美又色又爽又黄视频| 夜夜夜夜夜久久久久| 99国产极品粉嫩在线观看| 欧美日韩一级在线毛片| 亚洲 国产 在线| 在线观看日韩欧美| 亚洲 国产 在线| 操出白浆在线播放| 亚洲av成人不卡在线观看播放网| 热99re8久久精品国产| 日韩免费av在线播放| 麻豆成人午夜福利视频| 日韩欧美 国产精品| 精品欧美国产一区二区三| 午夜福利在线观看免费完整高清在 | 欧美最黄视频在线播放免费| 此物有八面人人有两片| 别揉我奶头~嗯~啊~动态视频| 全区人妻精品视频| 男女视频在线观看网站免费| 国产黄片美女视频| 午夜福利在线观看吧| 亚洲av二区三区四区| 色播亚洲综合网| 久久精品国产综合久久久| 美女黄网站色视频| 老司机午夜福利在线观看视频| 亚洲专区中文字幕在线| 免费在线观看成人毛片| 精品人妻1区二区| x7x7x7水蜜桃| 国产真实伦视频高清在线观看 | 国产精品嫩草影院av在线观看 | 精品一区二区三区av网在线观看| 久久精品影院6| 99热精品在线国产| 九九热线精品视视频播放| 真人做人爱边吃奶动态| 一个人免费在线观看的高清视频| 一本综合久久免费| 成人午夜高清在线视频| 成年版毛片免费区| www.www免费av| 中文亚洲av片在线观看爽| 日韩高清综合在线| 亚洲电影在线观看av| 给我免费播放毛片高清在线观看| 亚洲国产高清在线一区二区三| 性色avwww在线观看| 精品不卡国产一区二区三区| a在线观看视频网站| 亚洲成人免费电影在线观看| 国产淫片久久久久久久久 | 亚洲欧美日韩无卡精品| 国产精品久久久久久久久免 | 99精品欧美一区二区三区四区| 夜夜躁狠狠躁天天躁| АⅤ资源中文在线天堂| 中文字幕av在线有码专区| 亚洲av电影在线进入| 国产午夜精品论理片| 51午夜福利影视在线观看| 日韩国内少妇激情av| 欧美激情在线99| 婷婷精品国产亚洲av| 一区二区三区免费毛片| 国产一级毛片七仙女欲春2| 亚洲成av人片免费观看| 精品人妻偷拍中文字幕| 日本撒尿小便嘘嘘汇集6| 无遮挡黄片免费观看| 国产v大片淫在线免费观看| 色噜噜av男人的天堂激情| svipshipincom国产片| 久久精品国产亚洲av香蕉五月| 天堂影院成人在线观看| www.999成人在线观看| 日韩大尺度精品在线看网址| 精品久久久久久久毛片微露脸| 国产亚洲精品一区二区www| 日本免费a在线| 国产精品一区二区三区四区免费观看 | 国产aⅴ精品一区二区三区波| 免费av不卡在线播放| or卡值多少钱| 国产av一区在线观看免费| 村上凉子中文字幕在线| 麻豆国产av国片精品| 国产亚洲精品一区二区www| 亚洲精品色激情综合| 精品午夜福利视频在线观看一区| 97碰自拍视频| 亚洲五月婷婷丁香| 亚洲国产精品合色在线| 久久久精品欧美日韩精品| 久久久精品大字幕| 热99在线观看视频| 又爽又黄无遮挡网站| 国产免费一级a男人的天堂| 两个人视频免费观看高清| 免费av观看视频| 国产精品一区二区三区四区久久| 亚洲天堂国产精品一区在线| 91av网一区二区| 蜜桃久久精品国产亚洲av| 宅男免费午夜| 一区二区三区国产精品乱码| 欧美成人性av电影在线观看| 国产亚洲精品久久久久久毛片| 亚洲av免费高清在线观看| 欧美+亚洲+日韩+国产| 又爽又黄无遮挡网站| 热99在线观看视频| 一级毛片女人18水好多| 国内毛片毛片毛片毛片毛片| h日本视频在线播放| 亚洲成av人片免费观看| 日本五十路高清| 精品日产1卡2卡| 久久久成人免费电影| 中文字幕久久专区| 听说在线观看完整版免费高清| 老司机深夜福利视频在线观看| 亚洲aⅴ乱码一区二区在线播放| 小说图片视频综合网站| 18禁黄网站禁片午夜丰满| 日韩成人在线观看一区二区三区| www.熟女人妻精品国产| 亚洲精品国产精品久久久不卡| 一本综合久久免费| 中国美女看黄片| 国产高清视频在线播放一区| 在线免费观看不下载黄p国产 | 在线天堂最新版资源| av女优亚洲男人天堂| 日本与韩国留学比较| 国产久久久一区二区三区| 久久久久久久久久黄片| 亚洲国产日韩欧美精品在线观看 | 午夜福利在线观看免费完整高清在 | 9191精品国产免费久久| 国产精品久久久久久人妻精品电影| 亚洲电影在线观看av| x7x7x7水蜜桃| 国产伦精品一区二区三区视频9 | 色吧在线观看| 一区二区三区国产精品乱码| 91久久精品国产一区二区成人 | 一区二区三区免费毛片| 午夜福利高清视频| 老司机午夜福利在线观看视频| 五月伊人婷婷丁香| xxx96com| 国产精品女同一区二区软件 | 亚洲va日本ⅴa欧美va伊人久久| 最后的刺客免费高清国语| 中文字幕av成人在线电影| 国产精品久久久久久久久免 | 国产国拍精品亚洲av在线观看 | 久久香蕉国产精品| 一个人看视频在线观看www免费 | 女人高潮潮喷娇喘18禁视频| www.熟女人妻精品国产| 久久久久久大精品| 国产精品99久久99久久久不卡| 一级作爱视频免费观看| eeuss影院久久| 夜夜躁狠狠躁天天躁| 香蕉丝袜av| 国产视频一区二区在线看| 亚洲人与动物交配视频| avwww免费| 成人一区二区视频在线观看| 国产高清有码在线观看视频| 麻豆久久精品国产亚洲av| 国产单亲对白刺激| 美女高潮喷水抽搐中文字幕| 亚洲一区二区三区不卡视频| 男插女下体视频免费在线播放| 波多野结衣高清无吗| 日本在线视频免费播放| 综合色av麻豆| 在线a可以看的网站| 精品熟女少妇八av免费久了| 一a级毛片在线观看| 99精品在免费线老司机午夜| 久久久久久久精品吃奶| 国产精品综合久久久久久久免费| 午夜福利成人在线免费观看| 国产成人福利小说| 国产精品国产高清国产av| 久久久国产成人免费| 亚洲精品在线观看二区| e午夜精品久久久久久久| 欧美av亚洲av综合av国产av| 最近最新中文字幕大全电影3| 黄色视频,在线免费观看| 国产黄片美女视频| 无人区码免费观看不卡| 国产精品久久久人人做人人爽| 久久久久久久亚洲中文字幕 | 老鸭窝网址在线观看| 日本免费一区二区三区高清不卡| 性欧美人与动物交配| 少妇人妻一区二区三区视频| 亚洲人成网站在线播| 国产av不卡久久| 成人高潮视频无遮挡免费网站| 欧美色视频一区免费| 3wmmmm亚洲av在线观看| 免费在线观看日本一区| 91久久精品国产一区二区成人 | 美女cb高潮喷水在线观看| 亚洲自拍偷在线| 久久久精品大字幕| 亚洲av熟女| 国产爱豆传媒在线观看| 亚洲不卡免费看| a级毛片a级免费在线| 一级黄色大片毛片| 成年免费大片在线观看| 日本精品一区二区三区蜜桃| 国内精品一区二区在线观看| 在线观看日韩欧美| 久久这里只有精品中国| 免费av观看视频| 男人的好看免费观看在线视频| 国产真实伦视频高清在线观看 | 久久久久精品国产欧美久久久| 舔av片在线| 亚洲精品在线美女| 欧美黄色淫秽网站| 国产成年人精品一区二区| 757午夜福利合集在线观看| 日韩有码中文字幕| 丰满的人妻完整版| 在线观看舔阴道视频| 久久九九热精品免费| 国产一区二区在线观看日韩 | 岛国视频午夜一区免费看| 国产精品乱码一区二三区的特点| 日韩亚洲欧美综合| 国内揄拍国产精品人妻在线| 国产亚洲精品综合一区在线观看| 禁无遮挡网站| 色噜噜av男人的天堂激情| 亚洲中文日韩欧美视频| 午夜免费激情av| 97人妻精品一区二区三区麻豆| 可以在线观看毛片的网站| 久久精品91无色码中文字幕| 亚洲真实伦在线观看| 高潮久久久久久久久久久不卡| 免费搜索国产男女视频| 亚洲真实伦在线观看| 午夜精品一区二区三区免费看| 亚洲aⅴ乱码一区二区在线播放| 欧美+亚洲+日韩+国产| 观看免费一级毛片| 成人午夜高清在线视频| 日韩欧美国产一区二区入口| 两个人看的免费小视频| 欧美zozozo另类| 波多野结衣高清无吗| 99热这里只有是精品50| 69av精品久久久久久| 婷婷精品国产亚洲av在线| 淫妇啪啪啪对白视频| 亚洲aⅴ乱码一区二区在线播放| 神马国产精品三级电影在线观看| 国产精品一区二区免费欧美| 91久久精品国产一区二区成人 | 成人av在线播放网站| 国产伦精品一区二区三区视频9 | 国产一区二区亚洲精品在线观看|