A black hole toy model with non-local and boundary modes from non-trivial boundary conditions

We study gauge theories between two parallel boundaries with non-trivial boundary conditions, which serve as a toy model for black hole background with two boundaries near the horizon and infinite, aiming for a better understanding of the Bekenstein–Hawking entropy. The new set of boundary conditions allows boundary modes and non-local modes that interplay between the two boundaries. Those boundary modes and Wilson lines stretched between the two boundaries are carefully analyzed and are confirmed as physical variables in the phase space. Along with bulk fluctuation modes and topological modes, the partition function and entropy of all physical modes are evaluated via Euclidean path integral. It is shown that there are transitions between the dominance of different modes as we vary the temperature. The boundary fluctuation modes whose entropy is proportional to the volume dominate at high temperatures, and the boundary-area scaled boundary modes and Wilson lines are the more important at low temperatures. At super-low temperatures, when all the fluctuation modes die off, we see the topological modes whose entropy is the logarithm of the length scales of the system. The boundary modes and non-local modes should have their counterparts in a black hole system with similar boundary conditions, which might provide important hints for black hole physics.


Introduction
Gauge theories with non-trivial boundary conditions are important aspects of theoretical physics and can be used to understand lots of interesting physical phenomena [1][2][3][4][5].The would-be gauge degrees of freedom, which are no longer pure gauge, can become physical modes due to boundary conditions, which were suggested to explain the micro-states of black holes [4][5][6][7][8][9][10][11].Moreover, most theoretical physicists tend to believe that the boundary degrees of freedom are vital for a better comprehension of the quantum effects of gravity, for example in the path integral formulation [12] and in AdS/CFT [13,14].
As suggested in [4][5][6][7][8][9][10][11], it is important to carefully study the boundary would-be gauge degrees of freedom.However, simply saying the boundary pure gauge configurations λ bdy are physical due to the boundary conditions and counting the corresponding entropy of those modes does not bring us anything.A relatively proper way to deal with those modes is to introduce a boundary current J µ coupled with the boundary gauge fields in the action [15].Functional integration of the current in the path integral naturally introduces an effective action for the boundary would-be gauge modes, which is more or less proportional to (∂ µ λ bdy ) 2 .Note that the above procedure is not gauge invariant and the non-gauge invariance of the key ingredient of the story.The drawbacks of the above procedure are that boundary conditions are less transparent and the introduction of the boundary current seems to be artificial.
We are aiming to better understand the physics related to would-be gauge degrees of freedom from a set of nontrivial boundary conditions, where boundary modes and Wilson lines stretched between the two boundaries are allowed.More specifically, we study gauge fields living between two parallel plates with non-trivial boundary conditions, shown in (11), and carefully separate different parts of contributions in the presence of those boundaries.The boundary condition we are interested in is the one where we allow residual degrees of freedom for the component of A µ that is perpendicular to the boundaries to exist.The canonical formulation is carefully studied such that dynamical modes in the phase space (or Hilbert space) and the measure in the path integral is clear.The bulk fluctuation modes, boundary modes, bulk Wilson lines that stretched between two boundaries, and other topological modes should be considered as physical degrees of freedom in the current setup.The Wilson lines stretched between two boundaries are defined in equation (28), which captures the difference of boundary conditions on two boundaries.Then, we evaluate the thermal partition function and entropy via the Euclidean path integral.
As the temperature of the system varies from high temperatures to super-low temperatures, different modes dominate at different temperatures, and we can say that there are phase transitions between different modes.The bulk fluctuation modes always dominate at high temperatures, whose entropy is proportional to the volume.The entropies of the boundary modes and Wilson lines are both proportional to the area of the boundaries, which should be useful in understanding the Bekenstein-Hawking entropy on black hole backgrounds.In super-low temperature limits, all the fluctuation modes play a less important role, and we can see interesting competition between the constant modes and topological modes.Those behaviors in the super-low temperature limits are supposed to provide some hints for the extremal limit on black hole backgrounds.A careful analysis of the black hole case is devoted to further studies.Note that we are considering a more general set of boundary conditions compared to previous studies [5][6][7][8][9][10].Rather than only considering boundary modes on specific boundaries, the presence of the boundary-stretched Wilson lines makes the structure of the phase space richer, which can also be important for black hole entropy and building connection with black hole soft hairs.
The relevant research can be an important aspect of understanding the microscopic interpretation of the Bekenstein-Hawking entropy.Like in the brick wall model built by 't Hooft [16], there can be two boundaries: one near the horizon and the other at infinity.It was shown that black hole micro-states can be understood by studying the fluctuations (like a scalar field) on such backgrounds with two boundaries.The flat case we are going to study in this paper can be regarded as a toy model of the black hole case with gauge theories being considered, where effects due to the presence of two boundaries and non-local effects are the key ingredients of our model.We leave the study of the non-local effects on a black hole background for future research.Moreover, the study of bulk gauge theory with two boundaries might be helpful to the recent progress on wormhole geometry [17,18] and the factorization puzzle [19][20][21][22][23][24].We hope the non-local effects we study here can help us better understand the puzzles in AdS/CFT and quantum gravity.
The paper is organized as follows.In section 2, we discuss the basic setup and interesting boundary conditions.The section 3 is devoted to studying the canonical formulation of the theory, where we focus on the symplectic form and phase space.We also discuss the relationship between the canonical formulation and path integral, which leads to the Euclidean path integral part in section 4. We carefully analyze the contributions from different modes via path integral in section 4. In section 5, we demonstrate the transition between the dominance of different modes.Section 6 is the conclusion section.The appendix A provides more details of the calculations in the main text.This is the first paper in a series of papers that study the non-trivial boundary conditions of gauge theories.In the second paper of the series [25], we study a similar set of boundary conditions on black hole background and are trying to understand the difference between the microstructures of finite-temperature and extremal black holes.

