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

    一个人免费在线观看的高清视频| 国产精品美女特级片免费视频播放器 | 精品一区二区三区视频在线观看免费| 成人18禁在线播放| 免费观看的影片在线观看| 黄片小视频在线播放| 1000部很黄的大片| 成人精品一区二区免费| 美女cb高潮喷水在线观看 | 国产精品久久电影中文字幕| 男女之事视频高清在线观看| 免费看a级黄色片| 757午夜福利合集在线观看| 亚洲国产精品久久男人天堂| 日本 欧美在线| 国产私拍福利视频在线观看| 国产午夜福利久久久久久| 国产成人av激情在线播放| 久久九九热精品免费| 波多野结衣高清作品| 好男人电影高清在线观看| 搡老岳熟女国产| 亚洲专区字幕在线| 欧洲精品卡2卡3卡4卡5卡区| 免费观看的影片在线观看| 午夜福利免费观看在线| 97人妻精品一区二区三区麻豆| 亚洲最大成人中文| 欧美成人免费av一区二区三区| 日韩高清综合在线| 久久中文看片网| 波多野结衣高清无吗| 国产精品乱码一区二三区的特点| 在线播放国产精品三级| 99精品久久久久人妻精品| 免费看美女性在线毛片视频| 在线观看一区二区三区| 国产亚洲av高清不卡| 国产av麻豆久久久久久久| 国产成人欧美在线观看| a在线观看视频网站| 色综合站精品国产| 午夜福利欧美成人| 不卡一级毛片| 欧美乱色亚洲激情| 午夜成年电影在线免费观看| 毛片女人毛片| 精品久久久久久久毛片微露脸| 国产真人三级小视频在线观看| 一区二区三区激情视频| 黑人操中国人逼视频| 成人18禁在线播放| 99久久国产精品久久久| 69av精品久久久久久| 国产成年人精品一区二区| 午夜免费成人在线视频| 麻豆成人av在线观看| 中文字幕av在线有码专区| 精品久久久久久成人av| 亚洲欧美日韩高清专用| 亚洲精品中文字幕一二三四区| 欧美黑人欧美精品刺激| 亚洲乱码一区二区免费版| 国产精品久久久久久久电影 | 国产高清三级在线| 国产97色在线日韩免费| 免费看十八禁软件| 亚洲性夜色夜夜综合| 91麻豆精品激情在线观看国产| 久久欧美精品欧美久久欧美| 亚洲欧美精品综合一区二区三区| 国产伦精品一区二区三区四那| 成人三级做爰电影| 在线观看日韩欧美| 亚洲国产精品成人综合色| 成年女人永久免费观看视频| 精品国产三级普通话版| 日韩免费av在线播放| 一个人看的www免费观看视频| 美女被艹到高潮喷水动态| 欧美中文日本在线观看视频| 欧美成狂野欧美在线观看| 搡老妇女老女人老熟妇| 日韩三级视频一区二区三区| 久久亚洲精品不卡| 欧美+亚洲+日韩+国产| 变态另类成人亚洲欧美熟女| 午夜日韩欧美国产| 日韩欧美国产在线观看| 又爽又黄无遮挡网站| 欧美日韩一级在线毛片| 男女视频在线观看网站免费| 国产三级黄色录像| 一进一出好大好爽视频| 村上凉子中文字幕在线| 国产高潮美女av| 久久国产精品影院| 国产亚洲精品一区二区www| 国产探花在线观看一区二区| 亚洲男人的天堂狠狠| 99热6这里只有精品| 久久性视频一级片| 亚洲欧美日韩卡通动漫| 亚洲无线在线观看| 岛国视频午夜一区免费看| 国产不卡一卡二| 午夜免费成人在线视频| 免费在线观看影片大全网站| 91在线观看av| 亚洲无线观看免费| 成人三级做爰电影| 97超级碰碰碰精品色视频在线观看| 嫁个100分男人电影在线观看| 听说在线观看完整版免费高清| 亚洲,欧美精品.| 少妇丰满av| 亚洲自拍偷在线| 999久久久国产精品视频| 全区人妻精品视频| 亚洲色图av天堂| 成人永久免费在线观看视频| 国产精品野战在线观看| 99热这里只有是精品50| 国产激情久久老熟女| 亚洲人成伊人成综合网2020| tocl精华| 丝袜人妻中文字幕| 亚洲国产精品成人综合色| 制服丝袜大香蕉在线| 久久精品综合一区二区三区| 一二三四社区在线视频社区8| 欧美午夜高清在线| 精品国产三级普通话版| 国产成人系列免费观看| 亚洲真实伦在线观看| 成人av一区二区三区在线看| 九九在线视频观看精品| 婷婷精品国产亚洲av| 日本熟妇午夜| 中文资源天堂在线| 99国产精品99久久久久| 九色国产91popny在线| 美女扒开内裤让男人捅视频| 很黄的视频免费| 热99re8久久精品国产| 桃红色精品国产亚洲av| 精品久久久久久久末码| 99国产精品一区二区蜜桃av| 国内毛片毛片毛片毛片毛片| 国产午夜精品久久久久久| 久久久久久久久久黄片| 1024香蕉在线观看| 久久99热这里只有精品18| 日本 av在线| 国产真人三级小视频在线观看| 日韩欧美精品v在线| 精品人妻1区二区| 国产爱豆传媒在线观看| 欧美三级亚洲精品| 在线观看66精品国产| 激情在线观看视频在线高清| 欧美日韩瑟瑟在线播放| 日本熟妇午夜| 亚洲国产精品合色在线| 成人高潮视频无遮挡免费网站| 两性午夜刺激爽爽歪歪视频在线观看| 久久久久免费精品人妻一区二区| 9191精品国产免费久久| 99久久99久久久精品蜜桃| 国产亚洲精品综合一区在线观看| 久久久久久久久中文| 亚洲午夜理论影院| 9191精品国产免费久久| 在线观看舔阴道视频| 香蕉久久夜色| 特级一级黄色大片| 国产精品99久久久久久久久| 好看av亚洲va欧美ⅴa在| 久久精品亚洲精品国产色婷小说| 欧美乱码精品一区二区三区| 国产亚洲精品av在线| 国产精品亚洲美女久久久| 搞女人的毛片| 亚洲国产精品合色在线| 欧美日韩精品网址| 三级男女做爰猛烈吃奶摸视频| 小说图片视频综合网站| 欧美三级亚洲精品| 悠悠久久av| 男人舔女人的私密视频| 神马国产精品三级电影在线观看| 免费看a级黄色片| 搡老岳熟女国产| av欧美777| 18禁国产床啪视频网站| 天堂动漫精品| 变态另类丝袜制服| a级毛片a级免费在线| 国产精品精品国产色婷婷| 成在线人永久免费视频| 俄罗斯特黄特色一大片| 亚洲七黄色美女视频| 久久亚洲真实| 免费看十八禁软件| 午夜精品在线福利| 在线观看免费午夜福利视频| 日本撒尿小便嘘嘘汇集6| 日本黄色片子视频| 日韩高清综合在线| 亚洲欧美日韩卡通动漫| 欧美zozozo另类| 日韩高清综合在线| 国产淫片久久久久久久久 | 国产成人av激情在线播放| 国产激情欧美一区二区| 亚洲午夜理论影院| 欧美日韩瑟瑟在线播放| 小说图片视频综合网站| 999精品在线视频| 日韩欧美国产一区二区入口| 亚洲国产精品sss在线观看| 久久精品aⅴ一区二区三区四区| 久久久久久久久久黄片| 亚洲国产欧美人成| 观看美女的网站| 18美女黄网站色大片免费观看| 1024手机看黄色片| 桃色一区二区三区在线观看| 成人亚洲精品av一区二区| 国产aⅴ精品一区二区三区波| 国产欧美日韩精品一区二区| 久久精品国产综合久久久| 757午夜福利合集在线观看| 全区人妻精品视频| www.自偷自拍.com| 美女免费视频网站| av欧美777| 天堂动漫精品| 国产又色又爽无遮挡免费看| av欧美777| 男人舔女人的私密视频| 久久久国产欧美日韩av| 久久久成人免费电影| 久久香蕉国产精品| 亚洲性夜色夜夜综合| 最近最新免费中文字幕在线| 99精品久久久久人妻精品| 少妇的丰满在线观看| 久久国产精品人妻蜜桃| 成人欧美大片| 不卡av一区二区三区| 黑人巨大精品欧美一区二区mp4| 午夜福利高清视频| 五月伊人婷婷丁香| 欧美黑人欧美精品刺激| 亚洲九九香蕉| 一本久久中文字幕| 在线看三级毛片| 99热这里只有精品一区 | 国语自产精品视频在线第100页| 黄色 视频免费看| 男女午夜视频在线观看| 婷婷六月久久综合丁香| 国产三级黄色录像| 99久久无色码亚洲精品果冻| 国产真实乱freesex| 人妻久久中文字幕网| 亚洲在线观看片| 天天一区二区日本电影三级| 欧美日韩黄片免| 97人妻精品一区二区三区麻豆| 岛国视频午夜一区免费看| 精品国产乱码久久久久久男人| 一本久久中文字幕| 美女高潮的动态| 两个人看的免费小视频| 又黄又粗又硬又大视频| 国内精品美女久久久久久| av国产免费在线观看| 久久精品影院6| 亚洲国产精品成人综合色| 欧美日韩福利视频一区二区| 欧美另类亚洲清纯唯美| 美女被艹到高潮喷水动态| 亚洲午夜精品一区,二区,三区| 日韩免费av在线播放| 久久精品国产亚洲av香蕉五月| 国产精品爽爽va在线观看网站| 午夜福利在线观看吧| xxxwww97欧美| 淫妇啪啪啪对白视频| 在线十欧美十亚洲十日本专区| 亚洲一区二区三区色噜噜| 99re在线观看精品视频| 色综合站精品国产| 性色avwww在线观看| 亚洲一区高清亚洲精品| 夜夜躁狠狠躁天天躁| 欧美不卡视频在线免费观看| 好男人电影高清在线观看| 国产欧美日韩一区二区三| 国产精品精品国产色婷婷| 久久九九热精品免费| 熟女少妇亚洲综合色aaa.| 中文字幕高清在线视频| 日本黄色视频三级网站网址| 美女大奶头视频| 亚洲无线在线观看| 国产精华一区二区三区| www国产在线视频色| 日日夜夜操网爽| 国产精品亚洲一级av第二区| 日韩欧美精品v在线| 国产一区二区激情短视频| 日韩高清综合在线| 亚洲欧美精品综合久久99| 亚洲av中文字字幕乱码综合| 女同久久另类99精品国产91| 亚洲欧美日韩高清专用| 亚洲av电影不卡..在线观看| 少妇的逼水好多| 亚洲中文字幕日韩| 日韩 欧美 亚洲 中文字幕| 噜噜噜噜噜久久久久久91| 成人无遮挡网站| 91麻豆av在线| 搡老熟女国产l中国老女人| 国产高潮美女av| 一进一出好大好爽视频| 国产成人精品久久二区二区91| xxxwww97欧美| 99国产精品99久久久久| 国产99白浆流出| 免费av毛片视频| 91久久精品国产一区二区成人 | 久久久久久久午夜电影| 免费观看的影片在线观看| www.精华液| 欧美日韩瑟瑟在线播放| 男女午夜视频在线观看| 亚洲国产欧美人成| 亚洲欧美日韩高清专用| 亚洲天堂国产精品一区在线| www日本在线高清视频| 人人妻,人人澡人人爽秒播| 国产成人福利小说| 一本一本综合久久| 观看免费一级毛片| av天堂中文字幕网| 欧美在线黄色| 观看美女的网站| 97超级碰碰碰精品色视频在线观看| 99久久精品一区二区三区| 不卡一级毛片| 欧美黄色片欧美黄色片| 又大又爽又粗| 99久久成人亚洲精品观看| 老熟妇乱子伦视频在线观看| 18禁黄网站禁片免费观看直播| 俄罗斯特黄特色一大片| 在线免费观看不下载黄p国产 | 日日夜夜操网爽| 麻豆久久精品国产亚洲av| 宅男免费午夜| 十八禁人妻一区二区| 丰满人妻熟妇乱又伦精品不卡| 曰老女人黄片| avwww免费| 免费人成视频x8x8入口观看| 99久久综合精品五月天人人| 国产免费男女视频| 国产三级在线视频| 久久久精品欧美日韩精品| 午夜激情福利司机影院| 亚洲av中文字字幕乱码综合| 欧美一区二区精品小视频在线| 国产成人福利小说| 欧美日韩一级在线毛片| 成人欧美大片| 成人国产一区最新在线观看| 国产亚洲欧美在线一区二区| 欧美日韩福利视频一区二区| 国产精品野战在线观看| 国产乱人视频| 国产亚洲欧美在线一区二区| 18禁黄网站禁片午夜丰满| 国产激情欧美一区二区| 亚洲专区中文字幕在线| 国产成人福利小说| 亚洲人与动物交配视频| 午夜日韩欧美国产| 亚洲精品一区av在线观看| 欧美性猛交╳xxx乱大交人| 亚洲一区高清亚洲精品| 亚洲va日本ⅴa欧美va伊人久久| 啦啦啦免费观看视频1| 亚洲aⅴ乱码一区二区在线播放| 亚洲国产中文字幕在线视频| 亚洲欧美精品综合一区二区三区| 久久精品国产亚洲av香蕉五月| 中国美女看黄片| 真实男女啪啪啪动态图| 三级国产精品欧美在线观看 | 久久久久九九精品影院| 91老司机精品| 欧美一区二区精品小视频在线| 久久久国产精品麻豆| xxxwww97欧美| 在线十欧美十亚洲十日本专区| avwww免费| 18禁裸乳无遮挡免费网站照片| 国产在线精品亚洲第一网站| 12—13女人毛片做爰片一| 亚洲无线在线观看| 成人鲁丝片一二三区免费| 757午夜福利合集在线观看| 精品国产三级普通话版| 十八禁人妻一区二区| 国产探花在线观看一区二区| 又粗又爽又猛毛片免费看| 亚洲性夜色夜夜综合| 女人被狂操c到高潮| 亚洲av五月六月丁香网| 国产精品一及| 午夜精品久久久久久毛片777| 色综合婷婷激情| 美女免费视频网站| 国模一区二区三区四区视频 | 亚洲成人中文字幕在线播放| 欧美黑人巨大hd| xxx96com| 亚洲国产高清在线一区二区三| 亚洲精品粉嫩美女一区| 精品无人区乱码1区二区| 欧美绝顶高潮抽搐喷水| 国产精品av视频在线免费观看| 免费观看的影片在线观看| 真人做人爱边吃奶动态| 亚洲国产精品999在线| 国产精品久久视频播放| 一级毛片女人18水好多| 久久人人精品亚洲av| 亚洲av熟女| 久久久久久久久免费视频了| 欧美+亚洲+日韩+国产| 精品免费久久久久久久清纯| 国产黄a三级三级三级人| 在线观看免费午夜福利视频| 可以在线观看毛片的网站| 午夜福利高清视频| 天天添夜夜摸| 欧美丝袜亚洲另类 | 少妇人妻一区二区三区视频| 国产精品av久久久久免费| 最新在线观看一区二区三区| 亚洲专区字幕在线| 成人三级黄色视频| 在线观看午夜福利视频| 日日摸夜夜添夜夜添小说| www.自偷自拍.com| 亚洲国产精品成人综合色| xxx96com| 国产三级在线视频| 这个男人来自地球电影免费观看| 免费一级毛片在线播放高清视频| 精品人妻1区二区| 中国美女看黄片| 成人精品一区二区免费| 在线观看一区二区三区| 欧美高清成人免费视频www| 亚洲人成网站在线播放欧美日韩| 亚洲,欧美精品.| 中文字幕精品亚洲无线码一区| 欧美不卡视频在线免费观看| 亚洲精品色激情综合| 夜夜夜夜夜久久久久| 黑人巨大精品欧美一区二区mp4| 久久久久久久午夜电影| 精品99又大又爽又粗少妇毛片 | 国产亚洲av高清不卡| 90打野战视频偷拍视频| 欧美日韩亚洲国产一区二区在线观看| 人妻夜夜爽99麻豆av| 少妇裸体淫交视频免费看高清| 久久热在线av| 少妇人妻一区二区三区视频| 国产亚洲av嫩草精品影院| 亚洲欧美日韩东京热| 村上凉子中文字幕在线| 久久久久国内视频| 欧美乱码精品一区二区三区| 国产私拍福利视频在线观看| 国产成人一区二区三区免费视频网站| 日韩精品中文字幕看吧| 长腿黑丝高跟| 狂野欧美激情性xxxx| 亚洲电影在线观看av| 免费大片18禁| 一级毛片女人18水好多| 色噜噜av男人的天堂激情| 久99久视频精品免费| 婷婷精品国产亚洲av| 欧美性猛交黑人性爽| 欧洲精品卡2卡3卡4卡5卡区| 此物有八面人人有两片| 久久久久九九精品影院| 老司机深夜福利视频在线观看| 亚洲国产欧洲综合997久久,| 老司机福利观看| 亚洲 国产 在线| 性色avwww在线观看| 特级一级黄色大片| 九色成人免费人妻av| svipshipincom国产片| 女人高潮潮喷娇喘18禁视频| 久久九九热精品免费| 国产毛片a区久久久久| 综合色av麻豆| 亚洲在线自拍视频| 最新中文字幕久久久久 | 精品国产乱子伦一区二区三区| 国产一区在线观看成人免费| 热99re8久久精品国产| 亚洲 欧美 日韩 在线 免费| 欧美激情在线99| 久久精品国产综合久久久| 五月玫瑰六月丁香| 亚洲精品乱码久久久v下载方式 | 国产高潮美女av| 在线观看午夜福利视频| 亚洲欧美日韩东京热| 淫秽高清视频在线观看| 亚洲在线观看片| 久久99热这里只有精品18| 久久久国产成人精品二区| 欧美一区二区精品小视频在线| 这个男人来自地球电影免费观看| 日韩国内少妇激情av| 综合色av麻豆| 国产探花在线观看一区二区| 他把我摸到了高潮在线观看| 色综合欧美亚洲国产小说| 人人妻,人人澡人人爽秒播| 真人一进一出gif抽搐免费| 变态另类成人亚洲欧美熟女| 日本a在线网址| 韩国av一区二区三区四区| 国产一区二区在线av高清观看| 最近最新中文字幕大全电影3| 午夜a级毛片| 久久精品国产99精品国产亚洲性色| 男人舔女人的私密视频| tocl精华| 夜夜看夜夜爽夜夜摸| 一个人观看的视频www高清免费观看 | www.熟女人妻精品国产| 欧美在线黄色| av片东京热男人的天堂| 一区二区三区高清视频在线| 日本一二三区视频观看| av天堂中文字幕网| 啪啪无遮挡十八禁网站| 国产成+人综合+亚洲专区| or卡值多少钱| 窝窝影院91人妻| 国产欧美日韩精品亚洲av| 男女那种视频在线观看| 久久天堂一区二区三区四区| 国产成人啪精品午夜网站| 欧美在线黄色| 亚洲第一欧美日韩一区二区三区| 午夜影院日韩av| av片东京热男人的天堂| 国产精品久久久久久亚洲av鲁大| 久99久视频精品免费| 久久这里只有精品19| 麻豆久久精品国产亚洲av| 久久九九热精品免费| 国产精品久久久久久久电影 | 亚洲av熟女| 成人高潮视频无遮挡免费网站| 日本免费a在线| 特级一级黄色大片| 少妇裸体淫交视频免费看高清| 91九色精品人成在线观看| 中文字幕熟女人妻在线| 在线国产一区二区在线| 午夜福利免费观看在线| 波多野结衣巨乳人妻| 国产精品九九99| 国产激情偷乱视频一区二区| 国产精品九九99| 亚洲av中文字字幕乱码综合| 欧美日韩中文字幕国产精品一区二区三区| 成人一区二区视频在线观看| 国产黄a三级三级三级人| 日日夜夜操网爽| 琪琪午夜伦伦电影理论片6080| 天堂av国产一区二区熟女人妻| 琪琪午夜伦伦电影理论片6080| 麻豆一二三区av精品| 露出奶头的视频| 搡老熟女国产l中国老女人| 国内毛片毛片毛片毛片毛片| 熟妇人妻久久中文字幕3abv| 国产黄色小视频在线观看| 超碰成人久久| cao死你这个sao货| 国产精品一区二区精品视频观看| 午夜福利在线观看吧| 久久香蕉精品热| 久久性视频一级片| 黄色 视频免费看| 热99re8久久精品国产| 亚洲精品一区av在线观看|