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

    Hardware for multi-superconducting qubit control and readout*

    2021-11-23 07:24:14ZhanWang王戰(zhàn)HaiYu于海RongliLiu劉榮利XiaoMa馬驍XueyiGuo郭學儀ZhongchengXiang相忠誠PengtaoSong宋鵬濤LuhongSu蘇鷺紅YirongJin金貽榮andDongningZheng鄭東寧
    Chinese Physics B 2021年11期
    關(guān)鍵詞:東寧

    Zhan Wang(王戰(zhàn)) Hai Yu(于海) Rongli Liu(劉榮利) Xiao Ma(馬驍) Xueyi Guo(郭學儀)Zhongcheng Xiang(相忠誠) Pengtao Song(宋鵬濤) Luhong Su(蘇鷺紅)Yirong Jin(金貽榮) and Dongning Zheng(鄭東寧)

    1Institute of Physics,Chinese Academy of Sciences,Beijing 100190,China

    2Beijing Academy of Quantum Information Sciences,Beijing 100193,China

    3Beijing Liuhe Lianchuang Technology Co.,Ltd.,Beijing 100011,China

    4University of Chinese Academy of Sciences,Beijing 100049,China

    5Songshan Lake Materials Laboratory,Dongguan 523808,China

    Keywords: superconducting qubit, dispersive readout, arbitrary-waveform generator (AWG), analog-digital converter(ADC)

    1. Introduction

    Superconducting qubits as a promising solution for quantum computing is developing towards intermediate scale processors that consist of tens of qubits.[1,2]To characterize and perform research with these processors, the ability of control and readout of multiple qubits is essential. Furthermore, it is suggested that to realize a logic qubit a large number of physical qubits is needed,according to the error rate achieved currently. In this regard, not only control and readout of many qubits is essential, but also the used hardware requires good scalability, high synchronization and low latency. Moreover,in order to achieve fault-tolerant quantum computing, active quantum error correction of qubits is inevitable.[2-4]Some progress has been made in this direction.[5-7]For fault-tolerant quantum computing, high fidelity gate operations that meet the expected fault-tolerant threshold must be achieved and the feedback control time should be less than one percent of the decoherence time.[11]This means high accuracy in qubit control and low feedback latency are required.[8]

    There have been a number of reports on the literature to describe multi-qubit control and readout systems.[9-12]

    In this paper,we report our efforts in developing an electronic hardware system for the purpose of control and readout of multiple superconducting transmon qubits.We demonstrate the functionality of the hardware system by applying it to control and readout a 10-qubits device.

    2. Basics of control and readout of superconducting transmon qubits

    Superconducting qubits are designed to have energy level spacing around 5 GHz. To control the qubit state,microwave pulses with appropriate frequency, amplitude and phase are needed.The standard method is to define the pulse shape using anIQmixer. In-phase and quarture control signals are generated by AWGs. With theIandQsignals,the microwave signal from a microwave source can be modulated into a pulsed signal with desired frequency, phase and shape. For energy level tunable qubits, one may apply a fast or a slow DC bias to change the qubit energy level spacing. Thus, in order to control one qubit, three AWG channels (two forIQcontrol and one for fast bias control) and one DC bias channel are needed. For multi-qubit control,the number of AWG and DC bias channels should increase accordingly.

    The readout of transmon qubits is realized by the dispersive readout scheme.[14]The underlying physics is the coupling interaction between a resonator and a two-level atom(or qubit), being well understood in quantum optics. When the frequency detuning?, difference between the resonator frequencyωrand the qubit frequencyωq, is much larger than the coupling rateg, the system is in the so-called dispersive coupling regime. In this case,the transition of the qubit from the|0〉state to the|1〉state would lead to the change of the resonant frequency. In other words, detecting the qubit state is realized by observing the shift in the resonant frequency of a readout resonator interacting with the qubit. Apart from its nondestructive nature, another advantage of the dispersive readout scheme is that it can be easily adapted into multiplex readout form. In a multi-qubit device,each qubit is coupled to a readout resonator. The readout resonators are also coupled to a microwave transmission line and the resonant frequency of each resonator is made slightly different from others. In such an arrangement,one can detect a qubit by measuring the amplitude or phase change of the corresponding resonator.

    The electronic hardware for one readout channel includes a microwave source,two AWGs and a data acquisition board.The AWG signals mix with microwave signal via anIQmixer to form readout pulses that contain frequency values corresponding to the resonant frequency of each readout resonator.The readout pulses are then sent to the readout transmission line on the quantum chip. The signals from the chip are amplified and demodulated by the same local microwave signal via anotherIQmixer. The demodulatedIQsignals are acquired by the data acquisition board and processed either by on-board microprocessor or a computer. Usually,one readout transmission line can be used for detecting about ten qubits.

    In Fig. 2, schematics of experimental setup for multiqubit system is shown. The dilution refrigerator provides a milliKelvin low temperature environment. The control and readout hardware units are connected to the device via microwave coaxial cables, with attenuators and filters being inserted at different low temperature stages for each line to reduce external electromagnetic interference.

    From the above description, we can see that multichannel AWGs and DC bias sources together with data acquisition board are essential for the control and readout of multiqubit quantum chip.

    3. Design consideration

    Scalability To achieve good scalability of the system,we intentionally modularize the specific readout, control and other hardware functions. This means we can easily expand the system for more qubits. In our modularized hardware system,there are core board,clock board,readout board,control board and bias board.Each type of boards provide corresponding functions of qubit control and readout. We consider two kinds of scalability requirements. One is scalability within the chassis. If a new module (such as microwave AWG module)is needed, we can directly insert the corresponding modular board to the chassis, without considering its power supply,clock, trigger and other issues. Another one is the scalability between the chassis. If we need more and more readout and control hardware, we can link multiple chassis together through clock and trigger connections,which can ensure various chassis work together orderly.

    Synchronization High synchronization between different control channels is required to ensure accurate manipulation of qubit quantum state, and thus to reduce error accumulation. In our system,different module boards in a chassis are controlled by the same clock board, which distributes its high-frequency clock and trigger signals to other boards. In addition, the system can accept external reference clock and trigger signals to ensure the synchronization of different chassis.

    Latency The latency in qubit control and readout is mainly limited by three factors. The first one is the data transmission time between the master computer(usually a PC)and the slave computer as well as the time between the slave computer and FPGA. The following two methods are adopted to reduce this data transmission time. (i) Using direct memory access (DMA) transfer technology to improve the communication efficiency between the master and the slave computers.(ii)Instead of transferring all data to the master computer for data demodulation process,data acquired by the ADC are demodulated by the on-board FPGA of the readout module.

    The second limiting factor is due to the board control time. In other words, there is a response time for boards to store and execute after receiving a string of commands.When more than one boards are in use,we implement an asynchronous parallel control strategy in the program to reduce the module waiting time.

    The third limiting factor of latency is related to the feedback control.To realize quantum feedback,the received signal from qubits should be processed first and then send commands to the AWGs to generate control signals.In order to reduce the time delay, we integrate the AWGs and ADCs on one board and thus they are controlled by the same FPGA.The time required for ADC data acquisition to data processing, and then to AWG output can be optimized within one FPGA clock cycle.

    4. Details of hardware

    Figure 1 shows the photograph and schematic diagram of our hardware system. The system is mounted in a VPX-6U chassis. VPX motherboard in the chassis provides the power supply,communication hardware and RF connectors for each utility board.Up to 10 different types of boards can be inserted into this chassis. As shown in Fig.1(b),all boards communicate with the core board by peripheral component interconnect express (PCIE) bus. In the whole hardware system, the core board acts as the slave computer,which can be remotely controlled by the master computer. The function of clock board is to provide standard clock and trigger signals for other boards.

    The core board is a microcomputer with an Intel i75850EQ processor. All boards communicate with the core board through the bus switch (as shown in Fig. 1(b)) by the PCIE protocol. The core board communicates with the external measurement system through gigabit ethernet via User Datagram Protocol(UDP).

    The clock board is composed of an FPGA chip(XILINX XC7K325T),a low-noise broadband RF phase-locked chip(TI LMX2592)and other auxiliary hardware. The clock board accepts an external 10 MHz clock and an external trigger for synchronization between chassis. The function of the clock board is to provide standard clocks and different adjustable trigger signals for other boards. For example, the clock board provides a 1.25 GHz clock to the ADC chip,a 2.5 GHz clock to the digital-analog-converter(DAC)chips,and a low-frequency clock to the DC board.

    Fig. 1. (a) Photo graph of the hardware system. This chassis is a VPX-6U chassis. (b)Schematic diagram of the hardware system.

    Fig.2. Schematic experimental setup for superconducting multi-qubit systems.

    The clock board is composed of an FPGA chip(XILINX XC7K325T),a low-noise broadband RF phase-locked chip(TI LMX2592)and other auxiliary hardware. The clock board accepts an external 10 MHz clock and an external trigger for synchronization between chassis. The function of the clock board is to provide standard clocks and different adjustable trigger signals for other boards. For example, the clock board provides a 1.25 GHz clock to the ADC chip,a 2.5 GHz clock to the digital-analog-converter(DAC)chips,and a low-frequency clock to the DC board.

    The primary function of the readout board is to readout the quantum state of qubits. The readout board consists of two ADCs and two AWGs. As described in Section 2, the AWGs are used for the modulation of readout microwave signals and the ADCs are used for data acquisition. In addition,this board can also perform fast quantum feedback control.The FPGA on the readout board can simultaneously control the ADC and DAC to achieve fast feedback.The readout board is composed of an FPGA chip (XILINX XC7V690T), two DAC chips (LINEAR LTC2000IY-16), two ADC chips (TI ADC12D1800RF) and other peripheral auxiliary hardware.The DAC chips,as the same as used in the control board,function as two AWGs. The signal acquisition is performed by using two ADC chips with±1 V input voltage range, 12-bit resolution and 1.25 GHz sampling rate.

    The control boards are actually AWG boards. Each control board is composed of an FPGA chip (XILINX XC7V690T-2FF1761I), 6 DAC chips (LINEAR LTC2000IY-16) with 16-bit resolution and 2.5 GHz sampling rate, and other auxiliary chips. It is able to generate 6 channels of independent AWG signals with 16 bit resolution,600 MHz bandwidth and±1 V output voltage range. As shown in Fig. 2,the AWG signals are used to control qubit by modulating amplitude,frequency and phase of a microwave signal via anIQmixer. The modulated microwave pulses control the operation of the qubit in theXYdirection. Additionally, the AWG signal can also be used to generate a fast-changing waveform to directly control the qubitZbias to realize the gate of the qubit in theZdirection.

    The bias board is composed of an FPGA chip (XILINX XC7k325t),16 DAC chips(TI DAC5682Z)and related auxiliary hardware. The board provides 16 channels±10 V of DC voltage output with 20-bit resolution. As shown in Fig.2,the function of the bias board is to provide a DC bias for changing qubit working frequency to find a suitable operation and readout point. Additionally,the board can also be used to control the microwave switch,which can further reduce unwanted microwave leakage to the qubits and improve the control fidelity.

    4.1. Waveform generation

    In order to provide flexibility in operating the system at the different stages of experiments, we design two different waveform generation modes,namely,the direct mode and the sequence mode. The two waveform generation modes are shown in Fig.3.

    In the process of the direct mode,the DMA transfers data to FPGA,and FPGA loads it to AWG.Finally,the AWG generates waveforms. Although this method may cost more time for data transmission between components, it is simple and flexible to use during experiment preparation and test stage.

    Fig.3. Logic diagram for waveform generation. Two output modes of AWG are designed.

    In the sequence mode,we construct 64 memory units inside the FPGA.Each of the units consists of waveform information related to the id number, play settings and data. The information stored in the memory units can be recalled, arranged and compiled into a waveform in accordance with the requirement of quantum operations in experiments. For instance, a qubit gate operation waveform can be stored inside a sequence memory unit and the arrangement and compile of these sequence is equivalent to the operations of the entire experimental Hamiltonian. In this way, once the waveforms of the various gates are determined,we can control the entire operations with predetermined combination of the gates.

    4.2. Multi-qubit demodulation

    According to Section 2, theIQdata can be collected and demodulated by readout board to get information of qubit state.One may transmit all theIQraw data to master computer and perform data process,or use the FPGA to do it on-board.Obviously, the latter can reduce the latency of data transmission and processing, being beneficial to quantum feedback function.

    In quantum experiment,the data size is huge and the network communication between the master and slave computer is a time-consuming task. Realizing demodulation in FPGA can drastically decrease the time of data transmission. For instance, for one readout, 4096 16-bit data could be captured 2000 times by ADC.After the FPGA demodulation,only one 16-bit qubit state data needs to be transmitted and the whole data volume is reduced by 4096×2000 times, which greatly reduce the data transmission time.

    The demodulation algorithm in FPGA can be traced back to Fourier transform

    Here,the summation is over the time domain of data acquisition.

    As shown in Fig. 4(a), the demodulatedIandQdata of the|0〉and|1〉states can be displayed on theI-Qplane. The blue and orange dots are taken when the qubit is set at|0〉and|1〉states, respectively. Due to decoherence and noise, some of the orange dots are located in the blue cloud,and some blue dots appear in the blue cloud.

    Fig.4. Schematic diagram of qubit state rotation in IQ plane,before(a)and after(b)rotation.

    Therefore, the qubit state can be distinguished in FPGA by following operations.Firstly,we demodulate qubitIQdata,mark it on theI-Qplane and calculate the angle between the central axis of the two states and in-phase axis. According to the angle, we rotate the data so that the central axis of two clouds is parallel to the in-phase axis, as shown in Fig. 4(b).We can set a threshold between the two clouds for deciding the qubit state. After obtaining the rotate angle and the best threshold, these parameters can be used to distinguish subsequent qubit state. The corresponding algorithm steps in FPGA are shown in Fig.5 schematically.

    For readout of a multi-qubit device,each qubit can be demodulated by its own frequency. When measuringkqubits at the same time, the measured state is one of 2kpossible states.We can express the state using a binary number(such as 1010110100 for ten qubits)and assign it as the address of the registers in the FPGA.When this binary number appears once,the register will count once. Therefore, if we repeat readout measurementNtimes,and the same binary number appearsntimes,the register will countntimes. Then,the probability of the state isn/N.In this way,the information of the multi-qubit state can be obtained in the FPGA.

    Fig.5. The multi-qubits demodulation logic diagram in FPGA.

    4.3. Quantum feedback

    Active feedback control of qubits is an essential step to achieve quantum error correction. The low latency is needed in qubit feedback control.[11]Considering low latency,using FPGA as a parallel and fast processing method has more advantages than traditional using central processing unit(CPU).[2,13]

    In our system,we may use the readout board for fast feedback. The FPGA on the board controls both the ADCs and AWGs at the same time. Also, we use the real-time parallel demodulation algorithm in FPGA to reduce the calculation time as much as possible. In this demodulation algorithm,demodulation process(Eqs.(4)and(5))and the data acquisition process by ADC are performed simultaneously.

    Finally, we determine the total feedback latency related to the hardware by connecting AWGs output to the input of ADCs and measuring the time delay using the following procedure, as shown in Fig. 6. Initially, we let the AWGs generate a simulatedIQsignals that correspond to the|0〉or|1〉state of qubit. After receiving the data,the FPGA will capture the data from ADC and perform data demodulation process to distinguish if the|0〉or|1〉state information is received.Then,the FPGA will trigger AWGs to generate predetermined control pulses accordingly. The time domain corresponding to the feedback latency (plus 5 ns time delay due to 0.5 meter long coaxial cable connecting AWGs and ADCs)is indicated in Fig.6. The obtained feedback latency is about 178.4 ns.

    In Table 1, the latency for each step of feedback process is listed. The time on FPGA is determined by its clocks. One clock is 3.2 ns for the FPGA we use. The AWG output time of 4.4 ns is from the chip data sheet. The internal circuit delay is time for AWG signal propagating from AWG chip,through a differential amplifier and filters,to the output SMA connector on the chassis.

    Table 1. The parts of feedback logic latency.

    Fig.6. Quantum feedback control test data recorded on an oscilloscope. In this testing,the whole period time is 1.5μs,the readout state simulation time is 819.2 ns, the feedback latency is 178.4 ns and the trigger delay time is 65 ns.

    5. Hardware testing

    A series of tests have been conducted on our system, including hardware test and application test on qubit devices. A R&S?RTP084 oscilloscopes is used for AWG jitter test and a R&S?FSV3000 signal and spectrum analyzer is used for phase noise and spurious free dynamic range(SFDR)test.

    We test the jitter among two AWG channels by simultaneously outputting square signals to the oscilloscope and then recording the delays and standard deviations between these two AWG channels. After testing all six channels on one board,the standard deviations of jitter is about 8 ps.

    Fig.7. Phase noise test results of one AWG channel.

    The delays of AWG channels are affected by different length of outgoing line. There are three different length types(265 mm,225 mm,and 190 mm). We find that the delays between different AWG channels on the same control board are less than 1 ns and the delays of AWG channels between two different control boards are about 5-6 ns.

    In Fig.7,we show the phase noise measured on one AWG channel at seven different frequency values. The noise floor of the equipment is about?140 dBc/Hz. The phase noise with a frequency offset greater than 10 MHz tends to be consistent.Increasing the signal frequency yields downward movement of the phase noise curve.

    In the SFDR measurement,the harmonic noise represents the main spurious noise source, especially when the single frequency signal is used. As shown in Fig. 8, SFDR curves with and without harmonic noise are plotted,from 10 MHz to 450 MHz.

    Fig.8. Spurious free dynamic range of one AWG channel.

    6. Applications of the hardware in multi-qubit characterizations

    In order to verify the functionality in qubit control and readout for the hardware system we built, we carried out a number of characterization measurements on a home made 10-qubit device.

    Fig.9. The I-Q data for qubit-state differentiation. Blue and red dots are measured I-Q values when qubit is prepared in|0〉and|1〉state,respectively. A total of 5000 dots(repetitions)for each color.

    In Fig.9,we show the readout data of 10 qubits. As described in Section 2,these are dispersive readout scheme data with blue and red dots being taken when qubits are in the|0〉and|1〉state,respectively. In this case,the 10 qubits are readout simultaneously.

    The energy relaxation timeT1and the dephasing timeT*2results measured of one qubit are displayed in Fig. 10.The control pulses for these measurements are schematically shown in the insets.

    Fig.10. (a)T1 and(b)T*2 result of one qubit.

    We also performed square wave distortion calibration,with the hardware we built, for square wave signals used for control the energy level spacing. The qubits characterized are energy level tunable ones. In experiments, it is often required to change the energy level quickly by applying a square wave pulse to theZ-control line of one qubit. However, due to the presence of parasitic inductance and capacitance, an ideal square pulse is usually distorted when reaching to the chip, showing overshoot or undershoot near the rising edge and tailed falling edge.

    To correct the distortions, we need to determine the deformed square pulse firstly. This is carried out by using the following procedure. The corresponding qubit is biased to a relatively low frequencyfz(more than 1 GH below its maximum frequency) to improve its sensitivity to the pulse shape deformation. An amplitude fixed square pulse is applied to a qubit,and is followed aπpulse with frequencyfz. If there is no distortion, the qubit will be excited to the|1〉state by theπpulse. If there is small distortion after the square wave,the qubit will not be completely excited. By varying the compensation amplitude ?Z,we can find the value of full compensation when the maximum excitation probability is achieved. By changing the delay time ofπpulse, we can get a view of the distorted square wave response. With these response data,we can calculate the needed adjustments for the square pulses.

    In Fig. 11, we show the results for one qubit. As illustrated in Fig.11(a),the ideal square pulse(blue)shows distortion(yellow)due to non-ideal response function. After applying the calibrated pulse (green), the response shape reaching the qubit should be close to the ideal one. In Fig. 11(b), the data after the falling edge are shown. As shown in Fig. 8(b),we can measure the falling edge, which should be a straight line with zero voltage under ideal circumstances. The original uncalibrated data and the data after calibrations are shown.We repeat the pulse correction process three times,indicating as the first,second,and third-order calibration in the figure.

    Fig.11. (a)The Z pulse correction schematics and(b)correction results after the falling edge.

    Finally, we characterize single qubit gate fidelity using random benchmarking.[15]In this kind of measurement,a random sequence of gates belonging to the single qubit Clifford group is prepared and the experimental fidelity is obtained.This is set as the reference test. Then,to measure the fidelity of a specific gate,the gate is interleaved with m random Cliffords and perform the same measurement. By comparing with the reference, the average fidelitys of the gates are obtained.In Fig. 12, the representative reference gate sequence corresponding to 10 random Cliffords (separated by the dashed lines), and a recovery gate (in light blue region) are shown schematically. The measurement results displayed in Fig. 13 are for six gates of single qubit. For each gate, the operation time is 50 ns. The results show that averaged gate fidelity is better than 99.68%for the six gate we characterize.

    Fig.12. Representative reference gate sequency.

    Fig.13. Random bench marking result.

    7. Conclusion

    We have designed, built and tested a hardware system used for the control and readout of superconducting multiqubit devices. The system consists of multi-channel AWG boards, DC bias board, and readout boards. The AWG has 2.5 GHz sample rate and 16-bit voltage resolution. The system can be used for more and more energy level tunable qubit readout and control by connecting more and more boards and chassis. Efforts have been made to reduce the system latency and the total feedback latency related to the hardware is about 178.4 ns. The system has been used to characterize a 10-qubit device,demonstrating its usability in practical quantum experiments.

    猜你喜歡
    東寧
    Post-solitons and electron vortices generated by femtosecond intense laser interacting with uniform near-critical-density plasmas
    Fabrication of Josephson parameter amplifier and its application in squeezing vacuum fluctuations?
    Observation of the exceptional point in superconducting qubit with dissipation controlled by parametric modulation?
    寧波市海曙東寧工具有限公司
    An Analysis of the Difficulties and Learning Methods of English Grammar in Senior High Schools
    2019中俄文化節(jié) 相約東寧
    伙伴(2019年9期)2019-10-21 11:30:18
    Tunable coupling between Xmon qubit and coplanar waveguide resonator?
    Fabrication and characterization of ultra-low noise narrow and wide band Josephson parametric amplifiers?
    唯愛墜星海
    飛言情B(2017年3期)2017-03-28 21:22:09
    我們的作品
    色婷婷久久久亚洲欧美| 日韩中文字幕视频在线看片| 中国三级夫妇交换| 天堂8中文在线网| 国产精品久久久久久精品电影小说| 建设人人有责人人尽责人人享有的| 精品亚洲成国产av| 人体艺术视频欧美日本| 街头女战士在线观看网站| 国产精品欧美亚洲77777| 国产精品久久久久成人av| 久久精品国产鲁丝片午夜精品| 国产毛片在线视频| 人妻制服诱惑在线中文字幕| 久久精品国产亚洲av涩爱| 高清欧美精品videossex| 少妇猛男粗大的猛烈进出视频| 国产片特级美女逼逼视频| 国产精品熟女久久久久浪| 成年人午夜在线观看视频| 性色avwww在线观看| 少妇的逼好多水| 男女国产视频网站| 99久久综合免费| 我要看黄色一级片免费的| 国产一区二区三区av在线| 国产成人91sexporn| 成年人免费黄色播放视频| 国产成人freesex在线| 不卡视频在线观看欧美| av.在线天堂| 99精国产麻豆久久婷婷| av播播在线观看一区| 国产精品国产三级国产av玫瑰| 中文字幕最新亚洲高清| 亚洲精品,欧美精品| 久久午夜福利片| 国产免费一级a男人的天堂| 欧美亚洲 丝袜 人妻 在线| 欧美xxⅹ黑人| 纵有疾风起免费观看全集完整版| 亚洲经典国产精华液单| 欧美三级亚洲精品| 欧美日韩综合久久久久久| 岛国毛片在线播放| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲国产欧美日韩在线播放| 老女人水多毛片| 国产精品不卡视频一区二区| 国产免费现黄频在线看| 一级毛片 在线播放| 精品国产国语对白av| 亚洲精品乱码久久久v下载方式| 人妻制服诱惑在线中文字幕| 亚洲美女视频黄频| 特大巨黑吊av在线直播| 不卡视频在线观看欧美| 国产成人91sexporn| 日韩欧美精品免费久久| 精品视频人人做人人爽| 亚洲精品久久久久久婷婷小说| 黑人猛操日本美女一级片| 狂野欧美激情性bbbbbb| 亚洲性久久影院| 久久99精品国语久久久| 中文字幕av电影在线播放| 亚洲av不卡在线观看| 黄色一级大片看看| 久久久久久久久大av| 人妻人人澡人人爽人人| 大又大粗又爽又黄少妇毛片口| 亚洲国产精品一区二区三区在线| 亚洲欧洲日产国产| 丰满饥渴人妻一区二区三| 91成人精品电影| 国产老妇伦熟女老妇高清| 欧美日韩综合久久久久久| 香蕉精品网在线| 免费久久久久久久精品成人欧美视频 | 婷婷色麻豆天堂久久| 久久久久精品久久久久真实原创| 色吧在线观看| 色哟哟·www| 天美传媒精品一区二区| 日韩免费高清中文字幕av| 国产亚洲一区二区精品| tube8黄色片| 久久精品国产亚洲av涩爱| 满18在线观看网站| 3wmmmm亚洲av在线观看| 日本猛色少妇xxxxx猛交久久| 精品卡一卡二卡四卡免费| 黄片播放在线免费| 蜜桃国产av成人99| 国产有黄有色有爽视频| 亚洲国产欧美在线一区| 国产乱来视频区| 久久久久久久久久成人| 日本wwww免费看| 一级毛片我不卡| 国产一区有黄有色的免费视频| 三级国产精品欧美在线观看| 一级二级三级毛片免费看| 999精品在线视频| 国产欧美亚洲国产| 亚洲怡红院男人天堂| 一本久久精品| 日韩欧美一区视频在线观看| 亚洲熟女精品中文字幕| 91精品国产国语对白视频| 伊人久久国产一区二区| 天天躁夜夜躁狠狠久久av| 制服人妻中文乱码| 午夜精品国产一区二区电影| 熟女电影av网| 伦理电影免费视频| av线在线观看网站| 久久99一区二区三区| 夫妻午夜视频| 18禁动态无遮挡网站| .国产精品久久| 亚洲一级一片aⅴ在线观看| 五月天丁香电影| 日韩一区二区视频免费看| 久久久久久久久久久免费av| 欧美 日韩 精品 国产| 尾随美女入室| av在线app专区| 亚洲成人一二三区av| 美女中出高潮动态图| 成人亚洲精品一区在线观看| 91国产中文字幕| 亚洲内射少妇av| 性色avwww在线观看| 成年av动漫网址| 人妻一区二区av| 22中文网久久字幕| 五月玫瑰六月丁香| 精品卡一卡二卡四卡免费| 特大巨黑吊av在线直播| 热99久久久久精品小说推荐| 老司机影院毛片| 亚洲精品日韩在线中文字幕| 99热国产这里只有精品6| 久久99热这里只频精品6学生| 另类精品久久| 在线播放无遮挡| 免费看光身美女| 少妇被粗大猛烈的视频| 国产精品一区www在线观看| 欧美激情 高清一区二区三区| 精品国产乱码久久久久久小说| 91久久精品国产一区二区成人| 国产亚洲最大av| 中文精品一卡2卡3卡4更新| av视频免费观看在线观看| 久久久久久久久久成人| 精品人妻偷拍中文字幕| 有码 亚洲区| 亚洲精品日本国产第一区| 涩涩av久久男人的天堂| 十八禁网站网址无遮挡| 在线天堂最新版资源| 欧美人与性动交α欧美精品济南到 | 一级片'在线观看视频| 午夜免费鲁丝| 飞空精品影院首页| 国产亚洲精品第一综合不卡 | 色哟哟·www| 国产老妇伦熟女老妇高清| 女性被躁到高潮视频| 91成人精品电影| 大陆偷拍与自拍| 一区二区av电影网| 亚洲色图综合在线观看| 日本黄色片子视频| 午夜免费男女啪啪视频观看| 国产精品 国内视频| 国产亚洲午夜精品一区二区久久| 亚洲国产色片| 一区二区日韩欧美中文字幕 | 熟女电影av网| 五月玫瑰六月丁香| 国产精品99久久99久久久不卡 | 国产精品一二三区在线看| av在线app专区| 国产成人精品福利久久| 国产精品成人在线| 伦理电影免费视频| 国产亚洲一区二区精品| 高清视频免费观看一区二区| 亚洲一区二区三区欧美精品| 桃花免费在线播放| 美女福利国产在线| 成人手机av| 亚洲一区二区三区欧美精品| 午夜久久久在线观看| 99久久综合免费| 曰老女人黄片| 亚洲精品国产av蜜桃| 日本黄色日本黄色录像| 精品视频人人做人人爽| av.在线天堂| 亚州av有码| 亚洲av日韩在线播放| av.在线天堂| 91精品国产九色| 国产精品欧美亚洲77777| 少妇丰满av| 亚洲情色 制服丝袜| 丰满饥渴人妻一区二区三| 国产一区有黄有色的免费视频| 欧美日本中文国产一区发布| 欧美日韩视频精品一区| 久久97久久精品| 午夜精品国产一区二区电影| 国产白丝娇喘喷水9色精品| 日本vs欧美在线观看视频| av卡一久久| 精品人妻一区二区三区麻豆| 亚洲av中文av极速乱| 免费观看性生交大片5| 日本欧美视频一区| 亚洲美女视频黄频| 久久久精品区二区三区| 久久久亚洲精品成人影院| 国产成人免费无遮挡视频| 中文字幕久久专区| 色吧在线观看| 丝袜在线中文字幕| 国产永久视频网站| 一级毛片我不卡| 人人妻人人爽人人添夜夜欢视频| 一级片'在线观看视频| 久久久久久久国产电影| 中文欧美无线码| 女性被躁到高潮视频| 精品酒店卫生间| 激情五月婷婷亚洲| 亚洲av成人精品一区久久| 成人手机av| 色婷婷久久久亚洲欧美| 女人久久www免费人成看片| 黑人猛操日本美女一级片| 新久久久久国产一级毛片| 大香蕉久久成人网| 一本久久精品| av有码第一页| 人人澡人人妻人| 国产一区有黄有色的免费视频| 五月伊人婷婷丁香| a级毛色黄片| 热99久久久久精品小说推荐| 日韩中文字幕视频在线看片| 亚洲av综合色区一区| 精品酒店卫生间| 九九久久精品国产亚洲av麻豆| 日韩中字成人| 极品少妇高潮喷水抽搐| 在线观看一区二区三区激情| 成人免费观看视频高清| 久久热精品热| 国产伦精品一区二区三区视频9| 午夜老司机福利剧场| 欧美一级a爱片免费观看看| av有码第一页| 91国产中文字幕| 观看美女的网站| 99久久人妻综合| 大话2 男鬼变身卡| 久久ye,这里只有精品| 国产无遮挡羞羞视频在线观看| 欧美激情极品国产一区二区三区 | 少妇被粗大猛烈的视频| 国语对白做爰xxxⅹ性视频网站| 亚洲三级黄色毛片| 亚洲av国产av综合av卡| 国产高清有码在线观看视频| 欧美xxxx性猛交bbbb| 国产毛片在线视频| 国产极品粉嫩免费观看在线 | 亚洲av综合色区一区| 如何舔出高潮| 国产乱人偷精品视频| 欧美xxⅹ黑人| www.av在线官网国产| 99久久精品国产国产毛片| av国产久精品久网站免费入址| 久久国产亚洲av麻豆专区| 日本黄色片子视频| 99re6热这里在线精品视频| 欧美精品亚洲一区二区| h视频一区二区三区| 如日韩欧美国产精品一区二区三区 | 日本wwww免费看| 女性被躁到高潮视频| 国产毛片在线视频| 99热网站在线观看| a 毛片基地| 99re6热这里在线精品视频| 日本vs欧美在线观看视频| a 毛片基地| 黄色一级大片看看| 午夜激情久久久久久久| 久久精品久久久久久久性| 99久久精品一区二区三区| 国产免费一级a男人的天堂| 久久人人爽人人爽人人片va| 日本91视频免费播放| 欧美日韩精品成人综合77777| 麻豆精品久久久久久蜜桃| 大话2 男鬼变身卡| 国内精品宾馆在线| 欧美三级亚洲精品| 久久97久久精品| 黄色欧美视频在线观看| 十八禁网站网址无遮挡| 久久精品久久久久久久性| 午夜福利,免费看| 99热网站在线观看| 一二三四中文在线观看免费高清| 亚洲av不卡在线观看| 久久国产精品大桥未久av| 日韩精品有码人妻一区| 国语对白做爰xxxⅹ性视频网站| av专区在线播放| av播播在线观看一区| 性色avwww在线观看| 国产午夜精品一二区理论片| 日韩不卡一区二区三区视频在线| a级毛片黄视频| 99九九线精品视频在线观看视频| 激情五月婷婷亚洲| 日韩 亚洲 欧美在线| 一区二区日韩欧美中文字幕 | 观看美女的网站| 2018国产大陆天天弄谢| 毛片一级片免费看久久久久| 久久久久久久大尺度免费视频| 亚洲国产精品999| 久久精品久久久久久噜噜老黄| 欧美少妇被猛烈插入视频| 纯流量卡能插随身wifi吗| 国产熟女欧美一区二区| 最新的欧美精品一区二区| 久久久欧美国产精品| 老女人水多毛片| 街头女战士在线观看网站| 91在线精品国自产拍蜜月| 国产精品一区二区在线不卡| 嘟嘟电影网在线观看| 在线亚洲精品国产二区图片欧美 | 亚洲怡红院男人天堂| 免费观看无遮挡的男女| 九九爱精品视频在线观看| av.在线天堂| 天堂8中文在线网| 熟妇人妻不卡中文字幕| 纵有疾风起免费观看全集完整版| 国产精品女同一区二区软件| 精品人妻一区二区三区麻豆| 日韩一本色道免费dvd| 99久久综合免费| 精品人妻在线不人妻| 日韩不卡一区二区三区视频在线| 最近2019中文字幕mv第一页| av有码第一页| av一本久久久久| 18禁观看日本| 国产成人精品无人区| 丰满迷人的少妇在线观看| 另类精品久久| 老女人水多毛片| 中文乱码字字幕精品一区二区三区| 少妇 在线观看| 亚洲成人一二三区av| 日韩欧美一区视频在线观看| 久久久久精品久久久久真实原创| 国产免费视频播放在线视频| 涩涩av久久男人的天堂| 久久久久久久久久久久大奶| 一本色道久久久久久精品综合| 99热全是精品| 亚洲人成77777在线视频| 精品午夜福利在线看| 免费看av在线观看网站| 边亲边吃奶的免费视频| 国产视频首页在线观看| 亚洲精品456在线播放app| 午夜免费男女啪啪视频观看| 国产国语露脸激情在线看| 99久久综合免费| av在线app专区| 欧美日韩视频精品一区| 人妻少妇偷人精品九色| 伦理电影大哥的女人| 日本猛色少妇xxxxx猛交久久| 亚洲不卡免费看| 亚洲美女黄色视频免费看| 久久久久精品性色| 久久久久久久久久久免费av| 一级毛片电影观看| 高清在线视频一区二区三区| 看十八女毛片水多多多| 99精国产麻豆久久婷婷| 久久久久视频综合| 另类精品久久| 永久网站在线| 国产精品三级大全| 中文字幕精品免费在线观看视频 | 日韩一区二区视频免费看| 成人毛片a级毛片在线播放| 国产精品欧美亚洲77777| 超碰97精品在线观看| 成年美女黄网站色视频大全免费 | 日本91视频免费播放| 婷婷色麻豆天堂久久| 18禁在线播放成人免费| 亚洲第一av免费看| 亚洲精品亚洲一区二区| 久久午夜综合久久蜜桃| 免费黄色在线免费观看| 国产免费福利视频在线观看| 一本大道久久a久久精品| 亚洲丝袜综合中文字幕| 天天操日日干夜夜撸| 男女国产视频网站| 丝瓜视频免费看黄片| 欧美三级亚洲精品| 一级毛片 在线播放| 下体分泌物呈黄色| 老司机影院毛片| 国产成人av激情在线播放 | 九九爱精品视频在线观看| 日本免费在线观看一区| 夫妻午夜视频| 精品亚洲成a人片在线观看| 中文欧美无线码| 中文乱码字字幕精品一区二区三区| 国产精品一区二区三区四区免费观看| 有码 亚洲区| 五月玫瑰六月丁香| √禁漫天堂资源中文www| tube8黄色片| 亚洲精品乱码久久久久久按摩| 亚洲综合精品二区| 成人综合一区亚洲| 汤姆久久久久久久影院中文字幕| 亚洲精品日本国产第一区| 国产精品久久久久久久电影| 美女内射精品一级片tv| 黄色一级大片看看| 精品少妇久久久久久888优播| 欧美激情 高清一区二区三区| 亚洲精品国产色婷婷电影| av专区在线播放| 亚洲欧洲日产国产| 日韩精品免费视频一区二区三区 | 久久这里有精品视频免费| 精品少妇久久久久久888优播| 国产淫语在线视频| 国产男人的电影天堂91| 中文字幕人妻熟人妻熟丝袜美| 成人国产麻豆网| 中文字幕制服av| 边亲边吃奶的免费视频| 国产精品嫩草影院av在线观看| 美女脱内裤让男人舔精品视频| 久久午夜福利片| 久久久久久久久久久免费av| 卡戴珊不雅视频在线播放| 国产在线免费精品| 亚洲欧美日韩卡通动漫| 国产精品99久久久久久久久| 青春草亚洲视频在线观看| 国产精品国产av在线观看| 国产男人的电影天堂91| 亚洲av成人精品一区久久| 菩萨蛮人人尽说江南好唐韦庄| 飞空精品影院首页| 亚洲色图 男人天堂 中文字幕 | av播播在线观看一区| 我要看黄色一级片免费的| 久久亚洲国产成人精品v| 日韩大片免费观看网站| 中文字幕免费在线视频6| 99热国产这里只有精品6| 亚洲欧美日韩另类电影网站| 亚洲av中文av极速乱| 国产老妇伦熟女老妇高清| 久久人妻熟女aⅴ| 性色avwww在线观看| 黄色配什么色好看| 丰满乱子伦码专区| 黄色欧美视频在线观看| 国产欧美另类精品又又久久亚洲欧美| 亚洲国产精品国产精品| 又大又黄又爽视频免费| 一级黄片播放器| 视频在线观看一区二区三区| 久久国内精品自在自线图片| 久久影院123| av专区在线播放| 精品久久蜜臀av无| 国产深夜福利视频在线观看| 母亲3免费完整高清在线观看 | 久久久久久久大尺度免费视频| 免费高清在线观看视频在线观看| 韩国高清视频一区二区三区| 国产精品国产三级专区第一集| 秋霞在线观看毛片| 国产精品秋霞免费鲁丝片| 亚洲精品久久成人aⅴ小说 | 亚洲精品一二三| 久久午夜福利片| 亚洲激情五月婷婷啪啪| 国产深夜福利视频在线观看| 91久久精品国产一区二区三区| 国产免费现黄频在线看| 最近的中文字幕免费完整| 日韩伦理黄色片| 精品久久久噜噜| 久久综合国产亚洲精品| 久久人妻熟女aⅴ| 免费看av在线观看网站| 精品午夜福利在线看| 久久精品国产亚洲av涩爱| 日日摸夜夜添夜夜爱| 大香蕉久久网| 美女cb高潮喷水在线观看| 最新的欧美精品一区二区| 在线亚洲精品国产二区图片欧美 | 99久久人妻综合| 亚洲欧洲日产国产| 国产精品一区二区在线观看99| 51国产日韩欧美| 成人国产麻豆网| 成人午夜精彩视频在线观看| 好男人视频免费观看在线| 亚洲av国产av综合av卡| 国产亚洲精品第一综合不卡 | 亚洲精品中文字幕在线视频| 亚洲av福利一区| 看免费成人av毛片| 简卡轻食公司| 视频区图区小说| 考比视频在线观看| 好男人视频免费观看在线| 2022亚洲国产成人精品| 久久久国产欧美日韩av| 日韩成人av中文字幕在线观看| 寂寞人妻少妇视频99o| 国产精品久久久久久精品古装| 777米奇影视久久| 国产探花极品一区二区| 男人爽女人下面视频在线观看| 日本午夜av视频| 男人添女人高潮全过程视频| 91精品国产国语对白视频| 久久免费观看电影| 日日啪夜夜爽| 中国国产av一级| 你懂的网址亚洲精品在线观看| 亚洲成人一二三区av| 国产男人的电影天堂91| 亚洲中文av在线| 国产黄色视频一区二区在线观看| 亚洲伊人久久精品综合| 亚洲国产精品999| 免费日韩欧美在线观看| 免费播放大片免费观看视频在线观看| 少妇的逼水好多| 午夜老司机福利剧场| 成人手机av| 亚洲欧美色中文字幕在线| 日韩强制内射视频| 91精品三级在线观看| 久久午夜福利片| 亚洲av欧美aⅴ国产| 亚洲人成77777在线视频| 午夜91福利影院| 自线自在国产av| 国产精品三级大全| 岛国毛片在线播放| 一本一本综合久久| 满18在线观看网站| 免费看av在线观看网站| 边亲边吃奶的免费视频| 99热全是精品| 国产极品粉嫩免费观看在线 | 人体艺术视频欧美日本| 亚洲人与动物交配视频| 日韩精品有码人妻一区| 日韩电影二区| 制服人妻中文乱码| 波野结衣二区三区在线| 国产精品99久久99久久久不卡 | 国产亚洲一区二区精品| 亚洲av.av天堂| 91精品一卡2卡3卡4卡| 亚洲成人手机| 亚洲一区二区三区欧美精品| 自线自在国产av| 黑人猛操日本美女一级片| 国产成人精品福利久久| 精品久久久精品久久久| 日韩av不卡免费在线播放| 欧美成人精品欧美一级黄| 在线观看三级黄色| 亚洲久久久国产精品| videossex国产| 交换朋友夫妻互换小说| 国产精品久久久久久精品古装| 一本—道久久a久久精品蜜桃钙片| 亚洲精品一二三|