Boundary conditions
This section studies the boundary conditions for a U(1) gauge theory living between flat parallel plates.The situation we are mainly interested in is shown in figure 1, where we have a Maxwell field theory living between two parallel boundaries on the left-and right-hand side.The original Maxwell theory has the action where M is the 4-dimensional manifold, and e 2 is a dimensionless coupling constant.The 4-dimensional box has coordinate system x µ = (x a , r) = (t, x 2 , x 3 , r), where r is the radius direction, and x a are the directions along the boundaries.To perform the finite temperature field theory calculations, we Wick rotate the time direction t → −iτ such that τ becomes the Euclidean time with a periodicity β, i.e. the inverse temperature of the spacetime.If all the gauge fields die off near the boundary, this is just the blackbody radiation with two polarisation degrees of freedom after gauge fixing.However, interesting phenomena start to show up when we release the boundary conditions, meaning that there are extra boundary degrees of freedom allowed due to the nontrivial boundary conditions.Let us suppose the boundaries shown in figure 1 are labeled by r = r α , with the left and right boundaries located at r (l) = 0 or r (r) = L. L is the distance between the two plates.Then, in order to have a well-defined Hilbert space and variation principle, there would be some constraints for the boundary conditions.To see those constraints, let us first look at the variation of the action (1), which can be written as The on-shell variation of the action can be written as For the boundaries with normal vector n µ ∂ µ = ∂ r , as shown in figure 1, the variation can be written as To have a well-defined variation principle without adding any Gibbons-Hawking-like term, we have two obvious choices: We will detailedly discuss (5) later.For now, let us look at the other boundary components, i.e.
Note that the variation principle has no constraint on the components shown in (6), and thus can be arbitrary at first sight.The physically motivated metallic Casimir boundary condition [26][27][28][29][30][31][32][33], that requires the magnetic field normal to the boundary B ⊥ ∂M and electric field tangential to the boundary E ∂M to vanish.This is equivalent to asking Essentially, we are interested in the boundary conditions that respect the physically motivated boundary condition (7), while allowing δA r ∂M to be arbitrary.Now, let us discuss the boundary conditions (5) needed for a well-defined variation principle: For the first choice, we can have Neumann-like boundary condition Combined with (7), we have F µν ∂M = 0.The situation doesn't seem super interesting, because the only possible allowed boundary configurations are the pure gauge configurations.For the second choice F ra ∂M can take arbitrary values.There can be interesting physics when we allow δA r ∂M = 0, .
This is the situation we are going to investigate in the current paper.
So the boundary conditions we are mainly interested in are the ones that consistent with where f (x a ) can have local dependence of x a .A a ∂M are fixed configurations on boundaries, which don't need to be summed over in path integral.Meanwhile, there is no constraint on A r at the boundary, and correspondingly, F ra ∂M can be arbitrary.The Hilbert space is well-defined with the fixed boundary configurations A a ∂M , and we need to sum over different boundary configurations of A r in path integral.
Let us see what are the configurations that respect (11).For the A a components, the boundary configurations are fixed and we can let those fields be zero at the boundaries Note that the above setting also respects the Casimir boundary condition (7).A r ∂M can take different configurations at the left and right boundaries.Besides, the boundary configurations can fluctuate and have arbitrary x a dependence.So we need to separate the bulk and boundary configurations carefully.A r can be separated as follows where φ(x a ) are the configurations that makes Âr r=0 = 0, i.e.A r (r, x a ) r=0 = φ(x a )/L.L is added mainly for the dimension-counting reason.To write the fields in a more concordant way, we can also use Âa = A a to denote the bulk configurations that vanish on boundaries.Note that with the above separation, Âr does not equal zero at the boundary r = L, and we can also further decompose Âr into two parts.The part that satisfies L 0 dr Âr = 0 will not be our main concern here, and the part that satisfies L 0 dr Âr = 0 will be an important ingredient in our later calculations.One can also decompose Âr into the part that vanishes on both sides and the part that captures the difference between the two boundaries.Integrating the part that vanishes on both boundaries from one boundary to the other gives out zero, and L 0 dr Âr more or less captures the difference of A r between the two boundaries.The reason why we don't choose this decomposition is that from the bulk point of view, we need bulk on-shell configurations corresponding to those boundary degrees of freedom.The corresponding bulk configurations need to be on-shell such that we have a good separation of different modes in the action.What's more, as we will see in the next section, L 0 dr Âr has a good physical interpretation.Note that as already motivated in the introduction, we allow non-gauge invariance in the system to capture the possible would-be gauge modes that are wildly interested in gravitational systems.It is okay to add the would-be gauge modes that respect the boundary condition into the system as far as there are confirmed as physical through canonical analysis.So, before we actually integrate all the degrees of freedom in the Euclidean path integral, let us first analyze the theory's canonical phase space and see which degrees of freedom are physical.Only the physical degrees of freedom are supposed to be integrated over in the path integral.This is what we are going to do in the next section.

