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

    Some mathematical aspects of Anderson localization:boundary effect,multimodality,and bifurcation

    2022-11-11 07:53:06ChenJiaZiqiLiuandZhiminZhang
    Communications in Theoretical Physics 2022年11期

    Chen Jia,Ziqi Liu and Zhimin Zhang,3

    1 Applied and Computational Mathematics Division,Beijing Computational Science Research Center,Beijing 100193,China

    2 Department of Mathematical Sciences,Tsinghua University,Beijing 100084,China

    3 Department of Mathematics,Wayne State University,Detroit,Michigan 48202,United States of America

    Abstract Anderson localization is a famous wave phenomenon that describes the absence of diffusion of waves in a disordered medium.Here we generalize the landscape theory of Anderson localization to general elliptic operators and complex boundary conditions using a probabilistic approach,and further investigate some mathematical aspects of Anderson localization that are rarely discussed before.First,we observe that under the Neumann boundary condition,the low energy quantum states are localized on the boundary of the domain with high probability.We provide a detailed explanation of this phenomenon using the concept of extended subregions and obtain an analytical expression of this probability in the one-dimensional case.Second,we find that the quantum states may be localized in multiple different subregions with high probability in the one-dimensional case and we derive an explicit expression of this probability for various boundary conditions.Finally,we examine a bifurcation phenomenon of the localization subregion as the strength of disorder varies.The critical threshold of bifurcation is analytically computed based on a toy model and the dependence of the critical threshold on model parameters is analyzed.

    Keywords:landscape,spectrum,eigenvalue,eigenmode,eigenfunction,elliptic operator,Schr?dinger operator,confinement

    1.Introduction

    Anderson localization is a general wave phenomenon that applies to the transport of electromagnetic waves,quantum waves,spin waves,acoustic waves,etc[1–9].In particular,it is responsible for the metal-insulator transition in disordered alloys.Due to its importance and universality,it has fascinated scientists and mathematicians in many different areas and produced a huge body of literature in the past 60 years[10–12].Since Anderson’s seminal work[1],it is known that a Schr?dinger operator with a lattice potential can produce a strongly localized quantum state,provided that the degree of randomness(disorder)in the lattice is sufficiently large.This phenomenon has now been experimentally demonstrated in optic and electromagnetic systems[13–15].

    One of the most puzzling aspects of Anderson localization is the strong spatial confinement of quantum states of the system,i.e.the ability to maintain standing waves in a very restricted subregion of the domain,with their amplitudes decaying exponentially at long range even in the absence of confining force or potential.Recently,Filoche and Mayboroda[16]proposed a universal mechanism for Anderson localization.This mechanism reveals that in any vibrating system,there exists a hidden landscape of localization that partitions the original domain into many weakly coupled subregions.The boundaries of these subregions correspond to the valley lines of the landscape.The height of the landscape along its valley lines determines the strength of coupling between the subregions.This theory allows one to predict the localization behavior,including the geometry of the confining subregions,the energy of the quantum states,and the critical energy above which one can expect fully delocalized(conducting)states to appear.Moreover,by constructing an Agmon distance,it has been shown[17–19]that the inverse of the localization landscape can be interpreted as an effective confining potential that is responsible for the exponential decay of the localized quantum states even far from its main confining subregion.

    Despite more than half a century of research,many issues on the exact mechanism of Anderson localization still remain open[10–12].In the present paper,we consider Anderson localization of quantum states of a Schr?dinger equation with a disordered lattice potential under Neumann and Robin boundary conditions.Using a probabilistic approach that represents the solution of a second-order elliptic equation in terms of a reflecting diffusion,we extend the concepts of localization landscape and its valley lines to general nonsymmetric second-order elliptic operators and arbitrary boundary conditions,and we prove that the amplitude of any vibrational eigenmode can be controlled effectively by the landscape.The probabilistic approach is also applied to study the limit localization behavior of quantum states when the strength of disorder tends to infinity.

    More importantly,we investigate some novel mathematical aspects of Anderson localization that are seldom discussed in the literature.First,we observe that under the Neumann and Robin boundary conditions,the low energy quantum states are localized on the boundary of the domain with relative high probability.Moreover,we find that compared to the one-dimensional case,the eigenmodes of a twodimensional system are more likely to be localized on the boundary.These boundary effects are analyzed systematically,which is the first main contribution of the present paper.In addition,we find that the quantum states of a onedimensional system are often localized in multiple different subregions simultaneously and thus exhibit multiple peaks.We study the multimodal phenomenon in detail,which is the second main contribution of this paper.

    Finally,we find that under different strengths of disorder,the quantum states may be localized in completely different subregions.When the disorder is large,the first eigenmode tends to be localized in the largest low potential subregion with strong symmetry.However,when the disorder is small,the first eigenmode tends to be localized in the union of some low potential subregions with some potential barriers between them.As a result,the localization behavior of the system may undergo a bifurcation as the strength of disorder varies.We study such bifurcation phenomenon in detail based on a toy model,which is the third main contribution of this article.

    The structure of the present paper is organized as follows.In section 2,we extend the concepts of localization landscape and valley lines to general non-symmetric elliptic operators and general boundary conditions,prove the key Filoche–Mayborodaca inequality using a probabilistic approach,and study the limit localization behavior when the strength of disorder is large.In section 3,we examine the boundary effect of Anderson localization under Neumann and Robin boundary conditions,and clarify the reason why the low energy quantum states in high dimensions are more likely to be localized on the boundary of the domain using the concept of extended subregions.In section 4,we reveal why the eigenmodes may be localized in multiple different subregions simultaneously and compute the probability of this multimodal phenomenon in the one-dimensional case for various boundary conditions.In section 5,we investigate the bifurcation of the confining subregion as the strength of disorder varies based on a toy model,reveal the essence of the bifurcation phenomenon,and obtain the equation satisfied by the critical threshold of bifurcation.We conclude in section 6.

    2.Localization landscape and limit behavior

    2.1.Model

    Anderson localization is a famous wave phenomenon.Since Anderson’s seminal work[1],it was known that the low energy quantum states of a Schr?dinger operator with a disordered lattice potential are strongly localized within a small subregion,provided that the degree of randomness in the lattice is sufficiently large.

    Here we consider high-dimensional Anderson localization for the quantum states of the following stationary Schr?dinger equation with a Dirichlet,Neumann,or Robin boundary condition(Anderson’s original work[1]considers a tight-binding model which can be viewed as the discretization of the continuous model studied here):

    where H=-△+KV is the Hamiltonian with V=V(x)being a disordered potential,λ >0 is an energy level(eigenvalue),u=u(x)is the associated quantum state(eigenmode),and K >0 is the strength of disorder.The domain Ω=(0,1)dis the d-dimensional unit hypercube with each side being divided uniformly into N intervals.In this way,the hypercube Ω is divided into Ndsmaller hypercubes of the same size.In each small hypercube Ωk,the potential V is a constant with its value being sampled from a given probability distribution,which is often chosen as the Bernoulli distribution

    or the uniform distribution

    The values of V in these small hypercubes are independent of each other.The boundary condition of the model is rather general with g ≥0 being a given constant,h=h(x)≥0 being a given function on ?Ω,and n=n(x)being the unit outward pointing normal vector field on ?Ω.If g=0 and h=1,then it reduces to the Dirichlet boundary condition; if g=1 and h=0,then it reduces to the Neumann boundary condition; if g=1 and h ≠0,then it is the Robin boundary condition.

    which means that there is no net probability current across the boundary.Hence imposing the Neumann boundary condition can confine the system within the domain without an infinite potential well,as opposed to the Dirichlet boundary condition.This is the case of perfect reflection on the boundary(see section 5.2 of[20]for a detailed explanation).The Neumann boundary condition turns out to be very useful in the R-matrix theory of scattering[21,22]and the theory of quantum graphs[23].The Robin boundary condition is also considered here since it builds a bridge between the Dirichlet and Neumann boundary conditions.

    In fact,the above model can be generalized to a more complicated model as follows:

    2.2.Localization landscape

    A recent theory[16]has shown that under the Dirichlet boundary condition(g=0 and h=1),the spatial location of the localized quantum states of the eigenvalue problem(1)can be predicted by the solution of an associated Dirichlet problem

    where the solution w=w(x)is called the localization landscape.In fact,the theory in[16]was developed for symmetric elliptic operators and the Dirichlet boundary condition using the technique of the Green function.The landscape theory was further developed in[26]using a probabilistic approach and generalized in[19]to the Neumann boundary condition.Here we extend the concept of localization landscape to general non-symmetric elliptic operators and more complicated boundary conditions using a different probabilistic approach.

    Note that the boundary value problem(8)is obtained by setting the right-hand side of the eigenvalue problem(3)to be 1.Similarly,the landscape w also has a probabilistic representation using reflecting diffusion,which is given by(see appendix A for the proof)

    The probabilistic representations(7)and(9)are closely related.If we normalize the eigenmode u such that‖u‖∞=1,then we obtain the Filoche–Mayboroda inequality

    This inequality shows that when the energy level λ is not large,the eigenmode u must be small at those points where the landscape w is small.Hence the low energy eigenmodes must be localized to those subregions where the landscape is large.To see this,we illustrate the graphs of the landscape w and the first four eigenmodes u/λ for a one-dimensional problem when K is large(figure 1(a)).It can be seen that the inequality(10)indeed provides an accurate upper bound for the eigenmodes.Furthermore,the spatial locations of the confining subregions of these eigenmodes are perfectly predicted by the peak positions of the landscape.Recent numerical and theoretical studies[19,28]have shown that the order of the eigenmodes is closely related to the height of the peaks of the landscape.This phenomenon was observed in our simulations in figure 1(a),where the first four eigenmodes are localized around the highest four peaks of the landscape.

    2.3.Valley lines

    In the one-dimensional case,the local minima of the landscape w divide the domain Ω into many subintervals(figure 1(a)).Such subintervals give possible spatial locations where the eigenmodes can be localized.In the two-dimensional case,the confining subregions are determined by the valley lines of the landscape[16].In fact,the valley lines are defined as the lines of steepest descent,starting from the saddle points of the landscape and going to its local minima.In this way,the valley lines separate the domain Ω into many subregions.Along the valley lines,the values of w are expected to be small.From the inequality(10),the values of u are also small along the valley lines and thus the eigenmodes must be localized within the subregions enclosed by these lines.Following[16,17],we use the watershed algorithm proposed in[29]to compute the valley lines numerically.In higher dimensions,valley lines will become hypersurfaces and the results are similar.

    Figures 1(b),(c)illustrate the landscape,eigenmodes,and valley lines for a two-dimensional problem under the Neumann boundary condition.It is clear that the valley lines of the landscape separate the domain into many subregions and each eigenmode is exactly located in one of these subregions.Under the Dirichlet boundary condition,the eigenmodes must be located in the interior of the domain since they vanish on the boundary[16].However,under the Neumann boundary condition,some eigenmodes may be localized on the boundary(see the first,second,and fourth ones in figure 1(c)).The boundary values of these eigenmodes are large and their local maxima may even appear on the boundary.This phenomenon will never appear for the Dirichlet boundary condition.

    Here the eigenvalue problem(3)and the boundary value problem(8)are solved numerically by using the spectral element method instead of the classical finite difference method(FDM)or finite element method(FEM).In the spectral element method,the solution of the associated partial differential equation(PDE)is approximated by piecewise high-order polynomials and the Legendre polynomials are used as a basis to approximate the solution in each small hypercube Ωk.The degree of the Legendre polynomials is chosen to be 10 in the one-dimensional case and 6 in the two-dimensional case.In fact,classical FDM or FEM only have second-order convergence with respect to the mesh size,while the spectral element method can achieve exponential convergence with respect to the degree of the polynomials.Therefore,the spectral element method can obtain higher numerical accuracy with less computational costs than the other two methods.

    2.4.Limit behavior for large K

    Next,we focus on the limit behavior of the eigenmodes when the strength K of disorder is very large.For simplicity,we assume that the potential V restricted to each small hypercube Ωkis sampled from the Bernoulli distribution given in(2),which can only take the value of 0 or 1.To proceed,let

    denote the set of points at which the potential vanishes.We first consider the behavior of the eigenmodes outside D.Let x ∈Ω be the initial position of the reflecting diffusion Xtwhich solves the Skorokhod equation(6).When x ?D,we have V(Xs)=1 when s is small,which implies that

    It thus follows from(7)and(9)and the dominated convergence theorem that

    This shows that the landscape and eigenmodes must vanish outside D in the limit of K →∞.In other words,the eigenmodes can only be localized in the region where the potential attains its minimum.This coincides with the intuition that the quantum states tend to be localized in the region with low potential energy.

    We next focus on the behavior of the eigenmodes inside D.Clearly,D can be decomposed as the disjoint union of several connected components(subregions),i.e.D=D1∪D2∪…∪DM.Since the potential V vanishes in D,in the limit of K →∞,the eigenmode u must be the solution to the following local eigenvalue problem in each subregion Dk:

    Therefore,when K is very large,the spectrum of the Hamiltonian H=-L+KV is composed of the local eigenvalues of the operator -L in each subregion Dk.If an eigenvalue of H coincides with one of the local eigenvalues of -L in Dk,then the corresponding eigenmode will be localized in Dk; conversely,if an eigenvalue of H coincides with neither one of the local eigenvalues of -L in Dk,then the corresponding eigenmode will not be localized in Dk.Hence the limits in(11)enable us to decompose the eigenvalues problem(3)with random potential in the whole domain Ω into some local eigenvalue problems(12)with zero potential in the subregions Dk.When K is not very large,the landscape and eigenmodes are small but not zero on ?Dk?Ω and thus the domain Ω can be divided into the many weakly coupled subregions[16].

    To gain a deeper insight,we depict the potential,landscape,and eigenmodes for a one-dimensional problem under different values of K(figure 2(a)).When K is large,the landscape and eigenmodes are only localized in the region where the potential vanishes.However,this is not the case when K is relatively small.Figure 2(b)illustrates the potential and valley lines for a two-dimensional problem under different values of K.When K is small,the confining subregions enclosed by the valley lines may include many connected components of D.However,for sufficiently large K,the connected components of D are exactly separated by valley lines.

    3.Boundary effect

    In the previous discussion,we develop our theory for the general model(3).In the following,we only focus on the simpler model(1).For simplicity,we assume that the potential V restricted to each small hypercube Ωkis Bernoulli distributed.

    3.1.Localization on the boundary

    We have seen that under the Neumann or Robin boundary condition,the low energy quantum states are very likely to be localized on the boundary(see the second eigenmode in figure 1(a)and the first,second,and fourth eigenmodes in figure 1(c)).Such a phenomenon will never take place for the Dirichlet boundary condition.When an eigenmode is localized on the boundary,its boundary value must be very large and the local maxima of the eigenmode may even appear on the boundary.In what follows,each eigenmode is always normalized such that ‖u(x)‖∞=1.To better characterize the boundary effect,we introduce the following definition.In the one-dimensional case,we say that an eigenmode u is localized on the boundary if it satisfies

    In the two-dimensional case,similarly,we say that an eigenmode u is localized on the boundary if it satisfies

    and we say that it is localized in the corner if it satisfies

    Since the potential V is stochastic,the eigenmodes may or may not be localized on the boundary.For simplicity,we only consider the first eigenmode.The probability for the first eigenmode to be localized on the boundary is denoted by Pband the probability for the first eigenmode to be localized in the corner is denoted by Pc.These two probabilities are collectively referred to as boundary probabilities.Clearly,the boundary probabilities are both 0 for the Dirichlet boundary condition.However,for the Neumann or Robin boundary condition,these probabilities are strictly positive with their values depending on the parameters h,K,and p,where the function h is chosen to be a constant for simplicity.Figures 3(a)–(c)illustrate the boundary probabilities as h,K,and p vary under the Robin boundary condition.Here the boundary probabilities Pband Pcare computed numerically using the Monte Carlo method.Specifically,we generate the random potential 1000 times,solve the corresponding eigenvalue problems,and then count the frequency that the first eigenmode appears on the boundary or in the corner.The spectral element method enables us to solve the eigenvalue problem thousands of times in an acceptable time.It can be seen that the boundary probabilities decrease with the three parameters in both the one-and two-dimensional cases.When h is small,the first eigenmode is localized on the boundary with a relatively high probability.In particular,the probability Pbis exceptionally large in the two-dimensional case.When h ?1,the Robin boundary condition reduces to the Dirichlet one and thus all boundary probabilities tend to zero(figure 3(a)).This explains why the boundary probabilities decrease as h increases.

    We next focus on the dependence of the boundary probabilities on K.When K=h=0,the problem(1)reduces to the Laplacian eigenvalue problem with the Neumann boundary condition and thus the first eigenmode is a constant.According to the definition,a constant eigenmode is localized on the boundary(since its maximal boundary value is larger than 0.5).As K increases,the eigenmode may be localized in the interior of the domain and thus the boundary probabilities decrease.When K ?1,compared to the one-dimensional case,we find that the eigenmodes for a two-dimensional problem are more likely to be localized on the boundary(figure 3(b)).The effect of p on the boundary probabilities is much more complicated and will be discussed later.

    3.2.An explanation of the boundary effect

    Based on the simulations in figures 3(a)–(c),we have made two crucial observations:(i)when h is small,the eigenmodes are localized on the boundary with relatively high probability;(ii)when K is large,compared to the one-dimensional case,the eigenmodes for a two-dimensional problem are more likely to be localized on the boundary.Here we provide a detailed explanation for these two observations when h=0(Neumann boundary condition)and K ?1(strong disorder).In fact,the boundary effect can be explained intuitively using the energy functional theory.To see this,recall that for the Dirichlet and Neumann boundary conditions,the eigenmodes of(1)are actually the critical points of the energy functional[20]

    Hence the energy functional perspective also guides us to focus on the local eigenvalue problem in each subregion Dk.Intuitively,for the Dirichlet boundary condition,the function u must tend to zero on the boundary,which leads to a large|?u|and thus a large kinetic energy Jk(u)when Dkappears on the boundary.For the Neumann boundary condition,however,the function u does not have to be zero on the boundary,which makes the kinetic energy Jk(u)much smaller.This explains why the boundary effect is significant for the Neumann boundary condition.

    We next examine the boundary effect in more detail.For simplicity,we first focus on the one-dimensional case.For large K,we have shown that the eigenmode of the problem(1)restricted to each subregion Dkis approximately the solution to the local problem(12),which is a Laplacian eigenvalue problem.The first eigenmode must be localized in the subregion Dkwhere the local eigenvalue is the smallest.Suppose that Dk=[xk,xk+L],where L is the length of Dk.WhenDk∩?Ω = ?,i.e.xk>0 and xk+L <1,the local problem(12)can be rewritten as

    With the aid of the concept of extended subregions,we are able to transform a local eigenvalue problem with a mixed boundary condition into an equivalent eigenvalue problem with the Dirichlet boundary condition in an extended subregion.When h is small and K is large,the first eigenmode will be localized in the extended subregion with the smallest Dirichlet eigenvalue.In particular,in the one-dimensional case,the first eigenmode will be localized in the longest extended subregion.If the longest two extended subregions are the same in length,then the first eigenmode may be localized in these two subregions simultaneously.This phenomenon will be discussed in detail later.

    Note that in the one-dimensional case,there are at most two subregions around the boundary and the extended subregion is twice as long as the original one.However,in the two-dimensional case,the proportions of subregions that need to be extended are much higher(figure 3(e)).Moreover,if a subregion appears in the corner,then its area should be‘magnified’ four times.These two facts explain why compared to the one-dimensional case,the eigenmodes for a twodimensional problem have a much higher probability to be localized on the boundary.

    3.3.Dependence of the boundary probabilities on p

    In the previous discussion,we have explained why the boundary probabilities decrease as h and K increase.Here we focus on the dependence of the boundary probabilities on p.For simplicity,we only focus on the one-dimensional case for small h and large K.

    In the one-dimensional case,the domain Ω=(0,1)is divided uniformly into N subintervals of the same length 1/N.In each subinterval,the potential V can only take the value of 0 or 1 with p being the probability of V=0 and q=1-p being the probability of V=1.Then the domain Ω can be decomposed into many subregions with V=0 and V=1 alternatively(see the white and gray regions in figure 2(a)).Intuitively,when N is large,the lengths of these subregions are approximately independent and the number of subintervals included in each subregion with V=0(V=1)approximately follows a geometric distribution with parameter p(q).For small h and large K,we have shown that the first eigenmode must be localized in the longest extended subregion with V=0.Therefore,the probability that the first eigenmode is localized on the boundary is approximately equal to the probability that the longest extended subregion with V=0 appears on the boundary,which is given by(see appendix C for the proof)

    This formula reveals a complex relationship between the boundary probability Pbon the parameter p in the onedimensional case.

    To test our theory,we compare the theoretical prediction given in(21)with numerical simulations where the boundary probability Pbis computed by generating the random potential 1000 times and then solving the associated eigenvalue problem(1)(figure 3(f)).Interestingly,while the analytical expression is very complicated,it is in perfect agreement with simulation results.Moreover,we find that Pbdecreases with p.When p ≈0.5,we have Pb≈0.26,which means that the first eigenmode has a one in four chance of being localized on the boundary.In the two-dimensional case,the dependence of the boundary probabilities Pbor Pcon the parameter p is similar(figure 3(c)),but it is very difficult to obtain their analytical expressions.

    3.4.Boundary effects for other types of random potentials

    Thus far,the boundary effects were examined when the lattice potential V within each hypercube is randomly sampled from a Bernoulli distribution.From figure 1,it is clear that the boundary effect is also significant when the random potential V is sampled from a uniform distribution.A natural question is whether the boundary effect is also significant for other types of random potentials.

    To answer this,we compare the boundary effects for four types of random potentials sampled from the Bernoulli distribution,the normal distribution,the gamma distribution,and the uniform distribution,respectively,as shown in figure 4.Here the strength of disorder is chosen to be K=104.For each distribution type,we illustrate the boundary probabilities Pband Pcin the one- and two-dimensional cases for random potentials with the same mean μ=1/2 and three different standard deviations,where the standard deviations σ within each hypercube are chosen to beand σ=μ/3,respectively(see the four columns of figure 4).It is clear that the boundary effect is significant for all types of random potentials.Interestingly,we find that the boundary effect is insensitive to the standard deviation of the random potential within each hypercube.For the Bernoulli and gamma distributions,the boundary probabilities Pband Pcbecome slightly higher as the standard deviation decreases,i.e.as the distribution becomes more concentrated.However,the normal and uniform distributions give rise to the opposite effect,i.e.the boundary probabilities become slightly lower with the decrease of the standard deviation.Furthermore,we find that whether the low energy quantum states are localized on the boundary is remarkably affected by the distribution type.The random potentials sampled from the normal distribution lead to a much weaker boundary effect than those sampled from the other three types of distributions.

    4.Multimodality

    When K ?1,we have shown that the original system(1)can be decomposed into many local systems in smaller subregions D1,D2,...,DM.The spectrum of the original system is composed of the eigenvalues of the local systems.If an eigenvalue of the original system coincides with one of the local eigenvalues in the subregion Dk,then the associated eigenmode will be localized in Dk.Along this line,if multiple subregionsshare a common local eigenvalue,then the corresponding eigenmode may be localized in these subregions simultaneously and thus have multiple peaks.This phenomenon will be referred to as multimodality in this paper.Clearly,multimodality takes place when the system has multiple eigenvalues which are approximately equal.Figure 5(a)illustrates multimodality for a one-dimensional problem.In the one-dimensional case,we have shown that the first eigenmode is localized in the longest extended subregion when K ?1.If the longest two extended subregions are the same in length,then the first eigenmode may be localized in these two subregions simultaneously.In higher dimensions,any two subregions rarely have the same shape and thus rarely share a common eigenvalue.Hence multimodality in general will not occur in the high-dimensional case.Next,we only focus on multimodality in the one-dimensional case.

    As mentioned earlier,the domain Ω=(0,1)is divided uniformly into N subintervals of the same length 1/N.The number of subintervals included in each subregion with V=0(V=1)approximately follows a geometric distribution with parameter p(q).When K ?1,the first eigenmode must be localized in the longest extended subregion with V=0 for both the Dirichlet and Neumann boundary conditions.Thus the probability that the first eigenmode displays multimodality is approximately equal to the probability that there is more than one longest extended subregion.For the Dirichlet boundary condition,the probability of multimodality is approximately given by(see appendix C for the proof)

    For the Neumann boundary condition,the probability of multimodality is much more complicated and is approximately given by

    where[x]represents the largest integer less than or equal to x.

    To test our theory,we compare the analytical expressions given in(22)and(23)with numerical simulations for both the Dirichlet and Neumann boundary conditions(figures 5(b),(c)).It can be seen that the theoretical prediction is in excellent agreement with simulation results.For both boundary conditions,the probability of multimodality increases with the parameter p for large K.When p=0.5,the probability of multimodality is 0.28 for the Dirichlet boundary condition and 0.25 for the Neumann boundary conditions,which means that the first eigenmode has a one in four chance of displaying multimodality.The relatively large values of PDand PNalso suggest that multimodality is a common phenomenon in the one-dimensional case.

    5.Bifurcation of localization subregions

    5.1.Bifurcation analysis in a toy model

    From the simulations in figure 2(a),we can see that for different values of K,the eigenmodes may be localized in completely different subregions.When K ?1,the first eigenmode is localized in the extended subregion with the smallest local eigenvalue.In particular,in the one-dimensional case,the first eigenmode is localized in the longest extended subregion with V=0.When K is relatively small,however,the first eigenmode may be localized in somewhere else.This suggests that as K increases,the localization subregion for a given eigenmode may change and a bifurcation phenomenon may take place.Note that in figure 2(a),while the first eigenmode is not localized in the longest extended subregion with V=0 when K is small,it is localized in the union of some long subregions with V=0 that are separated by some short barriers with V=1.The union of these subregions is even longer than the longest extended subregion.In other words,when K is large,the localization behavior of the system is very sensitive to the short potential barriers,while it is insensitive to these barriers when K is small.We emphasize that while figure 2 only shows a bifurcation for the first eigenmode,similar phenomena can be also observed for other eigenmodes.

    To gain a deeper insight into the bifurcation phenomenon,we consider a one-dimensional toy model with periodic boundary conditions(figure 6(a)).In this model,the potential V takes the values of 0 and 1 alternately in the intervals with lengths L2/2,L1,L2,L3,L4,L3,and L2/2 respectively,i.e.

    The periodic boundary condition guarantees that the system has stronger symmetry and thus the theory will become much simpler.The potential V is composed of two ‘wells’ with lengths L1and 2L3+L4respectively,while there is a short barrier with length L4in the middle of the right well.The distance between the two wells is L2.For convenience,the left and right wells are denoted by W1=[x1,x2]and W2=[x3,x5],respectively.

    We next investigate the localization of the quantum states of the stationary Schr?dinger equation

    Here we require that the left well W1is the longest subinterval with V=0 and the right well W2is the union of two relatively long subintervals with V=0 that are separated by a very short subinterval with V=1.The total length of the two relatively long subintervals should be greater than the length of the longest subinterval so that bifurcation can occur.To ensure that the two wells do not interfere with each other,we further require the distance between the two wells to be large enough.To be more specific,the parameters L1,L2,L3,and L4should satisfy:(i)L1>L3,which grantees that W1is the longest subinterval;(ii)L1<2L3,which grantees that W2is longer than W1;(iii)L4<L3/2,which guarantees that the barrier in W2is short enough;(iv)L2>L1+2L3+L4,which grantees that the two wells are far enough;(v)L1+2L2+2L3+L4=1,which guarantees that the length of the whole interval is 1.

    For the toy model,we illustrate the first eigenmode of the system S given in(24)under different values of K(figure 6(b)).Since the right well is longer than the left one and the barrier in the right well is short,when K is small,the first eigenmode is mainly localized in the right well.Since the left well is the longest subinterval with V=0,when K is large,the first eigenmode is mainly localized in the left well.With the increase of K,the left peak of the eigenmode becomes higher and the right peak becomes lower.At a critical threshold of K,the two peaks are equal in height and bifurcation takes place—the confining subregion of the eigenmode transitions from the right well to the left one.For the first eigenmode u,we introduce

    which represents the height of the left peak relative to the total height of the two peaks.In general,the relative height F of the left peak for the first eigenmode increases with K(figure 6(c)).The critical threshold of K is defined as the value at which F=0.5.At the critical value K=Kc,the relative height F of the left peak changes sharply from a low to a high value,corresponding to the occurrence of bifurcation.

    When L2is large,the two wells are far apart and they have little influence on each other.Hence we can decompose the system S into two subsystems S1and S2with S1corresponding to the left peak of the eigenmode and with S2corresponding to the right peak.Specifically,we introduce two new potentials V1and V2,where V1is obtained from V by setting the values in the right well to 1 and V2is obtained from V by setting the values in the left well to 1(figure 6(a)),i.e.

    Then the original system S can be decomposed into two subsystems S1and S2,which are defined as

    Here the subsystem S1has potential V1and the subsystem S2has potential V2.Intuitively,for the first two eigenmodes of the system S,one is localized in the left well W1and the other is localized in the right well W2.The eigenmode localized in W1(W2)can be approximated by the first eigenmode of subsystem S1(S2).Figure 6(d)illustrates the first two eigenmodes of the original system and the first eigenmodes of the two subsystems.It can be seen that the first eigenmode of each subsystem almost coincides with one of the first two eigenmodes of the original system.Moreover,we illustrate the first two eigenvalues of the original system and the first eigenvalues of the two subsystems under different values of K in figure 6(e).Clearly,the first eigenvalues of the two subsystems serve as good approximations to the first two eigenvalues of the original system.

    We now provide a detailed explanation for the occurrence of bifurcation.From figure 6(e),the first eigenvalue of each subsystem increases with K.However,the eigenvalues of both subsystems increase with K at different rates(see two lines of different slopes in figure 6(e)).When K is small,the eigenvalue of subsystem S2is smaller than that of subsystem S1,and thus the first eigenmode of the original system S is localized in the right well W2.When K is large,the eigenvalue of S2exceeds that of S1,and thus the first eigenmode of S is localized in the left well W1.At the critical threshold K=Kc,the localization subregion of the first eigenmode transitions from W2to W1.In analogy to multimodality,bifurcation occurs when the first eigenvalues of both subsystems are approximately equal,i.e.when the first and the second eigenvalues of the original system are approximately equal.This critical eigenvalue is denoted by λ=λc.

    5.2.Computation of the critical threshold

    Clearly,the critical threshold of K plays a crucial role in the bifurcation phenomenon.Numerically,in order to compute the critical value,we need to solve the original system under many different values of K,which is very time-consuming.Here we provide an analytical theory of the critical threshold and the associated critical eigenvalue.

    To this end,we first simplify the subsystems S1and S2.For each subsystem Si,due to the periodic boundary condition,we can shift the potential Vialong the interval Ω=(0,1)without changing the eigenvalues.In particular,we move the well Wito the middle of the interval(figure 6(a))and the shifted potentials for the two subsystems are defined as

    5.3.Dependence of the critical threshold on model parameters

    We have seen that the critical threshold Kcand critical eigenvalue λcsatisfy the system of algebraic equations(28).Clearly,the value of Kcdepends on the parameters L1,L2,L3,and L4.However,the relationship between them is very complicated.Here we investigate their relationship numerically and reveal the key factors of bifurcation.

    Since L1+2L2+2L3+L4=1,there four parameters only have three degrees of freedom.For simplicity,we define

    where P1represents the total length of the two wells relative to the length of the whole interval,P2represents the length of the left well relative to the total length of the two wells,and P3represents the length of the barrier relative to the length of the right well.The relationship between Kcand P1,P2,and P3is complicated.However,we can obtain an approximate relationship between them using our analytical theory.Figures 7(a)–(c)illustrate the critical threshold Kcas P1,P2,and P3vary.Here we generate different values of P1,P2,or P3randomly from a large swathe of the parameter space and keep the other two parameters invariant.Interestingly,we find that Kcdepends on P1in a power law form as

    To check this,we plotlog(Kc)as a function oflog(P1)in figure 7(a)and find that there is an excellent linear relationship between them with an exceptionally high R2and a slope of -2.Similarly,we find that Kcdepends on P2in an exponential form as(figure 7(b))

    and Kcdepends on P3in a power law form as(figure 7(c))

    These results indicate that Kcdecreases with P1and P3at a power law speed and decreases with P2at an exponential speed.In particular,short potential wells(small P1),a short left well relative to the right well(small P2),and a short barrier in the right well(small P3)are more capable of producing a large critical value of K.The exponential decay of Kcwith P2suggests that the relative length of the left well affects the critical threshold more severely than the other two factors.

    6.Conclusions and discussion

    In the present paper,we focused on some new mathematical aspects of Anderson localization for the quantum states of a Schr?dinger equation with a disordered lattice potential(or more generally,the eigenmodes of a second-order elliptic operator with a random potential).In previous papers[16],the concepts of localization landscape and valley lines have been proposed for symmetric elliptic operators and the Dirichlet boundary condition.Here we generalized these concepts to non-symmetric elliptic operators and general boundary conditions using the probabilistic representation of the solution of a second-order elliptic equation by reflecting diffusions.We showed that the confining subregions of low energy quantum states are perfectly predicted by the peak positions of the landscape and are enclosed by the valley lines of the landscape.

    We also studied the localization behavior of the quantum states when the strength K of disorder is large.In the limit of K →∞,we demonstrated that the eigenmodes must be localized in the region with the lowest potential.When the potential is sampled from the Bernoulli distribution,this region can be decomposed into many connected components(subregions),which are exactly separated by the valley lines of the landscape.Then the original system can be decomposed into many uncoupled local systems in these subregions.We further showed that the spectrum of the original system is composed of the eigenvalues of the local systems.If an eigenvalue of the original system coincides with one of the local eigenvalues in a particular subregion,then the corresponding eigenmode will be localized in that subregion.In particular,if multiple subregions share a common local eigenvalue,then the corresponding eigenmode will have multiple peaks and thus display multimodality.In the onedimensional case,we derived the analytical expression for the probability that the first eigenmode displays multimodality when K is large under the Dirichlet or Neumann boundary condition,which was then used to study the influence of the parameter p on multimodality.

    For Neumann and Robin boundary conditions,we found that(i)compared to the interior of the domain,the eigenmodes tend to be localized on the boundary;(ii)compared to the one-dimensional case,the eigenmodes for a two-dimensional problem are more likely be localized on the boundary,especially when h is small and K is large.To explain this phenomenon,we proposed the concept of extended subregions,which are obtained from the original ones through mirror symmetry along the boundary,and showed that the low energy eigenmodes are more likely to be localized in large extended subregions with strong symmetry.Since(i)the extended subregion should in general be doubled when the original subregion is located on the boundary and(ii)the extended subregion,in general,has stronger symmetry than the original ones,and the eigenmodes tend to be localized on the boundary.Moreover,since(i)the proportion of subregions that need to extend is much higher in the two-dimensional case and(ii)the extended subregion should be magnified four times when the original subregion is located in the corner,the eigenmodes for a twodimensional problem are more likely to be localized on the boundary.In the one-dimensional case,we also derived an explicit expression for the probability that the first eigenmode is localized on the boundary,which was then used to examine the influence of the parameter p on the boundary effect.

    When K is large,the first eigenmode is localized in the extended subregion with the smallest local eigenvalue.In particular,in the one-dimensional case,the first eigenmode is localized in the longest extended subregion with V=0 for both the Dirichlet and Neumann boundary conditions.However,when K is small,the eigenmode is often localized in another region which is the union of some subregions with V=0 that are separated by some short barriers with V=1.At the critical value of K,the localization behavior of the system changes severely and bifurcation occurs.To gain a deeper insight into bifurcation,we consider a one-dimensional toy model whose potential has two wells:the left one is the longest subinterval with V=0 and the right one is the union of two relative long subintervals with V=0 separated by a short barrier with V=1 in between.For the toy model,we showed that bifurcation occurs when its first and second eigenvalues are approximately equal.Based on this property,we further obtained the equations satisfied by the critical threshold K=Kcand the critical eigenvalue λ=λc.The analytical theory is then used to study the dependence of bifurcation on model parameters and reveal the key factors of bifurcation.

    Due to the complexity of Anderson localization,we mainly develop the analytical theory of the boundary effect,multimodality,and bifurcation in the one-dimensional case when the lattice potential is randomly sampled from a Bernoulli distribution.In the study of bifurcation,we consider only a toy model whose potential has two wells with the right well having a short barrier in the middle.We anticipate that our results can be generalized to high dimensions and more general random potentials,and further stimulate work in this area.We also hope that the landscape theory and the probabilistic approach can be applied to high-frequency localization behavior[30,31]and eigenvalue problems in other disciplines such as biology and chemistry[32–34].

    Acknowledgments

    The authors acknowledge the support from National Natural Science Foundation of China with grants No.11871092,No.12131005,and NSAF U1930402.

    Appendix A.Proof of the Filoche–Mayboroda Inequality

    Appendix B.Equivalent system in the extended subregion

    Appendix C.Boundary effect and multimodality

    In the one-dimensional case,the domain Ω=(0,1)is divided uniformly into N subintervals of the same length 1/N.In each subinterval,the potential V is sampled from the Bernoulli distribution with parameter p.The values of V in these subintervals are independent of each other.Suppose that the values of V in the previous two subintervals are 1 and 0,respectively.Let X denote the number of the successive subintervals in which V takes the value of 0.Clearly,the event of X=n means that V=0 in the following n-1 subintervals and V=1 in the nth subinterval.Due to the independence of V in these subintervals,when N is sufficiently large,the random variable X approximately follows the geometric distribution with parameter p,i.e.

    where q=1-p.Similarly,let Y denote the number of successive subintervals in which V takes the value 1.Then when N is sufficiently large,the random variable X approximately follows the geometric distribution with parameter q,i.e.

    In this way,the domain Ω=(0,1)can be decomposed into many subregions with V=0 and V=1 alternatively.Two successive subregions are collectively called a period.Clearly,the mean length of each period is given by

    Therefore,the mean number of periods in the domain is roughly given by M=Npq.Clearly,the domain is composed of subregions with V=0 and V=1 alternatively.Let X1,X2,...,XMdenote the numbers of subintervals included in the successive subregions with V=0.When N is sufficiently large,these random variables are approximately independent.

    C.1.Calculation of boundary probability for the Neumann boundary condition

    Under the Neumann boundary condition,the first eigenmode is localized in the longest extended subregion with V=0 when K is large.Therefore,the event that the first eigenmode is localized on the boundary is equivalent to the event that the longest extended subregion with V=0 appears on the boundary.This event can be classified into the following three situations.

    First,we consider the situation where V(0)=V(1)=0,whose probability is q2.In this case,the lengths of the extended subregions with V=0 are given by 2X1,X2,… ,XM-1,2XM.Then the probability that the longest extended subregion appears on the boundary can be computed explicitly as

    Second,we consider the situation where V(0)=0 and V(1)=1,whose probability is pq.In this case,the lengths of the extended subregions with V=0 are given by 2X1,X2,… ,XM.Then the probability that the longest extended subregion appears on the boundary can be computed explicitly as

    Third,we consider the situation where V(0)=1 and V(1)=0,whose probability is pq.In this case,the lengths of the extended subregions with V=0 are given by X1,X2,… ,2XM.In analogy to the second case,the probability that the longest extended subregion appears on the boundary is given by

    Finally,using the total probability formula,the probability that the first eigenmode is localized on the boundary is given by

    which gives(21)in the main text.

    C.2.Calculation of the multimodal probability for the Dirichlet boundary condition

    Under the Dirichlet boundary condition,the first eigenmode is localized in the longest subregion with V=0 when K is large.If there is a unique longest subregion with V=0,then the first eigenmode is unimodal; if there are two or more longest subregions with V=0,then the first eigenmode is multimodal.Therefore,the event that the first eigenmode is unimodal is equivalent to the event that there is a unique longest subregion with V=0.This probability can be computed explicitly as

    which gives(22)in the main text.

    C.3.Calculation of the multimodal probability for the Neumann boundary condition

    Under the Neumann boundary condition,the first eigenmode is localized in the longest extended subregion with V=0 when K is large.Therefore,the event that the first eigenmode is unimodal is equivalent to the event that there is a unique longest extended subregion with V=0.This event can be classified into the following four situations.

    First,we consider the situation where V(0)=V(1)=0,whose probability is q2.In this case,the lengths of extended subregions with V=0 are given by 2X1,X2,… ,XM-1,2XM.Then the probability that there is a unique longest extended subregion is given by

    where[x]represents the largest integer less than or equal to x.Second,we consider the situation where V(0)=0 and V(1)=1,whose probability is pq.In this case,the length of extended subregions with V=0 are given by 2X1,X2,… ,XM.Then the probability that there is a unique longest extended subregion is

    Third,we consider the situation where V(0)=1 and V(1)=0,whose probability is pq.In this case,the probability that there is a unique longest extended subregion is given by p3=p2.Finally,we consider the situation in which V(0)=V(1)=1,whose probability is p2.In this case,the length of the extended subregions with V=0 are given by X1,X2,… ,XM.Then the probability that there is a unique longest extended subregion is given by p4=1-pD,where pDis the multimodal probability for the Dirichlet boundary condition.In summary,the probability that the first eigenmode is unimodal is given by

    which gives(23)in the main text.

    Appendix D.Critical threshold of bifurcation

    In the main text,we have shown that bifurcation occurs when the first eigenvalues of the two subsystemsandare equal.

    D.1.Eigenvalues of the subsystem ~S1

    D.2.Eigenvalues of the subsystem 2

    夜夜爽天天搞| 亚洲av中文字字幕乱码综合| 亚洲专区国产一区二区| 日本五十路高清| 美女高潮喷水抽搐中文字幕| 亚洲成av人片免费观看| 看免费av毛片| 国产免费男女视频| 无限看片的www在线观看| 国产午夜精品久久久久久| 国产精品免费一区二区三区在线| 男女床上黄色一级片免费看| 免费观看人在逋| av视频在线观看入口| 午夜福利视频1000在线观看| 老司机午夜十八禁免费视频| 精品午夜福利视频在线观看一区| x7x7x7水蜜桃| 国产单亲对白刺激| 国产又黄又爽又无遮挡在线| 韩国av一区二区三区四区| 手机成人av网站| 精品少妇一区二区三区视频日本电影| 亚洲色图av天堂| 久久精品91蜜桃| 18禁观看日本| 美女大奶头视频| 国产高清激情床上av| 亚洲男人天堂网一区| 亚洲无线在线观看| 亚洲黑人精品在线| 久久午夜亚洲精品久久| bbb黄色大片| 99久久综合精品五月天人人| 日韩欧美国产在线观看| 久久 成人 亚洲| 亚洲专区字幕在线| 最好的美女福利视频网| 性色av乱码一区二区三区2| 首页视频小说图片口味搜索| 成人国产综合亚洲| 香蕉av资源在线| 亚洲精华国产精华精| 欧美极品一区二区三区四区| 一级黄色大片毛片| 五月玫瑰六月丁香| 日本一区二区免费在线视频| 91九色精品人成在线观看| 嫩草影视91久久| 久久久久久久午夜电影| 女人被狂操c到高潮| 国产午夜精品论理片| 午夜福利在线在线| 中文字幕熟女人妻在线| 性欧美人与动物交配| 欧美日韩国产亚洲二区| 成人三级黄色视频| 亚洲aⅴ乱码一区二区在线播放 | 五月伊人婷婷丁香| 国产高清有码在线观看视频 | 久久久久久免费高清国产稀缺| 久久中文看片网| 99riav亚洲国产免费| 中文字幕最新亚洲高清| 欧美日韩黄片免| 五月玫瑰六月丁香| 黄色片一级片一级黄色片| 日日摸夜夜添夜夜添小说| 免费在线观看影片大全网站| 午夜福利在线在线| 精品少妇一区二区三区视频日本电影| 色尼玛亚洲综合影院| 精品日产1卡2卡| 免费在线观看日本一区| 亚洲一区中文字幕在线| 日韩欧美 国产精品| 啦啦啦免费观看视频1| 国产精品一区二区三区四区久久| 国产精品av久久久久免费| 高清毛片免费观看视频网站| e午夜精品久久久久久久| 最近最新中文字幕大全电影3| 亚洲专区字幕在线| 99久久久亚洲精品蜜臀av| 欧美黑人精品巨大| 欧美成人一区二区免费高清观看 | 久久 成人 亚洲| 日韩精品青青久久久久久| 一个人免费在线观看的高清视频| 女人被狂操c到高潮| 在线十欧美十亚洲十日本专区| 国产一区二区在线观看日韩 | 精品国产美女av久久久久小说| 久久久久国内视频| 国产激情偷乱视频一区二区| 精华霜和精华液先用哪个| 99国产精品一区二区蜜桃av| 夜夜躁狠狠躁天天躁| 每晚都被弄得嗷嗷叫到高潮| 国产黄色小视频在线观看| 亚洲一区高清亚洲精品| 国产精品免费视频内射| 久久香蕉国产精品| 欧美一级毛片孕妇| 亚洲中文字幕日韩| a级毛片a级免费在线| 两人在一起打扑克的视频| 不卡av一区二区三区| 一a级毛片在线观看| 真人做人爱边吃奶动态| 久久人人精品亚洲av| 亚洲一码二码三码区别大吗| xxx96com| 亚洲全国av大片| 免费看美女性在线毛片视频| 精品久久久久久成人av| 亚洲精品在线美女| 久久中文字幕人妻熟女| 欧美日韩福利视频一区二区| 午夜视频精品福利| 亚洲av成人不卡在线观看播放网| 国产精品久久久av美女十八| 少妇裸体淫交视频免费看高清 | 国产精品 欧美亚洲| 欧美性长视频在线观看| 男人舔奶头视频| 亚洲精品美女久久久久99蜜臀| 午夜福利在线观看吧| 桃色一区二区三区在线观看| 成人三级做爰电影| 可以在线观看毛片的网站| 岛国在线观看网站| 亚洲人成网站在线播放欧美日韩| 一区二区三区高清视频在线| 午夜福利欧美成人| www.999成人在线观看| 精品免费久久久久久久清纯| 亚洲欧美日韩高清专用| 欧美日韩亚洲国产一区二区在线观看| 国产探花在线观看一区二区| 变态另类成人亚洲欧美熟女| 久久久国产欧美日韩av| 九色成人免费人妻av| 亚洲精品久久国产高清桃花| 欧美成人免费av一区二区三区| 女人爽到高潮嗷嗷叫在线视频| 成人av一区二区三区在线看| 美女黄网站色视频| 中亚洲国语对白在线视频| 女人被狂操c到高潮| 黑人操中国人逼视频| 999久久久国产精品视频| 老汉色av国产亚洲站长工具| 午夜福利在线观看吧| 一级毛片女人18水好多| 午夜福利欧美成人| bbb黄色大片| 在线十欧美十亚洲十日本专区| 国产私拍福利视频在线观看| aaaaa片日本免费| 婷婷亚洲欧美| 深夜精品福利| 亚洲欧美激情综合另类| 精品国产亚洲在线| 性色av乱码一区二区三区2| 午夜久久久久精精品| 少妇粗大呻吟视频| 男女之事视频高清在线观看| 欧美最黄视频在线播放免费| 午夜a级毛片| 欧美乱妇无乱码| 亚洲av片天天在线观看| 国产私拍福利视频在线观看| 性欧美人与动物交配| 日本成人三级电影网站| 90打野战视频偷拍视频| svipshipincom国产片| 两个人的视频大全免费| 欧美黑人巨大hd| 看免费av毛片| 精华霜和精华液先用哪个| 男插女下体视频免费在线播放| 亚洲中文日韩欧美视频| 99热只有精品国产| 首页视频小说图片口味搜索| 欧美一区二区精品小视频在线| 免费观看人在逋| 中文在线观看免费www的网站 | 一卡2卡三卡四卡精品乱码亚洲| 老汉色∧v一级毛片| 极品教师在线免费播放| svipshipincom国产片| 18禁观看日本| 亚洲国产欧美一区二区综合| 国产伦一二天堂av在线观看| 老汉色av国产亚洲站长工具| 国内久久婷婷六月综合欲色啪| 久久精品人妻少妇| 一夜夜www| 在线看三级毛片| 九九热线精品视视频播放| 变态另类丝袜制服| 亚洲黑人精品在线| 亚洲欧美精品综合久久99| 国产伦在线观看视频一区| 免费看a级黄色片| 国产精品自产拍在线观看55亚洲| 午夜激情福利司机影院| 精品高清国产在线一区| 国产爱豆传媒在线观看 | 成熟少妇高潮喷水视频| 欧美绝顶高潮抽搐喷水| 国产在线观看jvid| 久久精品国产亚洲av高清一级| 亚洲av电影不卡..在线观看| 亚洲欧美精品综合一区二区三区| 欧美午夜高清在线| 很黄的视频免费| 亚洲欧美精品综合一区二区三区| 黄频高清免费视频| 欧美久久黑人一区二区| 欧美av亚洲av综合av国产av| 香蕉国产在线看| 国产免费男女视频| 中文字幕熟女人妻在线| 欧美黄色淫秽网站| 亚洲精品av麻豆狂野| 国产精品久久电影中文字幕| 精品高清国产在线一区| 两个人的视频大全免费| 国产精品精品国产色婷婷| 亚洲一区高清亚洲精品| 免费无遮挡裸体视频| 国产高清视频在线观看网站| e午夜精品久久久久久久| 成人国产综合亚洲| 欧美午夜高清在线| 一个人观看的视频www高清免费观看 | 久久午夜亚洲精品久久| 91老司机精品| 久99久视频精品免费| 国产又黄又爽又无遮挡在线| 国产精品99久久99久久久不卡| 久久中文看片网| 免费在线观看黄色视频的| 又黄又粗又硬又大视频| 国产精品久久久av美女十八| 狠狠狠狠99中文字幕| 校园春色视频在线观看| 国产乱人伦免费视频| 亚洲,欧美精品.| 国产欧美日韩一区二区精品| 亚洲国产欧美网| 此物有八面人人有两片| 久久久国产精品麻豆| 精品一区二区三区四区五区乱码| 久久久久久久久中文| 无遮挡黄片免费观看| 免费高清视频大片| 欧美乱妇无乱码| 天堂av国产一区二区熟女人妻 | 亚洲熟妇中文字幕五十中出| 亚洲欧美日韩高清专用| 亚洲国产欧美网| 美女午夜性视频免费| 亚洲精品久久国产高清桃花| 国产人伦9x9x在线观看| 精品人妻1区二区| 18禁美女被吸乳视频| 99热6这里只有精品| 免费无遮挡裸体视频| 久久99热这里只有精品18| 亚洲五月天丁香| 人人妻人人看人人澡| 两个人免费观看高清视频| 国产麻豆成人av免费视频| 曰老女人黄片| 在线免费观看的www视频| 国产av又大| 天天一区二区日本电影三级| 国产不卡一卡二| 欧美成人一区二区免费高清观看 | 亚洲精品国产一区二区精华液| www日本在线高清视频| 宅男免费午夜| 精品欧美国产一区二区三| 熟妇人妻久久中文字幕3abv| 美女午夜性视频免费| 亚洲国产欧美一区二区综合| 精品午夜福利视频在线观看一区| 亚洲欧美日韩东京热| 亚洲av五月六月丁香网| 国内揄拍国产精品人妻在线| 免费在线观看亚洲国产| 老熟妇乱子伦视频在线观看| АⅤ资源中文在线天堂| 亚洲一区高清亚洲精品| 欧美日韩精品网址| 国产精品九九99| 成年女人毛片免费观看观看9| 午夜激情av网站| 欧美+亚洲+日韩+国产| 精品免费久久久久久久清纯| 欧美色欧美亚洲另类二区| 麻豆一二三区av精品| 久久久国产欧美日韩av| 国产亚洲av嫩草精品影院| 国产激情偷乱视频一区二区| 18禁黄网站禁片免费观看直播| 久久久久九九精品影院| 91九色精品人成在线观看| 国产精品久久久人人做人人爽| 久久午夜亚洲精品久久| 韩国av一区二区三区四区| 欧美av亚洲av综合av国产av| 久久久久久大精品| 亚洲欧洲精品一区二区精品久久久| 亚洲午夜精品一区,二区,三区| 国产男靠女视频免费网站| 狠狠狠狠99中文字幕| 成人精品一区二区免费| 午夜福利18| 中国美女看黄片| 久久久久国产一级毛片高清牌| 大型黄色视频在线免费观看| 日本a在线网址| 级片在线观看| 91字幕亚洲| 国产精品久久久久久人妻精品电影| 99国产精品99久久久久| 搡老岳熟女国产| 一边摸一边抽搐一进一小说| 欧美最黄视频在线播放免费| 草草在线视频免费看| 免费在线观看成人毛片| 亚洲国产高清在线一区二区三| 美女高潮喷水抽搐中文字幕| 亚洲avbb在线观看| 亚洲精品中文字幕在线视频| 国产一区二区三区在线臀色熟女| 亚洲avbb在线观看| 一进一出抽搐gif免费好疼| 欧美zozozo另类| 欧美高清成人免费视频www| 99精品欧美一区二区三区四区| 欧美中文日本在线观看视频| 99国产精品一区二区三区| 色综合亚洲欧美另类图片| 欧美人与性动交α欧美精品济南到| 国产免费男女视频| 老汉色av国产亚洲站长工具| 夜夜夜夜夜久久久久| 久热爱精品视频在线9| 日韩 欧美 亚洲 中文字幕| 午夜精品在线福利| 成年版毛片免费区| 国产伦人伦偷精品视频| 亚洲av片天天在线观看| 久久久水蜜桃国产精品网| 九色国产91popny在线| 男人舔女人下体高潮全视频| 成人国语在线视频| 伊人久久大香线蕉亚洲五| 欧美成狂野欧美在线观看| 国产精品98久久久久久宅男小说| 51午夜福利影视在线观看| 欧美乱码精品一区二区三区| 变态另类丝袜制服| 黄色a级毛片大全视频| 色综合亚洲欧美另类图片| 美女扒开内裤让男人捅视频| 亚洲乱码一区二区免费版| 99久久国产精品久久久| 动漫黄色视频在线观看| www.www免费av| 亚洲精品一卡2卡三卡4卡5卡| 亚洲av熟女| 天天一区二区日本电影三级| 叶爱在线成人免费视频播放| 看黄色毛片网站| АⅤ资源中文在线天堂| 久久久久性生活片| 小说图片视频综合网站| 十八禁网站免费在线| 老司机午夜福利在线观看视频| 妹子高潮喷水视频| 欧美成人一区二区免费高清观看 | 曰老女人黄片| 国产精品免费一区二区三区在线| 亚洲一区高清亚洲精品| 国产高清视频在线播放一区| 亚洲五月天丁香| 9191精品国产免费久久| 日本精品一区二区三区蜜桃| 最新在线观看一区二区三区| 精品福利观看| 91av网站免费观看| 黄色a级毛片大全视频| avwww免费| 黄片大片在线免费观看| 啦啦啦韩国在线观看视频| 久久国产乱子伦精品免费另类| 成年人黄色毛片网站| 国产精品一区二区三区四区免费观看 | 韩国av一区二区三区四区| 亚洲激情在线av| 精品欧美一区二区三区在线| 精品国产超薄肉色丝袜足j| 在线观看美女被高潮喷水网站 | 久久午夜综合久久蜜桃| 精华霜和精华液先用哪个| 日本撒尿小便嘘嘘汇集6| 老汉色∧v一级毛片| 久久久久九九精品影院| 色综合婷婷激情| 深夜精品福利| 久久伊人香网站| 18禁观看日本| 99国产精品一区二区蜜桃av| 午夜久久久久精精品| 亚洲国产精品sss在线观看| 一本一本综合久久| 激情在线观看视频在线高清| 舔av片在线| 天堂av国产一区二区熟女人妻 | 岛国在线免费视频观看| 久久 成人 亚洲| 欧美黑人精品巨大| 欧美色欧美亚洲另类二区| 免费在线观看影片大全网站| 亚洲精华国产精华精| 欧美日韩国产亚洲二区| 亚洲国产欧美人成| 黄色丝袜av网址大全| 亚洲自拍偷在线| 操出白浆在线播放| 国产99白浆流出| 中国美女看黄片| 色老头精品视频在线观看| 国产亚洲欧美在线一区二区| 亚洲男人的天堂狠狠| 成人国产一区最新在线观看| 禁无遮挡网站| 日本三级黄在线观看| 一本综合久久免费| 琪琪午夜伦伦电影理论片6080| 淫妇啪啪啪对白视频| 国产成人系列免费观看| 国产精品久久视频播放| 久久久久免费精品人妻一区二区| 99热只有精品国产| 一卡2卡三卡四卡精品乱码亚洲| 好男人电影高清在线观看| 久久久久国内视频| 岛国在线观看网站| 国产精品日韩av在线免费观看| 国产av一区在线观看免费| 18禁黄网站禁片免费观看直播| 亚洲人成伊人成综合网2020| 不卡一级毛片| 中国美女看黄片| 日韩精品中文字幕看吧| 黄色视频,在线免费观看| 国产1区2区3区精品| 宅男免费午夜| 欧美另类亚洲清纯唯美| 一级a爱片免费观看的视频| 在线看三级毛片| 亚洲精品在线观看二区| 在线观看日韩欧美| 精品国产亚洲在线| 欧美成狂野欧美在线观看| 又紧又爽又黄一区二区| 69av精品久久久久久| 亚洲一区中文字幕在线| 99国产综合亚洲精品| 黄色毛片三级朝国网站| 美女高潮喷水抽搐中文字幕| 久久人妻av系列| 天天躁狠狠躁夜夜躁狠狠躁| 国产成人av教育| 欧美一区二区精品小视频在线| 午夜精品在线福利| 免费无遮挡裸体视频| av视频在线观看入口| 99久久99久久久精品蜜桃| 麻豆国产97在线/欧美 | 国产伦在线观看视频一区| 九九热线精品视视频播放| 两个人看的免费小视频| 久久中文看片网| 制服人妻中文乱码| 免费电影在线观看免费观看| 成人亚洲精品av一区二区| 美女大奶头视频| 国产成人欧美在线观看| 国产一区二区在线观看日韩 | 在线视频色国产色| 国产视频内射| 99国产极品粉嫩在线观看| 18禁黄网站禁片午夜丰满| 午夜两性在线视频| 欧美在线一区亚洲| 国产午夜精品论理片| 嫩草影院精品99| 中文字幕av在线有码专区| 欧美性猛交黑人性爽| 成人欧美大片| 长腿黑丝高跟| 看免费av毛片| АⅤ资源中文在线天堂| 两个人看的免费小视频| 久久精品国产清高在天天线| 久久人妻av系列| 久久久久久久久久黄片| 欧美中文综合在线视频| 久久精品aⅴ一区二区三区四区| 国产日本99.免费观看| 日韩免费av在线播放| 欧美人与性动交α欧美精品济南到| 人妻夜夜爽99麻豆av| 久9热在线精品视频| 嫩草影视91久久| av片东京热男人的天堂| 2021天堂中文幕一二区在线观| 亚洲人成电影免费在线| 日韩中文字幕欧美一区二区| 桃红色精品国产亚洲av| 国产精品久久久久久久电影 | 又粗又爽又猛毛片免费看| 在线观看免费午夜福利视频| 麻豆久久精品国产亚洲av| 国产亚洲av嫩草精品影院| 12—13女人毛片做爰片一| 在线观看日韩欧美| 此物有八面人人有两片| 两性夫妻黄色片| 黄色 视频免费看| 制服诱惑二区| 我的老师免费观看完整版| 国产熟女午夜一区二区三区| 午夜视频精品福利| 中文字幕精品亚洲无线码一区| 又粗又爽又猛毛片免费看| 免费在线观看黄色视频的| 国产男靠女视频免费网站| 精品欧美一区二区三区在线| 日韩 欧美 亚洲 中文字幕| 在线播放国产精品三级| 欧美黄色片欧美黄色片| 亚洲在线自拍视频| 免费在线观看完整版高清| 精品国产乱子伦一区二区三区| www.999成人在线观看| 久久 成人 亚洲| a级毛片a级免费在线| 久久精品国产综合久久久| 中文字幕熟女人妻在线| 国产伦在线观看视频一区| 性色av乱码一区二区三区2| 久久精品国产综合久久久| 欧美性猛交╳xxx乱大交人| 国产亚洲精品一区二区www| 一夜夜www| 国产精品1区2区在线观看.| 久久精品成人免费网站| 成人特级黄色片久久久久久久| 久久久国产成人免费| 国产精品久久久久久亚洲av鲁大| 久久精品国产99精品国产亚洲性色| 日韩精品免费视频一区二区三区| 欧美国产日韩亚洲一区| 国产午夜精品久久久久久| 国产精品1区2区在线观看.| 日本撒尿小便嘘嘘汇集6| 岛国视频午夜一区免费看| 757午夜福利合集在线观看| 一级作爱视频免费观看| 国产av一区二区精品久久| 女人爽到高潮嗷嗷叫在线视频| 一本一本综合久久| 亚洲精品久久成人aⅴ小说| 亚洲中文日韩欧美视频| 黄片大片在线免费观看| 啦啦啦免费观看视频1| 真人做人爱边吃奶动态| 久久久久九九精品影院| 啦啦啦免费观看视频1| 波多野结衣巨乳人妻| 99riav亚洲国产免费| 丁香欧美五月| 成人精品一区二区免费| 黄片大片在线免费观看| 日韩 欧美 亚洲 中文字幕| 日韩欧美在线乱码| 久久午夜亚洲精品久久| 亚洲片人在线观看| 香蕉久久夜色| 午夜精品久久久久久毛片777| 国产精品久久久人人做人人爽| 亚洲专区字幕在线| 久99久视频精品免费| 18禁国产床啪视频网站| 岛国在线观看网站| 老司机午夜十八禁免费视频| 天堂动漫精品| 成年版毛片免费区| 久久久久久亚洲精品国产蜜桃av| 久久中文看片网| 亚洲午夜精品一区,二区,三区| 日本撒尿小便嘘嘘汇集6| 久久精品影院6| 一二三四在线观看免费中文在| 国产区一区二久久| 久久久水蜜桃国产精品网|