## Entanglement production in the dynamical Casimir effect at parametric resonance

## Ivan Romualdo ∗

Departamento de F´ ısica -ICEx, Universidade Federal de Minas Gerais, CP 702, 30161-970, Belo Horizonte -MG, Brazil and Instituto de F´ ısica-Universidade de S˜o Paulo, CP 66318, 05315970- S˜o Paulo-SP, Brazil a a

## Lucas Hackl †

Max Planck Institute of Quantum Optics, Hans-Kopfermann-Str. 1, D-85748 Garching, Germany and Munich Center for Quantum Science and Technology, Schellingstraße 4, D-80799 M¨nchen, Germany u

Nelson Yokomizo ‡

Departamento de F´ ısica -ICEx, Universidade Federal de Minas Gerais, CP 702, 30161-970, Belo Horizonte -MG, Brazil

The particles produced from the vacuum in the dynamical Casimir effect are highly entangled. In order to quantify the correlations generated by the process of vacuum decay induced by moving mirrors, we study the entanglement evolution in the dynamical Casimir effect by computing the timedependent R´ enyi and von Neumann entanglement entropy analytically in arbitrary dimensions. We consider the system at parametric resonance, where the effect is enhanced. We find that, in (1 + 1) dimensions, the entropies grow logarithmically for large times, S A ( τ ) ∼ 1 2 log( τ ), while in higher dimensions ( n +1) the growth is linear, S A ( ) t ∼ λτ where λ can be identified with the Lyapunov exponent of a classical instability in the system. In (1 + 1) dimensions, strong interactions among field modes prevent the parametric resonance from manifesting as a Lyapunov instability, leading to a sublinear entropy growth associated with a constant rate of particle production in the resonant mode. Interestingly, the logarithmic growth comes with a pre-factor with 1 / 2 which cannot occur in time-periodic systems with finitely many degrees of freedom and is thus a special property of bosonic field theories.

## I. INTRODUCTION

The dynamical Casimir effect (DCE) describes the creation of particles from the vacuum by moving mirrors [1, 2]. Many of its basic features, as particle production and the dynamics of quantum fields with time-dependent boundary conditions, closely resemble analogous effects in sophisticated phenomena as the Hawking effect, the Unruh effect, or particle production in an expanding Universe. As a result, the DCE has been extensively investigated as a simple model to explore physical aspects of these more complicated systems in a more manageable context [3-7]. In addition, the DCE has been observed experimentally in a superconducting circuit [8], opening an avenue for the investigation of processes of particle production from the vacuum in tabletop experiments . 1

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00835v1-figure-81cb3556-260a-40f4-ba72-f9f3e26be545.jpg)

Because the DCE effects are usually small in laboratory conditions, it is interesting to analyze it in a condition of resonance [9-12]. In this work, we focus on the case where a mode of the field is at parametric resonance with the mirror oscillation, following [10]. The general setup of the DCE at parametric resonance is depicted in

∗ ivanromualdo@usp.br

† lucas.hackl@mpq.mpg.de

‡ yokomizo@fisica.ufmg.br

1 In the experiment [8], the idealized time-dependent boundary conditions imposed by the moving mirror are implemented in practice by changing the properties of the superconducting circuit by the application of a magnetic flux.

FIG. 1. Basic scenario of the dynamical Casimir effect in (1+ 1) dimensions. We have a standing mirror at the origin ( x = 0) and a moving mirror ( x = L t ( )). The distance between the mirrors oscillates as L t ( ) = L 1 [1+ glyph[epsilon1] sin (2 πt/T )] with period T . Note that the mirror only moves during a finite time interval.

Fig. 1. A one dimensional cavity is bounded by two mirrors, one of which is static, while the other is allowed to move. If the moving mirror is set to oscillate harmonically with twice the frequency of some mode of the field, such a mode becomes resonant. This leads to a strong enhancement of the effect, with particles being continuously produced from the vacuum, providing a natural scenario for eventual further experimental explorations of the DCE [2].

Several works deal with the most immediate aspects of the phenomenon, such as the time evolution of particle production and the Casimir force generated between the

plates (for a review, see section II of [6]). However, a highly unexplored direction is the dynamics of the quantum correlations that emerge due to particle creation and the mixing of field modes induced by the moving mirror. In this work, we contribute to fill this gap by analitically determining the time evolution of entanglement measures for the DCE at parametric resonance.

We compute the R´ enyi and the entanglement entropy production from an initial vacuum state for the DCE at parametric resonance in arbitrary dimensions. We consider a massless scalar field in a hyperrectangular cavity where the field is required to vanish on the boundaries. For t ≤ 0, the cavity is kept fixed. For t &gt; 0, we allow the length of the cavity to vary in one direction for a finite amount of time by keeping one mirror fixed and setting the other to oscillate according to

$$L ( t ) = L _ { 1 } [ 1 + \epsilon \sin { ( 2 \pi t / T ) } ]. \quad \quad ( 1 ) \quad \_.$$

After a finite number of oscillations, we then let the mirror return to its starting position, and the cavity is kept fixed again. The frequency of the mirror oscillation is set to twice that of some cavity mode, giving rise to parametric resonance. We select the lowest energy mode to be the resonant, since this is the most likely scenario to be reproduced in a laboratory. The time evolution of particle production in such conditions was fully analyzed in [10], where, under reasonable approximations, the Bogoliubov transformation between the initial and final states was computed analytically for all times in both (1 + 1) and (3 + 1) dimensions.

In order to compute the entropies, we apply recently developed symplectic techniques for the description of the dynamics of entanglement for Gaussian states in linear bosonic systems [13, 14]. In this framework, the R´nyi e entropy is first computed as a volume associated with the subsystem in phase space. Then, according to a theorem derived in [14], the asymptotic behavior of the R´nyi e and entanglement entropies will converge. Moreover, in the presence of Lyapunov instabilities, the asymptotic production of entanglement entropy is completely determined by the Lyapunov exponents. We show that in (1 + 1) dimensions the entanglement entropy grows logarithmically for large times. In higher dimensions, the resonance is associated with a Lyapunov instability, leading to an asymptotic exponential particle production and linear growth of the entanglement entropy, with a production rate equal to the Lyapunov exponent of the system.

This paper is structured as follows. In Section II, we review the Hamiltonian formalism of the dynamical Casimir effect and the analytical solution for the case of parametric resonance obtained in [10]. In Section III, we discuss the necessary background on symplectic techniques for the dynamics of entanglement of Gaussian states. In Section IV, we compute analytically the time evolution of the R´nyi and entanglement entropies in the e dynamical Casimir effect at parametric resonance. Finally, we summarize and discuss our results in Section V.

## II. HAMILTONIAN FORMALISM FOR THE DYNAMICAL CASIMIR EFFECT

## A. The model

Let φ t, x ( ) be a neutral massless scalar field in ( d +1)dimensions described by the Lagrangian

$$L = - \frac { 1 } { 2 } \int d ^ { n } x \left ( \partial ^ { \mu } \phi \right ) \left ( \partial _ { \mu } \phi \right ), \text{ \quad \ \ } \text{(2)}$$

where we adopt the mostly plus convention for the Minkowski metric, η µν = diag( -, + , . . . , +). Requiring the action to be stationary, one finds that the field satisfies the Klein-Gordon equation

$$\frac { \partial ^ { 2 } \phi } { \partial t ^ { 2 } } - \nabla ^ { 2 } \phi = 0 \,.$$

The associated momentum is

$$\pi ( t, x ) = \partial _ { t } \phi ( t, x ) \,,$$

and the Hamiltonian is

$$H = \frac { 1 } { 2 } \int d ^ { n } x \, \left [ \pi ^ { 2 } + ( \nabla \phi ) ^ { 2 } \right ] \,. \quad \quad ( 5 )$$

We are interested in the situation where the field is confined to a finite spatial volume V ( ). t More specifically, we consider a field restricted to a time-dependent parallelepiped bounded in some direction x 1 to a finite interval x 1 ∈ [0 , L ( )], t where the first boundary is fixed and the second boundary is allowed to move arbitrarily. In the remaining directions, the boundaries are fixed, so that x i ∈ [0 , L i ], i = 2 , . . . , d , for x ∈ V ( ). t Dirichlet boundary conditions are imposed at ∂V ( ): t

$$\phi ( t, x ) | _ { \partial V ( t ) } = 0 \,.$$

In particular,

$$\phi ( t, x ^ { 1 } = 0 ) = \phi ( t, x ^ { 1 } = L ( t ) ) = 0 \,. \quad \ \ ( 7 )$$

If the field φ represents a component of the electromagnetic field, then these boundary conditions correspond to perfect mirrors, and the field is confined to a cavity bounded by perfectly reflecting plates. The boundary is assumed to move only for a finite period of time T , returning at the end to the starting position L 1 :

$$L ( t ) = L _ { 1 } \,, \quad \text{for $t\leq 0$ and $t\geq T$} \,. \quad \quad ( 8 )$$

In the quantum theory, the motion of the mirrors will lead to particle production, and under this restriction the number of produced particles remains finite.

Let S be the space of solutions of the Klein-Gordon equation (3) with boundary conditions (6). This space is equipped with an invariant bilinear form

$$( \phi, \psi ) = i \int _ { V ( t ) } d ^ { n } x \left [ \phi ^ { * } ( \partial _ { t } \psi ) - ( \partial _ { t } \phi ^ { * } ) \psi \right ], \quad ( 9 )$$

which is not positive definite. The space of solutions can be decomposed into a direct sum of subspaces of positive and negative energy, S = S + ⊕S -. The decomposition is required to be such that (i) positive-energy solutions have positive squared norm ( φ, φ ) and (ii) the complex conjugate of a positive-energy solution is a negative-energy solution, S -= S + . The decomposition is not unique, and distinct decompositions will correspond to distinct choices of vacuum in the quantum theory. The restriction of the invariant bilinear form to S + defines an inner product in this subspace, turning it into a Hilbert space H . We denote by { u n ( t, x ) } an orthonormal basis of H . The solutions u n are positive-energy normal modes for the Klein-Gordon equation.

The canonical quantization of the scalar field φ t, x ( ) and its associated momentum π t, x ( ) provides field operators satisfying the canonical comutation relations at equal times:

$$[ \hat { \phi } ( t, x ), \hat { \pi } ( t, x ^ { \prime } ) ] = i \delta ( x - x ^ { \prime } ) \,, \quad \ \ ( 1 0 ) \quad \ m a { \infty }$$

$$[ \hat { \phi } ( t, x ), \hat { \phi } ( t, x ^ { \prime } ) ] = [ \hat { \pi } ( t, x ), \hat { \pi } ( t, x ^ { \prime } ) ] = 0 \,, \quad ( 1 1 )$$

and the equation of motion (3). For each choice of normal modes u n , the field operator can be expanded as

$$\hat { \phi } ( t, x ) = \sum _ { n } [ u _ { n } ( t, x ) \hat { a } _ { n } + u _ { n } ^ { * } ( t, x ) \hat { a } _ { n } ^ { \dagger } ] \,, \quad ( 1 2 ) \quad \underset { \text{plet} \ell } { \underbrace { 1 } }$$

where the ˆ a n , a ˆ † n are creation and annihilation operators satisfying the usual commutation relations:

$$[ \hat { a } _ { m }, \hat { a } _ { n } ^ { \dagger } ] = \delta _ { m n } \,, \ \ [ \hat { a } _ { m }, \hat { a } _ { n } ] = [ \hat { a } _ { m } ^ { \dagger }, \hat { a } _ { n } ^ { \dagger } ] = 0 \,. \quad ( 1 3 )$$

The quantized canonical momentum operator is obtained by taking the time derivative of the field operator in the representation (12). Distinct choices of normal modes correspond to distinct representations of the quantum field. The vacuum state | 0 〉 is the state annihilated by all ˆ a n 's and depends on the choice of H = S + (although not on the specific basis { u n } chosen for H ).

We introduce two special representations, called the in -and out -representations. For t &lt; 0, the Hamiltonian is time-independent, and one can find an orthonormal basis of positive energy in -solutions for which the timedependence is factored out in this region:

$$u _ { n } ^ { ( i n ) } ( t, x ) = \psi ^ { n } ( t, x ) \, \in \, e ^ { - i \omega _ { n } t } \,, \ \text{ for } t \leq 0 \,. \quad ( 1 4 ) \, \quad \text{ trans}$$

Similarly, one can find an orthonormal basis of positive energy out -solutions satisfying

$$u _ { k } ^ { ( o u t ) } ( t, x ) = \psi _ { k } ( t, x ) \in e ^ { - i \omega _ { k } t } \,, \quad \text{for $t\geq T$} \,. \quad ( 1 5 )$$

Both sets B in = { ψ , ψ n n ∗ } and B out = { ψ , ψ k ∗ k } are linear bases of S , but define distinct decompositions into positive and negative energy subspaces. The field operator can be expanded in either set of normal modes:

$$\dot { \phi } ( t, x ) & = \sum _ { k } [ \psi _ { k } ( t, x ) \hat { a } _ { k } + \psi _ { k } ^ { * } ( t, x ) \hat { a } _ { k } ^ { \dagger } ] \,, \quad ( 1 6 ) \\ & = \sum _ { n } [ \psi ^ { n } ( t, x ) \hat { b } _ { n } + \psi ^ { n * } ( t, x ) \hat { b } _ { n } ^ { \dagger } ] \,, \quad ( 1 7 )$$