Canonical formulation
We are going to work out the phase space and the equipped symplectic form in canonical formulation for the theory with non-trivial boundary conditions (11).The canonical formulation is always related to a Cauchy surface where the Hilbert space is defined.For the flat parallel plates case, a Cauchy surface is shown in Fig. 2. The canonical formulation can help us to better understand what are the physical degrees of freedom in the phase space.The phase space can be represented by Γ , which is an even-dimensional manifold with coordinates x I = {q i , p j }, where q i and p j are the canonical coordinates and momenta.For field theories, we have infinite-dimensional phase spaces Γ .After canonical quantization, the phase space can be turned into the Hilbert space of the theory.
For U(1) gauge theory with trivial boundary conditions, we can decompose the gauge fields into temporal and spatial directions and rewrite the gauge fields A µ as we have the Lagrangian density written in terms of V and A i as The corresponding conjugate momenta of fields V and A i can be written as where we have denoted Π V as the momentum for V and Π i as the momenta for A i .So, for trivial boundary conditions, the phase space we start with is Gauge fixing conditions help us further get rid of the unphysical degrees of freedom in the phase space and we end up with the standard two polarizations of photon in Maxwell's theory.Note that it is natural to use temporal gauge A t = 0 in the canonical formulation.However, there can be boundary subtleties in the phase space when we have nontrivial boundary conditions, like (11).What are those boundary subtleties?This can be answered by turning to the symplectic form of the theory and working out the Poisson bracket between the fields to make those boundary subtleties more explicit.Moreover, an explicit phase space and symplectic form can help us figure out the canonical variables and momentum, which should be integrated over in the path integral.The phase space is equipped with a closed, non-degenerate symplectic two-form Ω, which is defined as Ω IJ is invertible, and the inverse Ω IJ is defined by Ω IK Ω KJ = δ I J .Now equipped with the symplectic form, the classical Poisson bracket between functionals F and G can be defined as Quantum commutators can be obtained by the canonical quantization procedure.
The symplectic form of a field theory can be directly worked out on a chosen Cauchy surface [34].Let us consider a field theory with a Lagrange density L[Ψ ], where Ψ denotes an arbitrary collection of fields.Taking the variation of L, we have The equation of motion E = 0 kills the first term in (20).The (pre)-symplectic potential Θ[Ψ, δΨ ] is a D − 1 form and can be integrated over the chosen (D −1)-dimensional Cauchy surface.The symplectic current ω can be defined as where δ 1 and δ 2 can be regarded as variations with respect to two different transformations.Integrating the symplectic current ω over the Cauchy surface Σ, we finally get the symplectic form Ω written as Note that the choice of the (D−1)-form ω[Ψ, δ 1 Ψ, δ 2 Ψ ] also depends on the Cauchy surface.Specifying to the temporal surface Σ t with normal vector n t , we can work out each component of ω.Now, let us be specific to the U(1) gauge theory with non-trivial boundary conditions at hand.The variation of the action is similar to (2).The difference is Θ, as well as the symplectic form, is defined on Cauchy surface Σ, as shown in Fig. 2. Thus, the symplectic form Ω Σ can be worked out as Specifying on the chosen Cauchy surface Σ t with normal vector n µ ∂ µ = ∂ t , the above expression can be written in components as The symplectic form is essential in defining the Hamiltonian dynamics; only the fields equipped with a nontrivial symplectic form with their conjugate momentum can be regarded as dynamical variables in the phase space.We have decomposed the gauge fields into several parts as we discussed at the end of the section 2. We are going to see what are the symplectic partners of all the configurations in the symplectic form.Now one can put the variation of fields into the symplectic form (24), and rewrite it as We can separate the Âi part with other parts, and the symplectic form reads as The first term in (26) gives us the usual Poisson bracket of Maxwell's theory.Integrating over r in the second term gives out where we have used the boundary condition δ Ât ∂Σ = 0. Ω bdy can be regarded as the symplectic form of the dynamical variables due to the presence of boundary condition (11).
The above symplectic form tells us what are the extra physical degrees of freedom in the system, besides phase space (17).The non-local part in (27) can be denoted as a separate variable.Inspired by that, we define the quantity W as which will be called Wilson lines 1 .Note that field W captures the difference of A r on the two boundaries.From 1 The actual Wilson lines stretched between the two boundaries can be written as now on, we will extract W modes out of Âr such that Âr equals zero at both boundaries.The boundary symplectic form, defined on a codimension-2 surface, can be written as (30) The cross-term between W and φ can be canceled by refining the fields.For example, shifting φ → φ + iW , the above symplectic form can be written as With the above redefinition of fields, the overall symplectic form ( 26) can be expressed as where Πi denotes the conjugate momentum of Âi .Note that, beside the conjugate momentum defined in (16), we further define the conjugate momentum of W and φ as The Poisson brackets can be derived as As discussed at the beginning of this section, we need to add the degrees of freedom related to the boundary subtleties back to the phase space and perform path integration over those configurations in the path integral.Those boundary subtleties are the degrees of freedom related to the boundary configurations of A r .By the detailed symplectic form analysis in this subsection, we have found the zero modes φ(x a ) that take zero value for the longitudinal momentum along r direction and Wilson lines stretched between the two boundaries W have non-trivial symplectic partners and Poisson brackets.Those are the degrees with P denoting the path ordering.φ is the boundary configuration of Ar.The definition of the Wilson line can be exactly matched if we choose a different separation of degrees of freedom as discussed at the end of section 2.
of freedom needed to be added back.So the actual phase space should be Further gauge fixing conditions would help us to get rid of the bulk gauge redundancy of { Πi , Âi } such that we are only left with two bulk polarizations.This isn't our main concern in this paper.The new ingredients are φ and W , where φ is the zero modes along the r direction of A r , and W is the Wilson lines stretched between those two boundaries.
As discussed at the beginning of the paper, the modes that appear in the canonical formula should be included in the path integral method of calculating the partition function.It can easily be seen that the partition function Z =tr e −βH with the Hamiltonian H, can be written as a path integral over τ ∈ [0, β], with the compact Euclidean time τ = it (see e.g.[35] for more details).With the relation between Hamiltonian formula tr e −βH and Euclidean path integral DΨ e −S E being clear, we can apply all the canonical analysis to the path integral formalism.
In the next section, we are going to include all the physical degrees of freedom into the Euclidean path integral to evaluate their contribution to entropy.Note that the canonical analysis gives us some hints about what degrees of freedom are physical.Although we can use the canonical results as input and use the relation between the canonical formula and the Euclidean path integral to work out the partition function, we will study the path integral carefully.We may use different gauge fixing conditions if they are more convenient.

