• <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
    我們的作品
    1024香蕉在线观看| 色综合站精品国产| 91精品三级在线观看| 他把我摸到了高潮在线观看| 国产又爽黄色视频| 在线视频色国产色| 国产野战对白在线观看| 午夜影院日韩av| 无遮挡黄片免费观看| 波多野结衣高清无吗| 神马国产精品三级电影在线观看 | 国产精品二区激情视频| a级片在线免费高清观看视频| 另类亚洲欧美激情| 91麻豆精品激情在线观看国产 | 丝袜在线中文字幕| 少妇粗大呻吟视频| 国产99白浆流出| 999精品在线视频| 欧美日韩黄片免| 精品国产超薄肉色丝袜足j| 中文亚洲av片在线观看爽| 亚洲五月色婷婷综合| 黄网站色视频无遮挡免费观看| 亚洲欧美激情在线| 激情视频va一区二区三区| 精品一区二区三区视频在线观看免费 | 欧美激情久久久久久爽电影 | 久久婷婷成人综合色麻豆| 久久久久久人人人人人| 男人舔女人的私密视频| 99久久人妻综合| 亚洲少妇的诱惑av| 乱人伦中国视频| 色在线成人网| 精品国产一区二区三区四区第35| av欧美777| 国产欧美日韩一区二区三| 在线免费观看的www视频| 国产一区二区三区综合在线观看| 久久香蕉国产精品| av有码第一页| 大陆偷拍与自拍| 欧美在线一区亚洲| 麻豆av在线久日| 久久午夜亚洲精品久久| 69精品国产乱码久久久| 18禁黄网站禁片午夜丰满| 女人精品久久久久毛片| 久久久国产一区二区| 天堂影院成人在线观看| 怎么达到女性高潮| 精品卡一卡二卡四卡免费| 日韩国内少妇激情av| 国产精品1区2区在线观看.| 亚洲中文字幕日韩| 黄色 视频免费看| 精品人妻1区二区| 国产xxxxx性猛交| 国产免费av片在线观看野外av| 欧美老熟妇乱子伦牲交| 天堂√8在线中文| 午夜福利在线免费观看网站| 一夜夜www| 午夜福利在线观看吧| 如日韩欧美国产精品一区二区三区| 三级毛片av免费| 中文字幕色久视频| 法律面前人人平等表现在哪些方面| 黑人欧美特级aaaaaa片| 久久国产亚洲av麻豆专区| 欧美乱色亚洲激情| 久热爱精品视频在线9| 欧美亚洲日本最大视频资源| 精品电影一区二区在线| 亚洲成人精品中文字幕电影 | 久久天躁狠狠躁夜夜2o2o| 成年人黄色毛片网站| 色婷婷av一区二区三区视频| 校园春色视频在线观看| xxx96com| 中文字幕高清在线视频| 亚洲va日本ⅴa欧美va伊人久久| 一级a爱视频在线免费观看| 黄网站色视频无遮挡免费观看| 自线自在国产av| 国产又爽黄色视频| 亚洲第一青青草原| 天天添夜夜摸| 麻豆av在线久日| 精品一区二区三区四区五区乱码| 大码成人一级视频| 最近最新中文字幕大全电影3 | 国产黄a三级三级三级人| 国产精品av久久久久免费| 中亚洲国语对白在线视频| 久久久国产成人免费| 在线av久久热| 正在播放国产对白刺激| 国产1区2区3区精品| 757午夜福利合集在线观看| 成人特级黄色片久久久久久久| 乱人伦中国视频| 黄色女人牲交| 久久香蕉国产精品| 国产成人精品久久二区二区91| 久久性视频一级片| 亚洲免费av在线视频| 男女高潮啪啪啪动态图| 日日爽夜夜爽网站| 欧美日韩瑟瑟在线播放| 亚洲 欧美 日韩 在线 免费| 久久久久亚洲av毛片大全| 国产av一区二区精品久久| 国产一区二区在线av高清观看| 欧美日韩亚洲综合一区二区三区_| 18美女黄网站色大片免费观看| 欧美人与性动交α欧美精品济南到| 人人妻人人澡人人看| 熟女少妇亚洲综合色aaa.| 怎么达到女性高潮| 久久精品亚洲熟妇少妇任你| 丰满饥渴人妻一区二区三| 久久人人精品亚洲av| 国产一卡二卡三卡精品| 大陆偷拍与自拍| 这个男人来自地球电影免费观看| 18禁黄网站禁片午夜丰满| 日韩欧美在线二视频| 亚洲 欧美一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲五月色婷婷综合| 女人高潮潮喷娇喘18禁视频| 国产精品久久久人人做人人爽| 午夜福利免费观看在线| 亚洲精华国产精华精| 国产精品秋霞免费鲁丝片| 在线播放国产精品三级| 亚洲av片天天在线观看| 亚洲国产欧美网| 色播在线永久视频| 午夜日韩欧美国产| 黑人操中国人逼视频| 亚洲五月色婷婷综合| 精品人妻1区二区| 一进一出好大好爽视频| 99在线人妻在线中文字幕| 最近最新中文字幕大全免费视频| 国产精华一区二区三区| 欧美丝袜亚洲另类 | 自线自在国产av| 午夜福利一区二区在线看| 成人免费观看视频高清| 日韩av在线大香蕉| 欧美色视频一区免费| 91麻豆av在线| 看免费av毛片| a级毛片黄视频| 亚洲成人久久性| 亚洲精品久久成人aⅴ小说| 90打野战视频偷拍视频| 最近最新中文字幕大全免费视频| 一区在线观看完整版| 久久久国产一区二区| 国产成人免费无遮挡视频| 国产精品免费视频内射| 国产av在哪里看| 欧美老熟妇乱子伦牲交| 成在线人永久免费视频| a在线观看视频网站| 性少妇av在线| 黄片播放在线免费| 欧美 亚洲 国产 日韩一| 黄色 视频免费看| 男女下面进入的视频免费午夜 | 熟女少妇亚洲综合色aaa.| 久久 成人 亚洲| 久久国产精品影院| 最好的美女福利视频网| 亚洲视频免费观看视频| 欧美日韩亚洲综合一区二区三区_| 午夜免费激情av| 亚洲第一欧美日韩一区二区三区| 久久人人爽av亚洲精品天堂| 亚洲专区字幕在线| ponron亚洲| 午夜a级毛片| 久久久久久久午夜电影 | 老司机在亚洲福利影院| 亚洲av成人不卡在线观看播放网| 国产精品香港三级国产av潘金莲| av有码第一页| 日韩高清综合在线| 国产97色在线日韩免费| 精品久久久久久久久久免费视频 | 欧美av亚洲av综合av国产av| 国产三级黄色录像| av超薄肉色丝袜交足视频| www国产在线视频色| netflix在线观看网站| 国产欧美日韩一区二区三区在线| 国产成年人精品一区二区 | 在线观看免费视频日本深夜| 精品国产一区二区三区四区第35| 99精国产麻豆久久婷婷| 国产精品美女特级片免费视频播放器 | 91成年电影在线观看| 真人一进一出gif抽搐免费| 精品第一国产精品| 国产精品1区2区在线观看.| 一a级毛片在线观看| 男女之事视频高清在线观看| 性少妇av在线| 亚洲国产欧美网| 久久欧美精品欧美久久欧美| 亚洲熟妇熟女久久| 男人操女人黄网站| 国产日韩一区二区三区精品不卡| 久久国产精品影院| 久热这里只有精品99| 交换朋友夫妻互换小说| 国产亚洲欧美精品永久| 日韩 欧美 亚洲 中文字幕| 丝袜美腿诱惑在线| 日韩一卡2卡3卡4卡2021年| 首页视频小说图片口味搜索| 欧美久久黑人一区二区| 一边摸一边抽搐一进一小说| 欧美另类亚洲清纯唯美| 久久影院123| 窝窝影院91人妻| 亚洲欧美一区二区三区久久| 水蜜桃什么品种好| 欧美亚洲日本最大视频资源| 欧美日本亚洲视频在线播放| 一进一出好大好爽视频| 亚洲av电影在线进入| 久久亚洲真实| 久久香蕉激情| 日韩一卡2卡3卡4卡2021年| 国产精品 欧美亚洲| 少妇 在线观看| 国产亚洲av高清不卡| 天天躁夜夜躁狠狠躁躁| 国产精品国产av在线观看| 亚洲熟妇熟女久久| 怎么达到女性高潮| 精品国产国语对白av| 精品日产1卡2卡| 亚洲av成人一区二区三| av中文乱码字幕在线| 叶爱在线成人免费视频播放| 国产一区二区三区在线臀色熟女 | 一进一出好大好爽视频| 最近最新免费中文字幕在线| 最近最新中文字幕大全电影3 | 亚洲第一av免费看| 大型av网站在线播放| 老司机午夜福利在线观看视频| 美国免费a级毛片| 老司机靠b影院| 亚洲久久久国产精品| 国产av又大| 久久 成人 亚洲| 波多野结衣高清无吗| 亚洲熟妇熟女久久| 欧美日韩一级在线毛片| 色综合婷婷激情| 亚洲自拍偷在线| av在线天堂中文字幕 | 人人妻人人添人人爽欧美一区卜| 国产真人三级小视频在线观看| 欧美日韩乱码在线| av国产精品久久久久影院| 国产xxxxx性猛交| 午夜福利在线观看吧| 国产人伦9x9x在线观看| 欧美在线黄色| 最近最新免费中文字幕在线| 窝窝影院91人妻| 亚洲全国av大片| 亚洲人成电影观看| 亚洲精品久久午夜乱码| 午夜免费成人在线视频| 国产精品久久久久久人妻精品电影| 99久久综合精品五月天人人| 久久久久国产精品人妻aⅴ院| 80岁老熟妇乱子伦牲交| 国产欧美日韩一区二区精品| 欧美乱码精品一区二区三区| 91精品三级在线观看| 国产成人av教育| 久久精品国产综合久久久| 亚洲精品一区av在线观看| 久久中文字幕一级| 女人被狂操c到高潮| 亚洲熟女毛片儿| 高清av免费在线| 久久国产精品人妻蜜桃| 国产欧美日韩一区二区三区在线| 中文字幕人妻丝袜制服| 日本三级黄在线观看| 看黄色毛片网站| 久久中文看片网| 成人手机av| 脱女人内裤的视频| 男男h啪啪无遮挡| 亚洲欧洲精品一区二区精品久久久| 免费在线观看黄色视频的| 久久久久国内视频| 在线观看日韩欧美| 别揉我奶头~嗯~啊~动态视频| 三级毛片av免费| 最近最新免费中文字幕在线| 国产av在哪里看| 岛国在线观看网站| 日日摸夜夜添夜夜添小说| 两个人看的免费小视频| 国产av一区在线观看免费| 级片在线观看| 91大片在线观看| 久久人人爽av亚洲精品天堂| 国产成人av激情在线播放| 9色porny在线观看| 真人做人爱边吃奶动态| 久久久久久久久中文| 我的亚洲天堂| 免费一级毛片在线播放高清视频 | 在线观看免费视频日本深夜| 成人亚洲精品一区在线观看| 国产精品久久电影中文字幕| 亚洲中文字幕日韩| 一边摸一边抽搐一进一小说| 亚洲在线自拍视频| 国产亚洲精品第一综合不卡| 视频在线观看一区二区三区| 国产乱人伦免费视频| 日本精品一区二区三区蜜桃| 亚洲欧美激情综合另类| 手机成人av网站| 香蕉丝袜av| 免费久久久久久久精品成人欧美视频| 久久久国产欧美日韩av| 女人精品久久久久毛片| 黄色 视频免费看| 日韩大码丰满熟妇| 精品国产国语对白av| 国产亚洲欧美在线一区二区| 国产激情久久老熟女| 不卡一级毛片| www.www免费av| 色综合婷婷激情| 国产一区在线观看成人免费| 久久草成人影院| 精品一区二区三区视频在线观看免费 | 大型黄色视频在线免费观看| 丝袜美腿诱惑在线| 精品久久久久久久久久免费视频 | 波多野结衣一区麻豆| 亚洲精品av麻豆狂野| 88av欧美| 欧美一级毛片孕妇| 在线观看一区二区三区激情| 真人一进一出gif抽搐免费| 亚洲五月色婷婷综合| 久久香蕉国产精品| 法律面前人人平等表现在哪些方面| 每晚都被弄得嗷嗷叫到高潮| 国产精品九九99| 窝窝影院91人妻| 亚洲第一av免费看| 69精品国产乱码久久久| 国产精品久久久久成人av| 日韩视频一区二区在线观看| 成人影院久久| 国产野战对白在线观看| 国产黄色免费在线视频| 午夜激情av网站| 99在线人妻在线中文字幕| 亚洲一区高清亚洲精品| 每晚都被弄得嗷嗷叫到高潮| 视频区图区小说| 欧美日韩av久久| 男女床上黄色一级片免费看| 亚洲国产中文字幕在线视频| 美女午夜性视频免费| 男人的好看免费观看在线视频 | 在线国产一区二区在线| 欧美日韩国产mv在线观看视频| 91精品国产国语对白视频| 色综合婷婷激情| 在线国产一区二区在线| 777久久人妻少妇嫩草av网站| 久久人人97超碰香蕉20202| 丝袜在线中文字幕| 天天影视国产精品| 国产深夜福利视频在线观看| 成人av一区二区三区在线看| 亚洲美女黄片视频| 精品人妻1区二区| 在线观看免费视频网站a站| 免费av毛片视频| 美女高潮到喷水免费观看| 麻豆av在线久日| av片东京热男人的天堂| 欧美久久黑人一区二区| 亚洲国产欧美网| 在线观看66精品国产| 欧美av亚洲av综合av国产av| 在线观看免费日韩欧美大片| 99国产精品一区二区三区| 另类亚洲欧美激情| 国产97色在线日韩免费| 精品一区二区三区视频在线观看免费 | 亚洲精品美女久久av网站| 首页视频小说图片口味搜索| 亚洲精品在线美女| 性色av乱码一区二区三区2| 欧美日韩瑟瑟在线播放| av有码第一页| 黑人猛操日本美女一级片| 国产精品久久久人人做人人爽| av福利片在线| 欧美av亚洲av综合av国产av| 亚洲男人的天堂狠狠| 99国产极品粉嫩在线观看| 国产精品秋霞免费鲁丝片| 亚洲精品成人av观看孕妇| 叶爱在线成人免费视频播放| 老熟妇乱子伦视频在线观看| 国产一区二区激情短视频| 欧美日韩亚洲国产一区二区在线观看| 婷婷精品国产亚洲av在线| 操美女的视频在线观看| 丝袜在线中文字幕| 亚洲精华国产精华精| 看片在线看免费视频| 久久精品人人爽人人爽视色| 人妻丰满熟妇av一区二区三区| 亚洲第一av免费看| 久久这里只有精品19| 波多野结衣av一区二区av| 男人操女人黄网站| 9热在线视频观看99| 90打野战视频偷拍视频| 色综合婷婷激情| 亚洲成av片中文字幕在线观看| 精品日产1卡2卡| 神马国产精品三级电影在线观看 | 操美女的视频在线观看| 激情视频va一区二区三区| av福利片在线| 精品少妇一区二区三区视频日本电影| 久久人妻熟女aⅴ| 成人永久免费在线观看视频| 五月开心婷婷网| 久久影院123| 日韩欧美国产一区二区入口| 国产亚洲欧美98| 女性被躁到高潮视频| 亚洲精品久久成人aⅴ小说| 国产欧美日韩综合在线一区二区| 不卡av一区二区三区| 国产精品久久久人人做人人爽| 日本wwww免费看| 巨乳人妻的诱惑在线观看| 在线观看免费日韩欧美大片| 高清欧美精品videossex| 日本黄色日本黄色录像| 伦理电影免费视频| 亚洲精品一区av在线观看| 午夜老司机福利片| 精品一区二区三区视频在线观看免费 | 日韩视频一区二区在线观看| 国产99久久九九免费精品| av有码第一页| 最近最新中文字幕大全电影3 | 操美女的视频在线观看| 欧美亚洲日本最大视频资源| 超色免费av| 纯流量卡能插随身wifi吗| 视频在线观看一区二区三区| 久久精品亚洲熟妇少妇任你| videosex国产| 欧美国产精品va在线观看不卡| 亚洲人成77777在线视频| 一级毛片女人18水好多| av视频免费观看在线观看| 欧美大码av| 男女下面进入的视频免费午夜 | 黄色a级毛片大全视频| 日本 av在线| 久久久久久亚洲精品国产蜜桃av| 亚洲片人在线观看| 亚洲欧美精品综合久久99| 嫩草影视91久久| 中国美女看黄片| 999久久久精品免费观看国产| 国产精品久久视频播放| 国产精品国产高清国产av| 亚洲精品美女久久av网站| 久久精品亚洲精品国产色婷小说| 亚洲精品国产精品久久久不卡| 正在播放国产对白刺激| 精品一区二区三卡| 久9热在线精品视频| 丝袜在线中文字幕| 两性夫妻黄色片| 夜夜看夜夜爽夜夜摸 | 欧美中文综合在线视频| 亚洲少妇的诱惑av| 久久天躁狠狠躁夜夜2o2o| ponron亚洲| 久久久国产成人免费| 亚洲国产精品999在线| 18禁黄网站禁片午夜丰满| 亚洲熟女毛片儿| 视频区图区小说| 丝袜在线中文字幕| 国产黄色免费在线视频| 韩国av一区二区三区四区| 中国美女看黄片| 精品福利观看| 久久人妻熟女aⅴ| 日本黄色视频三级网站网址| 久久天躁狠狠躁夜夜2o2o| 日韩国内少妇激情av| 青草久久国产| 夜夜看夜夜爽夜夜摸 | 午夜老司机福利片| 99国产精品一区二区三区| 成人手机av| 国产99白浆流出| 如日韩欧美国产精品一区二区三区| 久久久久国产精品人妻aⅴ院| 国产av又大| 999久久久精品免费观看国产| 中亚洲国语对白在线视频| 久久久久久人人人人人| 亚洲七黄色美女视频| 成人18禁高潮啪啪吃奶动态图| 免费在线观看影片大全网站| 亚洲av熟女| 日韩大尺度精品在线看网址 | 午夜免费激情av| 高潮久久久久久久久久久不卡| 老司机亚洲免费影院| 91字幕亚洲| 中亚洲国语对白在线视频| 国产高清videossex| 色在线成人网| 亚洲人成电影免费在线| 亚洲欧美日韩另类电影网站| 啪啪无遮挡十八禁网站| 国产亚洲精品综合一区在线观看 | videosex国产| 精品高清国产在线一区| 亚洲精品美女久久av网站| 正在播放国产对白刺激| 在线观看一区二区三区| 91麻豆av在线| 成人av一区二区三区在线看| 1024视频免费在线观看| 黑人操中国人逼视频| 亚洲精品粉嫩美女一区| 亚洲av成人一区二区三| videosex国产| 真人做人爱边吃奶动态| 一边摸一边抽搐一进一出视频| av天堂在线播放| 大型av网站在线播放| 国产成+人综合+亚洲专区| 十八禁网站免费在线| 亚洲精华国产精华精| 欧美日韩瑟瑟在线播放| 国产亚洲精品久久久久5区| 国产免费男女视频| 嫩草影院精品99| 欧美日韩国产mv在线观看视频| 久久久久久久午夜电影 | 欧美日本亚洲视频在线播放| 欧美丝袜亚洲另类 | 精品久久久久久,| 亚洲成国产人片在线观看| 亚洲国产看品久久| 欧美成人免费av一区二区三区| 乱人伦中国视频| 国产精品美女特级片免费视频播放器 | 国产精品九九99| 日韩精品中文字幕看吧| 麻豆国产av国片精品| 国产亚洲精品第一综合不卡| 精品久久久久久久久久免费视频 | 久久午夜综合久久蜜桃| 亚洲人成77777在线视频| 丰满饥渴人妻一区二区三| 亚洲一区高清亚洲精品| 精品乱码久久久久久99久播| 中文字幕精品免费在线观看视频| 欧美精品一区二区免费开放| 美女大奶头视频| 可以在线观看毛片的网站| 久久久久久久久久久久大奶| 久久人人97超碰香蕉20202| 久久人妻av系列| 中文字幕另类日韩欧美亚洲嫩草| 成人特级黄色片久久久久久久| 中出人妻视频一区二区| 日韩三级视频一区二区三区| 国产精品1区2区在线观看.| 久久久久久亚洲精品国产蜜桃av| 88av欧美| a级毛片黄视频| 人人妻人人澡人人看| 欧美午夜高清在线|