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

    Selected topics of quantum computing for nuclear physics?

    2021-03-11 08:31:44DanBoZhang張旦波HongxiXing邢宏喜HuiYan顏輝EnkeWang王恩科andShiLiangZhu朱詩(shī)亮
    Chinese Physics B 2021年2期

    Dan-Bo Zhang(張旦波), Hongxi Xing(邢宏喜), Hui Yan(顏輝),Enke Wang(王恩科), and Shi-Liang Zhu(朱詩(shī)亮),?

    1Guangdong Provincial Key Laboratory of Quantum Engineering and Quantum Materials,School of Physics and Telecommunication Engineering,South China Normal University,Guangzhou 510006,China

    2Guangdong–Hong Kong Joint Laboratory of Quantum Matter,Frontier Research Institute for Physics,South China Normal University,Guangzhou 510006,China

    3Guangdong Provincial Key Laboratory of Nuclear Science,Institute of Quantum Matter,South China Normal University,Guangzhou 510006,China

    4Guangdong–Hong Kong Joint Laboratory of Quantum Matter,Southern Nuclear Science Computing Center,South China Normal University,Guangzhou 510006,China

    Keywords: quantum computing, nuclear physics, quantum field theory, quantum simulation, quantum algorithm

    1. Introduction

    Understanding how elementary particles form nuclear matter is fundamentally important. While the underling basic physical theory, namely, quantum chromodynamics (QCD),can be formulated concisely,it is notoriously hard to solve.[1]This is because interactions between quarks and gluons follow very distinct behaviors: quarks are asymptotically free at the high-energy scale,while strong couplings exist among quarks and gluons at low-energy. In this regard, QCD can be nonperturbative,leading to a failure of the perturbative expansion of Feynman diagrams which have gained remarkable success in quantum electrodynamics (QED). Lattice QCD has been developed,[2–4]which can rely on computational methods such as quantum Monto Carlo[5,6]and tensor networks.[7–14]The former suffers from the fermion sign problem. The tensor network approach expresses many-body wavefunction in a compressed way which is free of sign problem,but it may have an exponential growth of complexity for evaluating physical observations. While those numeral efforts are still pushing the limits on what quantum many-body problems can be solved on classical computers,they will meet some intrinsic difficulties that are guaranteed by computational complexity, which claims that quantum many-body problems are NP-hard.[6]

    To overcome the intrinsic difficulty of simulating a quantum world with classical computers,Feynman proposed quantum hardwares as simulation platforms in 1983.[15]Since then,simulation of physical systems has been an incentive for building quantum computers. Basically, quantum simulation can be divided into simulating static properties and real-time evolution of a physical system.[16]The corresponding quantum algorithms have been developed at the early age of quantum computing based on ideal quantum computers. One remarkable quantum algorithm is quantum phase estimation that can calculate the eigenenergies and access the corresponding eigenstates for quantum systems,[17]which has been used for solving molecular energies.[18]Recent scaling-up of quantum processors enables us to simulate static properties and real-time evolution of a quantum system to a larger size with noisy qubits. In such near-term noisy intermediate quantum(NISQ)era,[19]quantum algorithms are designed suitable for near-term quantum devices. Variational quantum algorithms are remarkable as candidates for fully exploiting the power of NISQ quantum computers. It receives special interests as it provides a practical approach for solving quantum chemistry,[20–25]quantum many-body systems,[26–28]and many other applications.[29–31]

    Nuclear physics can be considered as quantum manybody physics on and below the subatomic scale, and it is of great advantages in using quantum computers for ultimate solutions. However, simulation of nuclear physics meets some special challenges that are different than those of quantum chemistry and quantum many-body physics in condensed matter. Firstly,nuclear physics is studied via quantum field theory on continuous space–time, possible with infinite dimensional degrees of freedom per volume, and they should be encoded efficiently into an finite number of qubits on lattice. Moreover, the underling theory for nuclear physics should respect gauge-invariance, and in the corresponding quantum simulation,a physical Hilbert space should be kept. Those two challenges not only call for a formulation of lattice gauge theory for the purpose of efficient quantum computing,[32,33]but also the issue of realization of reliable simulation on noisy quantum processors which becomes tremendous important.[32,34–36]

    Efforts of simulating nuclear physics often benefit from advances on quantum simulations of condensed matter or electronic structures of molecules. There are different approaches for tackling this problem, and we may loosely divide it into two types: analog and digital quantum simulations, which may be developed with different motivations. Analog quantum simulation stands for early attempts for simulating gauge lattice theory, when people have realized that cold atoms can simulate quantum many-body problems in condensed matter. Different platforms for quantum simulation of lattice gauge field have been proposed, including cold atoms,[37–48]trapped ions,[27,49–53]superconducting circuit,[54,55]and Rydberg atoms.[56–60]Existing reviews can refer to Refs. [45,61]for details, and related artificial gauge fields for ultracold atoms were reviewed in Refs.[62,63]. Those proposals stress on how to realize the dynamical gauge field with controllable atomic interactions. It is realized that all proposals would inevitably require very complicated setups,suggesting that simulation of gauge field and nuclear physics is in general difficult and rather resource consuming. To date, a building block for simulating a coupling of gauge-field with fermion has been implemented in cold-atom experiment.[48]Remarkably, the local gauge invariant has been observed experimentally in a 71-site Bose–Hubbard quantum simulator.[36]Those advances pave the way for realistic analog quantum simulation of gauge lattice theory.

    The approach of digital quantum simulation for gauge field theory and nuclear physics begins later but receives more attention recently.[27,32,33,51,64–78,78–90]The digital way is programmable by compiling all operators into basic quantum gates and thus is much more flexible for simulating different quantum systems.[17]In 2012, Preskill et al. proved an exponential quantum speed-up in simulation of scattering problem for scalar relativistic quantum field theory with self-interactions (φ4theory),[64]with time complexity scaling up polynomial in the number of particles, their energy,and the desired precision, remarkably in the non-perturbative regime where classical algorithms fail to work. Although elaborated on a concrete model, it puts a solid base for simulating high-energy physics with quantum advantages from the aspect of computational complexity. More quantum algorithms have been proposed to simulate lattice gauge field and nuclear physics on near-term quantum processors. While lots of techniques can inherit from those of quantum computational chemistry, one key difference is the requirement to deal with gauge field and gauge-invariance. A very successful example is that gauge degree of freedom for the 1+1D Schwinger model can be eliminated.[27,51,91]With this reduction, simulations of ground state and real-time evolution of 1+1D Schwinger model have been demonstrated on trappedion platforms.[27,51]However,in general such an elimination is impossible, and systematic formulations of lattice gauge theory suitable on quantum computers have been investigated,especially on how to reduce the desired quantum resource for gauge field which has much redundancy. The approach from lattice gauge theory,however,is still too resource demanding for simulating nuclear physics interested in experimental observations. Another practical strategy, mostly pushed by the community of high-energy and nuclear physics,is to incorporate quantum algorithms as subroutines and wisdom of classical methods is exploited. This permits valuable quantum resources concentrated on classical difficulty parts, in order to exploit near-term noisy quantum processors to solve complicated nuclear physics. Such a hybrid quantum–classical strategy has been used in studying parton distributions,[68]evolution of non-equilibrium thermal states such as quark–gluon plasma,[66]and so on.

    Although sufficient advances have been made, quantum simulation of high-energy nuclear physics is still a young field,and it is just beginning to bring more researchers into this field.In this review,we aim to summarize recent advances,give basic concepts on quantum computing and quantum algorithms,introduce several representative works,and propose future directions. The review will be organized as follows.

    In Section 2, we first discuss how nuclear physics problems can be mapped into formulas of lattice field theory expressed with qubits,which can be solved on a quantum computer. Specifically, we reveal that the way of treating gauge field may be the key ingredient. We then introduce quantum algorithms related for simulating quantum systems,including both static properties and real-time evolution.

    In Section 3,,we give some specific topics and representative examples on the applications of quantum computing for nuclear physics. We first elaborate on a prototypical example,the 1+1D Schwinger model. This model,although describing quantum electrodynamics (QED), shares lots of key features with QCD, such as confinement and a topological theta vacuum. Thus, this minimal model can be used as a testbed for simulating lattice gauge field. We also give examples on the hybrid quantum-classical approach which uses quantum algorithms as a subroutine, including examples of parton physics and evolution of non-equilibrium thermal states.

    Finally,in Section 4,we give outlooks and summaries.

    2. Nuclear physics on a quantum computer

    In this section, we first discuss how to map nuclear physics problems onto a quantum computer,by reformulating gauge quantum field in the language of qubits. Then, we introduce quantum algorithms for solving static and dynamical properties of nuclear physics.

    2.1. Map nuclear physics onto a quantum computer

    The underling theory for nuclear physics is quantum field theory describing fermionic matter coupled with bosonic gauge fields on a continuous space–time background. However, a quantum computer consists of arrays of qubits, and each qubit owns 2-dimensional Hilbert space. To simulate nuclear physics on a quantum computer, it demands for a mapping of the original physical degree of freedom onto qubits at the first stage.

    On a quantum computer, the basic unit to encode information is quantum bit, which can store a superposed state of|0〉and|1〉. In general,information is expressed as a quantum state of multiple qubtis. With several hundreds of qubits, a quantum computer can store a quantum state that is beyond the capacity of classical computers. The quantum state is manipulated on a quantum computer with unitary evolutions,which can be decoupled into a sequence of quantum gates from a universal set of basic quantum gates. Then, information or interested quantities can be extracted by repeatedly measuring the final quantum states.Further post-processing on a classical computer may be required. The physical process and its corresponding simulation on a quantum computer are illustrated in Fig.1.

    Nuclear physics includes problems related to steady states,such as structure of nucleons and phases of matter,and real-time evolution, such as scattering problem and evolution of nuclear matter. While preparation of steady states and simulating the real-time evolution have been standard techniques in quantum computing and variants of approaches have been developed,it is important to firstly convert nuclear physics into formulas that a quantum computer can handle with. Keep in mind that basic components of a quantum computer are qubits for encoding information, a universal set of single- and twoqubit quantum gates for manipulating information, and measurements for extracting information.To study nuclear physics on a quantum computer, those field operators should be converted into qubit operators,and the space–time should be discretized, and some symmetries should be considered, especially the local gauge invariance(general Gauss law). We discuss them separately.

    Fig.1. Illustration of solving nuclear physics on a quantum computer. The upper is a physical process of scattering in the experiment of Large Hadron Collide(from https://www.symmetrymagazine.org/article/lhc-scientists-detect-most-favored-higgs-decay),while the bottom depicts a corresponding simulation with a quantum computer.

    Encode fields into qubits on latticeQuantum fields usually are described on a continuous space–time and they can be bosonic (infinite dimension) or fermionic, while qubits on a quantum computer are defined on a lattice with each site occupied by one qubit. To encode the quantum state of a quantum field on a quantum computer,we should discretize the space as a lattice,encode each local bosonic field into a finite number of qubits by a truncation,and encode a local fermionic field into a sequence of qubits by nonlocal transformations. Remarkably,to reduce the computational resource yet retain the essential physics,the discretization of space and the cut-off bosonic field should respect some physical arguments(such a consideration also exists for solving lattice gauge field on a classical computer). In other words, the original nuclear physics problem may have an infinite degrees of freedom and it is important to map it into a task solvable on a quantum computer with a finite number of qubits under desired precision.[64,69]Another important ingredient in the encoding stage is how to express gauge field with qubits. It is important to enforce local gauge-invariance(general Gauss law)for simulating lattice gauge theory on a quantum computer,[14,32–34]which corresponds to keep the simulation in the physical Hilbert space.Notably,a promising digitization scheme for the gauge field is discrete subgroup approximations,[92,93]which reduce the required number of qubits sufficiently by encoding continuous gauge filed with discrete subgroups.

    Unitary evolutionPhysical processes in nuclear physics can be both unitary and nonunitary,while states on a quantum computer evolve in a unitary process. For an unitary evolution such as evolutions of the Hamiltonian,one can decompose the unitary operator into a sequence of basic single-and two-qubit quantum gates with product-formula,by cutting the evolution time into short time periods. As an nonunitary process can be embedded in an unitary process of a larger system, one can introduce some auxiliary qubits to implement the nonunitary evolution. A common applied algorithmic design is by linear-combination-of-unitaries, using ancillary qubits[94–99]or continuous-variables.[99–103]

    MeasurementsIn nuclear physics, we often need to get distributions of particles, such as in the scattering problems.On a quantum computer, however, measurements often refer to computational basic(σzmeasurements),which may not directly correspond to find particles in specified directions or momentums. One can synthesize such physical measurements by unitary evolutions followed by measurements on computational basis.

    Digital vs. analog quantum computerIt should be reminded of a difference between digital and analog quantum computers. Analog quantum computers are designed specifically for simulating quantum systems. Unlike a digital quantum computer that uses an universal set of quantum gates to construct all unitary process,an analog quantum computer can directly simulate a Hamiltonian and its real-time evolution,by engineering the desired interactions. This limits its capacity for simulating different systems, while on the other hand, it can be easier to implement. Representative physical systems for analog quantum computers include cold atoms,[62]trapped ions,and superconducting circuits.Nevertheless,analog quantum simulator can be programmable, enlarging its capacity of simulating a different target Hamiltonian with a controllable Hamiltonian on the physical platform using variational methods. For instance, ground state of the 1+1D Schwinger model has been simulated in the trapped-ion system with longinteraction between qubits.[27]

    Continuous-variable approachGiven the great challenge of digitizing the gauge field, is it possible to directly encode the gauge (bosonic) field into continuous variables?There is a proposal for simulating bosonic systems with continuous-variable (CV) quantum computing for the Bose–Hubbard model.[104]Moreover,the hybrid discrete-continuous variables approach of quantum computing seems very natural for simulating the gauge theory,[102,103,105–109]since it is a system of hybrid variable due to a coupling between fermion matter and gauge field. Notably,the current mainstream quantum computing platforms often have naturally existing continuous variables,such as motional modes in trapped ions[109–111]and cavity modes in superconducting circuits,[112,113]and their interplay with qubits have been exploited in quantum information processing. However, the continuous-variable quantum computing approach can be more sensitive to quantum noisy.One solution is to encode a CV mode into many CV modes and use non-Gaussian resource (such as a qubit encoded into a CV mode)to correct Gaussian error.[114]In this regard, the approach of encoding the gauge field directly into continuous variables deserves for exploration,but more efforts on the experiments are awaited to demonstrate the continuous variable quantum computing at the level of quantum algorithm.

    2.2. Quantum algorithms for simulation

    To exploit the power of quantum computers for simulating nuclear physics, quantum algorithms are indispensable.After having converting a nuclear physics problem on a quantum computer, the remaining task is to design a quantum algorithm to solve it efficiently. Although quantum algorithms can be counter-intuitive and require special efforts to learn in general,it can be very natural for physicists when using quantum algorithms for solving quantum systems. Those quantum algorithms may be classified as two kinds: preparing steady states of a quantum system and simulating the real-time evolution. In the following,we give a brief introduction.

    For an desired precision ε,the time complexity is O(t2/ε). It is interesting to note that while such a formula is simple, it is powerful even on the near-term quantum devices, and the performance can be comparable with more advanced quantum algorithms for Hamiltonian simulations.[115–117]Remarkably,it is found that time-complexity can be reduced to O(t/ε)using the product-formula when the Hamiltonian can be written as H =HA+HB, where local terms in HAor HBcommute to each other but[HA,HB]/=0. This is because there is a destructive error interference that errors will cancel in different short time periods.[118]This indicates that the product-formula may be more efficient than expected.

    Solving steady states, including eigenstates and thermal states,is important for understanding static properties of quantum systems, but can be formidable challenge on a classical computer.Such a basic task is comparatively harder on a quantum computer than simulation of real-time evolution,and lots of efforts have been devoted into it. The first textbook algorithm is quantum phase estimation,[17]where eigenvalues of a Hamiltonian are written onto ancillary qubits and then the ancillary qubits are measured to read the eigenvalue as well as the associated eigenstates. Concretely, consider H|un〉=En|un〉and the initial state is|ψ〉=∑ncn|un〉,then each eigenstate can be obtained with a probability|cn|2. Thus it is important that the target eigenstate should have a large weighting in the initial state. Quantum phase estimation can be understood as signal processing from the time domain to the frequency domain. It involves Hamiltonian evolution at different time periods, and then uses quantum Fourier transformation to find eigenvalues at frequency domain. Both Hamiltonian evolution and quantum Fourier transformation are resource-costing, and consequently,quantum phase estimation is not suitable on the nearterm quantum devices. Quantum adiabatic algorithm provides an approach to prepare the ground state of a Hamiltonian.[119]It starts from a Hamiltonian whose ground state is easy to prepare,and by adiabatically tuning the Hamiltonian into the target Hamiltonian,the final state will be the ground state for the target Hamiltonian.The performance of the quantum adiabatic algorithm relies on a gap between the first excited state and the ground state along the adiabatic path. Thus it requires a good design for the adiabatic path which can be nontrivial. In practice, the adiabatic evolution should be very slow, and may be hard to finish within the coherent time. Thus, quantum adiabatic algorithm,although can be universal,may not be suitable in the era of NISQ.

    Another important class of quantum algorithms is variational quantum eigensolver (VQE),[20,21,23,24,26,27,120–122]which is considered promising for fully exploiting the power of NISQ quantum devices. It can rely on a shallow quantum circuit and a moderate number of qubits for solving classical intractable problems. The quantum circuit is parameterized and the parameters can be obtained by minimizing the energy.The optimization is a hybrid quantum-classical one. In this sense, variational quantum algorithms give a new paradigm for algorithmic design: rather than directly designing a quantum algorithm,it trains a quantum algorithm for a given task,by optimizing a cost function. The VQE is designed for solving the ground state by minimizing the energy,and variants of VQE have been developed for obtaining excited states,[123–126]thermal states at finite temperature,[127–131]quantum imaginary time evolution,[132]and general quantum processes.[133]

    3. Applications for nuclear physics

    Although at an early stage, quantum computing has had a broad applications for nuclear physics problems, and it is beyond the scope of this review for a throughout investigation. Rather, we attend to use some prototypical examples to illustrate how nuclear physics problems can be solved on a quantum computer. With concrete examples, the basic concepts and procedures may be revealed.We start with the 1+1D Schwinger model, which is a prototypical model for simulating gauge field with a state-of-art realistic quantum computer. Other examples include interesting problems ranging from scattering, evolution of non-equilibrium thermal states,nuclear structure at both high-energy and low-energy, and so on. In addition,we introduce a work that gives the time complexity of solving scattering for a scalar quantum field,which shows quantum advantage.

    3.1. Simulation of real-time evolution and ground state of lattice gauge theory: the Schwinger model

    The capacity of quantum simulation of lattice gauge field implies that nuclear physics can be studied on a quantum computer.However,the required quantum resource is too demanding in the NISQ era,and it is more practical to start with some minimal models. The Schwinger model is such a prototypical model that has been demonstrated on current quantum processors experimentally, including both ground state[27]and realtime evolution.[51]Moreover,the Schwinger model can reveal some important features shared with QCD. We thus use this model to illustrate a work flow for studying nuclear physics on a quantum computer: how the original formula of a quantum field problem should be mapped into a lattice spin model,what quantum algorithm should be chosen to solve the lattice spin model, and how desired physical quantities can be accessed with measurements on a quantum computer, possible with post-processing.

    We introduce the Hamiltonian for the Schwinger model with fixing gauge A0(x)=0,

    By eliminating the gauge field, the Hamiltonian is a bit more complicated as there are long-range interactions, but is still feasible on current quantum processors. Remarkably,trapped-ion quantum computers are renowned for its excellent connectivity, which is naturally suitable for simulating quantum systems with long-range interactions, while other platforms may suffer from limited connectivity. The following two experiments on trapped ions present the state-of-art simulations for the lattice gauge field on the real-time evolution and ground state properties,respectively.

    Real-time evolutionMany nuclear physics phenomena involve time evolution and it is fundamentally important to simulate the real-time evolution of gauge theory. While existing classical methods often meet difficulties, simulation of real-time evolution can be implemented very naturally on a quantum computer, once a real-time evolution of a Hamiltonian is compiled into a sequence of quantum gates with Trotter decomposition. The Schwinger model was raised for explaining the mechanism of creating pairs of fermion–antifermion from the bare vacuum. Such a process can be simulated with a language of qubits. Firstly, occupation of a fermion (antifermion)is encoded as|1〉(|0〉)on the even(odd)site. Then,the bare vacuum has no fermion at all, which can be represented as|0101010101···〉. With time evolution of the Hamiltonian, pairs of fermion–antifermion will appear due to hopping terms in the Hamiltonian.To implement e?itHon a quantum computer,a Trotter formula can decompose the time evolution into a product of short-time evolution of small Hamiltonian term,as in Eq.(1). Proliferation of fermion–antifermion pairs after a period of t can be revealed by repeated measurements.

    In trapped ions,a qubit is encoded as two internal energy levels in an ion. The experiment uses four ions to demonstrate the creation of fermion–antifermion pairs on realistic trappedion platforms,and the result fits good with theoretical results,with both idea evolution and the case of using the Trotter decomposition. Although with only a few qubits, this work has paved the way for simulating dynamical systems of large-scale systems that basically follows the same scheme.

    Ground stateWhile real-time evolution for a system can be implemented directly on a quantum computer, solving its ground state or excited states requires more efforts on quantum algorithms. Variational quantum eigensolver provides a feasible scheme on the NISQ quantum processor. The key point is to train a quantum circuit to prepare a variational ground state |Ψ(θ)〉=U(θ)|Ψ0〉 for the target Hamiltonian H,where parameters are optimized by minimizing the energy E(θ)=〈Ψ(θ)|H|Ψ(θ)〉. For the Schwinger model,one challenge is to adopt a wavefunction ansatz with enough expressive power under limited quantum resource, e.g., the parametrized quantum circuit should be shallow and yet can describe complicated quantum correlation of the ground state due to the long-range interactions. The current digital quantum computer is limited to small sizes as the required depths of quantum circuit can increase quickly with the system size.On the other hand, analog quantum simulator can realize an unitary operator by directly letting an engineered Hamiltonian evolves, which may be decomposed into a long sequence of quantum gates. Reference[27]devises a variational quantum simulator that can simulate the target Hamiltonian with engineered different Hamiltonian.

    For simulating the Schwinger model with trapped-ions,an advantage is that the engineered Hamiltonian of trappedions naturally has long-range interactions,which can be written as

    The parameter vector θ is trained by minimizing the variational energy E(θ). The optimization can be nontrivial for a noisy quantum processor, and it adopts a global optimizer called DIRECT,[27]which divides the parameter space into cells for search. It is demonstrated experimentally that the trapped-ion simulator can variationally solve the ground state for up to 20 qubits for the Schwinger model, which is quite remarkable. Moreover,the experiment uses zero energy variance to self-verify that the obtained state is indeed the eigenstate,since an eigenstate can be characterized by zero energy variance. Although only shown with eight qubits in experiment, self-verifying can be vital for large-size quantum simulator that is beyond the computational capacity of classical computers: it requires the quantum computer itself to verify if it gives a result with enough accuracy. In addition,a quantum phase transition is revealed on the trapped-ion simulator by tuning the mass across mc≈?0.7 (set a=1 and g=1).This indicates that the variational quantum simulator can efficiently solve ground states for the Schwinger model for a large range of parameters, especially around the quantum critical point where quantum correlations are strong. The work opens a new direction for simulating large-size quantum systems exploiting analog quantum simulators with great flexibility by using the variational method.

    We comment that elimination of the gauge field completely is only possible for 1+1D U(1) gauge theory. Nevertheless,the idea of reducing quantum resources for representing gauge field is important and possible,since gauge field has a redundancy in description. Many efforts have been denoted for this mapping stage, including non-Abelian gauge field in 1+1D,[12,70]1+2D Abelian[72]and non-Abelian gauge fields,and so on.

    3.2. Nuclear structure

    One center motivation of nuclear physics is to explore the internal structure of nucleus,which relies on an increasing of energy scale of probes for resolving finer structures. Protons and neutrons are basic ingredients of a nucleus at low-energy scales,but a proton or a neutron itself is a collection of quarks and gluons at large energy scales. Colliders have been built for detecting those structures at different energy scales from scattering cross sections. Yet, numeral methods meet challenges that are common for quantum many-body systems,especially for solving the structure of hardrons,which in nature is non-perturbative. We review recent efforts for solving nuclear structure on a quantum computer,[68,135]for an atomic nucleus and an hadron respectively, which adopt quite different strategies.

    3.2.1. Binding energy of atomic nucleus

    From the aspect of low-energy nuclear physics, an atom consists of interacting protons and neutrons that are bounded together. Although lattice QCD can be applied,a better starting point is to use effective field theory (EFT) that protons and neutrons are relevant degrees of freedom. For light nucleus, pionless EFT provides systematically improvable approach for modeling nuclear interaction. For this bound-state problem,a second-quantization of the Hamiltonian can be obtained by using a proper local basis,e.g.,a common choice can be the harmonic oscillator basis. As a truncation of the local basis can be applied, the problem can be solved in a finitedimensional Hilbert space. Further, the second-quantization fermionic Hamiltonian can be mapped into a qubit Hamiltonian. Such a procedure is very familiar in quantum chemistry and quantum many-body problems in condensed matter. And one can immediately recognize that the problem is much like quantum chemistry,and it is expected that an exponential wall prevents classical methods when the system size is becoming large.

    We point out that solving the Hamiltonian in Eq. (6) is actually easy by exact diagonalization since it involves a tridiagonal matrix.Why bother to use a quantum computer?Moreover,a N×N matrix is easy to diagonalize on a classical computer for N ~100,but can be challenge for the current quantum devices. It should be reminded that the VQE approach becomes valuable when solving nuclear structure becomes a many-body problem, and a simple exact diagonalization will face a M ?M matrix(M grows exponentially with N). It is expected that a road map for extending the scheme for quantum many-body nuclear physics is required.

    3.2.2. Parton physics

    Solving internal structure of a hadron is quite a different story. The partonic structure of a hadron depends on the energy scale to see it. For instance, a proton is a bounded state of three valence quarks at low energy scale, and no individual quark has been observed in experiments due to color confinement. On the other hand, at higher energy scale in deep inelastic scattering (DIS) with large momentum transfer, the hadron shall be modeled as a collection of charged point-like constituents, namely, parton. The parton contributes to the cross section with momentum xpμ, where x is the momentum fraction of the proton carried by the parton,and pμis the momentum of the hadron. Understanding hadronic structure in terms of the parton distribution functions(PDFs)over xpμand its generalization is a research frontier.[136–139]

    The PDFs f(x)serve as an input for explaining scattering cross sections in DIS experiments. However, the parton distribution itself is hard to compute from ab initio methods due to its nonperturbative nature. There are different approaches from first principles, such as light-front Hamiltonian and lattice QCD,which still await for exascale supercomputer to verify. Here we introduce a recent work that points out an avenue on a near-term quantum computer.

    Reference[68]presents two different schemes for calculating the parton distribution on a quantum computer: direct calculation of PDFs using operator formula and by evaluation of the hadronic tensor. For illustration purpose, the two schemes are demonstrated with 1+1D Thirring model. The lattice Hamiltonian writes

    Similar to the Schwinger model, the bare vacuum is set at m →∞limit, and |1〉 (|0〉) on the even (odd) site represents occupation of a fermion(antifermion),respectively.

    The operator formula for a definition of PDF (without gauge link)is given by

    The second scheme exploits a connection between PDF and the hadronic tensor. The hadronic tensor for the ddimension theory is defined as correlations of currents Jμ,

    which can characterize the many-body wavefunction |P〉through correlations and consequently reveal the hadronic structure to some degree. The connection between PDF and Wμν(q)is established as(via collinear factorization)

    The remaining task is to evaluate Wμν(q)on the state|P〉,which can refer to the linear response method on a quantum computer. Still, the preparation of the state |P〉 is quite involved,and reference[68]pointed out that quantum adiabatic state preparation may be employed.

    The above proposal suggests a way for studying parton physics on a quantum computer, but a demonstration is still awaited. This requires to efficiently prepare |P〉 (possible with VQE),evaluating the hadronic tensor,and extracting PDF through global fitting to the results of hadronic tensor obtained from quantum computing. Another important question is which lattice Hamiltonian should be a good starting point. The light-front Hamiltonian of QCD looks promising as it already has provided a framework for ab initio calculation of proton structure. Remarkably,VQE based on the lightfront Hamiltonian has just been demonstrated for calculating the structure of pion.[78]

    3.3. Quantum advantage: Scattering in scalar quantum field theory

    Scattering is central to nuclear physics as it is almost the only available experimental method. Consequently, calculating scattering amplitude is what in theory needed to do for a comparison with experiments. Different machineries of computational methods have been developed but there are some fundamental challenges. The difficulty can be revealed through a minimal model of quantum field theory,the φ4theory, which is a scalar theory with quartic self-interactions.When the coupling approaches the phase transition point,the perturbative method becomes unreliable; and even in the weak-coupling regime,precision can be not controllable with the Feynman diagram calculation. Also,it is beyond the ability of lattice field theory which is good at calculating static properties, while scattering is a time-evolution problem. On the other hand,time-evolution can be implemented on a quantum computer and it is thus expected that scattering can be investigated efficiently.

    A technical challenge is that the number of degree of freedom per unit volume is infinite, and it is necessary to give a mapping that can assign qubits to physical degrees of freedom in a controllable way under desired accuracy. For instance,truncation of Fock space can be applied, since the number of particles is constrained by the energy of incoming particles.Then,given the number of qubits,preparation of incoming two particles as two wave packets, and their scattering, should be formulated in the language of quantum circuit. The quantum advantage is that the depth of the quantum circuit scales only polynomial with the number of particles,their energy,and the desired precision,for both weak and strong couplings.[64]

    The space is discretized as a d-dimensional L×···×L lattice Ω with lattice spacing a. Let us start with the lattice Hamiltonian

    Here φ(x) is the field operator on site x, π(x) is the conjugate field satisfying [φ(x),π(y)]= ia?dδx,y, ?ais a finitedifference operator, m0is the particle mass for the noninteracting theory H0corresponding to λ0=0.

    Two incoming particles should be expressed as a manybody initial state with two separated wave packets on the lattice. The initial state can be prepared in an adiabatic way,which involves two steps: firstly, two separated wave packets are produced on the vacuum of the non-interacting Hamiltonian H0, where the vacuum is a Gaussian state and can be constructed exactly;secondly,the interaction is turned on and the system evolves adiabatically with Hamiltonian parameterized by H(s(t))=s(t)H0+(1 ?s(t))H , where H(0)=H0,H(1)=H, and s(t) is the adiabatic path. Note that the wave packets will propagate and broaden since it is not an eigenstate of H0,additional backward evolution is required to undo unnecessary dynamical phases. The time complexity for the initial state procedure scales as O((1/ε)1+d/2). The system of the two incoming particles then evolves for a period that scattering occurs. Then, the interaction is adiabatic turn-off and the scattering result is sampled by measuring the number operator of momentum modes defined in the free theory. With scaling analysis via effective field theory, the algorithm shows a complexity polynomial with ε,e.g.,in 1D its O((1/ε)?1.5?o(1)). Remarkably, at strong coupling, the time complexity still scales polynomial with 1/(λc?λ0)(λcis the quantum phase transition point), the momentum of incoming particles p, the maximum kinematically allowed number of outgoing particles nout,and the precision ε. This is in contrast to the known classical algorithm scaling exponentially with 1/(λc?λ0)and 1/ε,thus showing an exponentially quantum speed-up.

    The merit of this work is to make clear the point of solving quantum field theory on a quantum computer with quantum advantage. A scheme to reach this goal has also been pointed out. Still, some subroutines may be improved to make the algorithm more implementable in the NISQ era,and there indeed are some following jobs with more detailed algorithms.[69]

    3.4. Non-equilibrium dynamics

    In physics we often need to study non-equilibrium quantum systems at finite temperature. For nuclear physics,quark–gluon plasma, which is produced in heavy-ion collision or in the expansion of the early Universe, is an outstanding example but still limited knowledge is known, as simulations with classical methods are very hard. Basically, studies of such quantum systems involve time evolution of the density matrix.As a quantum computer essentially manipulates pure quantum states with unitary operation, an extension to the density matrix needs additional treatment. One can either refer to a subsystem by tracing out the ancillary,or view the density matrix as an ensemble of pure quantum states with a classical distribution.

    In total, the hybrid quantum-classical approach adopts a strategy that assigns quantum resource for the real-time evolution part which is hard for the classical computer,while the classical part can utilize the state-of-art classical algorithm to solve steady states, since state preparation can be demanding for quantum resource on near-term quantum computers. A combination of both makes this approach promising on NISQ quantum processors for studying non-equilibrium dynamics.Remarkably,the standard lattice field theory techniques could be incorporated which can improve the DMQMC.[140]

    4. Outlook and summary

    At this stage, we have given an introductory review of quantum computing for nuclear physics, but a throughout investigation is still beyond the scope of this review, given recent rapid expanding of this field. We further briefly give an outlook that focuses on this field within the NISQ era. Although formulations of lattice gauge theory for quantum computing have been developed in a few different approaches,it is clear that there still is a considerable gap for simulating QCD on NISQ quantum processors with limited quantum resources,and it is anticipated that more efficient schemes for expressing degrees of freedom of gauge field and its coupling to matter field would be developed. On the other hand, incorporating quantum algorithms as subroutines into a hybrid quantum–classical algorithm can be a very promising approach for solving nuclear physics with practical results.

    Another observation is that there is still a lack of consensus on the goal; perhaps it is inspiring that a roadmap can be proposed for demonstrating quantum advantages for nuclear physics on a quantum computer, in accordance with the scaling-up of near-term quantum processors and error mitigation techniques. Moreover, software methodology is important, as it enables systematic developments and rapid innovations on both quantum algorithms and applications to nuclear physics. Compared with quantum computing for chemistry[141]and quantum machine learning,[142]opensource packages are still lacked for nuclear physics, and it awaits for such packages for making an attempt to study nuclear physics on a quantum computer more friendly.

    In summary, we have given a brief review on recent advances on quantum computing for nuclear physics. We have clarified two different approaches from analog and digital quantum simulators,and the review has focused on the latter.We have outlined how degrees of freedom of nuclear physics problems can be mapped into a formula of qubits that is solvable on a quantum computer,and the corresponding quantum algorithms for both static properties and real-time evolution have been shortly discussed. Concretely,we have given some examples from recent outstanding works, ranging from simulation of lattice gauge theory, nuclear structure as well as quantum advantage in terms of time complexity for scattering problem. Lastly, we have pointed out that lots of efforts are still required to make quantum computing a playground for investigating nuclear physics on near-term quantum devices.

    Acknowledgment

    We thanks very helpful feedback from Henry S.Lamm.

    国产1区2区3区精品| 婷婷六月久久综合丁香| 欧美在线一区亚洲| 色av中文字幕| 亚洲国产精品合色在线| 国产精品综合久久久久久久免费| 国产一区二区激情短视频| 视频在线观看一区二区三区| 国产激情久久老熟女| 久久精品亚洲精品国产色婷小说| 国产亚洲精品第一综合不卡| 免费女性裸体啪啪无遮挡网站| 热99re8久久精品国产| 精品国产乱子伦一区二区三区| 亚洲专区国产一区二区| 国产精品,欧美在线| 夜夜爽天天搞| 一本精品99久久精品77| 日本黄色视频三级网站网址| 国产主播在线观看一区二区| 69av精品久久久久久| 国产精品亚洲av一区麻豆| 久久 成人 亚洲| 妹子高潮喷水视频| 亚洲久久久国产精品| 日韩高清综合在线| 精品久久久久久久久久免费视频| 变态另类成人亚洲欧美熟女| 一进一出好大好爽视频| 高清毛片免费观看视频网站| 久久国产精品人妻蜜桃| 妹子高潮喷水视频| 欧美日韩中文字幕国产精品一区二区三区| 热99re8久久精品国产| 久久人妻av系列| 麻豆成人av在线观看| 国产精品久久久av美女十八| 国产成人精品无人区| 国产v大片淫在线免费观看| 国产亚洲精品一区二区www| 神马国产精品三级电影在线观看 | 日韩视频一区二区在线观看| 亚洲人成电影免费在线| 午夜久久久在线观看| 国产成人精品久久二区二区免费| 国产黄片美女视频| 亚洲 国产 在线| 国产主播在线观看一区二区| 国内揄拍国产精品人妻在线 | 久久久久国产精品人妻aⅴ院| 老汉色∧v一级毛片| 日韩 欧美 亚洲 中文字幕| 欧美黑人精品巨大| 少妇裸体淫交视频免费看高清 | 欧美成人性av电影在线观看| 成年免费大片在线观看| 老司机福利观看| 亚洲av成人一区二区三| 亚洲一区二区三区色噜噜| 男女下面进入的视频免费午夜 | 在线播放国产精品三级| 91老司机精品| 一边摸一边做爽爽视频免费| 中文字幕人成人乱码亚洲影| 欧美日韩亚洲综合一区二区三区_| 99在线视频只有这里精品首页| e午夜精品久久久久久久| 国产国语露脸激情在线看| 欧洲精品卡2卡3卡4卡5卡区| 国产精品电影一区二区三区| 国产成人啪精品午夜网站| 午夜久久久久精精品| 欧美激情久久久久久爽电影| 午夜福利18| 女生性感内裤真人,穿戴方法视频| 日本 欧美在线| bbb黄色大片| 欧美大码av| 亚洲天堂国产精品一区在线| 香蕉久久夜色| 国产高清视频在线播放一区| 日本 av在线| 性色av乱码一区二区三区2| 制服诱惑二区| 国产激情久久老熟女| 韩国av一区二区三区四区| 亚洲电影在线观看av| 日韩精品青青久久久久久| 国产91精品成人一区二区三区| 亚洲九九香蕉| 国产精品自产拍在线观看55亚洲| 久久久国产成人免费| 一卡2卡三卡四卡精品乱码亚洲| 成人欧美大片| 欧美国产精品va在线观看不卡| 成人特级黄色片久久久久久久| 亚洲国产毛片av蜜桃av| 久久婷婷人人爽人人干人人爱| 亚洲精品美女久久久久99蜜臀| 听说在线观看完整版免费高清| а√天堂www在线а√下载| 日本a在线网址| 俄罗斯特黄特色一大片| 露出奶头的视频| 久久久久久久午夜电影| 国产精品,欧美在线| 久热这里只有精品99| 搡老岳熟女国产| 欧美日本亚洲视频在线播放| 他把我摸到了高潮在线观看| 每晚都被弄得嗷嗷叫到高潮| 午夜免费成人在线视频| 亚洲在线自拍视频| 女人爽到高潮嗷嗷叫在线视频| 老司机靠b影院| 国产又色又爽无遮挡免费看| 午夜福利一区二区在线看| 啦啦啦 在线观看视频| 亚洲欧美一区二区三区黑人| a在线观看视频网站| 一级作爱视频免费观看| 性欧美人与动物交配| 一级毛片女人18水好多| 女性被躁到高潮视频| 久久久久久人人人人人| 亚洲国产精品999在线| 亚洲av成人av| 在线观看午夜福利视频| 国产日本99.免费观看| 色尼玛亚洲综合影院| 午夜日韩欧美国产| 精品少妇一区二区三区视频日本电影| 亚洲美女黄片视频| 国产欧美日韩一区二区三| 99在线人妻在线中文字幕| 亚洲熟妇熟女久久| 国产亚洲精品久久久久5区| 欧美成人一区二区免费高清观看 | 搡老熟女国产l中国老女人| 手机成人av网站| 每晚都被弄得嗷嗷叫到高潮| 人人妻人人澡欧美一区二区| 免费av毛片视频| 国产精品自产拍在线观看55亚洲| 一边摸一边抽搐一进一小说| 欧美人与性动交α欧美精品济南到| 无人区码免费观看不卡| 国产精品久久电影中文字幕| 中文字幕久久专区| 国产精品亚洲美女久久久| 国产一区二区三区视频了| 一级a爱片免费观看的视频| 久久亚洲精品不卡| 搡老妇女老女人老熟妇| 国产极品粉嫩免费观看在线| 欧美中文综合在线视频| 啦啦啦韩国在线观看视频| 欧美乱色亚洲激情| 一区二区三区激情视频| 亚洲美女黄片视频| 黄片小视频在线播放| 久久久久久大精品| 一个人观看的视频www高清免费观看 | 1024香蕉在线观看| 亚洲七黄色美女视频| 亚洲久久久国产精品| 性色av乱码一区二区三区2| 亚洲国产欧洲综合997久久, | 女警被强在线播放| 午夜免费激情av| 男人的好看免费观看在线视频 | 18美女黄网站色大片免费观看| 黄色片一级片一级黄色片| 亚洲五月婷婷丁香| 久久九九热精品免费| e午夜精品久久久久久久| 在线播放国产精品三级| 黑人操中国人逼视频| 免费无遮挡裸体视频| 成人av一区二区三区在线看| 亚洲国产看品久久| www.999成人在线观看| 国产亚洲精品久久久久5区| 亚洲五月婷婷丁香| 99热6这里只有精品| 高清在线国产一区| 窝窝影院91人妻| 老熟妇乱子伦视频在线观看| 熟女电影av网| 少妇粗大呻吟视频| 嫩草影视91久久| 国产男靠女视频免费网站| 婷婷精品国产亚洲av在线| 美女大奶头视频| 国产免费av片在线观看野外av| 午夜免费观看网址| 日韩三级视频一区二区三区| av有码第一页| 日韩av在线大香蕉| 久久午夜亚洲精品久久| 亚洲人成电影免费在线| 欧美成人免费av一区二区三区| 国产激情久久老熟女| 黑人操中国人逼视频| 久热爱精品视频在线9| 欧美三级亚洲精品| 国产精品久久久久久亚洲av鲁大| 一a级毛片在线观看| 90打野战视频偷拍视频| 又大又爽又粗| 国产97色在线日韩免费| 啦啦啦观看免费观看视频高清| 日韩国内少妇激情av| 国产亚洲精品久久久久久毛片| 亚洲国产毛片av蜜桃av| 露出奶头的视频| av超薄肉色丝袜交足视频| 精品欧美国产一区二区三| 精品一区二区三区四区五区乱码| 中文在线观看免费www的网站 | 精品欧美一区二区三区在线| 可以在线观看毛片的网站| 久久精品亚洲精品国产色婷小说| 精品一区二区三区视频在线观看免费| 亚洲成人精品中文字幕电影| АⅤ资源中文在线天堂| 久久久久九九精品影院| 欧美成人免费av一区二区三区| 99精品在免费线老司机午夜| 国产亚洲精品久久久久久毛片| 久久久国产欧美日韩av| 天堂√8在线中文| 人成视频在线观看免费观看| 中出人妻视频一区二区| 国产真实乱freesex| 精品日产1卡2卡| www国产在线视频色| 色综合欧美亚洲国产小说| 亚洲第一电影网av| √禁漫天堂资源中文www| 国产精品亚洲美女久久久| netflix在线观看网站| 侵犯人妻中文字幕一二三四区| 黄网站色视频无遮挡免费观看| 无遮挡黄片免费观看| 一二三四社区在线视频社区8| 性色av乱码一区二区三区2| 麻豆国产av国片精品| 午夜久久久在线观看| 国产伦一二天堂av在线观看| 日本免费一区二区三区高清不卡| 国产精品免费一区二区三区在线| 久久香蕉国产精品| 禁无遮挡网站| 亚洲第一欧美日韩一区二区三区| www.熟女人妻精品国产| 久久热在线av| 中文字幕久久专区| 亚洲va日本ⅴa欧美va伊人久久| 国产视频一区二区在线看| 男人舔奶头视频| 久久婷婷人人爽人人干人人爱| 成人亚洲精品一区在线观看| 亚洲精华国产精华精| 免费在线观看亚洲国产| 久久午夜综合久久蜜桃| 无限看片的www在线观看| 一进一出抽搐gif免费好疼| 日韩欧美一区视频在线观看| 黄色视频不卡| 在线观看免费日韩欧美大片| 中出人妻视频一区二区| av电影中文网址| 久久久久久人人人人人| 国产精品久久久久久亚洲av鲁大| 日本精品一区二区三区蜜桃| 国产精品 国内视频| 亚洲第一欧美日韩一区二区三区| 亚洲av片天天在线观看| ponron亚洲| 十八禁网站免费在线| 免费看十八禁软件| 九色国产91popny在线| 日本熟妇午夜| 搡老熟女国产l中国老女人| 久久久久久人人人人人| 男人舔女人下体高潮全视频| 成人亚洲精品一区在线观看| 99在线视频只有这里精品首页| 亚洲国产精品久久男人天堂| 国产亚洲精品av在线| 成人国产一区最新在线观看| 少妇被粗大的猛进出69影院| 黄频高清免费视频| 免费在线观看亚洲国产| 99精品在免费线老司机午夜| 999精品在线视频| АⅤ资源中文在线天堂| 成年人黄色毛片网站| 高清在线国产一区| 亚洲专区字幕在线| 国产av在哪里看| 淫秽高清视频在线观看| 两人在一起打扑克的视频| 久久精品国产99精品国产亚洲性色| 91麻豆av在线| 色播在线永久视频| 国内久久婷婷六月综合欲色啪| 精品高清国产在线一区| 9191精品国产免费久久| 在线国产一区二区在线| 91九色精品人成在线观看| 国产成+人综合+亚洲专区| 老司机靠b影院| 在线免费观看的www视频| 最近最新免费中文字幕在线| 日韩欧美在线二视频| 中文资源天堂在线| 亚洲精品一区av在线观看| 一个人免费在线观看的高清视频| 美女 人体艺术 gogo| 欧美不卡视频在线免费观看 | 12—13女人毛片做爰片一| 久久久久久久午夜电影| 别揉我奶头~嗯~啊~动态视频| 免费观看精品视频网站| 中文亚洲av片在线观看爽| 日韩三级视频一区二区三区| 最新在线观看一区二区三区| 久久 成人 亚洲| 色哟哟哟哟哟哟| 久久天堂一区二区三区四区| 免费搜索国产男女视频| 精品久久蜜臀av无| 黑人操中国人逼视频| 大型黄色视频在线免费观看| 国产欧美日韩精品亚洲av| 久久久久九九精品影院| 国产片内射在线| 国产精品98久久久久久宅男小说| 18禁黄网站禁片免费观看直播| 听说在线观看完整版免费高清| 99久久国产精品久久久| 成人欧美大片| 亚洲全国av大片| 俄罗斯特黄特色一大片| 欧美色欧美亚洲另类二区| 亚洲中文日韩欧美视频| 不卡一级毛片| 精品久久久久久久毛片微露脸| 精品无人区乱码1区二区| 女人爽到高潮嗷嗷叫在线视频| 国产精品久久久久久亚洲av鲁大| 亚洲精品色激情综合| 黄片大片在线免费观看| 九色国产91popny在线| 在线观看www视频免费| 国产精品免费一区二区三区在线| 香蕉国产在线看| 精品高清国产在线一区| 国产av在哪里看| 搡老妇女老女人老熟妇| 精品久久久久久久久久免费视频| 精品电影一区二区在线| 国产精品爽爽va在线观看网站 | 热99re8久久精品国产| 看免费av毛片| 男人舔女人的私密视频| 亚洲熟妇中文字幕五十中出| 一级作爱视频免费观看| 在线国产一区二区在线| 亚洲美女黄片视频| 精品国产乱码久久久久久男人| 99久久综合精品五月天人人| 国产乱人伦免费视频| 久久久久久久久久黄片| 麻豆国产av国片精品| 色播在线永久视频| 美女免费视频网站| 亚洲五月天丁香| 夜夜躁狠狠躁天天躁| 日本免费一区二区三区高清不卡| 俄罗斯特黄特色一大片| 亚洲人成伊人成综合网2020| 最近最新中文字幕大全电影3 | 成人特级黄色片久久久久久久| 白带黄色成豆腐渣| 免费在线观看影片大全网站| 亚洲精品中文字幕在线视频| 精品久久久久久,| 欧美日韩中文字幕国产精品一区二区三区| 国产精品久久视频播放| 成人三级做爰电影| av电影中文网址| 韩国精品一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| aaaaa片日本免费| 在线十欧美十亚洲十日本专区| 亚洲国产看品久久| 久久性视频一级片| 波多野结衣巨乳人妻| 成人三级黄色视频| 麻豆av在线久日| 啦啦啦观看免费观看视频高清| 免费看a级黄色片| av超薄肉色丝袜交足视频| 亚洲 欧美一区二区三区| 久久欧美精品欧美久久欧美| 琪琪午夜伦伦电影理论片6080| 国产日本99.免费观看| 法律面前人人平等表现在哪些方面| 国产黄片美女视频| 一二三四社区在线视频社区8| bbb黄色大片| 女人爽到高潮嗷嗷叫在线视频| 亚洲欧美日韩无卡精品| 亚洲黑人精品在线| 亚洲欧洲精品一区二区精品久久久| 黄网站色视频无遮挡免费观看| 国产成年人精品一区二区| 免费在线观看黄色视频的| 波多野结衣巨乳人妻| 美女大奶头视频| 久久精品成人免费网站| 淫妇啪啪啪对白视频| 亚洲av成人不卡在线观看播放网| 亚洲人成77777在线视频| 91麻豆av在线| 亚洲国产高清在线一区二区三 | 精品久久久久久久毛片微露脸| 国产单亲对白刺激| 国产亚洲av高清不卡| 午夜福利一区二区在线看| 久久中文看片网| 91av网站免费观看| 国内精品久久久久久久电影| 欧美绝顶高潮抽搐喷水| 他把我摸到了高潮在线观看| 亚洲成人久久爱视频| 波多野结衣巨乳人妻| 婷婷精品国产亚洲av在线| 亚洲一码二码三码区别大吗| 黑人欧美特级aaaaaa片| 国产精品影院久久| 欧美乱妇无乱码| 老司机深夜福利视频在线观看| 久久久久久国产a免费观看| 日本一区二区免费在线视频| 99热这里只有精品一区 | 大型黄色视频在线免费观看| 狠狠狠狠99中文字幕| 最新在线观看一区二区三区| 天天一区二区日本电影三级| 久久亚洲精品不卡| 这个男人来自地球电影免费观看| 我的亚洲天堂| 日韩免费av在线播放| 最近最新免费中文字幕在线| 正在播放国产对白刺激| 一边摸一边做爽爽视频免费| 亚洲精品美女久久久久99蜜臀| 国产成人影院久久av| 在线观看舔阴道视频| 日本 av在线| 日本一区二区免费在线视频| 国产精品乱码一区二三区的特点| 亚洲精品色激情综合| 91麻豆av在线| 老汉色∧v一级毛片| 91成年电影在线观看| 最近最新中文字幕大全电影3 | 宅男免费午夜| 老司机午夜十八禁免费视频| 午夜激情av网站| 美女高潮喷水抽搐中文字幕| 国产成人一区二区三区免费视频网站| 一级a爱视频在线免费观看| 国产精品日韩av在线免费观看| svipshipincom国产片| 久久欧美精品欧美久久欧美| 老鸭窝网址在线观看| 午夜免费鲁丝| 91成人精品电影| 欧美成人一区二区免费高清观看 | 日本 av在线| 黄色片一级片一级黄色片| 欧美激情久久久久久爽电影| 亚洲精品国产精品久久久不卡| 美女扒开内裤让男人捅视频| 国产爱豆传媒在线观看 | 亚洲欧美一区二区三区黑人| 黑人欧美特级aaaaaa片| 美女大奶头视频| 久99久视频精品免费| 侵犯人妻中文字幕一二三四区| 哪里可以看免费的av片| 国产亚洲欧美精品永久| 不卡av一区二区三区| 波多野结衣av一区二区av| 给我免费播放毛片高清在线观看| 免费观看精品视频网站| 日本 av在线| 久久这里只有精品19| 亚洲中文字幕一区二区三区有码在线看 | 久久精品成人免费网站| 看黄色毛片网站| 亚洲精品美女久久av网站| 老司机午夜十八禁免费视频| 999久久久精品免费观看国产| 色精品久久人妻99蜜桃| av电影中文网址| 亚洲男人天堂网一区| 国产v大片淫在线免费观看| 久久久久久人人人人人| 日本 欧美在线| 亚洲最大成人中文| 极品教师在线免费播放| 免费在线观看日本一区| 午夜a级毛片| 国产区一区二久久| 国产爱豆传媒在线观看 | 欧美日韩乱码在线| 日日摸夜夜添夜夜添小说| 黄色女人牲交| 国产97色在线日韩免费| 91字幕亚洲| 在线观看日韩欧美| 妹子高潮喷水视频| 亚洲精品在线观看二区| 香蕉国产在线看| 国产精品亚洲一级av第二区| 狠狠狠狠99中文字幕| 欧美绝顶高潮抽搐喷水| 草草在线视频免费看| 黄色片一级片一级黄色片| 亚洲电影在线观看av| 国内毛片毛片毛片毛片毛片| 男女那种视频在线观看| 中文字幕人成人乱码亚洲影| 亚洲成人久久爱视频| 久久精品91蜜桃| 欧美日韩乱码在线| 欧美一级a爱片免费观看看 | 午夜视频精品福利| 制服诱惑二区| 亚洲一区中文字幕在线| 日韩欧美免费精品| 人妻丰满熟妇av一区二区三区| 国产在线观看jvid| 久久久久久亚洲精品国产蜜桃av| 日韩精品青青久久久久久| 亚洲一卡2卡3卡4卡5卡精品中文| 白带黄色成豆腐渣| 日韩av在线大香蕉| 亚洲最大成人中文| 狠狠狠狠99中文字幕| 香蕉久久夜色| 久久中文字幕一级| 黑人操中国人逼视频| 一个人免费在线观看的高清视频| 我的亚洲天堂| 亚洲成人国产一区在线观看| 亚洲欧美精品综合久久99| 午夜福利在线在线| 一级作爱视频免费观看| 国产一级毛片七仙女欲春2 | 搡老熟女国产l中国老女人| 琪琪午夜伦伦电影理论片6080| 精品一区二区三区av网在线观看| 亚洲第一青青草原| 亚洲天堂国产精品一区在线| 国产精品久久久人人做人人爽| 欧美中文综合在线视频| 男女下面进入的视频免费午夜 | 一级毛片精品| 亚洲精品粉嫩美女一区| 国产激情欧美一区二区| 黄色视频不卡| 一级毛片女人18水好多| 国产欧美日韩精品亚洲av| 国内精品久久久久精免费| 欧美黑人欧美精品刺激| 我的亚洲天堂| 在线天堂中文资源库| 免费一级毛片在线播放高清视频| 天堂√8在线中文| 欧美最黄视频在线播放免费| 啦啦啦 在线观看视频| 日日夜夜操网爽| 19禁男女啪啪无遮挡网站| 男男h啪啪无遮挡| 免费看美女性在线毛片视频| 成人国产综合亚洲| 亚洲精品久久国产高清桃花| 亚洲精品在线美女| 亚洲色图 男人天堂 中文字幕| 亚洲一区二区三区色噜噜| 麻豆成人午夜福利视频| 自线自在国产av| 97超级碰碰碰精品色视频在线观看| 午夜影院日韩av| 久久精品国产99精品国产亚洲性色| 免费看a级黄色片| 国产成人欧美在线观看| 国产又爽黄色视频| 亚洲国产中文字幕在线视频| 国产免费av片在线观看野外av| 欧美大码av| 波多野结衣高清作品| 黄色毛片三级朝国网站| 久久久久久国产a免费观看| 亚洲午夜理论影院| 精品免费久久久久久久清纯| 欧美日韩精品网址| 午夜免费激情av|