Euclidean path integral
After the canonical analysis in the previous section, it is clear that the dynamical modes are the bulk fluctuation modes Âµ , zero modes along r direction φ, and the Wilson lines stretched between the two boundaries W .Those fields are the ingredients that need to be included in the Euclidean path integral.Note that there should only be two bulk physical polarizations for the fields Âµ after gauge fixing.Specially caution is needed when handling gauge fixing, and we will deal with bulk gauge fixing conditions after the physics are clear to avoid gauging too much or too little.
As discussed in the previous subsection, the partition function can be written as a Euclidean path integral with the Euclidean action S E written as The above is the original formula, and we are going to do some massage according to the hints from the canonical analysis.We can separate the x a directions with r in the action Again, we are going to separate the gauge fields into different parts The Euclidean action can be written in terms of those modes as The second part of the above expression contains the interesting ingredients of our story and can be denoted as φ(x a ) is not a function of radius direction r, so the integral over r can just pass through φ, and gives out an extra L in the first term.Noticing that Âa equal zero at the boundaries ∂M, integrating over r in the above effective action gives out Denoting W (x a ) = i dr Âr , the original action can be rewritten as We have several remarks regarding the Euclidean path integral and the effective action (45): The original path integral ( 38) is a path integral over A µ , while we have different variables in action (45).There is no difference between Âa and A a .The integral over A r component can be divided into several different pieces.The zero modes φ and the Wilson lines W are the parts capturing the boundary configurations of A r .The bulk modes Âr that satisfy the conditions Âr r=0 = 0 and L 0 dr Âr = 0 should be regarded as the bulk contribution.We can always gauge fix Âr to zero, which doesn't kill any important physics.One of the main purposes of the canonical analysis in the previous section is to make clear the measure of different modes in the path integral.From the symplectic form (32) and Poisson brackets (34)(35)(36), the measure can be easily determined.
Putting the effective action (45) in the path integral (38), we can first work out the Gaussian integral over φ, which gives out det(∂ 2 ) −1/2 in the partition function.The above procedure also gives out an effective action for W , which is the action for a 3-dimensional massless scalar field.The determinant det(∂ 2 ) −1/2 getting from integrating φ can be rewritten as path integral.With all the above arguments, the effective action can be expressed as One can also check the path integrals with actions ( 45) and (46) give out the same result.So we will use (46) as the effective action in evaluating the partition function later on.As a direct analogy of the canonical analysis, the cross term between fields W and φ in (45) can also be canceled by shifting φ → φ + iW , which would give out the same result as (46).We have been super-careful about the gauge fixing such that we didn't gauge fix any interesting physics.
For example, we can gauge fix part of the bulk fields Âµ later, but we always need to make sure φ and W are not gauged away.There might be other interesting modes that are needed to be added back.Let us suppose we are dealing with a compact U(1) gauge theory.In the Euclidean background, the map between the background time circle τ ∼ τ + β and the compact gauge parameter allows us to include some topological modes for component A τ .The fundamental group of S 1 is Z.So the modes can be expressed as Those modes correspond to the large gauge transformation and might be physical.Those modes respect the boundary condition (11), but not (12).So, we do not include the modes (47) in the current calculation because of the more strict boundary condition (12).Now, let us evaluate the Euclidean path integral.For the first term corresponding to a Maxwell theory with vanishing boundary conditions, we can denote the partition function as Z Â.The field φ and W are 3-dimensional scalar fields living on a surface with coordinate x a , which will be regarded as the boundary contribution.We can separate the path integral as bulk and boundary parts The main task left is to evaluate the bulk and boundary partition functions.We are going to discuss those different modes for the remainder of this section, and evaluate the partition function and demonstrate the possible phase transitions in the next section.

