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

    Thermodynamics and microstructures of Euler-Heisenberg black hole in a cavity

    2023-10-11 05:30:46QinYuQiXuandJunTao
    Communications in Theoretical Physics 2023年9期

    Qin Yu,Qi Xu and Jun Tao

    Center for Theoretical Physics,College of Physics,Sichuan University,Chengdu,610065,China

    Abstract Euler–Heisenberg black holes with quantum electrodynamics(QED) correction are embraced by a cavity in this paper,which serves as a boundary of the black hole spacetime and contributes to the equilibrium of the system.We explore the thermodynamic properties of the black hole,including the phase transitions and phase structures.The small/large black hole phase transition occurs for a negative QED parameter,while the reentrant phase transition can be observed for a small positive QED parameter.Then the thermodynamic geometry is investigated to diagnose microscopic interactions of black hole thermodynamic systems.For the reentrant phase transition,the small black holes are dominated by repulsion for the first-order coexistence curve,while the interaction between the small black hole molecules could be attractive or repulsive for the small/large black hole phase transition.

    Keywords: black hole thermodynamics,phase transition,microstructure

    1.Introduction

    The investigation of black hole thermodynamics has shown crucial importance to quantum gravity and general relativity,based on the pioneering works by Hawking,Bekenstein and others[1–4].Between the thermal anti-de Sitter(AdS)space and black holes,the Hawking-Page phase transition was observed in Schwarzschild-AdS black holes for the first time[5].Inspired by the AdS/CFT correspondence between gravity systems and the conformal field theory[6–8],the Hawking-Page phase transition can be interpreted as a confinement/deconfinement transition of gauge field[9].Thermodynamic properties and phase structures of black holes were then intensively researched and have gained significant attention in recent years.

    The relation between black hole thermodynamics and conformal field theory deserves further discussion.For a fixed temperature at infinity,black holes are thermodynamically unstable in asymptotically flat space[10].Enclosing black holes with a cavity is a useful way to give solutions of it.Considering a heat bath around the cavity,the surface gives a fixed temperature.Black holes can be thermodynamically stable then,as in the case of AdS space,the horizon of the black hole can be near the fixed temperature point [11].This method is useful to investigate holography in asymptotically flat and AdS space as well.The Schwarzschild black hole in a cavity was researched by York for the first time [12].Compared with the Schwarzschild-AdS black holes,there is an overall similarity between the phase structures of these two cases.The Reissner-Nordstr o¨m(RN) black holes in a cavity showed a Hawking-Page phase transition in the grand canonical ensemble [13] and a van der Waals phase transition in the canonical ensemble [14],and it is closely similar to the phase structures of an RN-AdS black hole.There are almost no differences in phase structures between Gauss-Bonnet black holes in a cavity and Gauss-Bonnet-AdS black holes as well [10,15].However,some differences exist,such as phase structures of Born-Infeld black holes[16–18].For example,a large/small/large black hole reentrant phase transition[19–21]could be observed for Born-Infeld-AdS black holes,while it can not occur for any Born-Infeld black hole in a cavity.For black holes in a cavity and with the AdS boundary,the similarity or difference of phase structures may have interesting implications for black hole thermodynamics and holography[22].These investigations have aroused our interest to explore black holes in a cavity.

    In addition,the statistical description of the black hole microstates is still not entirely clear.Following the pioneering work by Weinhold [23],to represent the thermodynamic fluctuation theory,Ruppeiner constructed a Riemannian thermodynamic entropy metric and developed a systematic method to calculate the Ricci curvature scalar R of the metric[24].The sign of a Ruppeiner invariant is relative to the interparticle interaction: R >0 denotes repulsion while R <0 denotes attraction,and R=0 indicates no interaction [25].Whereafter,the Ruppeiner approach was applied in various black holes [26–46].From the Ruppeiner geometry,for the black holes with reentrant phase transition,the small black hole behaves like a Bose gas.However,for the black holes with usual first-order phase transitions,it behaves like a quantum anyon gas [45].Thus the phase transition may be affected by the microstructure of black holes [46].These studies show good examples of combining the phase transition with the microstructure of black holes.

    To understand the thermodynamics and microstructure of black holes,we pay close attention to nonlinear electrodynamics(NLED)and its extensive research.After first being introduced by Born and Infeld [47],nonlinear electrodynamics can be considered in the study of black holes.The Euler–Heisenberg electrodynamics is a particular example of nonlinear electrodynamics as well [48].Schwinger reformulated the modified Lagrange function for constant fields within the quantum electrodynamics (QED) framework [49].It is of particular interest to combine Euler–Heisenberg electrodynamics with the study of black holes,such as the asymptotically flat black holes in Euler–Heisenberg electrodynamics [50].Recently,Euler–Heisenberg black holes in AdS space with QED effects were considered to investigate the influence on black hole phase transitions and Ruppeiner geometry in the extended phase space [51].Some other relevant works have been considered in the literature[52–62].Inspired by this research,we investigate the phase structures of Euler–Heisenberg black holes in a cavity and find some interesting phenomena,especially the reentrant phase transition.Furthermore,from a microscopic point of view,the interaction is repulsive between the small black hole molecules in the reentrant phase transition case.

    The structure of this paper is as follows.In section 2,we derive the thermodynamics of Euler–Heisenberg black holes in a cavity.In section 3,the influence of a QED parameter and charge on black hole phase transitions is analyzed.Then different phase structures and phase diagrams are investigated for positive and negative QED parameters,respectively.In section 4,the microstructure of the black hole is studied by Ruppeiner geometry.Finally,we summarize the findings and make some discussions in section 5.Here we set G=?=c=1 for simplicity in this paper.

    2.Thermodynamics

    The metric of the static spherically symmetric charged Euler–Heisenberg black hole can be written as [55]

    where

    The parameter M is the ADM mass which can be obtained by solving the equation f(rh)=0,while rhis the event horizon of the black hole.The parameter q is the electric charge and a represents the strength of the QED correction[54].The metric reduces to the RN black hole case when the QED parameter a is set to be zero.

    By taking the Euclidean time τ=it,this metric turns into the Euclidean form.On a spacetime manifold M with a timelike boundary?M at r=rB,the Euclidean action is given by[12],

    The second term is related to the Gibbons-Hawking-York surface term,where γ is the metric on the boundary,K is the trace of the extrinsic curvature,and K0is subtracted to ensure the Gibbons-Hawking-York term vanishes in flat spacetime.The Euclidean action for the nonlinear electromagnetic field can be expressed as

    whereL(s,p)is the NLED Lagrangian that depends on two independent nontrivial scalars,

    The Fμν=?μAν-?νAμ,Aμis the gauge potential,and ∈μνρσis the completely antisymmetric Lorentz tensor.The surface term can be obtained by varying IEM,

    The total Euclidean action I is the sum ofgI,IEM,Isurf.Following the derivations of Euclidean action for an NLEDcharged black hole inside a cavity with the radius rB[16],the total Euclidean action takes the form as

    whereS=is the entropy of this black hole,and T is the temperature on the cavity boundary.The thermal energy of an NLED black hole with the cavity radius rBtakes the form of[10]

    and we can obtain the ADM mass by extending the cavity radius to infinity,

    The Lagrangian of the Euler–Heisenberg NLED is [48,55]

    For a →0 we recover the Maxwell electrodynamicsL=s,which meets the condition in equation (7).The thermal energy of the Euler–Heisenberg black hole takes the same form as equation (9).Then the temperature of Euler–Heisenberg black hole systems can be written as [17]

    and the Hawking temperature Thof Euler–Heisenberg black holes is

    It is natural to establish the first law of thermodynamics

    The electric potential can be determined by taking the derivative of the thermal energy with respect to the charge,which can be expressed as

    The phase transition can be characterized as the physical process of transition between one phase of a medium and another,as a result of changes in external conditions,such as temperature or pressure.The phase transitions can be classified based on the continuity of the free energy.For example,if the free energy is discontinuous,it is a zerothorder phase transition [63].If the free energy is continuous but its first derivative is discontinuous,it is a first-order phase transition.The investigation of black hole phase transitions usually requires the analysis of critical points,which can be determined by the following equations

    However,the analytical solution of critical points is difficult to be derived from the temperature of the Euler–Heisenberg black hole in a cavity,which impels us to solve it numerically.The type of phase structure can be roughly determined by the number of critical points,and the critical point can be combined with free energy for analysis.The free energy of black holes in the fixed charge ensemble takes the form as

    Several phases exist if there is more than one value of the free energy for a given temperature.Among these phases,the phase with the lowest free energy is globally stable.Then phase transitions can be judged by investigating whether there are multiple branches with the lowest free energy.At the critical point,the first derivative of the free energy curve is continuous,while the second derivative is discontinuous.

    In the following sections,we will further study the phase transition of the Euler–Heisenberg black hole in a cavity with different QED parameters.The microstructure of this black hole will be discussed by the Ruppeiner geometry as well.

    3.Phase transition and phase structure

    In this section,we shall study the phase transition of an Euler–Heisenberg black hole in a cavity.The phase behaviors of black hole thermodynamic systems will be discussed for different QED parameters.

    Based on features of the free energy,we find that phase transitions depend on the QED parameter a and charge q with a fixed cavity radius.Setting rB=3 without loss of generality,we show the phase transition types for different q and a in figure 1.The white region in figure 1(a) corresponds to the system without phase transition.In the green region I,a single first-order phase transition occurs.The reentrant phase transition exists for a >0,which is shown as the red region II.The black hole(BH)jumps from a stable large BH phase to a small one through a zeroth-order phase transition,and finally returns to the large BH phase through a first-order phase transition.The blue dashed line is given by the unstable critical point of charge for a >0,which is not related to the occurrence of phase transitions.

    There is a maximum value of the QED parameter for the existence of phase transitions,which can be labeled asamax.It is similar to the result in the Euler–Heisenberg-AdS black holes,where[51].Fora>amax,there is no real number solution from equation (16),so phase transitions can not exist for any value of charge q.It is difficult to obtain the solution analytically for the Euler–Heisenberg black hole in a cavity,which induces us to solve it numerically.The maximum value of the QED parameter allowing for phase transitions is determined by the cavity radius,and we can obtainamax=2.0848for rB=3 here.

    For a fixed QED parameter,the BH phase undergoes different phase transitions as the charge increases.For a <0,a first-order phase transition happens,where the small BH phase turns into a large one as the temperature increases.For a >0,we take a=0.08 as an example in figure 1(b),which will be discussed in section 3.1 as well.The BH phase remains stable for a small charge 0 <q <q1.When the value of the charge increases to the range of the red region q1<q <q2,the reentrant large/small/large BH phase transition occurs.After entering the green region,the reentrant phase transition disappears for q2<q <qc2,and the single first-order phase transition takes place.For q >qc2,there is no phase transition and one can no longer distinguish black holes between the small and large branches.

    Fig.1.Regions of different phase transitions for q and a.Left Panel(a):The white region represents the system without phase transition(PT).The green region I corresponds to the case that a single first-order phase transition occurs.The reentrant phase transition happens in the red region II.The blue dashed line corresponds to the unstable critical point for a >0,which does not present the occurrence of phase transitions.Right Panel(b):The black dashed line a=0.08 intersects the blue dashed line at qc1,and intersects the boundary of regions at q1,q2,qc2,respectively.The cavity radius is set to be rB=3.

    Fig.2.Temperature curve of the Euler–Heisenberg black hole in a cavity versus horizon radius for different values of q. Left Panel(a):The parameter a=-1,and the charge q is equal to 0.3,0.5,0.6909 and 0.8,respectively;Right Panel(b):The parameter a=0.08,and the charge q is equal to 0.1,0.1839,0.3,0.7098 and 0.8,respectively.The cavity radius rB is set to be 3.

    Furthermore,the critical points of the Euler–Heisenberg black holes in a cavity can be investigated by analyzing the temperature curve.The temperature T versus the horizon radius rhfor negative and positive QED parameters is plotted in figure 2,where the cavity radius rBis set to be 3.For a=-1,we obtain a critical charge qc=0.6909 from equation (16);For a=0.08,we obtain two critical charges qc1=0.1839 and qc2=0.7098.

    There are two extrema on the curve for q <qcas shown in figure 2(a).In this situation,one can divide the curve into three branches.The two sides of the curve correspond to the small BH (left side) and large BH (right side),respectively,while the middle branch represents the intermediate one.The heat capacity can be given as,which indicates that the positive and negative slopes of the temperature curve correspond to stable and unstable branches,respectively.Furthermore,two extrema will get closer as the charge increases,and then coincide at the critical point (red point).For q >qc,there is no extremal point.

    There are three extrema when qc1<q <qc2in figure 2(b),and two of them with larger rhwill coincide at the critical point as the charge q increases.In this case,there are four BH branches,including an unstable smaller branch,a stable small branch,an unstable intermediate branch and a stable large branch,which is a signal of possible phase transitions.For q <qc1,there are one unstable BH branch and one stable branch,thus no phase transition occurs.There is no phase transition for q >qc2either.

    Following the discussions above,the Euler–Heisenberg black hole in a cavity undergoes different phase transitions for negative and positive QED parameters,and we will discuss these two cases in the following sections.

    3.1.a >0 case

    As aforementioned,there will be no phase transition in any case whena>amax.For a small positive QED parametera<amax,there are two critical points,as shown in figure 2(b).In this section,we take a=0.08 and rB=3 for example,then the value of qc1,q1,q2and qc2can be obtained,which is 0.1839,0.2805,0.2907 and 0.7098,respectively.The phase structures for different charges are depicted in figure 3.

    Fig.3.The free energy versus the temperature. Upper Left Panel(a):No phase transition occurs for qc1 <q <q1,where we set q=0.279.Upper Right Panel(b):The reentrant phase transition for q1 <q <q2,where we set q=0.285.Lower Left Panel(c):The firstorder phase transition for q2 <q <qc2,and we choose q=0.3 as an example. Lower Right Panel(d):No phase transition occurs for q >qc2,where small and large BH phases can not be distinguished.We set q=0.8 here.The dashed curves represent unstable branches,while the solid curves stand for stable or metastable BHs.We set a=0.08 and rB=3.

    In figure 3(a),the analysis of the global minimum free energy shows that there is only one stable BH phase for qc1<q <q1and no phase transition occurs,which is the same as the case for a small charge.Here we take q=0.279 for example.

    The reentrant phase transition,a typical phase behavior including a zeroth-order phase transition and a first-order phase transition,occurs for q1<q <q2as shown in figure 3(b).Here we set the charge q=0.285.The solid curves stand for stable or metastable BHs,which can be judged that the former owns the lowest free energy at the same temperature.The dashed curves represent unstable BH branches.The black dashed line corresponds to the temperature of the beginning of the stable small BH branch (blue solid curve),which is larger than the initial temperature of the stable large BH branch (green solid curve) in this case.The phase behavior here is a zeroth-order phase transition,where the black hole system jumps from the stable large BH phase to the stable small BH phase along the black dashed line.Consequently,it indicates sudden changes for not only the horizon radius but also the free energy of the black hole.In addition,by following the rule of the lowest free energy,the black hole system undergoes a first-order phase transition,namely the small/large BH phase transition.There is only a sudden change in the horizon radius of the black hole as a nondifferentiable point in the F-T diagram.Analyzing the phase behavior of stable BH branches,we summarize the process of reentrant phase transition in through a zeroth-order phase transition from stable large BH phase to stable small BH phase,and then it returns to stable large BH phase through a first-order phase transition,namely the reentrant large/small/large BH phase transition.

    For q2<q=0.3 <qc2,we find that there is a first-order phase transition from the stable small BH phase to the stable large BH phase,as shown in figure 3(c).The other branches with positive heat capacity are metastable phases,corresponding to superheated small BH and supercooled large BH.With the continuous increase of charge,the swallow tail will disappear and shrink to a critical point q=qc2.

    For q=0.8 >qc2in figure 3(d),there is no phase transition,and one can not distinguish the large BH phase and the small BH phase anymore.

    To investigate the process from the distinguishable BH phase to the indistinguishable one,we take Δrh=rhl-rhsversus the temperature and charge for a=0.08 and rB=3 in figure 4.The parameter rhlis the radius of the event horizon for the large BH branch and rhsis the small one.The Δrhrepresents the sudden change of the horizon radius when the black hole comes up with a first-order phase transition from a small BH phase to a large BH phase.As the temperature decreases or the charge increases,Δrhdecreases and eventually turns to zero at the critical point,where the small and large BH phases can not be distinguished anymore.

    Fig.4.The sudden change Δrh when the black hole comes up a first-order phase transition.Left Panel(a):As the temperature decreases,the Δrh decreases to zero.Right Panel(b):As the charge increases,the Δrh decreases to zero.The QED parameter and radius of the cavity are set to be a=0.08 and rB=3.

    Fig.5.Phase diagram for a >0 by taking the charge q as an analogy of the pressure. Left Panel(a):First-order coexistence curve (red solid) and zeroth-order phase transition curve (blue solid) in the q-T plane.The black curve separates the regions with and without black holes.The coexistence curve separates black holes into small and large BH phases,and ends at the critical point of the small/large BH phase transition. Right Panel(b):The phase structure in the q-rh plane.We choose the parameter a=0.08 and the radius of cavity rB=3.

    The phase structures of Euler–Heisenberg black hole in a cavity are shown in figure 5.The coexistence curve [64] is represented by the red solid curve in figure 5(a).The minimum temperature curve is marked by the black solid curve,which divides regions with or without BH phases.The coexistence curve and the minimum temperature curve are connected by the zeroth-order phase transition curve (blue solid).The small BH region is between the coexistence curve and zeroth-order phase transition curve.The large BH regions are upon the coexistence curve and below the zeroth-order phase transition curve.We show the phase structure in the qrhplane in figure 5(b).The red curve corresponds to the coexistence curve of the small BH phase and large BH phase,while the blue ones denote the zeroth-order phase transition curve.The small and large BH phases are on the left and right sides of the coexistence curve,respectively.The phase coexistence region is bounded by the coexistence curve,which includes the superheated small BH branch,supercooled large BH branch and unstable BH branches.It is noteworthy that below the right branch of zeroth-order phase transition curve,another region for large BH phases exists.

    3.2.a <0 case

    The Euler–Heisenberg black hole in a cavity undergoes a single phase transition for a negative QED parameter,which is different from the case for a small positive QED parameter.For the condition that a=-1 and rB=3,we can obtain qc=0.6909,and the free energy F versus the temperature T for different values of charge q can be shown in figure 6.

    The swallow tail behavior occurs when the charge q <qcin figure 6(a),indicating that there is a first-order phase transition.With the increase of charge,the swallow tail shrinks to a point at the critical charge q=qc.After that,there will be no phase transition for a larger q,where a smooth free energy function F(T)appears.Furthermore,we investigate the free energy versus the temperature for the charge q=0.3 in figure 6(b).By judging the heat capacity of the black hole,one can tell that the red dashed curve represents the unstable BH phase,while the others are stable ones.As discussed before,the phase transition happens at the intersection of the blue and green solid curves.Since the stable system’s free energy is always lower compared to the unstable one,we can conclude that the small/large BH phase transition happens between a small BH phase and a large BH phase.

    Fig.6. Left Panel(a):The free energy versus the temperature for q=0.3 (blue solid curve),q=0.5 (purple solid curve),q=0.6909(green dashed curve)and q=0.8(orange solid curve).Right Panel(b):The free energy versus the temperature for q <qc,where we choose q=0.3.The small/large BH phase transition occurs at the intersection of the blue and green solid curve,which represents the stable small BH phase and the stable large BH phase,respectively.The red dashed curve is the unstable intermediate BH branch.We set a=-1 and rB=3 here.

    Fig.7.Phase diagram for a <0 by taking charge q as an analogy of pressure. Left Panel(a):First-order coexistence curve in the q-T plane,which separates the black hole into small BH phase and large BH phase.The coexistence curve ends at the critical point of the small/large BH phase transition.Right Panel(b):The phase structure in the q-rh plane.With the increase of radius,small and large BHs coexist in the region bounded by the red curve and the axis.We choose the parameter a=-1 and the cavity radius rB=3.

    The phase structures of Euler–Heisenberg black holes in a cavity for a <0 can be shown in figure 7.The coexistence curve approximately starts at q=0,as shown in figure 7(a),and ends at the critical point where the first-order phase transition disappears.As a coexistence curve,it divides the black hole into small and large BH phases in the q-T phase diagram.In figure 7(b),the coexistence region of the small and large BHs occupies below the coexistence curve.

    There are some differences between the cases that a <0 and 0<a<amax.The swallow tail behavior in the F-T diagram shows that the black hole admits a first-order phase transition for a <0,such as in figure 6(b).This phase transition exists even for a small charge.However,phase transitions do not occur at first for 0<a<amaxas the charge increases.In addition,the difference in critical points for a <0 and a >0 may be analyzed from the perspective of topology [65,66],which is considered as a change of topological charge.

    4.Ruppeiner geometry and microstructure

    After analyzing the effect of the QED parameter for phase transitions,we can study the microstructure of the Euler–Heisenberg black hole in a cavity by using the Ruppeiner geometry [24].The Ruppeiner metric gμνis defined as the negative second derivative of the entropy with respect to the independent variables xμin the thermodynamic system [24],and is given by

    where the xμ=(U,Y) and U is the internal energy and Y is another extensive variable of the system.The Ricci scalar of the Ruppeiner metric is called the Ruppeiner invariant R,whose value can shed some light on the microstructure of black holes.In figure 2,we judge the thermal stability by Cqwhich is equivalent to the slope of T(rh),and the Ruppeiner invariant diverges at the divergent point of the specific heat capacity [33].Hence,we choose the thermodynamic coordinates xμ=(U,Φ) where we define U as a new conjugate potential of E.The internal energy is the Legendre transformation of E,yielding

    By using equation(18),we obtain the line element of the Ruppeiner geometry with fixed a and rBusing the same method in [67,68]

    By introducing the total derivative of U,

    the line element now depends on T,Φ,S and q.

    Combining the following relations

    The associated Ruppeiner invariant R can be reduced to

    4.1.a >0 case

    In a certain range of the charge q and positive QED parameter a,the system undergoes a reentrant phase transition.As above,we take a=0.08 and rB=3 as an example.It has concluded in section 3.1 that only the large BH is the globally thermal stable phase for 0 <q <q1,while a reentrant phase transition occurs for q1<q <q2,and then a single first-order phase transition replaces in q2<q <qc2,followed by the case in which the small BH and large BH can not be distinguished for q >qc2.

    We can choose the same values of the charge q as in section 3.1 to show the details of the microstructure of the Euler–Heisenberg black hole in a cavity.According to the phase structure of figure 3,the charges q=0.279,0.285,0.3 and 0.8 are chosen to present the behavior of the Ruppeiner invariant R in figure 8.

    48. Kind heaven help us now!: Gretel prays for heavenly assistance. Once again, these religious references were added by the Grimms.Return to place in story.

    In the first three panels of figure 8,all of them show that R is always negative for the stable large BH (green solid curve) and the unstable intermediate BH (orange dashed curve),indicating attraction between these BH molecules.But R can be negative or positive for the stable small BH (blue solid curve) and the unstable smaller BH (red dashed curve)depending on the temperature,so the interaction between the BH molecules can be repulsive or attractive.

    For qc1<q <q1,the globally stable phase is the large BH with the minimum free energy,so no phase transition happens.In this case,shown in figure 8(a),R diverges to infinity at the minimum temperature and is always negative for the large BH,demonstrating that the interaction between BH molecules is always attractive for the globally stable phase.

    A reentrant phase transition exists for q1<q <q2,which consists of a first-order phase transition and a zeroth-order phase transition.We choose q=0.285 for this case,and the first-order phase transition between the stable small BH phase and the stable large BH phase occurs for T=0.0776,corresponding to the intersection of the small BH and the large BH in figure 3(b).The Ruppeiner invariant R is positive for the small BH,but it is negative for the large BH when the first-order phase transition happens,as shown in figure 8(b).Hence,the type of interaction changes when the system undergoes this phase transition.|R| of the small BH is larger than that of the large BH when the first-order phase transition occurs,so the interaction of the small BH is more intensive than that of the large BH.At the point of the zeroth-order phase transition between the large BH and small BH with T=0.0734,corresponding to the black dashed curve in figure 3(b),R=+∞for the small BH,while it is negative and finite for the large BH,so it displays that the interaction between BH molecules mutates from weak attraction to intensive repulsion as the temperature increases.

    There is a single first-order phase transition for q2<q <qc2,and the microstructure of this case is shown in figure 8(c).The first-order phase transition occurs at T=0.0770 for q=0.3,corresponding to the intersection of the small BH and large BH in figure 3(c).R is positive for the small BH but negative for the large BH,and |R| of the small BH is larger than that of the large BH,which indicates that the type of the interaction mutates and the attraction of the small BH is more intensive than the repulsion of the large BH at the first-order phase transition point.

    Fig.8.The Ruppeiner invariant R versus T with a fixed a=0.08.Upper Left Panel(a):For qc1 <q <q1,we set q=0.279.Upper Right Panel(b):q=0.285,in the range of q1 <q <q2. Lower Left Panel(c):q=0.3 is chosen as an example for q2 <q <qc2.In these three panels,there are four black hole solutions,and the red/orange dashed curves represent unstable smaller BH and intermediate BH,while the blue/green solid curves represent stable small/large BH,and the black dashed curve corresponds to the temperature of the first-order phase transition. Lower Right Panel(d):q >qc2,where we set q=0.8,there is only one stable solution (blue solid curve).

    No phase transition occurs for q >qc2.We could set q=0.8 to show the microstructure of the stable BH phase in figure 8(d).It shows R=+∞at T=0,and R →0-with T →+∞,which means that the attraction between BH molecules becomes weaker with the increasing temperature.The behavior of the Ruppeiner invariant is like the quantum anyon gas [45].

    4.2.a <0 case

    One can explore the microstructure of the black hole systems for a <0 by the same method in the last subsection.For the case in section 3.2,a first-order phase transition exists for a small charge q.Choosing a=-1 and rB=3 as an example,the system undergoes a first-order phase transition for 0 <q <qc,while no phase transition occurs for q >qc.For a visual representation,we set q=0.3 and 0.8 to present the behavior of the Ruppeiner invariant R in figure 9.The blue and green solid curves represent the stable small BH phase and the stable large BH phase,respectively,and the red dashed curve indicates the behavior of the unstable intermediate BH phase.

    For the q <qccase shown in figure 9(a),R is always negative for the large BH and intermediate BH,so attraction dominates between these BH molecules.Additionally,R →0-with T →+∞for the large BH indicates the attraction becomes weaker with temperature increasing.However,R can be positive for the small BH with low temperature,so repulsion dominates between these BH molecules.Considering the globally stable phase,one has R →+∞with T →0 for the extremal black hole,and R is always finite for T ≠0.A first-order phase transition between the small BH and large BH happens at Tp=0.0759,corresponding to the intersection of the small BH and the large BH in figure 6(b).At T=Tp,R is negative for both the small BH and the large BH,but |R|of the small BH is larger than that of the large BH.Thus,the attraction of the small BH is more intensive than that of the large BH.

    For q >qc,the microstructure of the black hole systems is displayed in figure 9(b).R →+∞with T →0+,and R →0-with T →+∞.The interaction of the BH molecules is always attractive,like the Bose gas.Hence,in the case of no phase transition,the microstructure of the Euler–Heisenberg black hole in a cavity for a <0 is different from that of RN black hole in a cavity studied in [37],in which R is always finite in a globally stable phase.

    Fig.9.The Ruppeiner invariant R versus temperature T with fixed a=-1.Left Panel(a):For q <qc,there are three black hole solutions,and we choose q=0.3 for visualization.The red dashed curve represents the unstable intermediate BH,while the blue/green solid curve represents the stable small/large BH,and the black dashed curve corresponds to the temperature of the first-order phase transition. Right Panel(b):q >qc,where q=0.8 to show the details,there is only one solution and it is stable (blue solid curve).

    Fig.10.The behavior of the Ruppeiner invariant R along the first-order coexistence curve.Left Panel(a):a=0.08,R depends on q along the first-order coexistence curve in figure 5(a).Right Panel(b):a=-1,the first-order coexistence curve is in figure 7(a).The green curve represents the small BH,while the red curve represents the large BH.

    We can compare the microstructure of the small BH with that of the large BH when the first-order phase transition occurs between a >0 and a <0.In the case of a=0.08,as the first-order phase transition occurs,the small BH intermolecular interaction is repulsive.However,in the case of a=-1,attraction dominates the small BH when the firstorder phase transition occurs.According to these differences,the Ruppeiner invariant along the first-order coexistence curve is worth studying.

    As mentioned before,a reentrant phase transition exists for q1<q <q2,and the special range is shown in the inset of figure 10(a).It is obvious that repulsion dominates in this specific range as the reentrant phase transition occurs,and R diverges positively with q=q1.As the charge increases,only the first-order transition exists,while the dominant interaction becomes attractive from repulsive.A change from low to high correlation(the magnitude of R)among components is one of the signs of a phase transition(or vice versa)[46].For a given charge,when decreasing the temperature,the phase transits from the large BH to the small BH in the most range of q,so the |R| becomes larger and the correlation increases.Additionally,in the range of q1<q <q2,the difference of|R| between the small BH and the large BH is larger than other ranges of charge,as shown in figure 10(a).Hence,we deduce that the larger difference in the correlation leads to a reentrant phase transition,which is the same as the conclusion in [46].

    Only a first-order phase transition occurs in the case of a <0.The divergence points of the Ruppeiner invariant are three instead of two for the small BH along the first-order coexistence curve.The divergence points correspond to q=0,qdand qc,where qd=0.0145.In the range of qd<q <qc,the microscopic behavior of this case is similar to that in the case of a >0.The inset of figure 10(b) shows that the interaction between the small BH molecules is intensively attractive with 0 <q <qd.

    Therefore,for the reentrant phase transition,only repulsion exists between the BH molecules,but both repulsion and attraction can dominate with the presence of the usual small/large BH phase transition.The distinction between the cases of negative and positive QED parameters is the extra range dominated by the attraction in the case of a <0,as illustrated in the inset of figure 10(b).So we deduce that the attraction may prevent the BH from undergoing an extra zeroth-order phase transition when the temperature decreases,and the negative QED parameter may help the formation of the extra attraction region to obstruct the reentrant phase transition.

    5.Conclusion

    In this paper,we focus on the thermodynamic properties of the Euler–Heisenberg black hole with the QED correction in a cavity.The phase transitions were investigated for different charges and QED parameters,and the Ruppeiner geometry was analyzed to explore the microscopic properties.

    We derived the first law of black hole thermodynamics via the thermal energy,and the free energy of black holes was given naturally.According to the free energy,the phase transitions of Euler–Heisenberg black holes in a cavity were analyzed,which was related to different charges q and QED parameters a.Based on this analysis,we found that there is a maximum value of a allowing for phase transitions,which is determined by the cavity radius.Then we investigated the temperature and critical points of the black hole systems for a >0 and a <0,respectively.The phase transitions were investigated in detail.For a small positive QED parameter,there could be two critical points and four BH branches,including two stable branches and two unstable branches.As the charge increases,the black hole system for a fixed QED parameter undergoes a process from a stable BH phase to a reentrant phase transition,to a first-order phase transition,and thereafter no phase transition occurs.For a negative QED parameter,there is a small/large BH phase transition,and only one critical point exists.This first-order phase transition occurs for a small charge,where there are two stable BH branches and one unstable intermediate BH branch.There is no phase transition for the charge above the critical point.

    Furthermore,we explore the microstructure of this black hole by the Ruppeiner approach [24].As a preliminary step,taking (U,Φ) as the coordinates in the thermodynamic phase space and defining a conjugate potential as the internal energy,we calculate the Ricci scalar of the Ruppeiner geometry for the black hole.Along the first-order coexistence curve,for the small BH,only repulsion dominates between the BH molecules when the reentrant phase transition occurs,while the dominant interaction can be attractive or repulsive when the usual small/large BH phase transition exists.Then an attraction region in the small BH phase is missing for the positive QED parameter.It infers that the attraction may affect the reentrant phase transition,and the negative QED parameter makes the extra attraction region for the small BH phase,so the zeroth-order phase transition is obstructed by the attraction region.

    Following the investigation,it is inspiring to find that the phase structure of the Euler–Heisenberg black hole in a cavity is similar to the black hole with the AdS boundary[51].Then,one can note that different thermodynamic coordinates may result in various microstructures,such as the investigation of the Kerr-Newman black hole in[39].So we will focus on the correlation between the phase structures and various microstructures in our future works.

    Acknowledgments

    We are grateful to Yihe Cao,Yuchen Huang,Ningchen Bai,Aoyun He,Xin Jiang and Peng Wang for useful discussions and valuable comments.This work is supported by NSFC(Grant No.12047573 and 12275184).

    亚洲精品色激情综合| 国语对白做爰xxxⅹ性视频网站| 精品人妻熟女av久视频| 男女那种视频在线观看| 黄色配什么色好看| 中文欧美无线码| 国产成人免费观看mmmm| 欧美激情在线99| 亚洲欧美中文字幕日韩二区| 亚洲国产精品国产精品| 久久久久久久久久人人人人人人| 少妇丰满av| 久久久精品欧美日韩精品| 亚洲精品国产av成人精品| 午夜精品国产一区二区电影 | 99视频精品全部免费 在线| 深夜a级毛片| 一级毛片电影观看| 水蜜桃什么品种好| 亚洲丝袜综合中文字幕| 国产有黄有色有爽视频| 国产一区二区在线观看日韩| 免费观看无遮挡的男女| 高清毛片免费看| 国产毛片在线视频| 亚洲欧美日韩无卡精品| 免费大片18禁| 国语对白做爰xxxⅹ性视频网站| 99热6这里只有精品| 国产成人a区在线观看| 黄色日韩在线| 亚洲综合色惰| 最近手机中文字幕大全| 亚洲美女搞黄在线观看| av天堂中文字幕网| 亚洲伊人久久精品综合| 精品久久久精品久久久| 一级爰片在线观看| 我的老师免费观看完整版| 亚洲av.av天堂| 精品久久久久久久久av| 嫩草影院入口| 国产一级毛片在线| 久久久久久久亚洲中文字幕| 精品酒店卫生间| 精品人妻视频免费看| 亚洲精品国产色婷婷电影| 久久久久网色| 午夜视频国产福利| 亚洲精品久久久久久婷婷小说| 大片电影免费在线观看免费| 少妇 在线观看| 美女被艹到高潮喷水动态| 国产免费福利视频在线观看| 夫妻午夜视频| 在线免费观看不下载黄p国产| 联通29元200g的流量卡| 男人舔奶头视频| 欧美少妇被猛烈插入视频| 亚洲最大成人手机在线| 久久久久久久亚洲中文字幕| 日韩人妻高清精品专区| 成人漫画全彩无遮挡| 少妇人妻精品综合一区二区| a级一级毛片免费在线观看| 欧美三级亚洲精品| 久久久国产一区二区| 毛片女人毛片| 日韩一区二区视频免费看| 永久免费av网站大全| 啦啦啦在线观看免费高清www| 国语对白做爰xxxⅹ性视频网站| 高清在线视频一区二区三区| 亚洲国产高清在线一区二区三| 亚洲国产精品专区欧美| 99久国产av精品国产电影| 又黄又爽又刺激的免费视频.| 日本免费在线观看一区| 国产男人的电影天堂91| 国产 精品1| 99久久精品一区二区三区| 交换朋友夫妻互换小说| 中国国产av一级| 最近的中文字幕免费完整| 午夜视频国产福利| 国产真实伦视频高清在线观看| 偷拍熟女少妇极品色| 啦啦啦啦在线视频资源| 国产精品.久久久| 久久精品国产自在天天线| 国产av国产精品国产| 亚洲国产成人一精品久久久| 全区人妻精品视频| 日韩av免费高清视频| 亚洲最大成人手机在线| 熟女av电影| 亚洲真实伦在线观看| 一级二级三级毛片免费看| 性色av一级| 欧美激情国产日韩精品一区| 亚洲天堂av无毛| 亚洲欧美日韩卡通动漫| 男女那种视频在线观看| 午夜日本视频在线| 少妇裸体淫交视频免费看高清| 成人鲁丝片一二三区免费| 熟女av电影| 精品人妻视频免费看| 欧美激情在线99| av线在线观看网站| 边亲边吃奶的免费视频| 成人特级av手机在线观看| 少妇熟女欧美另类| 国产毛片a区久久久久| 在线观看一区二区三区| 18禁裸乳无遮挡免费网站照片| 一级爰片在线观看| 国产亚洲5aaaaa淫片| 深爱激情五月婷婷| 亚洲av在线观看美女高潮| 亚洲成人av在线免费| 我的女老师完整版在线观看| 丝袜脚勾引网站| 80岁老熟妇乱子伦牲交| 免费观看a级毛片全部| 91在线精品国自产拍蜜月| 久久久精品免费免费高清| 亚洲精品乱码久久久v下载方式| av天堂中文字幕网| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 少妇的逼水好多| 午夜福利在线在线| 亚洲婷婷狠狠爱综合网| 涩涩av久久男人的天堂| 边亲边吃奶的免费视频| 啦啦啦中文免费视频观看日本| 女的被弄到高潮叫床怎么办| 夜夜看夜夜爽夜夜摸| 亚洲精品色激情综合| 国产美女午夜福利| www.av在线官网国产| 日本wwww免费看| 国产免费一区二区三区四区乱码| 免费人成在线观看视频色| 在线看a的网站| 综合色av麻豆| 99热这里只有精品一区| 交换朋友夫妻互换小说| 免费观看无遮挡的男女| 又黄又爽又刺激的免费视频.| 中文天堂在线官网| 91在线精品国自产拍蜜月| 伊人久久国产一区二区| 国产乱人视频| 精品少妇久久久久久888优播| 久久久久久久国产电影| 九九在线视频观看精品| 我要看日韩黄色一级片| 亚洲不卡免费看| 亚洲电影在线观看av| 精品久久久久久久久av| 欧美高清性xxxxhd video| 久久久久久久久久人人人人人人| 国产视频首页在线观看| 高清毛片免费看| 亚洲av日韩在线播放| 男女国产视频网站| 中文字幕免费在线视频6| 美女脱内裤让男人舔精品视频| 一本色道久久久久久精品综合| 乱系列少妇在线播放| 亚洲性久久影院| 日韩精品有码人妻一区| 最近2019中文字幕mv第一页| 黄色日韩在线| 青春草亚洲视频在线观看| 香蕉精品网在线| 最近中文字幕2019免费版| 国产精品三级大全| 亚洲精品一区蜜桃| 免费人妻精品一区二区三区视频| 午夜91福利影院| 美女中出高潮动态图| 在线观看免费日韩欧美大片| 中文字幕精品免费在线观看视频| 亚洲国产精品999| 一区二区av电影网| 亚洲精品第二区| 国产一区二区激情短视频 | 欧美精品一区二区免费开放| 人妻人人澡人人爽人人| 蜜桃在线观看..| 亚洲av成人精品一二三区| 国产老妇伦熟女老妇高清| 欧美激情 高清一区二区三区| 国产欧美日韩综合在线一区二区| 亚洲精品第二区| 亚洲av日韩精品久久久久久密 | av一本久久久久| 悠悠久久av| 男女之事视频高清在线观看 | 国产1区2区3区精品| 波多野结衣av一区二区av| 国产黄色免费在线视频| av电影中文网址| 考比视频在线观看| 精品国产一区二区久久| 久热爱精品视频在线9| 亚洲av日韩精品久久久久久密 | 久久影院123| 久久久久久人妻| 在线看a的网站| 天美传媒精品一区二区| h视频一区二区三区| 观看av在线不卡| 久久精品亚洲熟妇少妇任你| 亚洲久久久国产精品| 日韩大片免费观看网站| 飞空精品影院首页| 婷婷色麻豆天堂久久| 午夜福利视频精品| 欧美黄色片欧美黄色片| 2018国产大陆天天弄谢| av线在线观看网站| 国产一区二区三区av在线| 最近手机中文字幕大全| 女人久久www免费人成看片| 午夜久久久在线观看| 激情视频va一区二区三区| 精品国产露脸久久av麻豆| 亚洲av电影在线观看一区二区三区| www.自偷自拍.com| 制服诱惑二区| 97精品久久久久久久久久精品| 性高湖久久久久久久久免费观看| 搡老乐熟女国产| 国产乱人偷精品视频| 男女边摸边吃奶| 少妇猛男粗大的猛烈进出视频| 久热这里只有精品99| 久久精品久久久久久噜噜老黄| 黄色视频在线播放观看不卡| av.在线天堂| 欧美久久黑人一区二区| 免费少妇av软件| 亚洲色图 男人天堂 中文字幕| 夜夜骑夜夜射夜夜干| 亚洲,欧美,日韩| av又黄又爽大尺度在线免费看| 老汉色∧v一级毛片| 亚洲第一区二区三区不卡| 91老司机精品| 午夜免费观看性视频| 日本wwww免费看| 亚洲久久久国产精品| 大香蕉久久成人网| 777久久人妻少妇嫩草av网站| 人人妻人人澡人人爽人人夜夜| 欧美国产精品一级二级三级| 搡老岳熟女国产| 男人操女人黄网站| 国产日韩欧美视频二区| 成人毛片60女人毛片免费| 啦啦啦视频在线资源免费观看| 80岁老熟妇乱子伦牲交| 又粗又硬又长又爽又黄的视频| 悠悠久久av| 久久久久久久久久久久大奶| 日本欧美国产在线视频| 色网站视频免费| 观看美女的网站| 精品国产一区二区三区四区第35| 色婷婷av一区二区三区视频| 狠狠精品人妻久久久久久综合| 男女之事视频高清在线观看 | 天天躁日日躁夜夜躁夜夜| 高清视频免费观看一区二区| 久久久久精品国产欧美久久久 | 王馨瑶露胸无遮挡在线观看| 岛国毛片在线播放| 久久影院123| 国产日韩欧美在线精品| 国产97色在线日韩免费| 午夜精品国产一区二区电影| 青草久久国产| 少妇人妻精品综合一区二区| 美女大奶头黄色视频| 丰满少妇做爰视频| 国产精品一区二区在线观看99| 精品第一国产精品| 免费看不卡的av| 亚洲精品乱久久久久久| 人人妻,人人澡人人爽秒播 | 国产精品人妻久久久影院| 在线亚洲精品国产二区图片欧美| 国产一区二区三区av在线| 97人妻天天添夜夜摸| 国产av码专区亚洲av| 日本爱情动作片www.在线观看| 日韩一区二区视频免费看| 成人国产av品久久久| 男女之事视频高清在线观看 | 欧美 亚洲 国产 日韩一| 精品酒店卫生间| 国产精品女同一区二区软件| 亚洲精品一二三| 午夜激情久久久久久久| av在线app专区| 婷婷色综合大香蕉| 亚洲成人一二三区av| 日韩av在线免费看完整版不卡| 久久久久人妻精品一区果冻| 亚洲精品久久成人aⅴ小说| 日本vs欧美在线观看视频| 少妇 在线观看| 满18在线观看网站| 亚洲国产欧美网| 欧美老熟妇乱子伦牲交| 国产高清国产精品国产三级| 少妇被粗大的猛进出69影院| svipshipincom国产片| 丝袜喷水一区| 亚洲欧美成人综合另类久久久| 另类精品久久| 在线观看免费日韩欧美大片| 黄网站色视频无遮挡免费观看| 国产精品久久久久久人妻精品电影 | 日韩中文字幕欧美一区二区 | 大片电影免费在线观看免费| 欧美在线黄色| 大码成人一级视频| 高清在线视频一区二区三区| 亚洲av成人精品一二三区| 国产成人欧美| 日本色播在线视频| 欧美成人精品欧美一级黄| 日韩大码丰满熟妇| 国产精品免费视频内射| 亚洲国产精品一区二区三区在线| 我要看黄色一级片免费的| 午夜福利在线免费观看网站| 丝袜美足系列| 色网站视频免费| 精品少妇内射三级| 精品人妻熟女毛片av久久网站| 精品一区在线观看国产| 成人午夜精彩视频在线观看| 欧美中文综合在线视频| 精品一区二区三卡| 色综合欧美亚洲国产小说| 国产亚洲av片在线观看秒播厂| 久久精品久久久久久噜噜老黄| 丝袜在线中文字幕| 看非洲黑人一级黄片| 亚洲欧美日韩另类电影网站| 十八禁人妻一区二区| 国产精品免费视频内射| 国产黄色免费在线视频| 麻豆精品久久久久久蜜桃| 一级片免费观看大全| 热99国产精品久久久久久7| 天天添夜夜摸| 日韩免费高清中文字幕av| 老司机影院毛片| 美女福利国产在线| 亚洲人成77777在线视频| 丁香六月欧美| 波野结衣二区三区在线| 亚洲国产最新在线播放| 叶爱在线成人免费视频播放| 高清欧美精品videossex| 国产av精品麻豆| 亚洲国产欧美网| 亚洲av日韩精品久久久久久密 | 一本久久精品| av.在线天堂| 高清欧美精品videossex| 你懂的网址亚洲精品在线观看| 日韩欧美精品免费久久| 亚洲人成77777在线视频| 男男h啪啪无遮挡| 久久久国产精品麻豆| 午夜福利视频精品| 国产免费现黄频在线看| 亚洲欧美一区二区三区久久| 在线天堂最新版资源| 美国免费a级毛片| 久久国产精品大桥未久av| netflix在线观看网站| 午夜福利免费观看在线| av国产久精品久网站免费入址| 午夜福利免费观看在线| 国语对白做爰xxxⅹ性视频网站| 美女脱内裤让男人舔精品视频| 一级毛片我不卡| 人体艺术视频欧美日本| 精品久久久精品久久久| 男女免费视频国产| 丝袜脚勾引网站| 国产精品 国内视频| 丝袜脚勾引网站| 久久久亚洲精品成人影院| 中文乱码字字幕精品一区二区三区| 满18在线观看网站| 久久久久国产精品人妻一区二区| av在线观看视频网站免费| 青春草国产在线视频| 国产极品天堂在线| 亚洲成人一二三区av| 波多野结衣一区麻豆| 亚洲av综合色区一区| 观看av在线不卡| 狠狠精品人妻久久久久久综合| 国产极品粉嫩免费观看在线| 人人妻,人人澡人人爽秒播 | 国产视频首页在线观看| 99re6热这里在线精品视频| av国产精品久久久久影院| 五月开心婷婷网| 大片免费播放器 马上看| 久久精品人人爽人人爽视色| 婷婷色综合www| 在线看a的网站| 欧美日韩av久久| 国产精品 国内视频| 亚洲精品成人av观看孕妇| 一区二区三区精品91| 国产精品久久久久久精品古装| 制服诱惑二区| 亚洲av成人不卡在线观看播放网 | 久久精品亚洲av国产电影网| 女的被弄到高潮叫床怎么办| 亚洲伊人久久精品综合| 精品一区二区三卡| 久久影院123| 精品人妻熟女毛片av久久网站| 久久 成人 亚洲| 免费高清在线观看日韩| 五月天丁香电影| 男女边摸边吃奶| 国产亚洲一区二区精品| 99香蕉大伊视频| 亚洲伊人色综图| 国产一区亚洲一区在线观看| 国产午夜精品一二区理论片| av天堂久久9| 久久精品国产亚洲av涩爱| 亚洲色图综合在线观看| 国产有黄有色有爽视频| 777久久人妻少妇嫩草av网站| 国产精品国产三级专区第一集| 国产一级毛片在线| 日日撸夜夜添| 国产淫语在线视频| 久久ye,这里只有精品| 久久精品国产亚洲av涩爱| 中文字幕高清在线视频| 免费女性裸体啪啪无遮挡网站| a 毛片基地| 国产黄色免费在线视频| 在线亚洲精品国产二区图片欧美| videos熟女内射| 在线观看免费午夜福利视频| 黑人巨大精品欧美一区二区蜜桃| 在线观看三级黄色| 欧美激情 高清一区二区三区| av在线播放精品| 欧美最新免费一区二区三区| 老司机影院成人| 久久天躁狠狠躁夜夜2o2o | 两个人免费观看高清视频| av天堂久久9| 黄色一级大片看看| av不卡在线播放| 视频在线观看一区二区三区| 纯流量卡能插随身wifi吗| 99久久精品国产亚洲精品| 精品国产一区二区三区四区第35| 久久热在线av| 久久国产精品大桥未久av| 成年人午夜在线观看视频| 欧美日韩成人在线一区二区| 日韩制服丝袜自拍偷拍| 大码成人一级视频| 免费在线观看完整版高清| 毛片一级片免费看久久久久| 成人毛片60女人毛片免费| 建设人人有责人人尽责人人享有的| 操出白浆在线播放| 久久国产精品男人的天堂亚洲| www.熟女人妻精品国产| 国产一区二区激情短视频 | 国精品久久久久久国模美| 欧美日韩综合久久久久久| 国产精品久久久久久精品电影小说| 人妻一区二区av| 在线观看免费高清a一片| 免费女性裸体啪啪无遮挡网站| 少妇人妻精品综合一区二区| 国产欧美日韩综合在线一区二区| 日本欧美视频一区| 亚洲av电影在线观看一区二区三区| 日韩熟女老妇一区二区性免费视频| 亚洲欧美一区二区三区久久| 国产女主播在线喷水免费视频网站| 激情视频va一区二区三区| 国产一区二区在线观看av| 一区福利在线观看| 国精品久久久久久国模美| 国产伦人伦偷精品视频| 天堂俺去俺来也www色官网| 精品一区二区三区四区五区乱码 | 尾随美女入室| 天天躁日日躁夜夜躁夜夜| 色精品久久人妻99蜜桃| 亚洲国产成人一精品久久久| 亚洲国产欧美在线一区| 99国产精品免费福利视频| 欧美乱码精品一区二区三区| 色精品久久人妻99蜜桃| 伊人久久大香线蕉亚洲五| 高清av免费在线| 男女床上黄色一级片免费看| 日韩免费高清中文字幕av| 汤姆久久久久久久影院中文字幕| 亚洲av在线观看美女高潮| √禁漫天堂资源中文www| 国产在线视频一区二区| 欧美国产精品va在线观看不卡| 国产精品偷伦视频观看了| 高清在线视频一区二区三区| 亚洲精品国产色婷婷电影| 亚洲av日韩精品久久久久久密 | 精品国产超薄肉色丝袜足j| 免费不卡黄色视频| 精品一品国产午夜福利视频| 亚洲精品av麻豆狂野| 成人毛片60女人毛片免费| 国产成人欧美| 又粗又硬又长又爽又黄的视频| 十八禁网站网址无遮挡| 国产精品嫩草影院av在线观看| 国产黄频视频在线观看| 少妇精品久久久久久久| 黑人欧美特级aaaaaa片| av.在线天堂| 啦啦啦啦在线视频资源| 亚洲精品乱久久久久久| 亚洲久久久国产精品| 亚洲四区av| 亚洲熟女精品中文字幕| 国产色婷婷99| 亚洲精品av麻豆狂野| 人成视频在线观看免费观看| 如何舔出高潮| 国产精品免费大片| 国产成人精品福利久久| 国产精品99久久99久久久不卡 | 亚洲国产看品久久| 啦啦啦视频在线资源免费观看| 欧美精品人与动牲交sv欧美| 久久人人97超碰香蕉20202| 熟女少妇亚洲综合色aaa.| 国产成人精品福利久久| 欧美精品亚洲一区二区| 久久这里只有精品19| 国产又色又爽无遮挡免| 一级毛片我不卡| 国产精品成人在线| av在线老鸭窝| 亚洲精品久久午夜乱码| 男女免费视频国产| 纯流量卡能插随身wifi吗| 夫妻午夜视频| 黑丝袜美女国产一区| 亚洲精品成人av观看孕妇| av片东京热男人的天堂| 人人妻人人爽人人添夜夜欢视频| 不卡av一区二区三区| 哪个播放器可以免费观看大片| 男女无遮挡免费网站观看| 久久国产亚洲av麻豆专区| 欧美精品人与动牲交sv欧美| 亚洲图色成人| 亚洲国产精品一区三区| 欧美xxⅹ黑人| 午夜精品国产一区二区电影| www.熟女人妻精品国产| 丝袜美腿诱惑在线| 国产一级毛片在线| 亚洲国产中文字幕在线视频| 人人妻人人澡人人爽人人夜夜| 亚洲欧美激情在线| 亚洲一级一片aⅴ在线观看| 老汉色av国产亚洲站长工具| 久久久久精品性色| 日韩一卡2卡3卡4卡2021年| 久久av网站| 色婷婷久久久亚洲欧美| 成人国产麻豆网| 亚洲成人国产一区在线观看 | 久久久国产欧美日韩av| 久久99热这里只频精品6学生| 中文天堂在线官网| 亚洲伊人色综图| 天堂8中文在线网| 最近最新中文字幕大全免费视频 | 国产在线免费精品| 秋霞在线观看毛片| 精品午夜福利在线看| 国产精品.久久久| 青春草国产在线视频| 啦啦啦在线观看免费高清www| 久久精品人人爽人人爽视色| 国产亚洲一区二区精品|