where the ˆ b n , b ˆ † n are annihilation and creation operators in the in -representation, and the ˆ a , a k ˆ † k are annihilation and creation operators in the out -representation. Expanding the in -modes in the basis of out -modes,

$$\psi ^ { n } = \sum _ { k } ( \alpha _ { n k } \psi _ { k } + \beta _ { n k } \psi _ { k } ^ { * } ) \,, \quad \quad ( 1 8 )$$

and substituting into Eq. (17), we find

$$\hat { a } _ { k } = \sum _ { n } ( \alpha _ { n k } \hat { b } _ { n } + \beta _ { n k } ^ { * } \hat { b } _ { n } ^ { \dagger } ) \,, \quad \quad ( 1 9 )$$

$$\hat { a } _ { k } ^ { \dagger } = \sum _ { n } ( \alpha _ { n k } ^ { * } \hat { b } _ { n } ^ { \dagger } + \beta _ { n k } \hat { b } _ { n } ) \,. \quad \quad ( 2 0 )$$

This is a Bogoliubov transformation describing the relation between the two representations. The constants α nk , β nk are the Bogoliubov coefficients of the transformation. The transformation is encoded in a symplectic matrix:

$$M = \begin{pmatrix} \alpha ^ { T } & \beta ^ { \dagger } \\ \beta ^ { T } & \alpha ^ { \dagger } \end{pmatrix} \,.$$

The evolution of any property of the quantum field due to the motion of the boundaries for t ∈ [0 , T ] is completely characterized by the Bogoliubov coefficients. In particular, the process of particle creation from the initial vacuum has a simple description. Let the initial state for t ≤ 0 be the in -vacuum annihilated by all ˆ b n :

$$\hat { b } _ { n } \left | 0, i n \right \rangle = 0 \,, \quad \forall n. \quad \quad \quad ( 2 2 )$$

Since we are in the Heisenberg representation, the state of the field remains unchanged. However, the number of particles N k in a mode k observed after the motion of the boundary is described by the out -operators:

$$\langle N _ { k } \rangle = \langle 0, i n | \, \hat { a } _ { k } ^ { \dagger } \hat { a } _ { k } \, | 0, i n \rangle = \sum _ { n } | \beta _ { n k } | ^ { 2 } \,. \quad ( 2 3 )$$

We see that the vacuum is unstable if the coefficients β nk do not vanish. In the next section we will discuss formulas for the R´nyi and entanglement entropies in the e out region in terms of the coefficients of the Bogoliubov transformation. Before that, we review the calculation of the Bogoliubov coefficients for the dynamical Casimir effect at parametric resonance performed in [10].

## B. Dynamical Casimir effect at resonance

Let n be a multi-index n = ( n , . . . , n 1 d ). The in -modes ψ n are solutions of the Klein-Gordon equation (3) that reduce to standing waves in the in -region:

$$\psi ^ { n } = \left [ \prod _ { i = 1 } ^ { d } \frac { 1 } { \sqrt { \pi n _ { i } } } \sin \left ( \frac { n _ { i } \pi } { L _ { i } } x ^ { i } \right ) \right ] e ^ { - i \omega _ { n } t } \,, \quad \text{for $t\leq 0$} \,,$$

where