Bulk fluctuation modes
First of all, let us evaluate the partition function for bulk fluctuation modes Z Â.We will use the Faddeev-Popov method [36] to evaluate the partition function, by inserting the following identity with gauge fixing condition G = ∂ µ Âµ − c(x).Following the standard gauge fixing procedure, in Feynman gauge, we eventually get where C and C are ghost fields.After gauge fixing, the final result is the partition function for two bosonic polarizations If we define the energy and momenta of the gauge fields as (ω, p r , p 2 , p 3 ), the logarithm of Z Â can be calculated by working out the determinantal operator One can further evaluated the partition function by taking different limits of the length scales in the theory.We are going to evaluate the logarithm of the partition function in section 5 when we discuss different temperature limits.Note that we can also use different gauge fixing conditions, like the axial gauge or temporal gauge.The gauge fixing condition does not make much difference for the fluctuation modes as far as we keep two physical polarizations in the final result.

Fluctuation modes of φ and W
Let us evaluate the partition function of fields φ and W here.The action for φ and W can be written as which is the action for two massless scalar fields living on the boundary.Denoting the area of the boundary as "Area", the fluctuation modes of field φ and W can be expanded as The coefficient is chosen such that φs and W s are dimensionless and thus the integrals over d φ in the path integral give out dimensionless quantities.With this mode expansion, the corresponding partition function can be expressed as the logarithm of which is Z F is the partition function for two 3-dimensional massless scalar fields.We can then calculate the free energy and entropy of those modes in different temperature limits, and compare it with the bulk fluctuation modes, which will be the task for the next section.

