• <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).

    久久精品91蜜桃| 身体一侧抽搐| 桃色一区二区三区在线观看| 色综合站精品国产| 中文字幕人成人乱码亚洲影| 免费无遮挡裸体视频| 91老司机精品| 国产欧美日韩一区二区精品| 免费看美女性在线毛片视频| 国产精品久久久久久精品电影| 黄色成人免费大全| x7x7x7水蜜桃| 日本a在线网址| 欧美黑人欧美精品刺激| 国产视频一区二区在线看| 国产高清videossex| 狂野欧美白嫩少妇大欣赏| 国模一区二区三区四区视频 | 色吧在线观看| www.自偷自拍.com| 一二三四在线观看免费中文在| 99热精品在线国产| 欧美丝袜亚洲另类 | 色噜噜av男人的天堂激情| 亚洲中文字幕一区二区三区有码在线看 | cao死你这个sao货| 99久久精品热视频| 日本五十路高清| 日韩欧美在线乱码| 久久午夜综合久久蜜桃| 久9热在线精品视频| 亚洲欧美日韩卡通动漫| 草草在线视频免费看| 精品国产超薄肉色丝袜足j| 日韩欧美在线乱码| 日韩欧美国产在线观看| 1024手机看黄色片| 51午夜福利影视在线观看| 国产精品日韩av在线免费观看| 女人高潮潮喷娇喘18禁视频| 亚洲一区二区三区不卡视频| 成人18禁在线播放| 亚洲一区高清亚洲精品| 男人的好看免费观看在线视频| 免费观看人在逋| 国产精品久久久久久精品电影| 亚洲熟女毛片儿| 两个人看的免费小视频| 美女高潮的动态| 两个人视频免费观看高清| 99精品欧美一区二区三区四区| 一区福利在线观看| 在线视频色国产色| 国产综合懂色| 精品久久久久久久末码| 久久精品综合一区二区三区| 色吧在线观看| 日韩欧美精品v在线| 婷婷精品国产亚洲av在线| 他把我摸到了高潮在线观看| 久久久国产精品麻豆| 欧美成人性av电影在线观看| 毛片女人毛片| 午夜激情福利司机影院| 精品久久久久久久久久久久久| 亚洲精品乱码久久久v下载方式 | 18禁裸乳无遮挡免费网站照片| 国产高清有码在线观看视频| 99riav亚洲国产免费| 麻豆av在线久日| 又大又爽又粗| 久久九九热精品免费| 久久久色成人| 亚洲五月天丁香| 免费人成视频x8x8入口观看| 国产精品自产拍在线观看55亚洲| 男女做爰动态图高潮gif福利片| 国产三级在线视频| 男女午夜视频在线观看| 欧美日本亚洲视频在线播放| 99国产极品粉嫩在线观看| 97碰自拍视频| 悠悠久久av| 欧美激情久久久久久爽电影| 黑人操中国人逼视频| 日日夜夜操网爽| 亚洲人成网站高清观看| 久久国产精品人妻蜜桃| 国产aⅴ精品一区二区三区波| 亚洲天堂国产精品一区在线| 久久久久精品国产欧美久久久| 黑人欧美特级aaaaaa片| 色综合欧美亚洲国产小说| 亚洲国产色片| 香蕉国产在线看| 成年人黄色毛片网站| 国产美女午夜福利| 午夜激情欧美在线| 亚洲色图 男人天堂 中文字幕| a级毛片在线看网站| 欧美成人一区二区免费高清观看 | 欧美国产日韩亚洲一区| 国产综合懂色| 精品久久久久久久久久免费视频| 日韩欧美国产在线观看| 久久久久久久久免费视频了| 天堂网av新在线| 国产精品 欧美亚洲| 国产又黄又爽又无遮挡在线| 国产精品一区二区免费欧美| 精品国内亚洲2022精品成人| 中文字幕久久专区| 一级毛片女人18水好多| 9191精品国产免费久久| 欧美一级a爱片免费观看看| 午夜久久久久精精品| 免费观看精品视频网站| 亚洲国产精品久久男人天堂| 国产精品久久电影中文字幕| 午夜成年电影在线免费观看| 精品99又大又爽又粗少妇毛片 | 午夜精品久久久久久毛片777| 91久久精品国产一区二区成人 | 麻豆久久精品国产亚洲av| 白带黄色成豆腐渣| 人人妻人人看人人澡| 国产精品野战在线观看| 可以在线观看毛片的网站| 精品久久久久久久人妻蜜臀av| 啦啦啦韩国在线观看视频| 丁香六月欧美| 国产伦一二天堂av在线观看| 99热精品在线国产| 小蜜桃在线观看免费完整版高清| 精品一区二区三区av网在线观看| 99热这里只有是精品50| 欧美最黄视频在线播放免费| 黄色丝袜av网址大全| 亚洲avbb在线观看| 免费看日本二区| 国产久久久一区二区三区| 亚洲欧美激情综合另类| 丝袜人妻中文字幕| 日日干狠狠操夜夜爽| 欧美黑人巨大hd| 久久中文字幕一级| 两个人看的免费小视频| 国产久久久一区二区三区| 99久久综合精品五月天人人| 特大巨黑吊av在线直播| 伦理电影免费视频| 免费观看精品视频网站| 最新美女视频免费是黄的| 夜夜夜夜夜久久久久| 免费在线观看亚洲国产| 午夜精品久久久久久毛片777| 女警被强在线播放| 国产亚洲精品久久久com| 国产亚洲精品久久久com| 一个人观看的视频www高清免费观看 | 人人妻人人澡欧美一区二区| 动漫黄色视频在线观看| 一进一出抽搐gif免费好疼| 三级国产精品欧美在线观看 | 国产精品免费一区二区三区在线| 国产真人三级小视频在线观看| 啦啦啦免费观看视频1| 91av网一区二区| 五月伊人婷婷丁香| 国产一区在线观看成人免费| 日韩精品青青久久久久久| 国产三级在线视频| 91av网一区二区| 久久婷婷人人爽人人干人人爱| 亚洲熟妇熟女久久| 搡老妇女老女人老熟妇| 高清毛片免费观看视频网站| 国产成年人精品一区二区| 精品福利观看| 熟妇人妻久久中文字幕3abv| 不卡一级毛片| 国产男靠女视频免费网站| 淫秽高清视频在线观看| 在线看三级毛片| 人人妻人人看人人澡| www.自偷自拍.com| 亚洲av日韩精品久久久久久密| 两性夫妻黄色片| 两个人看的免费小视频| www.自偷自拍.com| 天天躁日日操中文字幕| 精品国产乱子伦一区二区三区| 久久久久久人人人人人| 欧美日韩乱码在线| 国产成人啪精品午夜网站| 亚洲熟女毛片儿| 十八禁网站免费在线| 操出白浆在线播放| 熟妇人妻久久中文字幕3abv| 久久中文看片网| 国产亚洲av嫩草精品影院| 99精品在免费线老司机午夜| 成人三级黄色视频| 久久久久九九精品影院| netflix在线观看网站| av国产免费在线观看| 日本 欧美在线| 中文字幕久久专区| 日韩成人在线观看一区二区三区| 嫩草影院入口| 免费一级毛片在线播放高清视频| 精华霜和精华液先用哪个| 嫩草影院入口| 色综合婷婷激情| 91av网一区二区| 国模一区二区三区四区视频 | 国产精品一及| 热99re8久久精品国产| 亚洲专区字幕在线| 1024香蕉在线观看| 在线免费观看不下载黄p国产 | 天天躁日日操中文字幕| 亚洲av中文字字幕乱码综合| 国产一区二区在线观看日韩 | 可以在线观看的亚洲视频| 搡老熟女国产l中国老女人| 在线观看一区二区三区| 成年人黄色毛片网站| 欧美高清成人免费视频www| 日韩 欧美 亚洲 中文字幕| 黄片大片在线免费观看| 亚洲成av人片在线播放无| 巨乳人妻的诱惑在线观看| 在线观看午夜福利视频| 91在线观看av| 宅男免费午夜| 女人被狂操c到高潮| 一个人看的www免费观看视频| 国产成年人精品一区二区| 国产日本99.免费观看| 老熟妇乱子伦视频在线观看| 国产亚洲欧美在线一区二区| 午夜福利成人在线免费观看| 精品电影一区二区在线| 亚洲国产欧美一区二区综合| 可以在线观看毛片的网站| 日本成人三级电影网站| 久久精品aⅴ一区二区三区四区| 俄罗斯特黄特色一大片| 欧美色欧美亚洲另类二区| 国产激情欧美一区二区| 1000部很黄的大片| 亚洲av成人av| 午夜福利视频1000在线观看| 黄色丝袜av网址大全| 亚洲va日本ⅴa欧美va伊人久久| 97碰自拍视频| 最新中文字幕久久久久 | 国产午夜精品久久久久久| 麻豆av在线久日| 99热精品在线国产| 真人一进一出gif抽搐免费| e午夜精品久久久久久久| 国内精品久久久久精免费| 国产黄片美女视频| 每晚都被弄得嗷嗷叫到高潮| 网址你懂的国产日韩在线| 无限看片的www在线观看| 成人精品一区二区免费| 国产aⅴ精品一区二区三区波| 一级毛片高清免费大全| 99精品在免费线老司机午夜| 特级一级黄色大片| 熟女电影av网| 亚洲成人久久性| 又紧又爽又黄一区二区| 亚洲美女黄片视频| 国产三级黄色录像| 免费人成视频x8x8入口观看| 午夜精品久久久久久毛片777| 校园春色视频在线观看| 欧美高清成人免费视频www| 两个人的视频大全免费| 婷婷六月久久综合丁香| 一个人免费在线观看电影 | 国产亚洲av高清不卡| 久久这里只有精品中国| 国产久久久一区二区三区| 九色成人免费人妻av| 日韩高清综合在线| 国产成人aa在线观看| 美女被艹到高潮喷水动态| 亚洲中文字幕一区二区三区有码在线看 | 免费在线观看视频国产中文字幕亚洲| 日韩欧美 国产精品| 99re在线观看精品视频| 老司机午夜十八禁免费视频| 熟女电影av网| 这个男人来自地球电影免费观看| 亚洲午夜精品一区,二区,三区| 久久久久久国产a免费观看| www.熟女人妻精品国产| 国产不卡一卡二| 国产成人一区二区三区免费视频网站| 国产一级毛片七仙女欲春2| 国产成人精品久久二区二区91| 非洲黑人性xxxx精品又粗又长| 美女扒开内裤让男人捅视频| 欧美另类亚洲清纯唯美| 日韩大尺度精品在线看网址| 亚洲无线观看免费| 国产午夜福利久久久久久| 美女免费视频网站| 97人妻精品一区二区三区麻豆| 欧美日本亚洲视频在线播放| 少妇丰满av| 欧美一区二区国产精品久久精品| 精品国产乱码久久久久久男人| 日韩三级视频一区二区三区| a级毛片在线看网站| 久久久国产成人精品二区| 99久久99久久久精品蜜桃| 一区二区三区激情视频| 美女黄网站色视频| 俺也久久电影网| 亚洲天堂国产精品一区在线| av天堂中文字幕网| 男女做爰动态图高潮gif福利片| av视频在线观看入口| 中出人妻视频一区二区| 一个人看视频在线观看www免费 | 精品久久久久久久久久久久久| 亚洲av成人精品一区久久| 日韩欧美国产在线观看| 久久久久国产精品人妻aⅴ院| 国产极品精品免费视频能看的| 久久香蕉精品热| 国产精品av视频在线免费观看| 岛国在线免费视频观看| 97超级碰碰碰精品色视频在线观看| 国产探花在线观看一区二区| 999精品在线视频| 亚洲精品国产精品久久久不卡| 18禁裸乳无遮挡免费网站照片| 欧洲精品卡2卡3卡4卡5卡区| 国产精品一区二区精品视频观看| 91麻豆精品激情在线观看国产| 高清毛片免费观看视频网站| 久久热在线av| 日韩中文字幕欧美一区二区| 美女黄网站色视频| 真实男女啪啪啪动态图| 日韩欧美三级三区| 亚洲欧美日韩高清在线视频| 哪里可以看免费的av片| 免费观看精品视频网站| 18禁国产床啪视频网站| 精品久久久久久久久久免费视频| 97超视频在线观看视频| 亚洲欧美日韩高清在线视频| 久久久水蜜桃国产精品网| 成人鲁丝片一二三区免费| 不卡av一区二区三区| 免费看a级黄色片| 十八禁网站免费在线| 在线免费观看不下载黄p国产 | 日本精品一区二区三区蜜桃| 亚洲av第一区精品v没综合| 成人午夜高清在线视频| 成人一区二区视频在线观看| 色吧在线观看| 欧美3d第一页| 色哟哟哟哟哟哟| 国产高清videossex| 久久国产精品影院| 国产伦人伦偷精品视频| 看免费av毛片| 欧美又色又爽又黄视频| 国产成人精品无人区| 国产精品野战在线观看| 国产一区在线观看成人免费| 亚洲色图 男人天堂 中文字幕| 亚洲国产欧美网| 精品国产美女av久久久久小说| 99在线视频只有这里精品首页| 亚洲欧美日韩高清在线视频| 九九久久精品国产亚洲av麻豆 | 麻豆一二三区av精品| 宅男免费午夜| 色综合欧美亚洲国产小说| 亚洲精品在线观看二区| 国产真实乱freesex| 熟女电影av网| 精品久久久久久久毛片微露脸| 久久久精品大字幕| 国产一区二区在线观看日韩 | 亚洲色图 男人天堂 中文字幕| 国产亚洲av高清不卡| 久久久久久久精品吃奶| 国产亚洲精品一区二区www| 亚洲18禁久久av| 国产久久久一区二区三区| 18禁黄网站禁片免费观看直播| 成人性生交大片免费视频hd| www日本在线高清视频| 麻豆国产97在线/欧美| 黄色丝袜av网址大全| 国产乱人视频| 一个人免费在线观看电影 | 国产一区二区在线观看日韩 | 少妇裸体淫交视频免费看高清| 亚洲欧美日韩高清在线视频| 18禁黄网站禁片免费观看直播| 精品一区二区三区视频在线 | 国产视频一区二区在线看| 亚洲av电影在线进入| av黄色大香蕉| 亚洲精品一区av在线观看| 日本一二三区视频观看| 久久精品亚洲精品国产色婷小说| 很黄的视频免费| 99re在线观看精品视频| svipshipincom国产片| 黄片小视频在线播放| 一级毛片女人18水好多| 淫秽高清视频在线观看| 又黄又爽又免费观看的视频| 黄片大片在线免费观看| 亚洲专区中文字幕在线| 中出人妻视频一区二区| 免费av毛片视频| 成人18禁在线播放| 哪里可以看免费的av片| 久久久久免费精品人妻一区二区| 亚洲黑人精品在线| 亚洲专区国产一区二区| 一本一本综合久久| 国产三级中文精品| 亚洲一区高清亚洲精品| 在线永久观看黄色视频| 免费观看精品视频网站| 免费看美女性在线毛片视频| 老熟妇乱子伦视频在线观看| tocl精华| 神马国产精品三级电影在线观看| 成人三级做爰电影| 国产私拍福利视频在线观看| 真人做人爱边吃奶动态| 日本免费a在线| 色老头精品视频在线观看| 99riav亚洲国产免费| 三级国产精品欧美在线观看 | 夜夜看夜夜爽夜夜摸| 精品一区二区三区视频在线 | 我要搜黄色片| h日本视频在线播放| 十八禁人妻一区二区| 悠悠久久av| 很黄的视频免费| 午夜福利欧美成人| 国产aⅴ精品一区二区三区波| 床上黄色一级片| 91在线观看av| 搡老妇女老女人老熟妇| 成人精品一区二区免费| 19禁男女啪啪无遮挡网站| 成人欧美大片| 精华霜和精华液先用哪个| 成年人黄色毛片网站| 国内精品美女久久久久久| 一本久久中文字幕| 日韩中文字幕欧美一区二区| 男人的好看免费观看在线视频| 天天一区二区日本电影三级| 五月伊人婷婷丁香| 热99re8久久精品国产| 国产精品一区二区精品视频观看| 91久久精品国产一区二区成人 | 国内少妇人妻偷人精品xxx网站 | 亚洲,欧美精品.| 成人一区二区视频在线观看| 18禁美女被吸乳视频| 久久香蕉精品热| 美女 人体艺术 gogo| 亚洲国产精品sss在线观看| 国产久久久一区二区三区| 亚洲精品美女久久久久99蜜臀| 久久久久免费精品人妻一区二区| 欧美成人一区二区免费高清观看 | 最新在线观看一区二区三区| 亚洲激情在线av| 窝窝影院91人妻| 亚洲七黄色美女视频| 国语自产精品视频在线第100页| 高清毛片免费观看视频网站| av天堂在线播放| 亚洲avbb在线观看| 91av网一区二区| 在线播放国产精品三级| 日韩三级视频一区二区三区| 少妇熟女aⅴ在线视频| 成人无遮挡网站| 午夜视频精品福利| 免费看日本二区| 免费av毛片视频| 中文字幕高清在线视频| 欧美在线一区亚洲| 三级国产精品欧美在线观看 | 丰满人妻一区二区三区视频av | 岛国在线观看网站| 亚洲欧美一区二区三区黑人| www日本在线高清视频| 哪里可以看免费的av片| 久久天躁狠狠躁夜夜2o2o| 小蜜桃在线观看免费完整版高清| 美女高潮喷水抽搐中文字幕| 日韩欧美在线二视频| 丰满的人妻完整版| 亚洲第一电影网av| 国产一区二区在线观看日韩 | 国产伦人伦偷精品视频| 香蕉久久夜色| 女生性感内裤真人,穿戴方法视频| 淫妇啪啪啪对白视频| 88av欧美| 久久久精品大字幕| 国产成年人精品一区二区| 久久久精品大字幕| 一二三四社区在线视频社区8| 最新美女视频免费是黄的| 在线永久观看黄色视频| 黄色丝袜av网址大全| 午夜免费观看网址| 中文在线观看免费www的网站| 国产蜜桃级精品一区二区三区| 在线播放国产精品三级| 亚洲片人在线观看| 色av中文字幕| 午夜a级毛片| 熟女少妇亚洲综合色aaa.| 深夜精品福利| 国产精品99久久久久久久久| 午夜精品在线福利| 中文字幕人成人乱码亚洲影| 亚洲 欧美一区二区三区| 69av精品久久久久久| 老汉色∧v一级毛片| 最新在线观看一区二区三区| 日本精品一区二区三区蜜桃| 欧美日韩瑟瑟在线播放| 欧美色视频一区免费| 亚洲国产看品久久| 99热精品在线国产| 久久亚洲真实| 一进一出抽搐gif免费好疼| 亚洲精品456在线播放app | 搞女人的毛片| 给我免费播放毛片高清在线观看| 亚洲自拍偷在线| 欧美大码av| 噜噜噜噜噜久久久久久91| 99视频精品全部免费 在线 | 久久天堂一区二区三区四区| 午夜日韩欧美国产| www.www免费av| 亚洲成人精品中文字幕电影| 亚洲人成网站高清观看| 51午夜福利影视在线观看| 日韩欧美 国产精品| 日韩欧美在线二视频| or卡值多少钱| 成人特级黄色片久久久久久久| 亚洲精品456在线播放app | 欧美日韩瑟瑟在线播放| 熟女人妻精品中文字幕| 他把我摸到了高潮在线观看| 亚洲国产日韩欧美精品在线观看 | 精品久久久久久久久久免费视频| 少妇裸体淫交视频免费看高清| 人妻夜夜爽99麻豆av| 99热这里只有精品一区 | 一级毛片女人18水好多| 性色av乱码一区二区三区2| 网址你懂的国产日韩在线| 亚洲五月婷婷丁香| 老汉色∧v一级毛片| e午夜精品久久久久久久| 亚洲国产精品sss在线观看| 国产精品 欧美亚洲| 丰满人妻一区二区三区视频av | 美女扒开内裤让男人捅视频| 亚洲中文av在线| 亚洲人成网站在线播放欧美日韩| 成年免费大片在线观看| 哪里可以看免费的av片| 免费av不卡在线播放| 亚洲一区高清亚洲精品| 日韩三级视频一区二区三区| 小说图片视频综合网站| 美女高潮的动态| 国产午夜精品久久久久久| 精品不卡国产一区二区三区| 欧美中文日本在线观看视频| 久久精品影院6| 麻豆成人av在线观看| 欧美黄色淫秽网站| 国产午夜精品久久久久久| 黑人巨大精品欧美一区二区mp4| 夜夜夜夜夜久久久久| 九色成人免费人妻av| 国产精品亚洲一级av第二区| 99久国产av精品| 亚洲中文av在线| 亚洲国产高清在线一区二区三| 成人三级做爰电影|