$$\omega _ { n } = \pi \sqrt { \sum _ { i = 1 } ^ { d } \left ( \frac { n _ { i } } { L _ { i } } \right ) ^ { 2 } } \,. \quad \quad \ \ ( 2 5 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

While the mirror is moving, the modes evolve in a nontrivial way. At each time t , we can expand them in a basis of instantaneous Fourier modes that satisfy the imposed boundary conditions,

$$\psi ^ { n } = & \sum _ { k } \left [ \prod _ { i = 1 } ^ { d } \frac { 1 } { \sqrt { \pi n _ { i } } } \sin \left ( \frac { n _ { i } \pi } { L _ { i } } x ^ { i } \right ) \right ] \\ & \times \sqrt { \frac { L _ { 1 } } { L ( t ) } } \frac { 1 } { \sqrt { \pi n _ { 1 } } } \sin \left ( \frac { n _ { 1 } \pi } { L ( t ) } x ^ { 1 } \right ) Q _ { k } ^ { n } ( t ) \,. \quad ( 2 6 )$$

Substituting this expansion into the Klein-Gordon equation (3), one finds equations of motion for the Fourier amplitudes Q n k ( ), t which we will discuss in a moment. The Fourier amplitudes satisfy the initial conditions

$$Q _ { k } ^ { n } ( 0 ) = \delta _ { k n } \,, \quad \dot { Q } _ { k } ^ { n } ( 0 ) = - i \omega _ { n } \delta _ { k n } \,. \quad ( 2 7 )$$

In the out -region, the evolution is again trivial. The positive energy out -modes are defined as

$$\psi _ { k } = \left [ \prod _ { i = 1 } ^ { d } \frac { 1 } { \sqrt { \pi k _ { i } } } \sin \left ( \frac { k _ { i } \pi } { L _ { i } } x ^ { i } \right ) \right ] e ^ { - i \omega _ { k } t } \,, \quad \text{for $t\geq T$} \,, \quad \text{from} \\ \text{Tho in mod $a^{\prime}n$ odd} \, \text{in} \, \text{or} \, \text{out} \, \text{in} \, \text{or} \, \text{out} \, \text{in} \, \text{or} \, \text{out} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \,\text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{inor} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \text{in} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{or} \, \text{in} \, \text{in} \, \text{or} \, \text{in} \, \text{in} \, \text{or} \, \text{in} \, \text{in} \, \text{or} \, \text{in} \, \text{in} \, \text{or} \, \text{in} \, \text{in} \, \text{or} \, \text{in} \, \text{in} \, \text{or} \, \text{in} \, \text{in} \, \text{or} \, \text{in} \, \text{in} \, \text{or} \, \text{in} \, \text{in} \, \text{or} \, \text{in} \, \text{in} \,\text{or} \, \text{in} \, \text{in} \, \text{or} \, \text{in} \, \text{in} \, \text{or} \, \text{in} \, \text{in} \, \text{or} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \,\text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in}\, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text` \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \tau \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \tt{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in} \, \text{in}$$

The in -mode ψ n evolves into a superposition of out -modes of positive and negative energy, as written in Eq. (18), with constant Bogoliubov coefficients α, β determined by:

$$Q _ { k } ^ { n } ( t ) = \alpha _ { n k } e ^ { - i \omega _ { k } t } + \beta _ { n k } e ^ { i \omega _ { k } t } \,, \ \text{ for } t \geq T \,. \quad ( 2 9 )$$

Up to this point, we considered a generic motion L t ( ) of the boundary. Now let the moving boundary oscillate harmonically around its initial position. When the frequency of the boundary oscillation is twice that of some normal mode, the system is said to be at parametric resonance. Explicit formulas for the Bogoliubov coefficients can then be obtained under an approximation of slow variation of the coefficients, as first shown in [10], and we review the relevant results in the following subsections.

## 1. (1+1) dimensions

In one spatial dimension, all multi-indices reduce to integer numbers, and the expansion (26) of the in -modes in the instantaneous basis becomes

$$\psi ^ { n } = \sum _ { k } \sqrt { \frac { L _ { 1 } } { L ( t ) } } \frac { 1 } { \sqrt { \pi n } } \sin \left ( \frac { n \pi } { L ( t ) } x \right ) Q _ { k } ^ { n } ( t ) \,. \quad ( 3 0 )$$

The evolution equation for the Fourier amplitudes reads

$$( 2 5 ) \quad \ddot { Q } _ { k } ^ { n } + \Omega _ { k } ^ { 2 } ( t ) Q _ { k } ^ { n } & = 2 \lambda ( t ) \sum _ { j } g _ { k j } \dot { Q } _ { j } ^ { n } + \dot { \lambda } ( t ) \sum _ { j } g _ { k j } Q _ { j } ^ { n } \\ \underset { \underset { \text{assis} } { n o n - 1 } } & + \lambda ^ { 2 } ( t ) \sum _ { j, l } g _ { j k } g _ { j l } Q _ { l } ^ { n } \,,$$

where

$$\Omega _ { k } ( t ) = \frac { \pi k } { L ( t ) } \,, \quad \lambda ( t ) = \frac { \dot { L } ( t ) } { L ( t ) } \,, \quad \quad ( 3 2 )$$

and the coefficients g form an antisymmetric matrix with components glyph[negationslash]

$$g _ { j k } = ( - 1 ) ^ { k + j } \frac { 2 j k } { k ^ { 2 } - j ^ { 2 } } \,, \quad \text{for $j\neq k$} \,. \quad ( 3 3 )$$

Consider the case where the first normal mode is resonant, that is, let the moving boundary oscillate harmonically around its initial position

$$L ( t ) = L _ { 1 } \left [ 1 + \epsilon \sin ( 2 \omega _ { 1 } t ) \right ] \,, \text{ \quad \ \ } ( 3 4 )$$

with twice the frequency of the first normal mode,

$$\omega _ { 1 } = \pi / L _ { 1 } \,. \quad \quad \quad ( 3 5 )$$

The functions Ω ( ) k t and λ t ( ) can then be determined from Eqs. (34) and (32). The oscillation amplitude is assumed to be small compared to the initial width L 1 of the boundary, glyph[epsilon1] glyph[lessmuch] 1. To first order in glyph[epsilon1] , the parameters λ and ˙ λ in the evolution equation oscillate harmonically with frequency 2 ω 1 , and λ 2 ∼ glyph[epsilon1] 2 can be neglected.

One can look for solutions of the evolution equation (31) of the form

$$Q _ { k } ^ { n } ( t ) = \alpha _ { n k } ( t ) e ^ { - i \omega _ { k } t } + \beta _ { n k } ( t ) e ^ { i \omega _ { k } t } \,, \quad ( 3 6 )$$

where the Bogoliubov coefficients are now allowed to be time-dependent, reaching their asymptotic values α nk , β nk at t = T . The time-dependent coefficients are assumed to vary slowly in time: after substituting (36) into Eq. (31), terms proportional to ¨ α and ¨ β are neglected, and the Bogoliubov coefficients are considered approximately constant during one period of oscillation of any normal mode. Multiplying the resulting equation of motion by e iω k t or e -iω k t and averaging over intervals T k = 2 π/ω k , one obtains a set of coupled differential equations for α and β , respectively. These are solved in [10]. Let us gather the relevant results for our purposes.

The complete elliptic integrals of the first and second kind are defined, respectively, as

$$K ( \kappa ) = \int _ { 0 } ^ { \pi / 2 } d \alpha \frac { 1 } { \sqrt { 1 - \kappa ^ { 2 } \sin ^ { 2 } \alpha } } \, \quad \ \ ( 3 7 )$$

$$E ( \kappa ) = \int _ { 0 } ^ { \pi / 2 } d \alpha \sqrt { 1 - \kappa ^ { 2 } \sin ^ { 2 } \alpha } \,. \quad \quad ( 3 8 )$$

It is convenient to introduce the new time variable

$$\tau = \frac { 1 } { 2 } \epsilon \omega t \,, \quad \quad \quad ( 3 9 )$$

and define the quantities

$$\kappa = \sqrt { 1 - e ^ { - 8 \tau } } \,, \quad \tilde { \kappa } = \sqrt { 1 - \kappa ^ { 2 } } = e ^ { - 4 \tau } \,. \quad ( 4 0 ) \quad \text{and} \, \, \mathfrak { t }$$

The lowest Bogoliubov coefficients are then given by

$$\alpha _ { 1 1 } = \frac { 2 } { \pi } \frac { E ( \kappa ) + \tilde { \kappa } K ( \kappa ) } { \left [ 1 + \tilde { \kappa } \right ] } \,, \quad \quad \left ( 4 1 \right ) \, \begin{array} { c } \sigma _ { i } \\ \text{of } t \end{array} \right ]$$

$$\beta _ { 1 1 } = - \frac { 2 } { \pi } \frac { E ( \kappa ) - \tilde { \kappa } K ( \kappa ) } { 1 - \tilde { \kappa } } \,. \quad \quad ( 4 2 )$$

Coefficients with one larger odd index are obtained from the recurrence relations

$$\text{the recurrence relations} \\ \sqrt { 3 } \, \alpha _ { 3 1 } & = - \beta _ { 1 1 } - \dot { \alpha } _ { 1 1 } \,, \\ \sqrt { 3 } \, \beta _ { 3 1 } & = - \alpha _ { 1 1 } - \dot { \beta } _ { 1 1 } \,, \\ \sqrt { n ( n + 2 ) } \, \alpha _ { n + 2, 1 } & = \sqrt { n ( n - 2 ) } \, \alpha _ { n - 2, 1 } - \dot { \alpha } _ { n 1 } \,, \quad n \geq 3 \,, \quad \text{as d} \, \\ \sqrt { n ( n + 2 ) } \, \beta _ { n + 2, 1 } & = \sqrt { n ( n - 2 ) } \, \beta _ { n - 2, 1 } - \dot { \beta } _ { n 1 } \,, \quad n \geq 3 \,, \quad \text{d} \, > 1$$

where the dots represent derivatives with respect to τ , and the relations

$$\alpha _ { 1, 2 j + 1 } & = ( - 1 ) ^ { j } ( 2 j + 1 ) \alpha _ { 2 j + 1, 1 } \,, \\ \beta _ { 1, 2 j + 1 } & = ( - 1 ) ^ { j } ( 2 j + 1 ) \beta _ { 2 j + 1, 1 } \,.$$

All coefficients with some even index vanish. In order to compute the evolved state of the subsystem formed only by the resonant mode, it is enough to know the Bogoliubov coefficients with some index equal to 1, which are all determined by the relations above.

In our analysis of entanglement production, the following formulas for infinite sums of products of Bogoliubov coefficients will play an important role:

$$\sum _ { \underset { \sim } { n = 1 } } ^ { \infty } \alpha _ { n 1 } \dot { \alpha } _ { n 1 } = \sum _ { n = 1 } ^ { \infty } \beta _ { n 1 } \dot { \beta } _ { n 1 } = - \alpha _ { 1 1 } \beta _ { 1 1 } \,, \quad \quad ( 4 3 ) \quad \text{ with}$$

$$\sum _ { \substack { n = 1 \\ \_ \_ } } ^ { \dots - \cdot } \left ( \alpha _ { n 1 } \dot { \beta } _ { n 1 } + \dot { \alpha } _ { n 1 } \beta _ { n 1 } \right ) = - \left ( \alpha _ { 1 1 } ^ { 2 } + \beta _ { 1 1 } ^ { 2 } \right ) \,, \quad ( 4 4 )$$

$$\sum _ { k, n = 1 } ^ { \stackrel { \stackrel { \stackrel { \stackrel { \stackrel { \stackrel { \stackrel { \stackrel { \stackrel { \stackrel { \stackrel { \stackrel { \stackrel { \stackrel { \stackrel { \stackrel { \stackrel { \stackrel { \stackrel { \stackrel { \stackrel { \stackrel { \stackrel { \stackrel { \stackrel { \stackrel \stack { \stackrel { \stackrel { \stackrel { \stackrel { \stack { \stackrel { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \ \stack { \ \stack { \ \ \ \ \ \ \ \ \ \ } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } ) } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } \ } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } { } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } {. \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \ \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \\null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null\null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \ \` \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \` \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \null \` \` \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \. \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ & \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ | \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \  \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
 \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ } \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

$$2. \ \left ( d + 1 \right ) d i m e n s i o n s$$

In higher dimensions, the expansion of the in -modes in the instantaneous basis has the general form (26). Upon substitution of this expansion into the Klein-Gordon equation, one obtains equations of motion for the Fourier coefficients Q n k similar to (31). The difference is that j, k, n are now multi-indices, the frequencies of the normal modes are

$$\Omega _ { k } ( t ) = \pi \sqrt { \sum _ { i = 1 } ^ { d } \left ( \frac { k _ { i } } { L _ { i } ( t ) } \right ) ^ { 2 } } \,, \quad \quad ( 4 6 )$$

and the coefficients g jk have a more complicated form, which is not relevant for our purposes, except for the fact that they are still constant.

A mode r is set at resonance by letting the boundary oscillate with twice the time-independent frequency ω r of the mode. One may take, for instance,

$$L ( t ) = L _ { 1 } \left [ 1 - \epsilon \cos ( 2 \omega _ { r } t ) \right ] ) \text{ \quad \ \ } ( 4 7 )$$

with

$$\omega _ { r } = \pi \sqrt { \sum _ { i = 1 } ^ { d } \left ( \frac { k _ { i } } { L _ { i } } \right ) ^ { 2 } } \,,$$

as done in [10]. Note that, in contrast to the onedimensional case, the frequencies are not equidistant for d &gt; 1. To first order in glyph[epsilon1] , the parameters λ and ˙ λ in the evolution equation oscillate harmonically with frequency 2 ω r , and λ 2 can be neglected, as before.

One can look again for solutions of the form (36) under a slow-variation approximation. The resulting equations for the Q 's are quite different, however, as compared to the one-dimensional case. The reason for that is the following. When multiplied by e ± iω k t , to first order in glyph[epsilon1] , the right-hand-side of Eq. (31) becomes a sum of terms proportional to exp[ ( i ± ω j ± ω k ± 2 ω r )]. The integral over a period T k = 2 π/ω k of each such term vanishes unless the sum of frequencies in the exponential vanishes. But since they are not equally spaced, this never happens, and the equations for the distinct modes all decouple. One is left with an infinite set of equations for independent oscillators with time-dependent frequencies

$$\Omega _ { k } ( t ) \sim \omega _ { k } \left [ 1 + 2 \gamma \cos ( 2 \omega _ { r } t ) \right ] \,, \quad \quad ( 4 9 )$$

with

$$\gamma = \frac { \epsilon } { 2 } \frac { \pi ^ { 2 } ( k _ { 1 } / L _ { 1 } ) ^ { 2 } } { \omega _ { k } ^ { 2 } } \,. \quad \quad ( 5 0 )$$

glyph[negationslash]

The calculation of the Bogoliubov coefficients can be done independently for each mode. Put Q k = Q k k . Fourier amplitudes G n k with k = n vanish, since the modes are decoupled. We need to solve

$$\ddot { Q } _ { k } + \Omega _ { k } ^ { 2 } ( t ) Q _ { k } = 0 \,, \quad \quad \quad ( 5 1 )$$

the equation of a time-dependent oscillator with a harmonically oscillating frequency. One can look again for solutions of the form (36), which we now write as

$$Q _ { k } ( t ) = \alpha _ { k } ( t ) e ^ { - i \omega _ { k } t } + \beta _ { k } ( t ) e ^ { i \omega _ { k } t } \,, \quad ( 5 2 )$$

with α k = α kk and β k = β kk . Under a slow variation approximation, terms proportional to ¨ α k and ¨ β k are neglected. After averaging over the fast oscillations, one finds for the resonant mode

$$\dot { \alpha } _ { r } = - i \omega _ { r } \gamma \beta _ { r } \,, \quad \dot { \beta } _ { r } = i \omega _ { r } \gamma \alpha _ { r } \,, \quad \text{(5.3)}$$

glyph[negationslash]

while the evolution of the non-resonant modes is trivial in this approximation: α , β k k = constant, for k = . Initial r conditions corresponding to the in -vacuum are given by

$$\alpha _ { r } ( 0 ) = 1 \,, \quad \beta _ { r } ( 0 ) = 0 \,. \quad \quad ( 5 4 ) \quad \text{an}$$

Integrating Eq. (53) with these initial conditions, we obtain the desired Bogoliubov coefficients

$$\alpha _ { r } ( t ) & = \cosh ( \omega _ { r } \gamma t ) \,, \\ \beta _ { r } ( t ) & = i \sinh ( \omega _ { r } \gamma t ) \,. \quad \quad \quad ( 5 5 ) \quad W \epsilon$$

The time-dependent symplectic transformation associated with such Bogoliubov coefficients, which describes the evolution of the resonant mode, is given by

$$M _ { r } ( t ) = \begin{pmatrix} \cosh ( \omega _ { r } \gamma t ) & - i \sinh ( \omega _ { r } \gamma t ) \\ i \sinh ( \omega _ { r } \gamma t ) & \cosh ( \omega _ { r } \gamma t ) \end{pmatrix} \,, \quad ( 5 6 ) \quad \text{where}$$

and has the simple form

$$M _ { r } ( t ) = \exp ( t K _ { r } ) \,, \ \ K _ { r } = \begin{pmatrix} 0 & - i \omega _ { r } \gamma \\ i \omega _ { r } \gamma & 0 \end{pmatrix} \quad ( 5 7 ) \quad \text{to mi} \\ \text{state}$$

with a time-independent symplectic generator K r .

## III. ENTANGLEMENT DYNAMICS OF GAUSSIAN STATES

In this section, we review the basic properties of Gaussian states. We focus on systems with a finite number of bosonic degrees of freedom, which we can later take to infinity to recover a bosonic field theory. In particular, we review compact expressions for the entropy and the von Neumann entanglement entropy in terms of the covariance matrix of the Gaussian state. Our conventions closely follow [13-15], while other reviews include [16, 17].

## A. Gaussian states

Before considering the field theory case, we focus on a system with N degrees of freedom. Here, the classical phase space V glyph[similarequal] R 2 N is equipped with a symplectic form Ω ab , i.e. an antisymmetric and non-degenerate bilinear form. For the quantum theory, we choose the basis

$$\hat { \xi } ^ { a } \equiv ( \hat { b } _ { 1 }, \cdots, \hat { b } _ { N }, \hat { b } _ { 1 } ^ { \dagger }, \cdots, \hat { b } _ { N } ^ { \dagger } ) \, \text{ \quad \ } ( 5 8 )$$

of creation and annihilation operators. This basis also fixes the form of the symplectic form to get the correct commutation relations

$$[ \hat { \xi } ^ { a }, \hat { \xi } ^ { b } ] = i \Omega ^ { a b } \text{ \ with \ } \Omega \equiv - i \begin{pmatrix} 0 & 1 \\ - 1 & 0 \end{pmatrix} \,. \quad ( 5 9 ) \text{ \quad } A ^ { \text{amr} } _ { \text{well-} ] }$$

The basis is not Hermitian, which means that we can define the 2 N × 2 N matrix C , such that

$$\hat { \xi } ^ { \dagger a } = C ^ { a } _ { \ b } \hat { \xi } ^ { b } \ \text{ with } \ C \equiv \begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix} \,. \quad \ ( 6 0 )$$

Given a linear observable O = f ξ a ˆ , we need to require a f ∗ a C a b = f b for O to be a Hermitian operator.

We can now introduce a special class of states | ψ 〉 that are fully characterized by their displacement vector z a and their covariance matrix G ab , defined as

$$z ^ { a } = \left \langle \psi \right | \hat { \xi } ^ { a } \left | \psi \right \rangle \,,$$

$$G ^ { a b } = \langle \psi | \, \hat { \xi } ^ { a } \hat { \xi } ^ { b } + \hat { \xi } ^ { b } \hat { \xi } ^ { a } \, | \psi \rangle - 2 z ^ { a } z ^ { b } \,. \quad \ \ ( 6 2 )$$

We refer to a state | ψ 〉 as Gaussian state and label it as | G,z 〉 if its linear complex structure J satisfies the condition

$$J ^ { 2 } = - 1 \text{ \ with \ } J ^ { a } _ { \, b } = G ^ { a c } \Omega _ { c b } ^ { - 1 } \,, \text{ \quad \ } ( 6 3 )$$

where we introduced the inverse symplectic form Ω -1 ab with Ω ac Ω -1 cb = δ a b . The linear map J : V → V is called linear complex structure because it represents the imaginary unit on the classical phase space, i.e. it squares to minus identity. We can use directly J to encode the state | G,z 〉 as solution to the equation

$$\frac { 1 } { 2 } ( \delta ^ { a } _ { \, b } - \mathrm i J ^ { a } _ { \, b } ) ( \hat { \xi } ^ { b } - z ^ { b } ) \left | G, z \right \rangle = 0 \,. \quad \ \ ( 6 4 )$$

While complex structures have been used to describe vacua of quantum fields in curved spacetime [18, 19] for many years, only recently they became a useful tool in quantum information to parametrize both bosonic [1315, 20] and fermionic Gaussian states [21-24] in a unified manner. As an example, we can choose z a = 0 and the covariance matrix to be given by

$$\stackrel { \dots } { c o - } } _ { \text{ons} } \quad G _ { 0 } \equiv \begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix} \quad \Rightarrow \quad J _ { 0 } = G _ { 0 } \Omega ^ { - 1 } \equiv \i \begin{pmatrix} 1 & 0 \\ 0 & - 1 \end{pmatrix} \,, \quad ( 6 5 )$$

which gives rise to the operator-valued vector

$$\frac { 1 } { 2 } ( \delta ^ { a } _ { b } - \mathrm i J ^ { a } _ { b } ) ( \hat { \xi } ^ { b } - z ^ { b } ) \equiv ( \hat { b } _ { 1 }, \cdots, \hat { b } _ { N }, 0, \cdots, 0 ) \quad ( 6 6 )$$

that annihilates the state | G , 0 0 , 〉 i.e. we just confirmed that the Gaussian state | G , 0 0 〉 is nothing else than the vacuum associated to the annihilation operators ˆ . b i For other choices of G ab and z a , we will find that (64) gives rise to a Bogoliubov transformation of the form

$$\hat { a } _ { i } = \sum _ { j } ( \alpha _ { j i } \hat { b } _ { j } + \beta _ { j i } ^ { * } \hat { b } _ { j } ^ { \dagger } ) + z _ { i } \,, \quad \quad ( 6 7 )$$

where ˆ a i represent a new set of annihilation operators annihilating the state | G,z 〉 under consideration.

A defining property of bosonic Gaussian states is the well-known Wick's theorem, i.e. that higher order n -point

correlation functions can be efficiently computed from 1and 2-point correlation functions. To state Wick's theorem, it is useful to define the connected n -point correlation function of the state | G,z 〉 as

$$C _ { n } ^ { a _ { 1 } \cdots a _ { n } } = \langle G, z | \left ( \hat { \xi } ^ { a _ { 1 } } - z ^ { a _ { 1 } } \right ) \cdots \left ( \hat { \xi } ^ { a _ { n } } - z ^ { a _ { n } } \right ) | G, z \rangle \. \ ( 6 8 )$$

We can then use the connected 2-point function given by

$$C _ { 2 } ^ { a b } = \frac { 1 } { 2 } ( G ^ { a b } + \mathrm i \Omega ^ { a b } ) \,, \quad \quad ( 6 9 ) \quad \, ;$$

to state Wick's theorem as the equality

$$C _ { 2 n + 1 } ^ { a _ { 1 } \cdots a _ { 2 n + 1 } } = 0 \,, & & ( 7 0 ) \quad \text{In pr} \\ & \quad \text{own} \infty.$$

$$C _ { 2 n } ^ { a _ { 1 } \cdots a _ { 2 n } } & = \sum ( \text{all 2-contractions of $C_{2}$} ) \quad \begin{pmatrix} \exp \\ \text{tions} \\ & = C _ { 2 } ^ { a _ { 1 } a _ { 2 } } \cdots C _ { 2 } ^ { a _ { 2 n - 1 } a _ { 2 n } } + \cdots, \end{pmatrix} \quad ( 7 1 ) \quad \begin{pmatrix} \exp \\ \text{upses t} \\ \text{denc} \end{pmatrix}$$

i.e. the connected n -point correlation function for odd n vanishes, while for even n , it is given by a sum over products of C ab 2 .

For the rest of this paper, we will focus on Gaussian states | G,z 〉 with z a = 0 , i.e. they are centered at the origin of phase space. We will write | G 〉 = | G, 0 〉 to simplify notation.

## B. Quadratic Hamiltonians

The most general class of Hamiltonians that preserve the family of Gaussian states, i.e. for which e tH ˆ | G,z 〉 is again Gaussian, takes the form

$$\hat { H } = \frac { 1 } { 2 } h _ { a b } \hat { \xi } ^ { a } \hat { \xi } ^ { b } + f _ { a } \hat { \xi } ^ { a } \,. \quad \quad ( 7 2 ) \quad \underset { \hat { H } } { \text{$\upsilon$} }$$

If we only look at Gaussian states | G,z 〉 with z a = 0, we need to choose f a = 0 in our Hamiltonian to ensure that z a remains equal to zero under time evolution. We will therefore restrict ourselves to the class of quadratic Hamiltonians

$$\hat { H } = \frac { 1 } { 2 } h _ { a b } \hat { \xi } ^ { a } \hat { \xi } ^ { b } \,, \quad \quad \quad ( 7 3 ) \quad \text{ne}$$

where we can derive a condition on h ab to ensure that ˆ H is Hermitian. For ˆ ξ † a = C a b ˆ , we find the condition ξ b

$$h _ { a b } = ( C ^ { \intercal } ) _ { a } ^ { \, c } h _ { d c } C ^ { d } _ { \, b } \,. \quad \quad ( 7 4 ) \quad$$

Furthermore, we can require h ab to be symmetric, i.e. h ab = h ba , because the anti-symmetric part will only contribute a constant to the energy. Defining U t ( ) = e -i tH ˆ , we can use the well-known Baker-Campbell-Hausdorff identity to compute

$$U ^ { \dagger } ( t ) \hat { \xi } ^ { a } U ( t ) = M ( t ) ^ { a } _ { \, b } \hat { \xi } ^ { b } \,, \quad \quad ( 7 5 )$$

where the symplectic transformation is given by

$$M ( t ) = e ^ { t K } \,, \ \text{ with } K ^ { a } _ { \, b } = \Omega ^ { a c } h _ { c b } \,. \quad ( 7 6 ) \quad \text{$\chi$}$$

For the dynamical Casimir effects, we will need to consider quadratic Hamiltonians with explicit time dependence

$$\hat { H } ( t ) = \frac { 1 } { 2 } h ( t ) _ { a b } \hat { \xi } ^ { a } \hat { \xi } ^ { b } \,,$$

which leads to a time evolution determined by the timeordered exponentials

$$U ( t ) = \mathcal { T } e ^ { - i \int _ { 0 } ^ { t } \hat { H } ( t ^ { \prime } ) \, d t ^ { \prime } } \,,$$

$$M ( t ) = \mathcal { T } e ^ { \int _ { 0 } ^ { t } K ( t ^ { \prime } ) \, d t ^ { \prime } } \text{ with } K ( t ) ^ { a } _ { \, b } = \Omega ^ { a c } h ( t ) _ { c b }.$$

In practice, we will not compute M t ( ) as a time-ordered exponential, but rather by solving the underlying equations of motion explicitly. Using (75), which also continues to be valid for Hamiltonian with explicit time dependence, we can compute the time evolution of the covariance matrix to be given by

$$G ( t ) ^ { a b } = M ( t ) ^ { a } _ { \ c } \, G ^ { c d } \, M ^ { \tau } ( t ) _ { d } ^ { \, b } \,. \quad \quad ( 8 0 )$$

For our class of Gaussian states, we therefore find

$$\left | G _ { t } \right \rangle = U ( t ) \left | G _ { 0 } \right \rangle = \left | M ( t ) G _ { 0 } M ^ { \intercal } ( t ) \right \rangle \,. \quad ( 8 1 )$$

A special class of time-dependent Hamiltonians are Floquet Hamiltonians, i.e. time-periodic Hamiltonians with period T satisfying

$$\hat { H } ( t + T ) = \hat { H } ( t ) \,. \quad \quad \quad ( 8 2 )$$

Such Hamiltonians can be described stroboscopically by only looking at the evolved states at t = nτ with n ∈ Z . For this, it is sufficient to compute U τ ( ) which allows us to define a time-independent effective Hamiltonian ˆ H eff = 1 τ log U τ ( ). For quadratic Hamiltonians ˆ ( H t ), we find ˆ H eff = 1 2 h eff ab ˆ ˆ ξ a ξ b to be also quadratic with

$$h _ { a b } ^ { \text{eff} } = \frac { 1 } { T } \Omega _ { a c } ^ { - 1 } \left ( \log M ( T ) \right ) ^ { c } _ { b }, \text{ \quad \ \ } ( 8 3 )$$

where M τ ( ) can be evaluated as the time-ordered exponential (79) or by integrating the equations of motion.

## C. Entanglement measures

Given a pure Gaussian state | G 〉 and a choice of subsystem A ⊂ V (inducing a tensor product decomposition H = H A ⊗ H B ), we can quantify the amount of quantum correlations between A and its complement B using different entanglement measures. We are particularly interested in the von Neumann entropy S A ( | G 〉 ) and the R´ enyi entropy R A ( | G 〉 ) of the reduced state

$$\rho _ { A } = \text{Tr} _ { B } \left | G \right \rangle \left \langle G \right | \,,$$

where we trace out the degrees of freedom in B .

One can show that the mixed state ρ A is also Gaussian, i.e. it continues to satisfy Wick's theorem and can

be efficiently parametrized by the restriction of G ab to the subsystem A . Mathematically speaking, we have a decomposition of the classical phase space V as a direct sum V = A ⊕ B , as well as of its dual V ∗ = A ∗ ⊕ B ∗ , such that the bilinear form G : V ∗ × V ∗ → R can be meaningfully restricted to

↦

$$G _ { A } \colon A ^ { * } \times A ^ { * } \colon ( a _ { 1 }, a _ { 2 } ) \mapsto G ( a _ { 1 }, a _ { 2 } ) \,. \quad ( 8 5 ) \quad ^ { \text{ let } v } _ { \ u m o t }$$

In practice, we use the matrix representation of G with respect to a basis ( ˆ ξ a A , ξ ˆ a B )

$$\hat { \xi } _ { A } ^ { a } \equiv ( \hat { b } _ { 1 } ^ { A }, \cdots, \hat { b } _ { N ^ { A } }, \hat { b } _ { 1 } ^ { A \dagger }, \cdots, \hat { b } _ { N _ { A } } ^ { A \dagger } ) \,, \quad ( 8 6 ) \quad \stackrel { \dots } { H a n }$$

$$\hat { \xi } _ { B } ^ { a } \equiv ( \hat { b } _ { 1 } ^ { B }, \cdots, \hat { b } _ { N _ { B } } ^ { B }, \hat { b } _ { 1 } ^ { B \dagger }, \cdots, \hat { b } _ { N _ { B } } ^ { B \dagger } ) \,, \quad ( 8 7 )$$

with ˆ ξ a = ( ˆ ξ a A , ξ ˆ a B ) and N = N A + N B . Note that this basis is complex, which implies that also G correspond to the analytically continued bilinear form on the complexified phase space V C and its dual V ∗ C .

The key result is that ρ A is fully characterized by G A , which implies that any function of ρ A can be written purely in terms of G A . The von Neumann entropy S A ( | G 〉 ) and the R´nyi entropy e R ( n ) A ( | G 〉 ) of order n are

$$S _ { A } ( | G \rangle ) = - \text{Tr} \rho _ { A } \log \rho _ { A } \,, \quad \ \ ( 8 8 ) \quad ^ { \ w h } _ { \colon \dots \, }$$

$$R _ { A } ^ { ( n ) } ( | G ) ) = \frac { 1 } { 1 - n } \text{Tr} ( \rho _ { A } ^ { n } ) \,, \quad \quad ( 8 9 ) \quad \ a c t \, \\ M ^ { \tau }$$

and we define R A ( | G 〉 ) := R (2) A ( | G 〉 ) as the R´ enyi entropy, dropping the reference to order 2 in this case. To express these functions in terms of G A , it is useful to define the restricted linear complex structure

$$( J _ { A } ) ^ { a } _ { \, b } = - ( G _ { A } ) ^ { a c } ( \Omega _ { A } ) ^ { - 1 } _ { c b } \,, \quad \ \ ( 9 0 ) \quad \ m o.$$

where Ω A is the restriction of Ω to the subsystem A . An important difference between the full linear complex structure J (describing a pure Gaussian state) and the restricted linear complex structure J A (describing in general a mixed Gaussian state) lies in the fact that -J 2 A is not necessarily the identity-instead, we have the inequality

$$- J _ { A } ^ { 2 } \geq 1 \,, \quad \quad \ \ ( 9 1 ) \quad \ \frac { \nu } { \text{p} }$$

i.e. J 2 A has negative eigenvalues whose magnitude is at least 1, but possibly larger. The magnitude of these eigenvalues or, put differently, the failure of J A to square to minus identity provides a measure of entanglement between A and B in the state | G 〉 . The formulas for S A ( | G 〉 ) and R A ( | G 〉 ) are then given by [13, 25, 26]

$$R _ { A } ( | G \rangle ) = \frac { 1 } { 2 } \log \det J _ { A } = \frac { 1 } { 2 } \log \frac { \det G _ { A } } { \det \Omega _ { A } } \,, \quad ( 9 2 ) \quad \frac { \dashv \dashv. } { R \acute { e n y } }$$

$$S _ { A } ( | G ) ) = \overset { \sim } { \text{Tr} } \left ( \frac { 1 + i J _ { A } } { 2 } \log \left | \frac { 1 + i J _ { A } } { 2 } \right | \right ) \,, \quad ( 9 3 ) \\ \text{where the absolute value } | \cdot | \text{ in the second equation is } \quad w h e r$$

where the absolute value | · | in the second equation is meant in terms of eigenvalues. Note that in the basis from (58) leads to the standard form of Ω (and Ω A ) from (59), which leads to det Ω A = ( -1) N A .

## D. Theorems on entanglement asymptotics

For classically unstable quadratic systems and Gaussian initial states, there exist several theorems that describe the asymptotic behavior of the entanglement entropy and the R´ enyi entropy for large times. For this, let us introduce the concept of Lyapunov exponents and unstable quadratic systems .

Definition 1 (Lyapunov exponents) . Given a quadratic Hamiltonian system H t ( ) = 1 2 h t ( ) ab ξ a ξ b with linear Hamiltonian flow

$$M ( t ) = \mathcal { T } \exp \int _ { 0 } ^ { t } d t ^ { \prime } K ( t ^ { \prime } ) \text{ \quad \ \ } ( 9 4 )$$

(where K t ( ) a b = Ω ac h t ( ) cb ) on the classical phase space V , we define the Lyapunov exponent λ glyph[lscript] of a linear observable glyph[lscript] ∈ V ∗ as the limit

$$\lambda _ { \ell } = \lim _ { t \to \infty } \log \frac { \| M ^ { \intercal } ( t ) \ell \| } { \| \ell \| } \,, \quad \quad ( 9 5 )$$

where the norm ‖·‖ can be induced by any positive definite inner product on the dual phase space. Note that the action on the dual phase space is given by the transpose M ᵀ ( ) t a b = M t ( ) b a .

Definition 2 (Unstable quadratic system) . A classical quadratic Hamiltonian H t ( ) = 1 2 h t ( ) ab ξ a ξ b is called unstable if there exists at least one linear observable glyph[lscript] ∈ V ∗ with positive Lyapunov exponent λ glyph[lscript] &gt; 0 . There are at most 2 N distinct Lyapunov exponents and we can always find a symplectic basis D L = ( glyph[lscript] 1 , · · · , glyph[lscript] 2 N ) with associated Lyapunov exponents λ i := λ glyph[lscript] i satisfying

$$\lambda _ { 1 } \geq \dots \geq \lambda _ { N } \geq 0 \geq \lambda _ { N + 1 } \geq \dots \geq \lambda _ { 2 N } \quad ( 9 6 )$$

and λ 2 N +1 -i = -λ i , i.e. all Lyapunov exponents come in conjugate pairs.

Proof. The proof can be found in definition 2 of appendix A.2 in [14].

With these definitions at hand, we can now review the key theorems that we will use in subsequent sections to study dynamics of R´ enyi entropy (of order 2) and von Neumann entanglement entropy.

Theorem 1 (R´ enyi entropy as phase space volume) . The R´ enyi entropy (of order 2 ) R A ( | G 〉 ) = R (2) A ( | G 〉 ) is given by

$$R _ { A } ( | G \rangle ) = \log \text{Vol} _ { G } \mathcal { V } _ { A } \,, \text{ \quad \ \ } ( 9 7 )$$

where V A ⊂ A ∗ ⊂ V ∗ represents an arbitrary parallelepiped of symplectic volume 1 .

$$P r o o f. \text{ The proof can be found in Sec. 6.2 of } [ 1 4 ]. \quad \Box$$

Theorem 2 (Asymptotic entanglement production) . We consider an unstable and regular 2 Hamiltonian system with N degrees of freedom and quadratic Hamiltonian ˆ ( H t ) = 1 2 h t ( ) ab ˆ ˆ ξ a ξ b , whose Lyapunov exponents are

$$\lambda _ { 1 } \geq \dots \geq \lambda _ { N } \geq 0 \geq \lambda _ { N + 1 } \geq \dots \geq \lambda _ { 2 N } \,. \quad ( 9 8 )$$

An initial Gaussian state | G 0 〉 will evolve into | G t 〉 , where the time-dependent covariance matrix is given by

$$G _ { t } = M ( t ) G _ { 0 } M ^ { \intercal } ( t ) \,. \quad \quad ( 9 9 ) \quad \text{cor}$$

A system decomposition V = A ⊕ B of the classical phasespace V (with dim A = 2 N A and dim B = 2 N B ) induces a tensor product decomposition H = H A ⊗ H B that we can use to compute entanglement entropy S A ( | G t 〉 ) . The asymptotics of the R´nyi entropy e R A ( | G t 〉 ) and the entanglement entropy S A ( | G t 〉 ) is given by

$$R _ { A } ( | G _ { t } \rangle ) \sim S _ { A } ( | G _ { t } \rangle ) \sim \Lambda _ { A } \, t \,, \quad \text{with} \quad \Lambda _ { A } = \sum _ { i = 1 } ^ { 2 N _ { A } } \lambda _ { i } \,, \quad \text{from} \quad \cdot \\ ( 1 0 0 ) \quad \cdot$$

where we assume the subsystem to be generic, i.e. above statement applies to all subsystems except a set of measure zero.

Proof. A detailed proof can be found in Sec. 2.1 and 2.2 of [14].

Next, we will compare the rate of entanglement entropy production in the dynamical Casimir effect with these predictions.

## IV. ENTANGLEMENT PRODUCTION AT PARAMETRIC RESONANCE

In this section, we compute the time-dependent R´nyi e and entanglement entropies at resonance. In (1 + 1)-D, the normal modes are strongly coupled by the evolution and become highly entangled. We compute the entropies explicitly for the subsystem associated with the resonant mode. In higher dimensions, the modes decouple under the assumed approximations, and the normal modes are not entangled. In this case, we consider a generic subsystem that intersects nontrivially the resonant mode.

## A. (1+1) dimensions

Let the scalar field be initially in its vacuum state | 0 , in 〉 . We work in the Heisenberg representation, so that the state of the system is constant and the time

2 The word 'regular' is used to exclude some pathological cases discussed in Appendix A.3 of [14].

evolution is encoded in the observables. Annihilation and creation operators ˆ a k ( ) t , a ˆ † k ( ) for a mode t k at each time t are defined in terms of the time-dependent Bogoliubov coefficients as in Eq. (19) :

$$\hat { a } _ { k } ( t ) = \sum _ { n } \left [ \alpha _ { n k } ( t ) \, \hat { b } _ { n } + \beta _ { n k } ^ { * } ( t ) \, \hat { b } _ { n } ^ { \dagger } \right ] \,. \quad \ ( 1 0 1 )$$

Such operators are associated with the instantaneous basis at t . It is convenient to use the variable τ = glyph[epsilon1]ω t/ 2 to compute the evolution of the entropies.

The covariance matrix at time τ is given by

$$G ^ { a b } ( \tau ) = \langle 0, i n | \left ( \hat { \xi } ^ { a } \hat { \xi } ^ { b } + \hat { \xi } ^ { b } \hat { \xi } ^ { a } \right ) | 0, i n \rangle \, \quad \ \ ( 1 0 2 )$$

where ˆ ξ a = (ˆ a k ( τ ) , a ˆ † k ( τ )) in the complex representation. The restriction of the covariance matrix to a subsystem A defines the reduced covariance matrix G A associated with it. The R´ enyi entropy R A ( | G 〉 ) can be computed from G A as

$$R _ { A } = \frac { 1 } { 2 } \log ( - \det G _ { A } ) \,, \quad \quad ( 1 0 3 )$$

where we used Eq. (92) and the fact that det(Ω A ) = -1 for a single degree of freedom ( N A = 1). Similarly, the entanglement entropy is given by

$$S _ { A } ( | G ) ) = \text{Tr} \left ( \frac { \mathbb { 1 } + \mathrm i J _ { A } } { 2 } \log \left | \frac { \mathbb { 1 } + \mathrm i J _ { A } } { 2 } \right | \right ) \quad ( 1 0 4 ) \\ \text{With } J _ { A } = - G _ { A } \, \omega _ { A } ^ { - 1 } \text{ as explained in } E q. \ ( 9 3 ). \ \text{We will}$$

With J A = -G A ω -1 A as explained in Eq. (93). We will now compute these quantities.

Expressing ˆ ξ a in terms of creation and annihilation in -operators ˆ b † m and ˆ b m using Eq. (101), and representing the covariance matrix as

$$G _ { A } ( \tau ) = \begin{pmatrix} G _ { A } ^ { 1 1 } & G _ { A } ^ { 1 2 } \\ G _ { A } ^ { 2 1 } & G _ { A } ^ { 2 2 } \end{pmatrix} \,, \quad \quad ( 1 0 5 )$$

we find for the matrix components

$$\nu \max \, \infty \, \varphi \, \text{max} \, \text{comp.} \\ G _ { A } ^ { 1 1 } & = \sum _ { n } 2 \alpha _ { n 1 } ( \tau ) \beta _ { n 1 } ^ { * } ( \tau ) \\ G _ { A } ^ { 1 2 } = G _ { A } ^ { 2 1 } & = \sum _ { n } \alpha _ { n 1 } ( \tau ) \alpha _ { n 1 } ^ { * } ( \tau ) + \beta _ { n 1 } ( \tau ) \beta _ { n 1 } ^ { * } ( \tau ) \\ G _ { A } ^ { 2 2 } & = \sum _ { n } 2 \beta _ { n 1 } ( \tau ) \alpha _ { n 1 } ^ { * } ( \tau ) \,. \quad ( 1 0 6 ) \\ \text{from now on a will omit the arguments of the time}$$

From now on, we will omit the arguments of the timedependent Bogoliubov coefficients for conciseness.

Since α nm and β nm are real for any m and n , we have

$$G _ { A } ^ { 1 1 } = G _ { A } ^ { 2 2 } = \sum _ { n } 2 \alpha _ { n 1 } \beta _ { n 1 } \, \text{ \quad \ \ } ( 1 0 7 )$$

$$G _ { A } ^ { 1 2 } = G _ { A } ^ { 2 1 } = \sum _ { n } \left ( \alpha _ { n 1 } ^ { 2 } + \beta _ { n 1 } ^ { 2 } \right ) \,. \quad \ \ ( 1 0 8 )$$

Differentiating with respect to τ , these sums can be evaluated using Eqs. (43)-(45) and then integrated analytically.

Let us start with the term G 11 A . Taking its derivative with respect to τ yields

$$\frac { d G _ { A } ^ { 1 1 } } { d \tau } = \sum _ { n } 2 \left ( \dot { \alpha } _ { n 1 } \beta _ { n 1 } + \alpha _ { n 1 } \dot { \beta } _ { n 1 } \right ) \,. \quad ( 1 0 9 ) \quad \text{$G$} \quad \text{where}$$

Then, according to Eq. (44), we have

$$\frac { d G _ { A } ^ { 1 1 } } { d \tau } = - 2 \left ( \alpha _ { 1 1 } ^ { 2 } + \beta _ { 1 1 } ^ { 2 } \right ) \,. \quad \quad ( 1 1 0 )$$

Similarly, taking the derivative of the term G 12 A , we find:

$$\frac { d G _ { A } ^ { 1 2 } } { d \tau } = \sum _ { n } 2 \left ( \alpha _ { n 1 } \dot { \alpha } _ { n 1 } + \beta _ { n 1 } \dot { \beta } _ { n 1 } \right ) \,. \quad ( 1 1 1 ) \quad \text{where} \quad \text{con} \quad \text{the}$$

From Eq. (43), we conclude that

$$\frac { d G _ { A } ^ { 1 2 } } { d \tau } = - 4 \alpha _ { 1 1 } \beta _ { 1 1 } \,. \quad \quad ( 1 1 2 ) \quad \text{ in}$$

Since G 11 A = G 22 A and G 12 A = G 21 A , the only remaining step is to integrate Eqs. (110) and (112).

The explicit form of the Bogoliubov coefficients α 11 and β 11 is presented in Eq. (41). Inserting these formulas in Eq. (110) we find

$$\frac { d G _ { A } ^ { 1 1 } } { d \tau } & = - \frac { 1 6 \left \{ [ E ( \kappa ) - \tilde { \kappa } ^ { 2 } K ( \kappa ) ] ^ { 2 } + \tilde { \kappa } ^ { 2 } [ E ( \kappa ) - K ( \kappa ) ] ^ { 2 } \right \$$

where the complete elliptic integrals of the first and second kind, K κ ( ) and E κ ( ), are defined in Eqs. (38) and (37), and κ and ˜ are given by Eq. (40). κ

From the definitions of κ and ˜, we can easily check κ that

Substituting this relation in (113), we can integrate it analytically (see [10]), obtaining:

$$G _ { A } ^ { 1 1 } & = - \frac { 4 } { \pi ^ { 2 } \kappa ^ { 2 } } \left [ E ( \kappa ) - \tilde { \kappa } ^ { 2 } K ( \kappa ) \right ] \left [ K ( \kappa ) - E ( \kappa ) \right ] + C _ { 1 }, \\ \quad \quad \quad$$

where C 1 is a constant of integration to be determined from the initial conditions.

In a similar way, substituting (41) in Eq. (112), we find

$$\frac { d G _ { A } ^ { 1 2 } } { d \tau } = \frac { 1 6 } { \pi ^ { 2 } \kappa ^ { 2 } } [ E ^ { 2 } ( \kappa ) - \tilde { \kappa } ^ { 2 } K ^ { 2 } ( \kappa ) ]. \quad ( 1 1 6 )$$

This equation can again be integrated [10], leading to

$$G _ { A } ^ { 1 2 } = \frac { 4 } { \pi ^ { 2 } } E ( \kappa ) K ( \kappa ) + C _ { 2 } \,, \quad \quad ( 1 1 7 )$$

where C 2 is a constant to be determined from the initial conditions.

To fix the integration constants, we need to evaluate the expressions in Eqs. (115) and (115) at τ = 0 and compare them with the components of the initial covariance matrix G A (0). But at τ = 0, the Bogoliubov transformation is trivial, α mn (0) = δ mn , β mn (0) = 0. Hence, the initial covariance matrix is given by

$$G _ { A } ( 0 ) = \left ( \begin{array} { c c } 0 & 1 \\ 1 & 0 \end{array} \right ). \quad \quad \quad ( 1 1 8 )$$

Now, in the small-time limit τ glyph[lessmuch] 1, the elliptic integrals can be expanded as [27]

$$K ( \kappa ) = \frac { \pi } { 2 } \left \{ 1 + \frac { 1 } { 4 } \kappa ^ { 2 } + \frac { 9 } { 6 4 } \kappa ^ { 4 } + \cdots \right \} \,, \quad ( 1 1 9 )$$

$$E ( \kappa ) = \frac { \pi } { 2 } \left \{ 1 - \frac { 1 } { 4 } \kappa ^ { 2 } - \frac { 3 } { 6 4 } \kappa ^ { 4 } - \cdots \right \} \,. \quad \ \ ( 1 2 0 )$$

In addition, in this limit we see that κ → 0 and ˜ κ → 1. Keeping only terms up to first order in τ , we find the matrix elements

$$G _ { A } ^ { 1 1 } = G _ { A } ^ { 2 2 } \sim - 2 \tau + C _ { 1 } \,, \quad \quad ( 1 2 1 )$$

$$G _ { A } ^ { 1 2 } = G _ { A } ^ { 2 1 } \sim 1 + C _ { 2 } \,. \quad \quad \quad ( 1 2 2 )$$

$$d \tau = \frac { \kappa d \kappa } { 4 \tilde { \kappa } ^ { 2 } } \,. \quad \quad \quad ( 1 1 4 ) \quad \ t h$$

Hence, the integration constants are C 1 = C 2 = 0, and the matrix G A ( τ ) assumes the form

$$G _ { A } ( \tau ) = \left ( \begin{array} { c } - \frac { 4 } { \pi ^ { 2 } \kappa ^ { 2 } } \left [ E ( \kappa ) - \tilde { \kappa } ^ { 2 } K ( \kappa ) \right ] \left [ K ( \kappa ) - E ( \kappa ) \right ] & \frac { 4 } { \pi ^ { 2 } \kappa ^ { 2 } } E ( \kappa ) K ( \kappa ) \\ \frac { 4 } { \pi ^ { 2 } } E ( \kappa ) K ( \kappa ) & - \frac { 4 } { \pi ^ { 2 } \kappa ^ { 2 } } \left [ E ( \kappa ) - \tilde { \kappa } ^ { 2 } K ( \kappa ) \right ] \left [ K ( \kappa ) - E ( \kappa ) \right ] \end{array} \right ) \,.$$

Furthermore, we can analyze the long-time behavior of the R´ enyi entropy. According to [27], when τ glyph[greatermuch] 1, the leading terms of the asymptotic expansion of the elliptic integrals are

$$K ( \kappa ) \sim \log \frac { 4 } { \tilde { \kappa } } + \frac { 1 } { 4 } \left ( \log \frac { 4 } { \tilde { \kappa } } - 1 \right ) \tilde { \kappa } ^ { 2 } + \cdots \quad ( 1 2 4 )$$

$$E ( \kappa ) \sim 1 + \frac { 1 } { 2 } \left ( \log \frac { 4 } { \tilde { \kappa } } - \frac { 1 } { 2 } \right ) \tilde { \kappa } ^ { 2 } + \cdots \quad \quad ( 1 2 5 )$$

Therefore, the covariance matrix becomes

$$G _ { A } ( \tau ) = \frac { 4 } { \pi ^ { 2 } } \left ( \begin{array} { c c } 1 - \zeta & \zeta \\ \zeta & 1 - \zeta \end{array} \right ), \quad \ \ ( 1 2 6 )$$

where ζ ≡ log 4 + 4 τ . We find that, in the asymptotic limit, the R´nyi entropy is e

$$R _ { A } ( \tau ) \sim \frac { 1 } { 2 } \log \left ( \frac { 1 6 ( 8 \tau + \log 1 6 - 1 ) } { \pi ^ { 4 } } \right ) \quad \ ( 1 2 7 ) \quad \ ^ { \lambda } \ 1 ( \tau )$$

$$\sim \frac { 1 } { 2 } \log \frac { 1 2 8 } { \pi ^ { 4 } } + \frac { 1 } { 2 } \log \tau \,.$$

The R´ enyi entropy converges fast to its asymptotic behavior, as shown in Fig. 2. In particular, we see in Fig. 3 that the first asymptotic expansion from Eq. (127) approaches the exact solution exponentially fast, while the leading behavior from Eq. (128) keeps a finite offset.

For a single degree of freedom, the relation between the R´ enyi entropy R A and the entanglement entropy S A is given by

$$S _ { A } ( \tau ) = s \left ( e ^ { R _ { A } ( \tau ) } \right ) \,, \quad \quad ( 1 2 9 ) \quad \text{wh}$$

with

$$s ( x ) & = \left ( \frac { x + 1 } { 2 } \right ) \log \left ( \frac { x + 1 } { 2 } \right ) - \left ( \frac { x - 1 } { 2 } \right ) \log \left ( \frac { x - 1 } { 2 } \right ) \cdot \quad \text{any sys} \\ \intertext { r u \colon c........................................................................................................................................................................................................ } }$$

This function behaves as

$$s ( x ) \sim \log x + ( 1 - \log 2 ) \ \text{ as } \ x \to \infty \,. \quad ( 1 3 1 ) \quad \text{linea}$$

Hence, for large x , we have

$$S _ { A } ( \tau ) \sim R _ { A } ( \tau ) + ( 1 - \log 2 ) \ \text{ as } \ x \tau \to \infty \,. \ \ ( 1 3 2 ) \quad \text{ence} \, \Omega$$

Combining Eqs. (128) and (132), we obtain the asymptotic behavior of the entanglement entropy given by

$$S _ { A } ( \tau ) \sim 1 + \frac { 1 } { 2 } \log \frac { 3 2 } { \pi ^ { 4 } } + \frac { 1 } { 2 } \log \tau \,, \quad ( 1 3 3 ) \quad \text{effe} \,.$$

FIG. 2. Comparison of the exact time-dependent R´ enyi entropy R A ( τ ) of a resonant mode in (1 + 1) dimensions with its asymptotic expansions given in Eqs. (127) and (128).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00835v1-figure-9115655e-cc32-45b7-94c8-366fdf4c9ce4.jpg)

FIG. 3. Difference between the exact R´nyi entropy e R A ( τ ) of a resonant mode in (1 + 1) dimensions and the asymptotic expansions given in Eqs. (127) and (128).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00835v1-figure-dee482ce-8e76-408c-9d5d-755ec26eeebc.jpg)

where we find the constant to be 1 + 1 log 32 /π 2 ≈ 0 44. .

2 In short, we found that in one spatial dimension, the asymptotic growth of the R´ enyi and entanglement entropies is logarithmic. As discussed in Section III B, to any system with a periodic quadratic Hamiltonian there is an associated effective time-independent Hamiltonian that describes its stroboscopic evolution. Moreover, if such an effective Hamiltonian displays a Lyapunov instability, then the asymptotic growth of the entropy is linear, according to the general results summarized in Section III D. Therefore, we conclude that for d = 1 the DCE displays no Lyapunov instability, despite the presence of parametric resonance. As stressed in [10], the case of d = 1 is special in the sense that, after the averaging of the fast oscillations, the field modes remain strongly coupled. The intermode interactions compete with the effect of parametric resonance, keeping the production of particles in the resonant mode linear, in opposition

FIG. 4. Comparison of the time-dependent entanglement entropy S A ( τ ) of a resonant mode in (1 + 1) dimensions with its asymptotic expansion given in Eq. (133).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00835v1-figure-1ed36072-5169-43c2-99e6-2cf2b6fb5a04.jpg)

to what happens in higher dimensions, where it grows exponentially. Our result show that this suppression of particle creation is reflected in the absence of a Lyapunov instability for the stroboscopic evolution.

The possible asymptotic behaviors of the entanglement entropy for systems with time-independent quadratic Hamiltonians and a finite number of degrees of freedom was fully analyzed in [26]. In the presence of instabilities, the asymptotic growth is linear, in agreement with the previous results obtained in [14], discussed in Section III D. In the absence of instabilities, the entanglement entropy either grows logarithmically, for so-called metastable systems, or is bounded. Hence, the DCE in d = 1 behaves qualitatively as a metastable system, but there is an important difference. In [26], it was proved that the prefactor multiplying the logarithmic function must be an integer, i.e. S A ( ) t ∼ n log t , n ∈ N . In contrast, we found S A ( ) t ∼ 1 / 2 log t , according to Eq. (133). This provides evidence that the extension of the results of [26] to systems with an infinite number of degrees of freedom is nontrivial, with new classes of asymptotic behavior arising in this limit. In fact, in preliminary numerical investigations, we found that, by truncating the number of degrees of freedom in the one dimensional DCE, the entanglement entropy becomes bounded, with an initial regime of logarithmic growth. As the size of the truncated system is increased, the regime of logarithmic growth becomes longer, consistently with the fact that for the infinite system such a regime has an infinite duration. This mechanism illustrates how a new class of logarithmic asymptotic growth can arise in the nontrivial limit of infinite number of degrees of freedom and may be investigated elsewhere.

## B. (2+1) dimensions

We now proceed to the case of 2 dimensions. The frequencies of the normal modes are determined by Eq. (46), where k is now a set of two integers. We let the moving mirror oscillate with twice the frequency ω r = ω 11 of the first mode r = (1 1), ,

$$L ( t ) = L _ { 1 } [ 1 - \epsilon \cos ( 2 \omega _ { r } t ) ], \quad \quad ( 1 3 4 ) \quad \text{sub}$$

so that the resonant mode is r . The length of the cavity along the x 2 -direction is kept fixed. The resonant frequency is

$$\omega _ { r } = \pi \sqrt { \frac { 1 } { L _ { 1 } ^ { 2 } } + \frac { 1 } { L _ { 2 } ^ { 2 } } } \,. \quad \quad ( 1 3 5 ) \quad \ m c _ { \text{for} }$$

The R´ enyi entropy can be computed from Eq. (92) as before. The main difference is that, as discussed in section II B 2, the evolution of the system is now much simpler. Under the slow-variation approximation, all normal modes decouple, and the evolution is nontrivial only for the resonant mode. The problem is thus reduced to that of a one dimensional parametric oscillator at resonance. The nontrivial Bogoliubov coefficients are given by Eq. (55)

$$\alpha _ { 1 1 } ( t ) & = \cosh ( \omega _ { r } \gamma t ), \\ \beta _ { 1 1 } ( t ) & = - i \sinh ( \omega _ { r } \gamma t ) \,, \quad \quad ( 1 3 6 )$$

where γ , defined in (50), now reads

$$\gamma = \frac { \epsilon \pi ^ { 2 } } { 2 L _ { 1 } ^ { 2 } \omega _ { r } ^ { 2 } } = \frac { \epsilon } { 2 } \frac { L _ { 2 } ^ { 2 } } { L _ { 1 } ^ { 2 } + L _ { 2 } ^ { 2 } } \,. \quad \quad ( 1 3 7 )$$

We start again from the covariance matrix in order to compute the R´ enyi entropy. The sums in the formulas (106) for the matrix components for each mode are now easily computed, since only the first term of each sum is nonzero. For the resonant mode r = (1 1), we obtain ,

$$G _ { r } ( t ) = \left ( \begin{array} { c c } G _ { r } ^ { 1 1 } & G _ { r } ^ { 1 2 } \\ G _ { r } ^ { 2 1 } & G _ { r } ^ { 2 2 } \end{array} \right ) \,, \quad \quad ( 1 3 8 )$$

with

$$G _ { r } ^ { 1 1 } = - G _ { 1 } ^ { 1 1 } = 2 i \cosh ( \omega _ { r } \gamma t ) \,,$$

$$G _ { r } ^ { 1 2 } = G _ { 1 } ^ { 2 1 } = \cosh ^ { 2 } ( \omega _ { r } \gamma t ) + \sinh ^ { 2 } ( \omega _ { r } \gamma t ) \,. \quad ( 1 4 0 )$$

For the remaining modes, the covariance matrix keeps its initial form throughout the evolution, glyph[negationslash]

$$G _ { k } ( t ) = \left ( \begin{array} { c c } 0 & 1 \\ 1 & 0 \end{array} \right ) \,, \ \ k \neq ( 1, 1 ) \,. \quad \ \ ( 1 4 1 )$$

In order to compute the entanglement entropy via Eq. (93), we need to project the covariance matrix to a subspace of interest. Here, a special care must be taken so that a nontrivial evolution is obtained. Since distinct modes evolve independently, if we choose the subsystem to be a single mode, then the entropy will be zero. In order for the entropy to be nonzero, we should consider subspaces that intersect both the resonant mode and its complement.

As a simple example, we study a mixture of two modes. In order to define the subsystem A of interest, let us first introduce an extended subsystem E including the resonant mode and some other (arbitrary) mode s . The only mode with a nontrivial evolution is the resonant mode; for all others, the covariance matrix keeps its initial form. As a result, regardless of which mode s one chooses to build the subsystem E , the projection of the covariance matrix to E will always have the form

$$G _ { E } ( t ) = \left ( \begin{array} { c c c c } 2 i \cosh ( \omega _ { r } \gamma t ) & 0 & \cosh ^ { 2 } ( \omega _ { r } \gamma t ) + \sinh ^ { 2 } ( \omega _ { r } \gamma t ) & 0 \\ 0 & 0 & 0 & 1 \\ \cosh ^ { 2 } ( \omega _ { r } \gamma t ) + \sinh ^ { 2 } ( \omega _ { r } \gamma t ) & 0 & - 2 i \cosh ( \omega _ { r } \gamma t ) & 0 \\ 0 & 1 & 0 & 0 \end{array} \right ).$$

$$( 1 4 2 )$$

We now define a new set of annihilation operators

$$\hat { \tilde { a } } _ { 1 } = \frac { \hat { a } _ { 1 } + \hat { a } _ { s } } { \sqrt { 2 } }, \quad \hat { \tilde { a } } _ { s } = \frac { \hat { a } _ { 1 } - \hat { a } _ { s } } { \sqrt { 2 } } \quad \ \ ( 1 4 3 ) \quad \ \ a g \epsilon$$

and their conjugates ˜ ˆ a † 1 and ˆ . ˜ a † s This symplectic transformation defines a new basis on the subsystem E . The transformation matrix B that connects the new and the old basis is

As discussed in [14], in any system displaying Lyapunov instabilities, the asymptotic production of entropy for a generic subsystem takes place at a constant rate Λ A , which is determined by the Lyapunov exponents of the system. We will now show that Λ A = ω γ r in the present case, in agreement with Eq. (149).

$$B = \frac { 1 } { \sqrt { 2 } } \left ( \begin{array} { c c c c } 1 & 1 & 0 & 0 \\ 1 & - 1 & 0 & 0 \\ 0 & 0 & 1 & 1 \\ 0 & 0 & 1 & - 1 \end{array} \right ). \quad \begin{array} { c c c c } & \exp \\ & \text{of } f \\ ( 1 4 4 ) & \text{larg} \\ & \text{cha} \\ & \text{1} \end{array}$$

In order to mix the modes and find a nonzero entropy, we apply this transformation to G E

$$\tilde { G } _ { E } ( t ) = B G _ { E } B ^ { - 1 }, \quad \quad ( 1 4 5 ) \quad \underset { \iota \iota } { o u }$$

and then restrict to the subspace spanned by the first transformed mode, which gives the 2 × 2 restricted covariance matrix in the transformed basis, which we write as

$$\tilde { G } _ { A } ( t ) = \left ( \begin{array} { c c } \frac { 1 } { 2 } i \sinh ( 2 \omega _ { r } \gamma t ) & \cosh ^ { 2 } ( \omega _ { r } \gamma t ) \\ \cosh ^ { 2 } ( \omega _ { r } \gamma t ) & - \frac { 1 } { 2 } i \sinh ( 2 \omega _ { r } \gamma t ) \end{array} \right ).$$

It is now simple to compute the R´nyi entropy as e

$$R _ { A } ( t ) & = \frac { 1 } { 2 } \log \left ( - \det ( \tilde { G } _ { A } ) \right ) \\ & = \frac { 1 } { 2 } \log \left ( \cosh ^ { 2 } ( \omega _ { r } \gamma t ) \right ). \quad \quad ( 1 4 7 ) \quad \text{Cor} \\ & \text{asy} ]$$

In the limit of t glyph[greatermuch] 1, the R´nyi entropy becomes e

$$R _ { A } ( t ) \sim - \log 2 + \omega _ { r } \gamma t \,. \quad \quad ( 1 4 8 ) \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \ \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \ \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \,$$

From Eq. (132), the entanglement entropy has the same asymptotic behavior, differing only by an offset

$$S _ { A } ( t ) \sim 1 - 2 \log 2 + \omega _ { r } \gamma t \,. \quad \quad ( 1 4 9 )$$

Hence, entropy is produced at a constant rate for large times, with a rate determined by the frequency ω r of the resonant mode, the amplitude of oscillation glyph[epsilon1] of the mirror, and the dimensions L , L 1 2 of the cavity, via the parameter γ . We see in Fig. 5 that, for large t , Eq. (149) gives the same asymptotic behavior of the exact entanglement entropy computed from Eqs. (129) and (147).

This result can be interpreted as signaling the presence of an instability in the system, associated with the fast amplification of excitations in the resonant mode.

The asymptotic rate of entanglement entropy production Λ A for a generic subsystem with N A degrees of freedom is given by the sum of the largest 2 N A Lyapunov exponents [14]. Since we are dealing with a single degree of freedom, N A = 1 here, and Λ A is the sum of the two largest Lyapunov exponents. We will later give a precise characterization of such generic subsystems.

glyph[negationslash]

The Lyapunov exponents can be computed from the Hamiltonian flow M t ( ) of the system. In the special case where M t ( ) = exp( Kt ), they are simply the eigenvalues of the time-independent symplectic generator K . In our case, the generator decomposes into a series of 2 × 2 blocks K k associated with the independent modes k , as discussed in Section II B 2. The evolution of the nonresonant modes is trivial, K k = 0, for k = r , and the associated Lyapunov exponents vanish. The symplectic generator of the resonant mode is given in Eq. (57). The Lyapunov exponents are the eigenvalues

$$\lambda _ { 1 } ^ { ( r ) } = \omega _ { r } \gamma \,, \quad \lambda _ { 2 } ^ { ( r ) } = - \omega _ { r } \gamma \,. \quad \ \ ( 1 5 0 )$$

The two largest Lyapunov exponents of the system are then λ ( r ) 1 = ω γ r and 0, leading to

$$S _ { A } ( t ) \sim \Lambda _ { A } t \,, \quad \Lambda _ { A } = \omega _ { r } \gamma \,. \quad \ \ ( 1 5 1 )$$

Comparing Eqs. (149), (150) and (151), we see that the asymptotic rate of growth of the entanglement entropy is precisely the value of the positive Lyapunov exponent, which is also equal to Λ A .

The coefficent Λ A can also be computed directly from the equation of motion of the normal modes in the case of parametric resonance. A simple model of parametric resonance is provided by the Mathieu equation

$$\ddot { x } ( t ) + \omega ^ { 2 } ( t ) x ( t ) = 0 \,, \quad \quad \ \ ( 1 5 2 )$$

where

$$\omega ^ { 2 } ( t ) = \omega _ { 0 } ^ { 2 } + \alpha \cos ( \Omega t ) \text{ \quad \ \ } ( 1 5 3 )$$

is the time-dependent natural frequency, which oscillates around ω 0 with frequency Ω = 2 ω 0 and amplitude α . It is a well known fact that the Floquet exponents µ can be computed directly from this equation as

$$\mu = \frac { \alpha } { 4 \omega _ { 0 } } \,. \quad \quad \quad ( 1 5 4 )$$

FIG. 5. Comparison of the time-dependent entanglement entropy S A ( τ ) of a resonant mode in (2 + 1) dimensions with its asymptotic expansion given in Eq. (149).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00835v1-figure-d43c3100-d5e8-4579-8a71-06c04e456e64.jpg)

In addition, the Floquet exponent is the positive Lyapunov exponent for the stroboscopic evolution in this case. Therefore, in the dynamical Casimir effect, we can write the time-dependent frequency Ω ( ) of the resonant r t mode, given by Eq. (46), in the form (153) and compute the Lyapunov exponents directly from Eq. (154). With r = (1 1), we have ,

$$\Omega _ { r } ^ { 2 } = \pi ^ { 2 } \left [ \frac { 1 } { L _ { 1 } ^ { 2 } ( 1 - \epsilon \cos ( 2 \omega _ { r } t ) ) ^ { 2 } } + \frac { 1 } { L _ { 2 } ^ { 2 } } \right ] \,. \quad ( 1 5 5 )$$

To first order in glyph[epsilon1] ,

$$\Omega _ { r } ^ { 2 } = \omega _ { r } ^ { 2 } + 2 \epsilon \frac { \pi ^ { 2 } } { L _ { 1 } ^ { 2 } } \cos ( 2 \omega _ { r } t ) \,. \quad \quad ( 1 5 6 ) \quad \text{$\text{$th$} \\ \text{$gpu$}$$

Comparing with the Mathieu equation, we see that α = 2 glyph[epsilon1]π 2 /L 2 1 . Substituting into Eq. (154), setting ω 0 = ω r , and using Eq. (137), we find

$$\mu = \omega _ { r } \gamma = \Lambda _ { A } \,, \quad \quad \quad ( 1 5 7 )$$

which agrees with (150).

In the more explicit derivation of the entropy production rate leading to Eq. (149), we considered a special choice of subsystem A . The result remains valid in a much more general context, however. In fact, the arguments leading to the same result in Eq. (151) apply to a large class of generic finite-dimensional subsystems, which can be fully characterized from Theorem 2. In the next section, we provide a detailed characterization of such generic subsystems that applies to any spatial dimension d ≥ 2.

## C. Higher dimensions

For any spatial dimension d ≥ 2, the dynamical Casimir effect at resonance for the lowest energy mode has the same overall features in the slow-variation approximation, as discussed in Section II B 2. The time evolution is trivial for all modes, except for the resonant one, which behaves as a single oscillator at parametric resonance. The evolution of the resonant mode is described by the symplectic transformation (56), with the natural frequency ω r and the coefficient γ given by Eqs. (48) and (50), respectively. A Lyapunov instability is associated with the resonance. For a generic subsystem, entropy is then produced asymptotically at a constant rate Λ A = ω γ r , which is equal to the positive Lyapunov exponent, Λ A = λ ( r ) 1 . We now present a simple characterization of such generic subsystems and prove that they satisfy Λ A = ω γ r .

Let glyph[lscript] 1 , glyph[lscript] 2 ∈ V ∗ be the eigenvectors of the generator K ᵀ whose eigenvalues λ , λ 1 2 are the nonzero Lyapunov exponents λ ( r ) 1 = -λ ( r ) 2 ,

$$\ell _ { 1 } = \frac { 1 } { \sqrt { 2 } } \begin{bmatrix} 1 \\ i \end{bmatrix} \,, \quad \ell _ { 2 } = \frac { 1 } { \sqrt { 2 } } \begin{bmatrix} - 1 \\ i \end{bmatrix} \,. \quad ( 1 5 8 )$$

They form a canonical pair of V , Ω ab ( glyph[lscript] 1 ) a ( glyph[lscript] 2 ) b = 1. The remaining eigenvectors of K ᵀ have vanishing eigenvalues. We can construct a linear basis D L = { glyph[lscript] i } of V composed of canonically conjugate eigenvectors of K that includes glyph[lscript] 1 and glyph[lscript] 2 . The evolution of the basis vectors under the Hamiltonian flow has a simple form:

$$M ^ { \intercal } ( t ) \, \ell _ { 1 } & = e ^ { \omega _ { r } \gamma t } \ell _ { 1 } \,, \\ M ^ { \intercal } ( t ) \, \ell _ { 2 } & = e ^ { - \omega _ { r } \gamma t } \ell _ { 2 } \,, \\ M ^ { \intercal } ( t ) \, \ell _ { i } & = \ell _ { i } \,, \quad i > 2 \,.$$

This basis is adequate for analyzing the expansion of subsystems with time, which determines the asymptotic growth of the R´nyi and entanglement entropies, as dise cussed in Section III D.

Let A be a proper symplectic subspace of V representing a subsystem of interest. Its symplectic complement B is defined as

$$B = \left \{ v \in V \ | \ \Omega _ { a b } v ^ { a } u ^ { b } = 0, \forall u \in A \right \}. \quad \ \ ( 1 6 0 )$$

We have V = A ⊕ B , and any vector can be uniquely decomposed as v = v A + v B , with v A ∈ A and v B ∈ B . The symplectic projection of v to the subsystem A is the component v A in this decomposition. The decomposition V = A ⊕ B induces a unique decomposition V ∗ = A ∗ ⊕ B ∗ of the dual phase space based on the isomorphism induced by Ω -1 ab , i.e. we can define A ∗ = Ω { -1 ab v b | v a ∈ A } and B ∗ = { Ω -1 ab v b | v a ∈ B } . We call A (or, equivalently, A ∗ ) a generic subsystem if:

- (i) glyph[lscript] 2 has a nonzero symplectic projection on A ∗ ;
- (ii) glyph[lscript] 2 has a nonzero symplectic projection on B ∗ .

The condition (i) implies that there exists a vector θ 1 ∈ A ∗ such that glyph[negationslash]

$$\Omega ^ { a b } ( \theta _ { 1 } ) _ { a } ( \ell _ { 2 } ) _ { b } \neq 0 \,. \quad \quad ( 1 6 1 )$$

glyph[negationslash]

Expanding θ 1 in the basis D L , θ 1 = ∑ k ′ θ k ′ 1 glyph[lscript] k ′ , and exploring the fact that D L is a Darboux basis, we find that θ 1 ′ 1 = 0, that is, θ 1 has a nonzero expansion coefficient in the first basis vector of D L .

glyph[negationslash]

glyph[negationslash]

Now let D A = { θ r | r = 1 , . . . , 2 N A } be a Darboux basis of the subsystem A that includes θ 1 . Expanding its elements in the basis D L , we can write θ r = θ k ′ r glyph[lscript] k ′ . By a simple procedure of Gaussian elimination, we can construct a new linear basis ˜ D A = { ˜ θ r } of A such that ˜ θ 1 = θ 1 , and ˜ θ 1 r = 0, for r = 1. Moreover, it follows from condition (ii) that, for all r , ˜ θ k ′ r = 0 for some k ′ = 2, otherwise glyph[lscript] 2 would be an element of A and have a vanishing symplectic projection on B . In short, the only vector with a contribution from the expanding unstable mode is θ 1 , and all other modes have a contribution from some stable mode. As a result, for large times, we have glyph[negationslash]

$$& \| M ( t ) \, \theta _ { 1 } \| \, \in \, e ^ { \omega _ { r } \gamma t } \,, \quad \quad \ \ t r a r _ { \substack { i \geq 2 \,. \quad \\ \| M ( t ) \, \theta _ { i } \| \rightarrow \text{constant} \,, \ \ i \geq 2 \,. \quad \ \ ( 1 6 2 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

We are interested in computing how the volume of the unit cube with sides θ r evolves for large times, Vol( M ᵀ ( ) t D A ). Since the basis D A is related to ˜ D A by a fixed transformation, their volumes are related by a time-independent factor, Vol( M t ( ) D A ) ∝ Vol( M t ( ) ˜ D A ). But then, from Eq. (162),

$$\text{Vol} ( M ^ { \tau } ( t ) \mathcal { D } _ { A } ) \in e ^ { \omega _ { r } \gamma t } \,, \quad \quad ( 1 6 3 ) \quad ^ { \text{act} } _ { \text{tha} }$$

and we find for the subsystem exponent:

$$\Lambda _ { A } = \lim _ { t \to \infty } \frac { 1 } { t } \log \frac { \text{Vol} ( M ^ { \intercal } ( t ) \mathcal { D } _ { A } ) } { \text{Vol} ( \mathcal { D } _ { A } ) } \quad \ \ ( 1 6 4 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

$$= \omega _ { r } \gamma \,. \quad \text{ \quad \ \ } ( 1 6 5 ) \quad \text{ on} \, \widehat { \text{In} } \,.$$

This completes the proof that Λ A = ω γ r for a generic subsystem satisfying conditions (i) and (ii).

Note that in (2 + 1) dimensions and higher, the Lyapunov instability is restricted to the resonant mode, so that the condition (ii) is violated for the subsystem composed of the resonant mode only. This is the reason why we had to consider a subsystem mixing more than one mode to obtain a nontrivial behavior in Section IV B.

## V. DISCUSSION

We analyzed the production of R´ enyi and entanglement entropies for the dynamical Casimir effect at resonance in arbitrary dimensions. In our settings, a d -dimensional parallelepided has one of its perfectly reflecting boundaries oscillating harmonically, while all others are fixed. By setting the frequency of oscillation equal to twice that of some normal mode r , the system is set at parametric resonance and particles are continuously produced from the vacuum. We considered the regime where the oscillations of the mirror are small, with a relative amplitude glyph[epsilon1] = ∆ L/L glyph[lessmuch] 1. We combined analytical techniques introduced in [10] for the study of pair creation under these circumstances with general symplectic techniques applicable to the analysis of entanglement production in general time-dependent linear systems [13, 14] in order to provide a thorough description of the build up of correlations and entanglement at resonance.

The system has two characteristic time scales: the period of the resonant mode, ω -1 r , and the time scale ( ω glyph[epsilon1] r ) -1 for which the system starts to depart considerably from the initial configuration at the ground state, which depends on the amplitude of oscillation of the mirror. Since glyph[epsilon1] glyph[lessmuch] 1, these scales are well separated. By averaging over the fast oscillations of normal modes, the equations of motion for the mode amplitudes are simplified and can be solved exaclty, as shown in [10]. The cases of d = 1 and 3 are treated in [10], where a Heisenberg representation is used for d = 1 and a Schr¨dinger repreo sentation for d = 3. We extended the treatment for arbitrary dimensions, and showed that it can be implemented in the Heisenberg representation in any dimension.

The moving boundary induces a nontrivial evolution of the cavity modes, which in general includes two kinds of interactions: particle creation from the vacuum and mode coupling. The vacuum becomes unstable due to the time-dependent boundary conditions, and the created excitations can be transferred from one mode to another.

For any spatial dimension d ≥ 2, the intermode interactions are strongly suppressed. This is due to the fact that the energies required to excite nonresonant modes do not match differences between energy levels of the resonant mode, making such transitions unlikely. The amplitude of the resonant mode then increases exponentially, being continuously pumped by the parametric resonance, providing an example of a Floquet instability. In the approximation where the fast oscillations are integrated out, the mode amplitude grows monotonically, and the original Floquet instability then corresponds to a Lyapunov instability. We found that, for any subsystem that intersects both the resonant mode and its complement, R´ enyi and entanglement entropies are produced at a constant rate for large times, S A ∼ R A ∼ Λ A t , where the production rate is equal to the positive Lyapunov exponent of the system, as expected from general results of [14]. It can be explicitly written as Λ A = λ ( r ) 1 = ω γ r , with ω r and γ given by Eqs. (48) and (50).

An alternate technique for studying a Floquet instability in a system with periodicity τ consists of looking at its stroboscopic evolution. That is, if M t ( ) is its Hamiltonian flow, one focus on the discrete series of instants M nτ ( ), n ∈ Z . An effective, time-independent Hamiltonian can then be defined in terms of τ -1 log M τ ( ) such that at full periods its evolution agrees with that of the original system. This effective Hamiltonian displays a Lyapunov instability associated with the original Floquet instability. Production of entropy in systems with Floquet instabilities was studied in detail in [14] within the stroboscopic approach. Applying it to the dynamical Casimir effect at resonance for d ≥ 2, we found the same results for the entanglement entropy production as in the

approximation introduced in [10], where the fast oscillations are integrated out. Our result suggests that, for systems satisfying a slow-variation condition, M τ ( ) ∼ 1 , the technique of averaging out the fast oscillations provides a good approximation for computing the Lyapunov exponents of the effective Hamiltonian associated with the stroboscopic dynamics of the system.

The case of one spatial dimension displays special features that lead to a different asymptotic behavior for the R´ enyi and entanglement entropies. In constrast to what happens in higher dimensions, for d = 1, the energy required to excite any normal mode is a multiple of the difference glyph[planckover2pi1] ω r between energy levels of the lowest energy mode, which we set at resonance. All modes are then strongly coupled, and excitations produced at the resonant mode can transition to other modes. The energy injected in the field by the moving mirror thereby spreads throughout an infinite number of modes. As a result, two new effects take place. First, the number of particles in the resonant mode, instead of growing exponentially, increases only linearly with time, as first recognized in [10]. In addition, the resonant mode is strongly entangled with the nonresonant modes. We found that the associated R´nyi and entanglement entropies of the resoe nant mode display a logarithmic growth, S A ∼ 1 / 2 log τ , where τ = glyph[epsilon1]ω t/ 2 is a dimensionless time variable. It is interesting that, even in the presence of parametric resonance, the system does not display features expected for unstable systems, i.e., a linear production of entropy. Instead, its behavior is qualitatively similar to that of metastable Hamiltonians discussed in [26], except for the fact that the logarithmic function is here multiplied by a factor of 1 / 2, while for metastable systems the prefactor is always an integer.

The techniques here applied for the analysis of entropy production in the dynamical Casimir effect at resonance can also be applied to other mirror trajectories. For finite cavities, the adaptation of the procedure followed

- [1] G. T. Moore, Journal of Mathematical Physics 11 , 2679 (1970).
- [2] V. V. Dodonov, Phys. Scripta 82 , 038105 (2010).
- [3] P. C. W. Davies and S. A. Fulling, Proc. Roy. Soc. Lond. A356 , 237 (1977).
- [4] R. D. Carlitz and R. S. Willey, Phys. Rev. D36 , 2327 (1987).
- [5] I. H. Brevik, K. A. Milton, S. D. Odintsov, and K. E. Osetrin, Phys. Rev. D62 , 064005 (2000), arXiv:hepth/0003158 [hep-th].
- [6] V. Dodonov, Advances in Chemical Physics 119 , 309 (2001).
- [7] M. P. E. Lock and I. Fuentes, New J. Phys. 19 , 073005 (2017), arXiv:1607.05444 [quant-ph].
- [8] C. M. Wilson, G. Johansson, A. Pourkabirian, M. Simoen, J. R. Johansson, T. Duty, F. Nori, and P. Delsing, Nature 479 , 376 (2011).
- [9] V. Dodonov and A. Klimov, Physics Letters A 167 , 309

in this work should be straightforward. A natural question for further developments is how to implement our techniques for the dynamical Casimir effect with a single mirror. As is well known, one of the main motivations for the study of the dynamical Casimir effect is that it can mimic the Hawking effect [28]. The application of symplectic techniques for the study of entropy production in general field theories is discussed in [14], and we hope our approach can be extended along such lines to the analysis of the time evolution of the entanglement entropy in models simulating aspects of particle creation in the Hawking effect [4, 29].

## ACKNOWLEDGMENTS

We thank Eugenio Bianchi and Rodolfo R. Soldati for inspiring discussions during the development of this project. We thank Eugenio Bianchi for sharing his notes on Floquet instabilities and for many invaluable discussions during the conception and initial stages of this work. IR acknowledges support from CAPES and CNPq. LH thanks the members of the physics department at the Universidade Federal de Minas Gerais for their hospitality during his visit in the summer of 2018. LH acknowledges support through the Brazil-U.S. Exchange Program of the American Physical Society which made this visit possible. LH is funded by the the Max Planck Harvard Research Center for Quantum Optics and supported by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germanys Excellence Strategy EXC-2111 39081486. NY acknowledges support from CNPq, Brazil, under the grant PQ 306744/2018-0, and from the Programa Institucional de Aux´ ılio ` Pesquisa de Docentes Rec´ em-Contratados, a PRPq/UFMG.