Other interesting modes
There are some other interesting topological modes of W .The Wilson lines stretched between the two boundaries can be denoted as Because it's always inside of an exponential function, L 0 dr Âr is compact with periodicity 2π.The requirement that the Wilson lines are single-valued allows us to include the elements of fundamental group S 1 .In the Euclidean background, the background time circle τ ∼ τ + β allows the field W to wind around the S 1 circle and have some winding modes 2πnτ /β.Now the field W has compact constant modes and novel winding modes which are interesting to deal with.The constant modes contribution of W can always be written as The winding mode contribution can be written as Z w equals 1 when the coefficient Area 2e 2 βL is very large, since the mode with n = 0 dominants.When the coefficient is very small, we can change the sum into a Gaussian integral.The overall partition function of φ and W is the product of constant modes Z 0 , winding modes Z w , and fluctuation modes Z F discussed previously.
As a summary of the different modes and corresponding partition functions we have got, we can express the overall partition function as We have two bulk polarizations in Z Â, two collections of fluctuation modes from W and φ in ln Z F , constant modes Z 0 and winding modes Z w .The logarithm of the overall partition function can be expressed as We will directly evaluate different parts of (59) in the next section.

Transition between different phases
In this section, we evaluate and compare the partition function shown in (59) in different temperature limits.The partition function contains contributions from the bulk fluctuation modes Âµ , the zero modes φ, and the Wilson lines W .The detailed calculations of the partition function are included in Appendix A, such that we don't drown in the tasteless details.The main content of this section is devoted to the discussion of different behaviors and phase transitions.
There are three different dimensional length scales in the theory, the inverse temperature β, the distance between the two boundaries L, and the length scale of the boundary √ Area.The coupling constant e 2 is dimensionless.We are going to compare the inverse temperature β with other length scales in the theory and call them different temperature limits.In the different temperature limits, we can study different behaviors of the partition function.We are mainly interested in the following three different temperature limits.
The so-called high-temperature limit is the limit when we have β L √ Area.β is the smallest length scale in the system.In this temperature limit, the bulk fluctuation modes Âµ should be the most important contribution.The second temperature limit we are interested in is the low-temperature limit, where we have L β √ Area.The distance between the two boundaries is way smaller than the inverse temperature β, and all the high-frequency modes along the r direction will be gapped.The zero modes φ and the Wilson lines W start to play the most important role in this limit.The last case is the super-low temperature limit, where we have L √ Area β.The temperature is super low, and all the fluctuation partition functions that proportional to the temperature are disappeared.The logarithm contributions of fields W shown in the previous section become the most important ones.
Let us discuss those three temperature limits separately.The qualitative behavior of the entropy is illustrated in figure 4. The overall entropy is a summation of 3. A sketch of the entropy of fields φ and W with varying temperature.At high temperatures, the entropy scales as Area × T 2 .For lower temperatures, the entropy scales as the logarithm of temperature and coupling constant.The second picture is an enlarged version of the low-temperature region.The red dashed line is an auxiliary line showing ln T .As can be seen from the figure, the entropy goes to zero in the super-low temperature limit because the contributions coming from zero modes and winding modes cancel each other.different contributions, while figure 4 illustrates the contributions from different modes.The solid red curves show the dominant contributions.More details of the calculations can be found in Appendix A, the main results and behaviors are discussed below.

Case I: High temperature limit
In the high-temperature limit, nothing is special, and we expect to see the usual result of black body radiation in a box because the bulk fluctuation modes, whose entropy is proportional to T 3 , are the most important contribution.The partition function Z Â and entropy S Â of bulk modes Âµ are shown below which is exactly the blackbody radiation result in a flat box.In order to avoid confusion with action S, we use S Â to denote the corresponding entropy.When the tempera- Fig. 4. A sketch of the entropy of the whole system in different temperature limits.The actual entropy is the sum of different contributions, and the red line demonstrates the dominant contribution.There are two transitions of the dominants shown in the figure.The bulk fluctuation modes always dominate in the high-temperature limit, and the entropy scales as the volume multiplied by temperature cubed.For lower temperature becomes, the area contribution starts to dominate.At super-low temperatures, the fluctuation contribution is not important anymore, and the only contribution is from the constant modes and winding modes of field W .A more clear curve of the entropy near the origin is shown in the second panel of figure 3.
ture is high, the contributions from φ and W proportional to the area of the boundary are small compared with the bulk radiation.So in the high-temperature limit, the dominant contribution always comes from the bulk fluctuation modes Âµ , which scales as the volume between the two boundaries multiplied by the temperature cubed.

