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

    人人妻人人澡人人看| 久久ye,这里只有精品| 99视频精品全部免费 在线| 九色亚洲精品在线播放| 在线观看免费高清a一片| 日本与韩国留学比较| av在线老鸭窝| 免费人妻精品一区二区三区视频| 国产av一区二区精品久久| 美女视频免费永久观看网站| 亚洲国产最新在线播放| 日本色播在线视频| 黑丝袜美女国产一区| 亚洲少妇的诱惑av| av线在线观看网站| 一区二区三区免费毛片| 一区在线观看完整版| 精品酒店卫生间| 嘟嘟电影网在线观看| 久久 成人 亚洲| 一级毛片 在线播放| 熟女av电影| 日韩视频在线欧美| 最近2019中文字幕mv第一页| 国产视频首页在线观看| 成人毛片60女人毛片免费| 老司机影院成人| 日韩伦理黄色片| 亚洲av中文av极速乱| 国产男女超爽视频在线观看| 免费黄频网站在线观看国产| 免费观看无遮挡的男女| 黑人猛操日本美女一级片| 大码成人一级视频| 天天躁夜夜躁狠狠久久av| 日韩一区二区三区影片| 久久久久精品久久久久真实原创| 在线精品无人区一区二区三| 热99久久久久精品小说推荐| 欧美一级a爱片免费观看看| 中文字幕最新亚洲高清| 国产国语露脸激情在线看| 日本-黄色视频高清免费观看| 欧美+日韩+精品| 欧美变态另类bdsm刘玥| 极品少妇高潮喷水抽搐| 如何舔出高潮| 51国产日韩欧美| 夫妻午夜视频| 九九久久精品国产亚洲av麻豆| 最新的欧美精品一区二区| 日韩欧美一区视频在线观看| 高清毛片免费看| 午夜激情av网站| 水蜜桃什么品种好| 精品国产一区二区久久| 免费看光身美女| 国产av一区二区精品久久| 国产又色又爽无遮挡免| av在线app专区| 女性被躁到高潮视频| 人人妻人人添人人爽欧美一区卜| 亚洲精品色激情综合| 精品久久久久久久久亚洲| 99视频精品全部免费 在线| 亚洲精品一二三| 18禁观看日本| 国产黄色视频一区二区在线观看| 亚洲欧美清纯卡通| 亚洲欧洲精品一区二区精品久久久 | 久久久久久久亚洲中文字幕| 午夜激情av网站| 国产免费视频播放在线视频| 一边亲一边摸免费视频| 男的添女的下面高潮视频| 国产精品国产av在线观看| 久久久久久人妻| 国产精品99久久99久久久不卡 | 免费久久久久久久精品成人欧美视频 | 91在线精品国自产拍蜜月| 亚洲av.av天堂| 99热网站在线观看| 日本wwww免费看| 日日摸夜夜添夜夜添av毛片| 久久久久久久大尺度免费视频| 国产男人的电影天堂91| 最黄视频免费看| 少妇人妻久久综合中文| 国产极品粉嫩免费观看在线 | 国产一区二区三区av在线| 午夜福利网站1000一区二区三区| 少妇高潮的动态图| videos熟女内射| 一区二区三区免费毛片| 日韩中文字幕视频在线看片| 五月伊人婷婷丁香| 久久久久精品久久久久真实原创| 欧美日韩精品成人综合77777| 亚洲精品亚洲一区二区| 黄色一级大片看看| 秋霞伦理黄片| 久久久久久久大尺度免费视频| 国内精品宾馆在线| 我要看黄色一级片免费的| 大片免费播放器 马上看| 夫妻午夜视频| 免费久久久久久久精品成人欧美视频 | 亚洲婷婷狠狠爱综合网| av一本久久久久| 国产极品粉嫩免费观看在线 | 欧美性感艳星| 亚洲精品久久午夜乱码| 亚洲精品乱久久久久久| 久久女婷五月综合色啪小说| 蜜桃在线观看..| 91精品一卡2卡3卡4卡| 两个人免费观看高清视频| 狠狠精品人妻久久久久久综合| 国产精品99久久99久久久不卡 | 亚洲内射少妇av| 国产黄频视频在线观看| 五月开心婷婷网| 中国美白少妇内射xxxbb| av在线app专区| 中文字幕制服av| 我要看黄色一级片免费的| 91精品国产国语对白视频| 亚洲伊人久久精品综合| 午夜日本视频在线| 母亲3免费完整高清在线观看 | 亚洲精品,欧美精品| 久久久国产精品麻豆| 在线精品无人区一区二区三| 黑丝袜美女国产一区| 日韩av不卡免费在线播放| 亚洲欧美一区二区三区黑人 | 午夜免费鲁丝| 嫩草影院入口| 99热全是精品| 久久久久国产精品人妻一区二区| 视频中文字幕在线观看| 国产精品无大码| 插逼视频在线观看| 国产成人精品久久久久久| 国产永久视频网站| 亚洲国产精品专区欧美| 亚洲欧洲日产国产| 91久久精品电影网| 欧美日本中文国产一区发布| 久久av网站| 亚洲中文av在线| av视频免费观看在线观看| 人人妻人人爽人人添夜夜欢视频| 看十八女毛片水多多多| 人体艺术视频欧美日本| 久久鲁丝午夜福利片| 亚洲国产日韩一区二区| 久久精品国产亚洲av涩爱| 国产一区有黄有色的免费视频| 亚洲精品一区蜜桃| 午夜福利视频在线观看免费| 国产精品不卡视频一区二区| 晚上一个人看的免费电影| 亚洲成人手机| 2021少妇久久久久久久久久久| 亚洲欧美一区二区三区国产| 蜜桃国产av成人99| 亚洲精品成人av观看孕妇| 一本久久精品| 亚洲成人手机| 亚洲精品自拍成人| 一级毛片aaaaaa免费看小| 又黄又爽又刺激的免费视频.| 又黄又爽又刺激的免费视频.| 亚洲精品久久成人aⅴ小说 | 精品少妇黑人巨大在线播放| 久久久久国产网址| 亚洲人成网站在线观看播放| 午夜影院在线不卡| 黄色欧美视频在线观看| 国产色爽女视频免费观看| 久久精品国产a三级三级三级| 狠狠精品人妻久久久久久综合| 国产69精品久久久久777片| 日本黄色片子视频| 丝袜在线中文字幕| 国产欧美日韩一区二区三区在线 | 日本av免费视频播放| 国产精品偷伦视频观看了| 日韩精品有码人妻一区| 亚洲伊人久久精品综合| 日韩精品有码人妻一区| av免费在线看不卡| 久久久久久伊人网av| 精品一区二区免费观看| 亚洲综合色惰| 女人精品久久久久毛片| 中国国产av一级| 国产成人a∨麻豆精品| 五月天丁香电影| 日韩 亚洲 欧美在线| 国产69精品久久久久777片| 极品少妇高潮喷水抽搐| 赤兔流量卡办理| 肉色欧美久久久久久久蜜桃| 你懂的网址亚洲精品在线观看| 熟女av电影| 国产成人一区二区在线| 在线播放无遮挡| av有码第一页| 美女大奶头黄色视频| 亚洲av男天堂| 99久久中文字幕三级久久日本| 亚洲怡红院男人天堂| 久久精品久久久久久噜噜老黄| 日本vs欧美在线观看视频| 色婷婷av一区二区三区视频| 国产高清国产精品国产三级| 午夜91福利影院| 香蕉精品网在线| 国产精品国产av在线观看| 亚洲国产精品国产精品| 美女视频免费永久观看网站| 一级毛片aaaaaa免费看小| 国产精品不卡视频一区二区| 国产精品一区二区在线观看99| 久久人妻熟女aⅴ| 看非洲黑人一级黄片| 日韩三级伦理在线观看| 午夜久久久在线观看| av又黄又爽大尺度在线免费看| 成人国语在线视频| 欧美少妇被猛烈插入视频| 亚洲欧美清纯卡通| 亚洲精品乱码久久久v下载方式| 婷婷成人精品国产| 免费黄频网站在线观看国产| 成人综合一区亚洲| 亚洲性久久影院| 欧美三级亚洲精品| 国产一区二区三区av在线| 亚洲av成人精品一二三区| 五月玫瑰六月丁香| 久久久久久久精品精品| 国产精品女同一区二区软件| 国产成人a∨麻豆精品| 春色校园在线视频观看| 99热全是精品| 久久精品久久久久久久性| 欧美日韩亚洲高清精品| 国产在线一区二区三区精| 午夜福利视频在线观看免费| 国产不卡av网站在线观看| 成年人免费黄色播放视频| 欧美xxⅹ黑人| 中文天堂在线官网| 国产精品国产三级专区第一集| 国产精品一国产av| 亚洲精品乱久久久久久| 一级毛片电影观看| 激情五月婷婷亚洲| 国产免费视频播放在线视频| 成人二区视频| 久久久精品免费免费高清| 精品人妻熟女毛片av久久网站| 国产探花极品一区二区| 成年美女黄网站色视频大全免费 | 欧美日韩视频高清一区二区三区二| 最近中文字幕高清免费大全6| 最近2019中文字幕mv第一页| 亚洲美女搞黄在线观看| 久久久久久久亚洲中文字幕| 国产精品不卡视频一区二区| 国产精品国产三级国产av玫瑰| 亚洲av成人精品一二三区| 男女边吃奶边做爰视频| 久久人人爽人人片av| 国产免费又黄又爽又色| 肉色欧美久久久久久久蜜桃| 99久久中文字幕三级久久日本| 好男人视频免费观看在线| 纯流量卡能插随身wifi吗| 国产成人精品福利久久| 免费人妻精品一区二区三区视频| 狂野欧美激情性xxxx在线观看| 国产精品一二三区在线看| 日韩电影二区| 男女免费视频国产| 国产男女内射视频| 国产69精品久久久久777片| av有码第一页| 欧美精品国产亚洲| 丰满饥渴人妻一区二区三| 国产精品国产av在线观看| 日本vs欧美在线观看视频| 国产一级毛片在线| 久久久国产精品麻豆| 中文字幕最新亚洲高清| 亚洲人成网站在线观看播放| 亚洲av福利一区| 中文字幕免费在线视频6| 午夜日本视频在线| 五月天丁香电影| 五月开心婷婷网| 丰满迷人的少妇在线观看| 青春草视频在线免费观看| 成人国产麻豆网| 超碰97精品在线观看| 男女免费视频国产| 亚洲av免费高清在线观看| 夫妻性生交免费视频一级片| 亚洲国产精品一区三区| 91精品一卡2卡3卡4卡| 精品久久蜜臀av无| 亚洲欧美色中文字幕在线| 亚洲国产av影院在线观看| 久久精品久久久久久久性| 久久99一区二区三区| 亚洲综合精品二区| 免费av不卡在线播放| 国产高清有码在线观看视频| 欧美日本中文国产一区发布| 国产精品女同一区二区软件| 国产又色又爽无遮挡免| 一级a做视频免费观看| 国产深夜福利视频在线观看| 日本wwww免费看| 久久精品久久久久久久性| 久久99一区二区三区| 婷婷色综合www| 国产精品一区二区三区四区免费观看| 青春草亚洲视频在线观看| 国产免费又黄又爽又色| 亚洲一级一片aⅴ在线观看| 亚洲精品乱码久久久v下载方式| 免费观看a级毛片全部| 精品人妻熟女av久视频| 久久免费观看电影| 久久久久国产精品人妻一区二区| 亚洲精品视频女| 亚洲精品国产av成人精品| 欧美激情极品国产一区二区三区 | 精品一区二区三区视频在线| 成人毛片60女人毛片免费| 成人无遮挡网站| 一个人看视频在线观看www免费| 亚洲精品国产色婷婷电影| 午夜影院在线不卡| 日韩伦理黄色片| av国产久精品久网站免费入址| 97超碰精品成人国产| 国产欧美另类精品又又久久亚洲欧美| 丁香六月天网| 日韩中文字幕视频在线看片| 中文字幕免费在线视频6| 久久99精品国语久久久| 曰老女人黄片| 寂寞人妻少妇视频99o| 天堂8中文在线网| 久热这里只有精品99| 午夜久久久在线观看| 国产淫语在线视频| 大码成人一级视频| 久久狼人影院| 男女国产视频网站| 日韩一区二区视频免费看| 一边亲一边摸免费视频| 久久狼人影院| 777米奇影视久久| 亚洲国产精品国产精品| 极品人妻少妇av视频| 熟妇人妻不卡中文字幕| 欧美日韩精品成人综合77777| 亚洲国产欧美日韩在线播放| 99热国产这里只有精品6| 纯流量卡能插随身wifi吗| 日韩伦理黄色片| 嘟嘟电影网在线观看| 免费高清在线观看视频在线观看| 国产日韩欧美视频二区| 日本av免费视频播放| 婷婷成人精品国产| 国产极品粉嫩免费观看在线 | 一区二区三区精品91| a级毛色黄片| 日韩一本色道免费dvd| 最近的中文字幕免费完整| 亚洲av欧美aⅴ国产| 亚洲国产精品成人久久小说| 男的添女的下面高潮视频| 老司机影院成人| 免费观看在线日韩| 国产片特级美女逼逼视频| 国产有黄有色有爽视频| 国模一区二区三区四区视频| 成年美女黄网站色视频大全免费 | 亚洲美女黄色视频免费看| 日本黄色片子视频| 99热国产这里只有精品6| 亚洲精品乱码久久久v下载方式| 菩萨蛮人人尽说江南好唐韦庄| 中文字幕制服av| 美女主播在线视频| 精品99又大又爽又粗少妇毛片| 午夜av观看不卡| 18禁动态无遮挡网站| 国产 一区精品| 亚洲国产欧美日韩在线播放| 精品久久久久久电影网| 乱码一卡2卡4卡精品| 免费不卡的大黄色大毛片视频在线观看| 国产在线一区二区三区精| 一区二区av电影网| 又黄又爽又刺激的免费视频.| 啦啦啦在线观看免费高清www| 亚洲精品久久成人aⅴ小说 | 亚洲精品久久成人aⅴ小说 | 一区二区日韩欧美中文字幕 | 美女内射精品一级片tv| 人人妻人人澡人人爽人人夜夜| 最后的刺客免费高清国语| 你懂的网址亚洲精品在线观看| 亚洲综合精品二区| 男女免费视频国产| 亚洲精品成人av观看孕妇| freevideosex欧美| 亚洲成人一二三区av| 如日韩欧美国产精品一区二区三区 | 午夜视频国产福利| 久久久久久久精品精品| 久久久久久久久久人人人人人人| 插阴视频在线观看视频| 狠狠婷婷综合久久久久久88av| 亚洲四区av| 久久久亚洲精品成人影院| 国产精品99久久99久久久不卡 | 久久久国产欧美日韩av| 免费久久久久久久精品成人欧美视频 | 日本av手机在线免费观看| 中文字幕制服av| 精品久久久久久久久亚洲| 久热久热在线精品观看| 看非洲黑人一级黄片| 如何舔出高潮| 十分钟在线观看高清视频www| 亚洲三级黄色毛片| 18禁在线播放成人免费| 伦理电影大哥的女人| 一级毛片我不卡| 全区人妻精品视频| 欧美精品亚洲一区二区| 欧美精品一区二区免费开放| 亚洲图色成人| 精品久久久久久久久av| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 超色免费av| 亚洲av不卡在线观看| 日日撸夜夜添| 少妇被粗大的猛进出69影院 | 最近手机中文字幕大全| 在线看a的网站| 久久精品夜色国产| 黄片无遮挡物在线观看| 久久久久久人妻| 日韩中文字幕视频在线看片| 亚洲欧美成人精品一区二区| 啦啦啦中文免费视频观看日本| 五月开心婷婷网| 大香蕉97超碰在线| 一级,二级,三级黄色视频| 丝袜在线中文字幕| 中国三级夫妇交换| 考比视频在线观看| 91在线精品国自产拍蜜月| 99国产综合亚洲精品| 中文乱码字字幕精品一区二区三区| 亚洲精品久久久久久婷婷小说| 免费看av在线观看网站| 七月丁香在线播放| av专区在线播放| 日本av手机在线免费观看| 久久久久久久久久久免费av| 国语对白做爰xxxⅹ性视频网站| 寂寞人妻少妇视频99o| 午夜激情av网站| 国产 一区精品| 中文字幕人妻熟人妻熟丝袜美| 日本猛色少妇xxxxx猛交久久| 蜜桃在线观看..| 国产成人免费无遮挡视频| 亚洲第一区二区三区不卡| 韩国av在线不卡| 精品国产乱码久久久久久小说| 国产亚洲av片在线观看秒播厂| 97在线视频观看| 日日摸夜夜添夜夜添av毛片| 国产免费视频播放在线视频| 两个人的视频大全免费| 欧美成人精品欧美一级黄| 精品久久久久久电影网| 青春草视频在线免费观看| 久久久久人妻精品一区果冻| 亚洲成色77777| 一本大道久久a久久精品| 国产黄色免费在线视频| 中文字幕制服av| 最近最新中文字幕免费大全7| 一级毛片aaaaaa免费看小| 9色porny在线观看| 国产熟女欧美一区二区| 99热全是精品| 少妇丰满av| 精品一区二区三卡| 免费大片18禁| 肉色欧美久久久久久久蜜桃| 日本黄色片子视频| 亚洲精品乱久久久久久| 免费看不卡的av| 日日爽夜夜爽网站| 交换朋友夫妻互换小说| 纯流量卡能插随身wifi吗| 欧美精品国产亚洲| 麻豆成人av视频| videossex国产| 午夜日本视频在线| 99热这里只有精品一区| 免费av中文字幕在线| 99re6热这里在线精品视频| 午夜久久久在线观看| tube8黄色片| 久久久国产一区二区| 在线亚洲精品国产二区图片欧美 | freevideosex欧美| 老司机影院成人| 免费av不卡在线播放| 能在线免费看毛片的网站| 日韩欧美一区视频在线观看| 97在线视频观看| 久久狼人影院| 国国产精品蜜臀av免费| 日韩欧美一区视频在线观看| 国产视频首页在线观看| 中国美白少妇内射xxxbb| 黄色毛片三级朝国网站| 国产精品一区二区在线观看99| 99久久综合免费| 亚洲人成77777在线视频| 国产在视频线精品| 美女脱内裤让男人舔精品视频| 午夜91福利影院| 永久网站在线| 狂野欧美白嫩少妇大欣赏| 国产一区二区三区av在线| 久热久热在线精品观看| av免费观看日本| 女的被弄到高潮叫床怎么办| 制服诱惑二区| 中文字幕制服av| 最黄视频免费看| 亚洲图色成人| 爱豆传媒免费全集在线观看| 国产精品麻豆人妻色哟哟久久| 2018国产大陆天天弄谢| h视频一区二区三区| 青春草亚洲视频在线观看| 熟女av电影| 啦啦啦中文免费视频观看日本| 成年人免费黄色播放视频| 国产免费一区二区三区四区乱码| 欧美三级亚洲精品| 最近2019中文字幕mv第一页| 亚洲精品视频女| 人妻系列 视频| 国产一区二区三区av在线| 欧美 日韩 精品 国产| 又大又黄又爽视频免费| 久久精品国产a三级三级三级| xxx大片免费视频| 在线观看免费视频网站a站| av天堂久久9| www.av在线官网国产| 午夜免费鲁丝| 成年人免费黄色播放视频| 好男人视频免费观看在线| 免费高清在线观看日韩| 日本欧美国产在线视频| a级毛片在线看网站| 另类精品久久| 亚洲国产av新网站| 国产日韩欧美亚洲二区| 亚洲内射少妇av| 久久av网站| 美女内射精品一级片tv| 国产精品久久久久久久电影| 国产片特级美女逼逼视频| 国产精品 国内视频| 亚洲精品456在线播放app| 又大又黄又爽视频免费| 精品亚洲乱码少妇综合久久| 18禁裸乳无遮挡动漫免费视频| 一区二区日韩欧美中文字幕 | 国产一区二区在线观看av| 日韩电影二区| 国产视频首页在线观看| 在线精品无人区一区二区三| 好男人视频免费观看在线| 精品熟女少妇av免费看| 在线 av 中文字幕| 国产精品一国产av| 午夜视频国产福利| 我要看黄色一级片免费的| 日韩精品有码人妻一区| 国产毛片在线视频| 欧美精品亚洲一区二区| 18禁裸乳无遮挡动漫免费视频|