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

Gauge theories with non-trivial boundary conditions are important aspects of theoretical physics and can be used to understand lots of interesting physical phenomena Regge:1974zd ; Wadia:1976fa ; Gervais:1978kn ; Barnich:2019qex ; Barnich:2018zdg . 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 Carlip:1998wz ; Donnelly:2014fua ; Donnelly:2015hxa ; Donnelly:2016auv ; Barnich:2018zdg ; Barnich:2019qex ; Blommaert:2018oue ; Blommaert:2018rsf . 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 Gibbons:1976ue and in AdS/CFT Stanford:2017thb ; Mertens:2022ujr .

As suggested in Carlip:1998wz ; Donnelly:2014fua ; Donnelly:2015hxa ; Donnelly:2016auv ; Barnich:2018zdg ; Barnich:2019qex ; Blommaert:2018oue ; Blommaert:2018rsf , it is important to carefully study the boundary would-be gauge degrees of freedom. However, simply saying the boundary pure gauge configurations $\lambda_{\text{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^{\mu}$ coupled with the boundary gauge fields in the action Cheng:2020vzw . 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 $(\partial_{\mu}\lambda_{\text{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_{\mu}$ 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 Donnelly:2014fua ; Donnelly:2015hxa ; Donnelly:2016auv ; Barnich:2018zdg ; Blommaert:2018oue ; Blommaert:2018rsf . 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 Beken- stein-Hawking entropy. Like in the brick wall model built by ’t Hooft tHooft:1984kcu , 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 Penington:2019kki ; Almheiri:2019qdq and the factorization puzzle Harlow:2019yfa ; Blommaert2021 ; Saad:2021rcu ; Saad:2021uzi ; An:2023dmo ; Cheng:2022nra . 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 Cheng2302 , 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.

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

$$S=-\frac{1}{4e^{2}}\int_{\mathcal{M}}d^{4}x~{}F^{\mu\nu}F_{\mu\nu}\,,$$

(1)

where $\mathcal{M}$ is the 4-dimensional manifold, and $e^{2}$ is a dimensionless coupling constant. The 4-dimensional box has coordinate system $x^{\mu}=(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\to-i\tau$ such that $\tau$ becomes the Euclidean time with a periodicity $\beta$, 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_{\alpha}$, 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

$$\delta S=\frac{1}{e^{2}}\int_{\mathcal{M}}d\tau d^{3}x~{}\partial_{\mu}F^{\mu\nu}\delta A_{\nu}-\frac{1}{e^{2}}\int_{\partial\mathcal{M}}d\tau d^{2}x~{}n_{\mu}F^{\mu\nu}\delta A_{\nu}\,.$$

(2)

The on-shell variation of the action can be written as

$$\delta S=-\frac{1}{e^{2}}\int_{\partial\mathcal{M}}d\tau d^{2}x~{}n_{\mu}F^{\mu\nu}\delta A_{\nu}\,.$$

(3)

For the boundaries with normal vector $n^{\mu}\partial_{\mu}=\partial_{r}$, as shown in figure 1, the variation can be written as

$$\delta S=-\frac{1}{e^{2}}\int_{\partial\mathcal{M}}d\tau d^{2}x~{}F^{ra}\delta A_{a}\,.$$

(4)

To have a well-defined variation principle without adding any Gibbons-Hawking-like term, we have two obvious choices:

$$F^{ra}\big{|}_{\partial\mathcal{M}}=0\,,~{}~{}~{}\text{or}~{}~{}~{}~{}\delta A_{a}\big{|}_{\partial\mathcal{M}}=0\,.$$

(5)

We will detailedly discuss (5) later. For now, let us look at the other boundary components, i.e.

$$n_{\mu}{}^{*}F^{\mu\nu}\big{|}_{\partial\mathcal{M}}\,,~{}~{}~{}~{}~{}~{}\delta A_{r}\big{|}_{\partial\mathcal{M}}\,.$$

(6)

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 Milton2001 ; Bordag2009 ; Jaffe2005 ; Chernodub:2017gwe ; Chernodub:2017mhi ; Chernodub:2018pmt ; Chernodub:2022izt ; Chernodub:2023dok , that requires the magnetic field normal to the boundary $\textbf{B}_{\perp}\big{|}_{\partial\mathcal{M}}$ and electric field tangential to the boundary $\textbf{E}_{\parallel}\big{|}_{\partial\mathcal{M}}$ to vanish. This is equivalent to asking

$$n_{\mu}{}^{*}F^{\mu\nu}\big{|}_{\partial\mathcal{M}}=0\,.$$

(7)

Essentially, we are interested in the boundary conditions that respect the physically motivated boundary condition (7), while allowing $\delta A_{r}\big{|}_{\partial\mathcal{M}}$ to be arbitrary.

Now, let us discuss the boundary conditions (5) needed for a well-defined variation principle:

• $\diamond$

  For the first choice, we can have Neumann-like boundary condition

  $$F^{ra}\big{|}_{\partial\mathcal{M}}=0\,.$$

  (8)

  Combined with (7), we have $F^{\mu\nu}\big{|}_{\partial\mathcal{M}}=0$. The situation doesn’t seem super interesting, because the only possible allowed boundary configurations are the pure gauge configurations.

• $\diamond$

  For the second choice

  $$\delta A_{a}\big{|}_{\partial\mathcal{M}}=0,$$

  (9)

  $F^{ra}\big{|}_{\partial\mathcal{M}}$ can take arbitrary values. There can be interesting physics when we allow

  $$\delta A_{r}\big{|}_{\partial\mathcal{M}}\neq 0,.$$

  (10)

  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

$$\delta A_{a}\big{|}_{\partial\mathcal{M}}=0,~{}~{}~{}~{}~{}~{}~{}\delta A_{r}\big{|}_{\partial\mathcal{M}}=f(x^{a})\,.$$

(11)

where $f(x^{a})$ can have local dependence of $x^{a}$. $A_{a}\big{|}_{\partial\mathcal{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}\big{|}_{\partial\mathcal{M}}$ can be arbitrary. The Hilbert space is well-defined with the fixed boundary configurations $A_{a}\big{|}_{\partial\mathcal{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

$\displaystyle A_{a}\big{|}_{r=r_{\alpha}}=0\,.$

(12)

Note that the above setting also respects the Casimir boundary condition (7). $A_{r}\big{|}_{\partial\mathcal{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

$$A_{r}(x^{\mu})=\hat{A}_{r}(x^{\mu})+\frac{\phi(x^{a})}{L}\,,$$

(13)

where $\phi(x^{a})$ are the configurations that makes $\hat{A}_{r}\big{|}_{r=0}=0$, i.e. $A_{r}(r,x^{a})\Big{|}_{r=0}=\phi(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 $\hat{A}_{a}=A_{a}$ to denote the bulk configurations that vanish on boundaries. Note that with the above separation, $\hat{A}_{r}$ does not equal zero at the boundary $r=L$, and we can also further decompose $\hat{A}_{r}$ into two parts. The part that satisfies $\int_{0}^{L}dr~{}\hat{A}_{r}=0$ will not be our main concern here, and the part that satisfies $\int_{0}^{L}dr~{}\hat{A}_{r}\neq 0$ will be an important ingredient in our later calculations. One can also decompose $\hat{A}_{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 $\int_{0}^{L}dr~{}\hat{A}_{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, $\int_{0}^{L}dr~{}\hat{A}_{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.

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 $\Gamma$, 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 $\Gamma$. 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_{\mu}$ as

$$A_{\mu}=(-V,A_{i}),~{}~{}~{}~{}A^{\mu}=(V,A^{i})\,,$$

(14)

we have the Lagrangian density written in terms of $V$ and $A_{i}$ as

$$\mathcal{L}=\frac{1}{2e^{2}}(\dot{A}^{i}+\partial^{i}V)(\dot{A}_{i}+\partial_{i}V)-\frac{1}{2e^{2}}F^{ij}\partial_{i}A_{j}\,.$$

(15)

The corresponding conjugate momenta of fields $V$ and $A_{i}$ can be written as

$$\Pi_{V}=\frac{\partial\mathcal{L}}{\partial\dot{V}}\,~{}~{}~{}~{}\Pi^{i}=\frac{\partial\mathcal{L}}{\partial\dot{A}_{i}}=\frac{1}{e^{2}}(\dot{A}^{i}+\partial^{i}V)\,,$$

(16)

where we have denoted $\Pi_{V}$ as the momentum for $V$ and $\Pi^{i}$ as the momenta for $A_{i}$. So, for trivial boundary conditions, the phase space we start with is

$$\Gamma_{0}=\{V,A_{i},\Pi_{V},\Pi^{i}\}\,.$$

(17)

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 $\Omega$, which is defined as

$$\Omega=\frac{1}{2}\Omega_{IJ}\,\mathrm{d}x^{I}\wedge\,\mathrm{d}x^{J}\,.$$

(18)

$\Omega_{IJ}$ is invertible, and the inverse $\Omega^{IJ}$ is defined by $\Omega^{IK}\Omega_{KJ}=\delta^{I}_{J}$. Now equipped with the symplectic form, the classical Poisson bracket between functionals $F$ and $G$ can be defined as

$$\{F,G\}=\Omega^{IJ}\frac{\delta F}{\delta x^{I}}\frac{\delta G}{\delta x^{J}}\,.$$

(19)

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 Wald:1999wa . Let us consider a field theory with a Lagrange density $\mathcal{L}[\Psi]$, where $\Psi$ denotes an arbitrary collection of fields. Taking the variation of $\mathcal{L}$, we have

$$\delta\mathcal{L}=E\cdot\delta\Psi+\,\mathrm{d}\Theta\,.$$

(20)

The equation of motion $E=0$ kills the first term in (20). The (pre)-symplectic potential $\Theta[\Psi,\delta\Psi]$ is a $D-1$ form and can be integrated over the chosen $(D-1)$-dimensional Cauchy surface. The symplectic current $\omega$ can be defined as

$$\omega[\Psi,\delta_{1}\Psi,\delta_{2}\Psi]=\delta_{1}\Theta[\Psi,\delta_{2}\Psi]-\delta_{2}\Theta[\Psi,\delta_{1}\Psi]\,,$$

(21)

where $\delta_{1}$ and $\delta_{2}$ can be regarded as variations with respect to two different transformations. Integrating the symplectic current $\omega$ over the Cauchy surface $\Sigma$, we finally get the symplectic form $\Omega$ written as

$$\Omega[\Psi,\delta_{1}\Psi,\delta_{2}\Psi]=\int_{\Sigma}\omega[\Psi,\delta_{1}\Psi,\delta_{2}\Psi]\,.$$

(22)

Note that the choice of the $(D-1)$-form $\omega[\Psi,\delta_{1}\Psi,\delta_{2}\Psi]$ also depends on the Cauchy surface. Specifying to the temporal surface $\Sigma_{t}$ with normal vector $n^{t}$, we can work out each component of $\omega$.

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 $\Theta$, as well as the symplectic form, is defined on Cauchy surface $\Sigma$, as shown in Fig. 2. Thus, the symplectic form $\Omega_{\Sigma}$ can be worked out as

$\displaystyle\Omega_{\Sigma}=-\frac{1}{e^{2}}\int_{\Sigma}d^{3}x~{}n^{\mu}\delta F_{\mu\nu}\wedge\delta A^{\nu}\,.$

(23)

Specifying on the chosen Cauchy surface $\Sigma_{t}$ with normal vector $n^{\mu}\partial_{\mu}=\partial_{t}$, the above expression can be written in components as

$$\Omega_{\Sigma_{t}}=-\frac{1}{e^{2}}\int_{\Sigma_{t}}d^{3}x~{}\delta F_{ti}\wedge\delta A^{i}\,.$$

(24)

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

	$\displaystyle\Omega_{\Sigma_{t}}$	$\displaystyle=$	$\displaystyle-\frac{1}{e^{2}}\int_{\Sigma_{t}}d^{3}x~{}\Big{[}\delta(\hat{F}_{tr}+\frac{\dot{\phi}}{L})\wedge\delta(\hat{A}^{r}+\frac{\phi}{L})$		(25)
			$\displaystyle~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}+\delta\hat{F}_{t2}\wedge\delta\hat{A}^{2}+\delta\hat{F}_{t3}\wedge\delta\hat{A}^{3}\Big{]}\,.$

We can separate the $\hat{A}_{i}$ part with other parts, and the symplectic form reads as

	$\displaystyle\Omega_{\Sigma_{t}}$	$\displaystyle=$	$\displaystyle-\frac{1}{e^{2}}\int_{\Sigma_{t}}d^{3}x~{}\delta\hat{F}_{ti}\wedge\delta\hat{A}^{i}\,$
			$\displaystyle-\frac{1}{e^{2}}\int_{\Sigma_{t}}d^{2}xdr~{}\Big{[}\frac{\delta\dot{\phi}}{L}\wedge\delta\hat{A}^{r}+\frac{\delta\dot{\phi}}{L}\wedge\frac{\delta\phi}{L}+\delta\hat{F}_{tr}\wedge\frac{\delta\phi}{L}\Big{]}\,.$

The first term in (3) gives us the usual Poisson bracket of Maxwell’s theory. Integrating over $r$ in the second term gives out

	$\displaystyle\Omega_{\text{bdy}}$	$\displaystyle=$	$\displaystyle-\frac{1}{L\cdot e^{2}}\int d^{2}x~{}\Big{[}\delta\dot{\phi}\wedge\delta(\int_{0}^{L}dr\hat{A}^{r})+\delta\dot{\phi}\wedge\delta\phi$		(27)
			$\displaystyle~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}~{}+\delta(\int_{0}^{L}dr\dot{\hat{A}}_{r})\wedge\delta\phi\Big{]}\,,$

where we have used the boundary condition $\delta\hat{A}_{t}\big{|}_{\partial\Sigma}=0$. $\Omega_{\text{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

$$W=i\int_{0}^{L}dr~{}\hat{A}_{r}\,,$$

(28)

which will be called Wilson lines¹¹1 The actual Wilson lines stretched between the two boundaries can be written as $$\mathcal{W}\propto\mathcal{P}\exp\left[i(\int_{0}^{L}dr~{}\hat{A}_{r})+i\phi\right]\,,$$ (29) with $\mathcal{P}$ denoting the path ordering. $\phi$ is the boundary configuration of $A_{r}$. 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. . Note that field $W$ captures the difference of $A_{r}$ on the two boundaries. From now on, we will extract $W$ modes out of $\hat{A}_{r}$ such that $\hat{A}_{r}$ equals zero at both boundaries.

The boundary symplectic form, defined on a codimension-2 surface, can be written as

$$\Omega_{\text{bdy}}=-\frac{1}{e^{2}L}\int d^{2}x~{}\left[-i\delta\dot{\phi}\wedge\delta W-i\delta\dot{W}\wedge\delta\phi+\delta\dot{\phi}\wedge\delta\phi\right]\,.$$

(30)

The cross-term between $W$ and $\phi$ can be canceled by refining the fields. For example, shifting $\phi\to\phi+iW$, the above symplectic form can be written as

$$\Omega_{\text{bdy}}=-\frac{1}{e^{2}L}\int d^{2}x~{}\left[\delta\dot{\phi}\wedge\delta\phi+\delta\dot{W}\wedge\delta W\right]\,.$$

(31)

With the above redefinition of fields, the overall symplectic form (3) can be expressed as

	$\displaystyle\Omega_{\Sigma_{t}}$	$\displaystyle=$	$\displaystyle-\frac{1}{e^{2}}\int_{\Sigma_{t}}d^{3}x~{}\delta\hat{F}_{ti}\wedge\delta\hat{A}^{i}$		(32)
			$\displaystyle-\frac{1}{e^{2}L}\int d^{2}x~{}\left[\delta\dot{\phi}\wedge\delta\phi+\delta\dot{W}\wedge\delta W\right]\,$
		$\displaystyle=$	$\displaystyle-\frac{1}{e^{2}}\int_{\Sigma_{t}}d^{3}x~{}\delta\hat{\Pi}_{i}\wedge\delta\hat{A}^{i}$
			$\displaystyle-\frac{1}{e^{2}L}\int d^{2}x~{}\left[\delta\Pi_{W}\wedge\delta W+\delta\Pi_{\phi}\wedge\delta\phi\right]\,,$

where $\hat{\Pi}^{i}$ denotes the conjugate momentum of $\hat{A}_{i}$. Note that, beside the conjugate momentum defined in (16), we further define the conjugate momentum of $W$ and $\phi$ as

$$\Pi_{W}=\dot{W},~{}~{}~{}~{}\Pi_{\phi}=\dot{\phi}\,.$$

(33)

The Poisson brackets can be derived as

		$\displaystyle\frac{1}{e^{2}}[~{}\hat{\Pi}^{i}(r,x^{2},x^{3}),\hat{A}_{j}(r^{\prime},x^{\prime 2},x^{\prime 3})~{}]=i\delta^{i}_{j}~{}\delta(r-r^{\prime})~{}\delta^{2}(x-x^{\prime})\,,$		(34)
		$\displaystyle\frac{1}{e^{2}L}[~{}\Pi_{W}(x^{2},x^{3}),W(x^{\prime 2},x^{\prime 3})~{}]=i\delta^{2}(x-x^{\prime})\,,$		(35)
		$\displaystyle\frac{1}{e^{2}L}[~{}\Pi_{\phi}(x^{2},x^{3}),\phi(x^{\prime 2},x^{\prime 3})~{}]=i\delta^{2}(x-x^{\prime})\,.$		(36)

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 $\phi(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 of freedom needed to be added back. So the actual phase space should be

$$\Gamma=\left\{~{}\hat{\Pi}^{i},\hat{A}_{i},\Pi_{\phi},\phi,\Pi_{W},W~{}\right\}\,.$$

(37)

Further gauge fixing conditions would help us to get rid of the bulk gauge redundancy of $\{\hat{\Pi}^{i},\hat{A}_{i}\}$ such that we are only left with two bulk polarizations. This isn’t our main concern in this paper. The new ingredients are $\phi$ and $W$, where $\phi$ 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^{-\beta H}$ with the Hamiltonian $H$, can be written as a path integral over $\tau\in[0,\beta]$, with the compact Euclidean time $\tau=it$ (see e.g. Kapusta2011 for more details). With the relation between Hamiltonian formula tr $e^{-\beta H}$ and Euclidean path integral $\int\mathcal{D}\Psi~{}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.

After the canonical analysis in the previous section, it is clear that the dynamical modes are the bulk fluctuation modes $\hat{A}_{\mu}$, zero modes along $r$ direction $\phi$, 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 $\hat{A}_{\mu}$ 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

$$Z=\int\mathcal{D}A_{\mu}~{}e^{-S_{E}}\,,$$

(38)

with the Euclidean action $S_{E}$ written as

$$S_{E}=\frac{1}{4e^{2}}\int_{\mathcal{M}}d\tau d^{3}x~{}F^{\mu\nu}F_{\mu\nu}\,.$$

(39)

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

$\displaystyle S_{E}=\frac{1}{4e^{2}}\int_{\mathcal{M}}d\tau d^{3}x~{}F^{ab}F_{ab}+\frac{1}{2e^{2}}\int_{\mathcal{M}}d\tau d^{3}x~{}F^{ra}F_{ra}\,.$

(40)

Again, we are going to separate the gauge fields into different parts

$\displaystyle A_{a}=\hat{A}_{a}\,,~{}~{}~{}~{}~{}~{}A_{r}=\hat{A}_{r}+\frac{\phi(x^{a})}{L}\,.$

(41)

The Euclidean action can be written in terms of those modes as

	$\displaystyle S_{E}$	$\displaystyle=$	$\displaystyle\frac{1}{4e^{2}}\int_{\mathcal{M}}d\tau d^{3}x~{}\hat{F}^{\mu\nu}\hat{F}_{\mu\nu}$		(42)
			$\displaystyle+\frac{1}{2e^{2}}\int_{\mathcal{M}}d\tau d^{3}x~{}\left[-\frac{2}{L}~{}\hat{F}^{ra}\partial_{a}\phi+\frac{\partial^{a}\phi\partial_{a}\phi}{L^{2}}\right]\,.$

The second part of the above expression contains the interesting ingredients of our story and can be denoted as

$$S_{\text{bdy}}=\frac{1}{2e^{2}}\int_{\mathcal{M}}d\tau d^{2}xdr~{}\left[\frac{\partial^{a}\phi\partial_{a}\phi}{L^{2}}+\frac{2}{L}~{}\hat{F}^{ar}\partial_{a}\phi\right]\,.$$

(43)

$\phi(x^{a})$ is not a function of radius direction $r$, so the integral over $r$ can just pass through $\phi$, and gives out an extra $L$ in the first term. Noticing that $\hat{A}_{a}$ equal zero at the boundaries $\partial\mathcal{M}$, integrating over $r$ in the above effective action gives out

$\displaystyle S_{\text{bdy}}$

$\displaystyle=$

$\displaystyle\frac{1}{2e^{2}L}\int d\tau d^{2}x~{}\Big{[}{\partial^{a}\phi\partial_{a}\phi}-2i~{}\partial^{a}(i\int dr\hat{A}^{r})\partial_{a}\phi\Big{]}\,.$

Denoting $W(x^{a})=i\int dr\hat{A}_{r}$, the original action can be rewritten as

	$\displaystyle S_{E}$	$\displaystyle=$	$\displaystyle\frac{1}{4e^{2}}\int_{\mathcal{M}}d\tau d^{3}x~{}\hat{F}^{\mu\nu}\hat{F}_{\mu\nu}$		(45)
			$\displaystyle+\frac{1}{2e^{2}L}\int d\tau d^{2}x~{}\left[{\partial^{a}\phi\partial_{a}\phi}-2i~{}\partial^{a}W\partial_{a}\phi\right]\,.$

We have several remarks regarding the Euclidean path integral and the effective action (45):

• $\diamond$

  The original path integral (38) is a path integral over $A_{\mu}$, while we have different variables in action (45). There is no difference between $\hat{A}_{a}$ and $A_{a}$. The integral over $A_{r}$ component can be divided into several different pieces. The zero modes $\phi$ and the Wilson lines $W$ are the parts capturing the boundary configurations of $A_{r}$. The bulk modes $\hat{A}_{r}$ that satisfy the conditions $\hat{A}_{r}\big{|}_{r=0}=0$ and $\int_{0}^{L}dr\hat{A}_{r}=0$ should be regarded as the bulk contribution. We can always gauge fix $\hat{A}_{r}$ to zero, which doesn’t kill any important physics.

• $\diamond$

  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-36), the measure can be easily determined.

• $\diamond$

  Putting the effective action (45) in the path integral (38), we can first work out the Gaussian integral over $\phi$, which gives out $\det(\partial^{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(\partial^{2})^{-1/2}$ getting from integrating $\phi$ can be rewritten as path integral. With all the above arguments, the effective action can be expressed as

  	$\displaystyle S_{E}$	$\displaystyle=$	$\displaystyle\frac{1}{4e^{2}}\int_{\mathcal{M}}d\tau d^{3}x~{}\hat{F}^{\mu\nu}\hat{F}_{\mu\nu}$		(46)
  			$\displaystyle+\frac{1}{2e^{2}L}\int d\tau d^{2}x~{}\left[{\partial^{a}\phi\partial_{a}\phi}+\partial^{a}W\partial_{a}W\right]\,.$

  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 $\phi$ in (45) can also be canceled by shifting $\phi\to\phi+iW$, which would give out the same result as (46).

• $\diamond$

  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 $\hat{A}_{\mu}$ later, but we always need to make sure $\phi$ 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 $\tau\sim\tau+\beta$ and the compact gauge parameter allows us to include some topological modes for component $A_{\tau}$. The fundamental group of $S_{1}$ is Z. So the modes can be expressed as

  $$A_{\tau}\ni\frac{2\pi n}{\beta}\,,~{}~{}~{}~{}~{}n\in\textbf{Z}\,.$$

  (47)

  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_{\hat{A}}$. The field $\phi$ 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

	$\displaystyle Z$	$\displaystyle=$	$\displaystyle Z_{\hat{A}}\times\int\mathcal{D}\phi~{}\mathcal{D}W~{}$		(48)
		$\displaystyle\times$	$\displaystyle\exp\left[-\frac{1}{2e^{2}L}\int d\tau d^{2}x(\partial^{a}\phi\partial_{a}\phi+\partial^{a}W\partial_{a}W)\right]\,.$

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.

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 $\hat{A}_{\mu}$, the zero modes $\phi$, 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 $\beta$, the distance between the two boundaries $L$, and the length scale of the boundary $\sqrt{\text{Area}}$. The coupling constant $e^{2}$ is dimensionless. We are going to compare the inverse temperature $\beta$ 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.

• $\diamond$

  The so-called high-temperature limit is the limit when we have $\beta\ll L\ll\sqrt{\text{Area}}$. $\beta$ is the smallest length scale in the system. In this temperature limit, the bulk fluctuation modes $\hat{A}_{\mu}$ should be the most important contribution.

• $\diamond$

  The second temperature limit we are interested in is the low-temperature limit, where we have $L\ll\beta\ll\sqrt{\text{Area}}$. The distance between the two boundaries is way smaller than the inverse temperature $\beta$, and all the high-frequency modes along the $r$ direction will be gapped. The zero modes $\phi$ and the Wilson lines $W$ start to play the most important role in this limit.

• $\diamond$

  The last case is the super-low temperature limit, where we have $L\ll\sqrt{\text{Area}}\ll\beta$. 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 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.

For lower temperatures, when we have $L\ll\beta\ll\sqrt{\text{Area}}$, the situation starts to change. In this temperature range, finite $\beta$ 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 $\hat{A}_{a}$ are killed by the boundary conditions $\hat{A}_{a}\big{|}_{\partial\mathcal{M}}=0$, so we get no contribution from $\hat{A}_{a}$ components. Fortunately, the zero modes of $A_{r}$, i.e. $\phi$, is a survivor. Moreover, $W$ plays a similar role as $\phi$. The partition function and entropy for $\phi$ and $W$ can be obtained as

	$\displaystyle\ln Z$	$\displaystyle=$	$\displaystyle\frac{1}{2}\ln\frac{2\pi^{2}\text{Area}}{e^{2}\beta L}-\frac{1}{6\pi}{\beta\text{Area}}\cdot{\Lambda^{3}}+\frac{\zeta(3)}{\pi}\frac{\text{Area}}{\beta^{2}}+\ln Z_{w}\,,$
	$\displaystyle\mathcal{S}$	$\displaystyle=$	$\displaystyle(1-\beta\partial_{\beta})\ln Z$
		$\displaystyle=$	$\displaystyle\frac{3\zeta(3)}{\pi}{\text{Area}}\times T^{2}+\text{logarithm corrections}\,.$

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 $\phi$ 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.

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 $\frac{L\cdot\text{Area}}{\beta^{3}}\gg 1$. The bulk fluctuation modes are already frozen to death in the previous stage, and now it is the turn for the ones of $\phi$ and $W$. The logarithm contributions from the constant modes and winding modes can be written as

$\displaystyle\ln Z=\frac{1}{2}\ln\frac{2\pi^{2}\text{Area}}{e^{2}\beta L}+\ln Z_{w}\,.$

(64)

with

$\displaystyle\ln Z_{w}=\left\{\begin{matrix}0~{}\,;&&~{}~{}~{}~{}\frac{\text{Area}}{2e^{2}L\beta}\gg 1\\ -\frac{1}{2}\ln[\frac{\text{Area}}{\beta L}]+\ln e~{}.&&~{}~{}~{}~{}\frac{\text{Area}}{2e^{2}L\beta}\ll 1\end{matrix}\right.$

(65)

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}\ll 1$ such that

$$\frac{1}{e^{2}}\frac{\text{Area}}{\beta L}\ll 1\,,$$

(66)

the constant modes contribution of $\phi$ is canceled by the winding modes. However, in the strong coupling limit

$$\frac{1}{e^{2}}\frac{\text{Area}}{\beta L}\gg 1\,,$$

(67)

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}\times 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 $\phi$ and $W$ is proportional to $T^{2}\times\text{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.

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-dimension-one 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 $\hat{A}_{\mu}$, zero longitudinal momentum modes of $A_{r}$ which is $\phi(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

$$\mathcal{S}_{\hat{A}}\propto\text{Volume}\times T^{3}\,.$$

(68)

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 $\beta$ becomes small, and the zero modes $\phi$ and the Wilson lines $W$ that behave like boundary scalar fields give out the dominant contributions. The entropy of those modes scales as

$$\mathcal{S}_{\phi,W}\propto\text{Area}\times T^{2}\,.$$

(69)

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 non-local modes should have their counterparts in a black hole system with similar boundary conditions, and the non-local 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 Hawking:2016sgy ; He:2014laa ; Haco:2018ske ; Hawking:2016msc ; Strominger2018 ; Cheng:2020vzw ; Cheng:2022xyr ; Cheng:2022xgm , Barnich’s non-proper degrees of freedom Barnich:2018zdg ; Barnich:2019qex ; Alessio:2020lpk ; Aggarwal:2022rrp , and edge modes Donnelly:2015hxa ; Donnelly:2016auv ; Blommaert:2018oue ; Blommaert:2018rsf ; Balachandran:1994vi ; Kabat:1994vj ; Kabat:2010nm ; Seraj:2016jxi ; Geiller:2017xad ; Seraj:2017rzw ; Henneaux:2018gfi . More concrete connections between those boundary effects due to boundary conditions are worth further understanding.