Case II: Low temperature limit
For lower temperatures, when we have L β √ Area, the situation starts to change.In this temperature range, finite β means that 1/L is very big, and the energy needed to excite high-frequency modes along the r is super high.Thus, the modes along the r direction are gapped, and we are only left with zero modes along this direction.Even the zero modes of Âa are killed by the boundary conditions Âa ∂M = 0, so we get no contribution from Âa components.Fortunately, the zero modes of A r , i.e. φ, is a survivor.Moreover, W plays a similar role as φ.The partition function and entropy for φ and W can be obtained as The entropy of the thermal fluctuation modes along the boundary direction is proportional to the area of the plates times temperature squared.There are extra contributions from constant modes and winding modes if e 2 is not very large, which is proportional to the logarithm of temperature and the coupling constant.The logarithm contribution is mainly controlled by the coupling constant e 2 , which can surpass the area contribution for suitable e 2 .We will discuss the logarithm contribution later because those terms will be the most important contribution as the temperature is even lower.So in this temperature limit, the entropy of the system mainly comes from the fluctuation modes of φ and W , which is more or less proportional to the area of the plates times the temperature squared.
It is worth noticing that the situation here is similar to the Kuluza-Klein (K-K) reduction along the radius direction.The energy of the K-K tower is proportional to 1/R, where R is the length scales of extra dimensions.When R is very small, we can only see the zero modes of the K-K tower, and the effective theory is lower-dimensional.
Case III: Super-low temperature limit As the temperature becomes even lower, all the thermal fluctuation contributions proportional to the temperature will not survive.In the so-called super-low temperature limit, we have L•Area As can be seen from ( 65), the coupling constant 1/e 2 is a controller of constant modes.For weak coupling, where we have 1/e 2 1 such that the constant modes contribution of φ is canceled by the winding modes.However, in the strong coupling limit we can always see the contribution of the constant modes, which always scales as the logarithm of the coupling constant and temperature.The corresponding behavior of the entropy of constant modes and winding modes is shown in the second picture of figure 3.
As a summary of this section, let us qualitatively illustrate the basic behavior of the entropy corresponding to (60).The bulk fluctuation modes are the most important modes at super-high temperatures, whose entropy should be the blackbody radiation, i.e.T 3 × V .The new ingredient of our story is the modes due to the boundary condition.At high temperatures, the entropy of the fluctuation modes of φ and W is proportional to T 2 × Area.While as temperature goes lower and lower, the fluctuation modes play a less and less important role, and the zero modes and winding modes that contribute as the logarithm of the temperature and e 2 start to dominate.However, as the temperature becomes super low, the contribution from ln Z w cancels the one from zero modes, and the overall entropy goes to a constant.The entropy of the gauge theory with the given boundary condition is illustrated in figure 4.