(1992).

- [10] V. V. Dodonov and A. B. Klimov, Phys. Rev. A53 , 2664 (1996).
- [11] A. Lambrecht, M.-T. Jaekel, and S. Reynaud, Physical review letters 77 , 615 (1996).
- [12] G. Plunien, R. Sch¨tzhold, u and G. Soff, Physical review letters 84 , 1882 (2000).
- [13] E. Bianchi, L. Hackl, and N. Yokomizo, Phys. Rev. D92 , 085045 (2015), arXiv:1507.01567 [hep-th].
- [14] E. Bianchi, L. Hackl, and N. Yokomizo, JHEP 03 , 025 (2018), arXiv:1709.00427 [hep-th].
- [15] L. F. Hackl, 'Aspects of gaussian states: Entanglement, squeezing and complexity,' (2018).
- [16] J. Eisert and M. Plenio, International Journal of Quantum Information 1 , 479 (2003).
- [17] C. Weedbrook, S. Pirandola, R. Garc´ ıa-Patr´n, o N. J. Cerf, T. C. Ralph, J. H. Shapiro, and S. Lloyd, Rev. Mod. Phys. 84 , 621 (2012).

- [18] A. Ashtekar and A. Magnon, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences 346 , 375 (1975).
- [19] R. M. Wald, Quantum field theory in curved spacetime and black hole thermodynamics (University of Chicago Press, 1994).
- [20] L. Hackl and R. H. Jonsson, arXiv preprint arXiv:1904.06246 (2019).
- [21] L. Vidmar, L. Hackl, E. Bianchi, and M. Rigol, Physical review letters 119 , 020601 (2017).
- [22] L. Vidmar, L. Hackl, E. Bianchi, and M. Rigol, Physical review letters 121 , 220602 (2018).
- [23] L. Hackl and R. C. Myers, Journal of High Energy Physics 2018 , 139 (2018).
- [24] L. Hackl, L. Vidmar, M. Rigol, and E. Bianchi, Physical Review B 99 , 075123 (2019).
- [25] R. D. Sorkin, in Tenth International Conference on General Relativity and Gravitation (held Padova, 4-9 July, 1983), Contributed Papers , Vol. 2 (1983) pp. 734-736.
- [26] L. Hackl, E. Bianchi, R. Modak, and M. Rigol, Phys. Rev. A97 , 032321 (2018), arXiv:1710.04279 [hep-th].
- [27] I. S. Gradshteyn and I. M. Ryzhik, Table of integrals, series, and products (Academic press, 2014).
- [28] P. C. W. Davies and S. A. Fulling, Proc. Roy. Soc. Lond. A348 , 393 (1976).
- [29] M. R. R. Good, P. R. Anderson, and C. R. Evans, Phys. Rev. D88 , 025023 (2013), arXiv:1303.6756 [gr-qc].
Ivan Romualdo Lucas Hackl Nelson Yokomizo
Entanglement production in the dynamical Casimir effect at parametric resonance
The particles produced from the vacuum in the dynamical Casimir effect are
highly entangled. In order to quantify the correlations generated by the
process of vacuum decay induced by moving mirrors, we study the entanglement
evolution in the dynamical Casimir effect by computing the time-dependent
R\'enyi and von Neumann entanglement entropy analytically in arbitrary
dimensions. We consider the system at parametric resonance, where the effect is
enhanced. We find that, in (1+1) dimensions, the entropies grow logarithmically
for large times, $S_A(\tau)\sim\frac{1}{2}\log(\tau)$, while in higher
dimensions (n+1) the growth is linear, $S_A(t)\sim \lambda\,\tau$ where
$\lambda$ can be identified with the Lyapunov exponent of a classical
instability in the system. In $(1+1)$ dimensions, strong interactions among
field modes prevent the parametric resonance from manifesting as a Lyapunov
instability, leading to a sublinear entropy growth associated with a constant
rate of particle production in the resonant mode. Interestingly, the
logarithmic growth comes with a pre-factor with $1/2$ which cannot occur in
time-periodic systems with finitely many degrees of freedom and is thus a
special property of bosonic field theories.