Conclusion
We analyze the partition function of the U(1) gauge field living between two parallel boundaries with boundary condition (11) in this paper.The canonical analysis helps us understand what are the dynamical variables in the phase space (or Hilbert space).We also get the measure of different fields by working out the symplectic form of the theory with the given boundary conditions.As shown in figure 1, the radius coordinate is labeled by r and the transverse coordinates are x a .Besides the edge modes due to the boundary condition, there are also non-local modes due to the physics interplay between the two boundaries.Those modes are non-local effects, like the Wilson lines stretched between the boundaries, but behave like co-dimensionone fields which are pretty similar to the boundary edge modes.So, the physical modes of the theory at hand contain four different parts: bulk fluctuation modes Âµ , zero longitudinal momentum modes of A r which is φ(x a ), boundary stretched Wilson lines W (x a ), constant modes and winding modes.
Putting all of the above modes into the Euclidean path integral, we can work out the partition function of the theory that contains contributions from the four parts.The bulk fluctuation modes always play the dominant role at very high temperatures, whose entropy scales as the volume of the bulk multiplied by temperature curbed The modes arising because of the boundary conditions become more and more important as the temperature becomes lower.For lower temperatures, the ratio between the distance of the two boundaries L and inverse temperature β becomes small, and the zero modes φ and the Wilson lines W that behave like boundary scalar fields give out the dominant contributions.The entropy of those modes scales as As the temperature becomes super-low, no fluctuation mode can be seen, and we are left with some constant modes and topological modes contributions.The entropy of those modes is approximately the logarithm of the coupling constant and the temperature.The qualitative behavior of the entropy of different modes is shown in figures 3 and 4.
The flat parallel plates case is supposed to serve as a good toy model for the more general situation in curved spacetime.We would like to see if a similar phenomenon also shows up in the black hole background (or even for the wormhole background).The boundary modes and nonlocal modes should have their counterparts in a black hole system with similar boundary conditions, and the nonlocal modes might be understood as soft hair of the system.The phase transitions suggest that there might be similar phenomena for black holes.If not, the special properties of the gravitational system are then encoded in the difference.We leave the related issues for future research.
At the end of the paper, let us briefly comment on the first set of boundary conditions, shown in (8).It's not ridiculous to add all configurations that respect the boundary condition, including the boundary gauge modes.For the U(1) gauge theory, there are bulk on-shell configurations that have one-to-one correspondences with those boundary configurations, and the entropy of those modes can be counted.Moreover, the non-gauge invariance of the boundary condition we are mainly interested in this paper naturally captures that kind of physics.Those boundary pure gauge modes are soft modes because of the vanishing Hamiltonian for the would-be gauge modes.Moreover, it was suggested there are possible connections between those modes and soft hair degrees of freedom [15,[37][38][39][40][41][42][43]], Barnich's non-proper degrees of freedom [4,5,44,45], and edge modes [7][8][9][10][46][47][48][49][50][51][52].More concrete connections between those boundary effects due to boundary conditions are worth further understanding.
Note that the first part involving UV cutoff Λ is a constant in the free energy because the logarithm of the partition function is linear in β.Therefore, the entropy of those modes can be written as So when e 2 is large enough, the contribution of the constant modes ln Z 0 can be canceled.Nevertheless, this does not matter because the constant modes contribution of φ is always much smaller than the fluctuation modes contributions in this temperature limit.The statistical entropy of the fluctuation modes of φ and W can be computed as Compared with the volume contribution, all the area and logarithm contributions are not going to be important.As shown in equation ( 62), the most important contribution always comes from the bulk fluctuation modes Âµ , which scales as the volume times the temperature cubed.

A.2 Low temperature limit
High temperature is boring because we can only see the bulk fluctuation modes.As the temperature goes lower, when we have L β √ Area, interesting phenomena due to the boundary condition start to show up.
First, let us look at the bulk fluctuation modes ln Z Â.In this limit, the distance between the two plates is very small compared to the inverse temperature β.Assuming finite temperature, we have ω m = 2πm/β.Small L implies the high-frequency modes along the r direction are gapped, and we would only see zero modes along the r direction.Âµ vanish on the boundary, so the zero modes of Âµ along the r direction is killed by the boundary conditions.The only surviving zero modes are the zero modes of A r namely φ, which will be discussed separately.In the low-temperature limit and also in the super-low temperature limit, we will never see any contribution from bulk fluctuation modes anymore.So we can conclude that the entropy from the bulk fluctuation modes is S Â = 0 . (85) However, the zero modes φ survived and the partition functions for φ and W are not changed.We still have ln Z Area βL 1, the contribution of the constant modes ln Z 0 is canceled by the winding modes contribution, we only left with fluctuation modes contribution.The overall entropy is The entropy of the system now scales the area times the temperature squared.

A.3 Super-low temperature limit
As the temperature becomes even lower, we have L √ Area β, which is the low-temperature limit.In this temperature limit, not only the contribution from Âµ can be ignored, the fluctuation modes of φ and W are not important at all.As for the constant modes and winding modes, we have which cancels the constant modes and the overall entropy tends to be a small constant.

Fig. 1 .
Fig. 1.U(1) gauge theory living between two parallel boundaries.The orange surface is a Cauchy surface with constant time.

Fig. 2 .
Fig. 2. The bulk manifold M, the boundaries ∂M, Cauchy surface Σt, and the boundaries of the Cauchy surface ∂Σ are all illustrated here.

β 3 1 .
The bulk fluctuation modes are already frozen to death in the previous stage, and now it is the turn for the ones of φ and W .The logarithm contributions from the constant modes and winding modes can be written as ln Z

)β 2 2 n=0 = 1 ,e 2 βL 1 ,
For a similar reason, one can show that the partition function for φ(x a ) and W can be written as ln Z + ln Z w .(78)Thewinding modes contribution Z w shown above depends on the value of coupling constant e 2 .(80) thus ln Z w = 0.However, when e 2 is big enough such that 1 Area (81) the coefficient inside of the exponential function is very small, and we can change the sum into an integral.Thus we have Z w ≈ dn e

1 ,
we have ln Z w = 0.The overall entropy can be written as Area β 2 + ln Z w .Area βL 1 , we have ln Z w = 0, and the overall partition function can be written as