Non-Equilibrium Quantum Brain Dynamics: Super-Radiance and Equilibration in 2 + 1 Dimensions

We derive time evolution equations, namely the Schrödinger-like equations and the Klein–Gordon equations for coherent fields and the Kadanoff–Baym (KB) equations for quantum fluctuations, in quantum electrodynamics (QED) with electric dipoles in 2+1 dimensions. Next we introduce a kinetic entropy current based on the KB equations in the first order of the gradient expansion. We show the H-theorem for the leading-order self-energy in the coupling expansion (the Hartree–Fock approximation). We show conserved energy in the spatially homogeneous systems in the time evolution. We derive aspects of the super-radiance and the equilibration in our single Lagrangian. Our analysis can be applied to quantum brain dynamics, that is QED, with water electric dipoles. The total energy consumption to maintain super-radiant states in microtubules seems to be within the energy consumption to maintain the ordered systems in a brain.


Introduction
Numerous attempts to understand memory in a brain have been made over one hundred years starting at the end of 19th century. Nevertheless, the concrete mechanism of memory still remains an open question in conventional neuroscience [1][2][3]. Conventional neuroscience is based on classical mechanics with neurons connected by synapses. However, we still cannot answer how limited connections between neurons describe mass excitations in a brain in classical neuron doctrine.
Quantum field theory (QFT) of the brain or quantum brain dynamics (QBD), is one of the hypotheses expected to describe the mechanism of memory in the brain [4][5][6]. Experimentally, several properties of memory, namely the diversity, the long-term but imperfect stability and nonlocality (Memory is diffused and non-localized in several domains in a brain. It does not disappear due to the destruction in a particular local domain. The term 'nonlocality' does not indicate nonlocality in entanglement in quantum mechanics.), are suggested in [7][8][9]. The QBD can describe these properties by adopting infinitely physically or unitarily inequivalent vacua in QFT, distinguished from quantum mechanics which cannot describe unitarily inequivalence. Unitarily inequivalence represents the emergence of the diversity of phases and allows the possibility of spontaneous symmetry breaking (SSB) [10][11][12][13]. The vacua or the ground states appearing in SSB describe the stability of the states. Furthermore, the QFT can describe both microscopic degrees of freedom and macroscopic matter [10]. To describe stored information, we can adopt the macroscopic ordered states in QFT with SSB involving

The Two-Particle-Irreducible Effective Action and Time Evolution Equations
We begin with the following Lagrangian density to describe quantum electrodynamics (QED) with electric dipoles in 2 + 1 dimensions in the background field method [52][53][54][55], where A is the background coherent photon fields, a is the quantum fluctuations of photon fields, F µν [A] = ∂ µ A ν − ∂ ν A µ is the field strength, the α 1 is a gauge fixing parameter, the m is the mass of a dipole, the I is the moment of inertia, u i = (cos θ, sin θ) is the direction of dipoles and 2ed e is the absolute value of dipole vector. The variable θ represents the degrees of freedom of rotation of dipoles in 2 + 1 dimensions. The dipole-photon interaction term −2ed e Ψ * u i ΨF 0i [A + a] has the similar form to that in [27]. We shall expand the electric dipole fields Ψ and Ψ * by the angular momentum and consider only the ground state and the first excited states in energy-levels. Then we can write them as, in 2 + 1 dimensions. (In 3 + 1 dimensions, we might expand Ψ and Ψ * by spherical harmonics.) We can rewrite the terms in the above Lagrangian as, We also write the dipole-photon interaction term with electric fields F 0i = −E i by, with the direction of dipoles u i = (cos θ, sin θ). Next, we show two-particle-irreducible (2PI) effective action [47][48][49] for electric dipole fields and photon fields. Starting with the above Lagrangian density, we write the generating functional with the gauge fixing condition for quantum fluctuation, gauge fixing :a 0 = 0, (7) and perform the Legendre transformations. Then we arrive at, where the C represents the Keldysh contour [50,51] shown in Figure 1, the spatial dimension d = 2, the bar represents the expectation value · with the density matrix. The 3 × 3 matrix i∆ −1 0 (x, y) is defined as follows, for −1, 0 and 1, and the iD −1 0,ij (x, y) is written by, where i and j run over spatial components 1, · · ·, d = 2 in 2 + 1 dimensions. The 3 × 3 matrix ∆(x, y) is, where ∆ −10 (x, y) = T C δψ −1 (x)δψ * 0 (y) with time-ordered product T C in the closed-time path contour. The Green's function of dipole fields ∆ −10 (x, y) is also written by the 2 × 2 matrix ∆ ab −10 (x, y) with a, b = 1, 2 in the contour. The Green's function for photon fields D ij (x, y) represents, Figure 1. Closed-time path contour C. The label "1" represents the path from t 0 to ∞ and the label "2" represents the path from ∞ to t 0 .
Finally we write time evolution equations for coherent fields and quantum fluctuations. The 2PI effective action satisfies the following equations, due to the Legendre transformation of the generating functional. Equation (13) is written by, with iΣ ≡ − 1 2 δΓ 2 δ∆ . The matrix of self-energy Σ can be written by diagonal elements, since we can neglect the off-diagonal elements which are higher order of the coupling expansion. Equation (17) represents the Kadanoff-Baym equations for electric dipole fields in the two-energy-level approximation in 2 + 1 dimensions. Similarly, the Kadanoff-Baym equation for photon fields in Equation (14) is written by, with iΠ ≡ − δΓ 2 δD . Equation (15) is given by, with, Equation (20) represents the Klein-Gordon equations for spatial dimensions i = 1 and 2. Equation (16) is written by, and their complex conjugates. They are Schrödinger-like equations for coherent dipole fields. Equations (23) and (24) and their complex conjugates give the following probability conservation, We shall define J 0 (x) as, Then since we can use where the time dependent term in the time integral might be interpreted as an initial charge, but it is set to be zero. This equation represents the Poisson equation for scalar potential A 0 given by ∇ 2 A 0 = ∇ · µ with the vector of dipole moments −µ on the right-hand side in Equation (26). (Since the Fourier transformedÃ 0 (q) is written byÃ 0 (q) ∝ (q iμ i )/q 2 with µ i =μ i δ(r), the electric field E j = −∇ j A 0 (r) is proportional to q e iq·r q j q iμ i q 2 . If we can also apply the analysis in this section to the case in 3

Kinetic Entropy Current in the Kadanoff-Baym Equations and the H-Theorem
In this section, we derive a kinetic entropy current from the Kadanoff-Baym equations with first order approximation of the gradient expansion and show the H-theorem for the leading-order approximations in the coupling expansion based on [56][57][58]. The analysis in this section is similar to that in open systems (the central region connected to the left and the right region) [59]. Since (−1, 1) and (1, −1) components in i∆ −1 0 (x, y) in Equation (9) are zero, the same procedures to rewrite the Kadanoff-Baym equations as those in open systems [59][60][61][62][63] can be adopted. We set t 0 → −∞.
First, we shall write the Kadanoff-Baym equations in Equation (17) for each components. By multiplying the matrix ∆ from the right in Equation (17) and taking the (0, 0) component, we can write it as, where the (0, 0) component of the matrix . By taking (α, 0) component, we can write it as, It is convenient to introduce the Green's functions ∆ g,αα as, Then by using Equations (29) and (30), we can write ∆ α0 as, Equation (31) means the propagation from y to x with zero angular momentum, change of angular momentum at w and the propagation from w to x with angular momentum α = ±1. By using Equation (31), we can rewrite Equation (28) as, The second term on the left-hand side in Equation (32) represents the propagation from y to w with zero angular momentum, the change of the angular momentum to α = ±1 at w due to the coherent electric fields, the propagation from w to x and the change of the angular momentum from α = ±1 to zero due to the coherent electric fields. In a similar way to φ 4 theory in open systems [59], we can derive, where we have used, The (α, α) components of the Kadanoff-Baym equations are written by, and, where we have used Equations (31) and (34). Next, we shall perform the Fourier transformation ( d(x − y)e ip·(x−y) ) with the relative coordinate x − y of the (0, 0) and (α, α) components of the Kadanoff-Baym equations. We use the 2 × 2 matrix notation in the closed-time path with a, b, c, d = 1, 2. Equations (32) and (33) are transformed as, i∆ ac 00 (X, p) where X = x+y 2 , σ z = diag(1, −1), and the U αα (X, p) is the Fourier transformation, with the definition of ζ and |E|, and, The • is expanded by the derivative of X [64][65][66][67] as, with the definition of the Poisson bracket, We find that the U αα represents the change of momenta of dipoles as shown in Figure 2a.
In a similar way to [59], in the 0th and the first order in the gradient expansion in Equations (37) and (38), we can derive the following retarded Green's function, with the retarded parts (the subscript 'R') ∆ 00,R = i(∆ 11 00 − ∆ 12 00 ), Σ 00,R = i(Σ 11 00 − Σ 12 00 ) and U αα,R = i(U 11 αα − U 12 αα ). By taking the imaginary part of the retarded Green's function ∆ 00,R (X, p), we can derive the spectral function ρ 00 = i(∆ 21 00 − ∆ 12 00 ) = 2iIm∆ 00,R (X, p) which represents the information of dispersion relations. Similarly, the (α, α) components of the Kadanoff-Baym equations are written as, and, where, and, We can also write for ∆ cb g,αα (X, p) as, In the 0th and the first order in the gradient expansion in Equations (46) and (47), we can derive, with Here we have used the solution in the 0th and the first order in the gradient expansion in Equations (50) and (51) given by, with Σ αα,R = i(Σ 11 αα − Σ 12 αα ). The derivation is the same as [59]. The imaginary part of the retarded Green's function ∆ αα,R (X, p) multiplied by 2i represents the spectral function ρ αα = i(∆ 21 αα − ∆ 12 αα ) = 2iIm∆ αα,R (X, p) which represents the information of dispersion relations. In addition, the Kadanoff-Baym equations for photons in Equation (19) are written by, with, Next we shall derive the self-energy in the leading-order (LO) of the coupling expansion in Equation (6). The (a, b) = (1, 2) and (2, 1) component of i Γ 2 2 are given by, where t represents the transposition. It is convenient to rewrite, where T and L represent the transverse and the longitudinal part, respectively. The LO self-energy By Fourier-transforming with the relative coordinate x − y and multiplying δ ij − The second equation is due to the spatial dimension d = 2. Similarly, we arrive at, The Fourier transformation of the LO self-energy iΣ 12 00 (x, y) = − 1 2 δΓ 2,LO δ∆ 21 00 (y,x) is, Similarly, This self-energy is shown in Figure 2b. Similarly we can derive, and, Finally we derive a kinetic entropy current in the first order approximation in the gradient expansion and show the H-theorem in the LO approximation in the coupling expansion. By taking a difference of Equations (32) and (33), we arrive at, We use the Kadanoff-Baym Ansatz ∆ 12 where we just rewrite the (1, 2) and the (2, 1) components with the spectral parts ρ 00 , Σ 00,ρ and U αα,ρ and distribution functions f 00 , γ 00 and γ U,αα . The distribution functions f 00 , γ 00 and γ U,αα approach the Bose-Einstein distributions near equilibrium states. In the first order approximation in the gradient expansion in Equation (69) for (a, b) = (1, 2) and (2, 1), we can derive, (Rewrite (a, b) = (1, 2) and (2, 1) components in Equation (69), then we can show the collision We shall multiply ln i∆ 12 00 ρ 00 in (a, b) = (1, 2) component in Equation (69) and ln i∆ 21 00 ρ 00 in (2, 1) component in Equation (69), take the difference of them and integrate with p . By the use of Equation (70), we arrive at, with the definition of entropy current s µ matter,00 for (0, 0) component, We can derive the Boltzmann entropy p [(1 + n) ln(1 + n) − n ln n] with the number density n(X, p) in the quasi-particle limit ImU αα,R = ImΣ 00,R → 0 in the same way as in [58]. Similarly, we can derive a kinetic entropy current for (αα) components. >From Equations (46) and (47), we can derive We use the Kadanoff-Baym Ansatz ∆ 12 for distribution functions f αα , γ αα and γ V,αα by writing the (a, b) = (1, 2) and (2, 1) components in the Kadanoff-Baym equations (74). We can also show, from Equations (50) and (51). We shall multiply ln (74), take the difference of them and integrate with p . By using Equations (75) and (76), we arrive at, with the definitions of entropy current s µ matter,αα for (αα) components, In this derivation, we have used the same way as that in open systems in [59]. We can also derive the following equations for the Kadanoff-Baym equations for photons with the Kadanoff-Baym Ansatz D 21 (1 + f L ) and D 12 L = ρ L i f L with distribution functions f T and f L and spectral functions ρ T and ρ L , with the entropy current for photons, As a result, the total entropy current s µ = s µ matter,00 + ∑ α s µ matter,αα + s µ photon satisfies, where we have used the inequality (x − y) ln x y ≥ 0 for real variables x and y with x > 0 and y > 0. The equality is satisfied in with γ g,αα ∼ f αα in first order in the gradient expansion. We have shown the H-theorem in the LO approximation in the coupling expansion and in the first order approximation in the gradient expansion. There is no violation in the second law in thermodynamics in the dynamics.

Time Evolution Equations in Spatially Homogeneous Systems and Conserved Energy
In this section, we write time evolution equations in spatially homogeneous systems and show a concrete form of the conserved energy density.
It is convenient to introduce the statistical functions F 00 = . The variables of these functions are (X 0 , p 0 , p) with the center-of-mass coordinate X 0 = x 0 +y 0 2 and p given by the Fourier transformation with the relative coordinate x − y in variables (x, y) in Green's functions and self-energy in Section 2. The statistical functions and parts are real at any time when we start with real statistical functions at initial time. The spectral functions are given by taking the difference of (2, 1) and (1, 2) components multiplied by i, namely ρ 00 = i(∆ 21 00 − ∆ 12 00 ). They represent the information of which states can be occupied by particles in (p 0 , p) (dispersion relations). The spectral parts in self-energy are given by taking the difference of (2, 1) and (1, 2) components multiplied by i (and written by the subscript ρ), namely ∆ g,αα,ρ = i(∆ 21 g,αα − ∆ 12 g,αα ), Σ 00,ρ = i(Σ 21 00 − Σ 12 00 ) and so on. The spectral functions and parts are pure imaginary at any time when we start with pure imaginary spectral functions at initial time. We can use the real statistical parts labeled by the subscripts F and the pure imaginary spectral parts labeled by the subscript ρ in self-energy in the time evolution. We use the subscript 'R', 'F' and 'ρ' to represent the retarded, statistical and spectral parts in self-energy, respectively. The Kadanoff-Baym equation for the statistical and spectral functions are given by, ReU αα,R , F 00 + Re∆ 00,R , Σ 00,F − ∑ α U αα,F ReU αα,R , ρ 00 + Re∆ 00,R , Σ 00,ρ − ∑ α U αα,ρ = 0, (83) and longitudinal parts given by changing the label T to L in Equations (88) and (89). We can use Equation (69) in the previous section to derive Equations (82) and (83), for example. We can write, U αα,F (X, p) = (ed e ) 2 E(X) 2 ∆ g,αα,F (p + α∂ζ), U αα,ρ (X, p) = (ed e ) 2 E(X) 2 ∆ g,αα,ρ (p + α∂ζ), V αα,F (X, p) = (ed e ) 2 E(X) 2 F 00 (p − α∂ζ), V αα,ρ (X, p) = (ed e ) 2 E(X) 2 ρ 00 (p − α∂ζ).
In case we start with initial condition E 2 (X 0 = 0) = 0, ∂ 0 E 2 (X 0 = 0) = 0 and symmetric Green's functions for α → −α in spatially homogeneous systems, we can use ∂ζ = 0 in the above equations at any time. We can write the self-energy as, where we have omitted the label of the center-of-mass cordinate X in Green's functions and self-energy. We find that the Π T,F (k) = Π L,F (k) are symmetric (Π T,F (−k) = Π T,F (k)) under k → −k and that Π T,ρ = Π L,ρ are anti-symmetric (Π T,ρ (−k) = −Π T,ρ (k)) under k → −k, for any Green's functions for dipole fields. When we prepare initial conditions with symmetric F T,L and anti-symmetric ρ T,L for photons, we can derive symmetric F T,L and anti-symmetric ρ T,L at any time. In addition, since Π(k)'s are proportional to (k 0 ) 2 , there is no mass gap for incoherent photons for the leading-order self-energy in the coupling expansion. The velocity of gapless modes of incoherent photons will decrease when we increase the density of dipoles in this theory. Finally, we show the energy density E tot . In the spatially homogeneous system in the 2 + 1 dimensions, we can derive ∂E tot ∂X 0 = 0 with the energy density given by, where we have used the KB equations in this section, the Klein-Gordon Equation (20) and the Schödinger-like Equations (23) and (24) in Section 2. The first term represents the contribution of nonzero angular momenta for coherent dipole fields. The second term represents the contribution by electric fields E i = ∂ 0 A i . The third and the fourth terms represent the contribution by quantum fluctuations for dipoles and photons, respectively. When the temperature is nonzero T = 0 at equilibrium states and the spectral width in the spectral functions is small enough, statistical functions which are proportional to the Bose-Einstein distributions 1 e p 0 /T −1 give temperature-dependent terms mT 2 for dipole fields and ∝ T 3 for photon fields in 2 + 1 dimensions. The fifth term represents the potential energy in processes in Figure 2a. The sixth, seventh and eighth terms represent the potential energy in processes in Figure 2b. The coefficients in the sixth and seventh terms are not 1 3 but 1. While the factor 1 might look like a contradiction with the preceding research in [68,69] which suggest that the factor 1 3 appears in the interaction with 3-point-vertex, the factor 1 appears due to time derivative (∂ 0 ) 2 in self-energy for dipole fields and photon fields.

Super-Radiance
In this section, we show the super-radiance in time evolution equations for coherent fields with the rotating wave approximations neglecting non-resonant terms and quantum fluctuations. We have used the derivations in [70,71] for background coherent fields.
We shall consider only k 0 = 1 2I in this section and we expand the electric field E 1 and the transition rateψ 0ψ * 1 as, We consider the following case, Neglect non-resonant terms like e ±2ik 0 x 0 and quantum fluctuations (Green's functions ∆ 01 and ∆ 10 ) (the rotating wave approximation). Then from Equations (99)-(101), we arrive at the Maxwell-Bloch equations, We assume that , Z and R are independent of the spatial coordinate of the x 2 direction. We shall change → i in the above equations and assume real functions R = R * and = * . Then we can write, We find the conservation law with the definition B 2 ≡ 2R 2 + Z 2 , The relation ∂B ∂x 0 = 0 represents the probability conservation since we can rewrite B 2 = 2|ψ 1 | 2 + |ψ 0 | 2 2 by Equation (103) and Z ≡ 2|ψ 1 | 2 − |ψ 0 | 2 . We also find the following conservation law, which represents the energy conservation. By this relation, we might be able to estimate the maximum energy density of electric fields, in case there is no external energy supply. We derive the following solutions in Equations (108)-(110), with ∂θ ∂x 0 = √ 2ed e and the constant B in a similar way to [71]. The θ(x 0 ) swings around the position θ = π with the frequency Ω = ed e √ k 0 B in case we start with initial conditions at around θ 0 ∼ π (|ψ 1 | 2 = 0), since we can rewrite Equation (108) as The B is the order of the number density of dipoles. We introduce the damping term 1 L for the release of radiation and the propagation length L in Equation (108). We can write, In κ = 1 L time derivative, we can neglect the first term in the above equations, then The solution is, and, with τ R = κ (ed e ) 2 k 0 B and τ 0 = −τ R ln(tan θ 0 2 ). The τ R ∝ 1/B ∼ 1/N with the number of dipoles N represents the relaxation time of electric fields in the super-radiance. When N dipoles decay within time scales 1/N, the intensity of electric fields becomes the order N 2 (super-radiant decay with correlation among dipoles), not N (spontaneous decay without correlation among dipoles).

Higgs Mechanism and Tachyonic Instability
In this section, we rewrite time evolution equations for coherent fields with only real functions. We assume the spatially homogeneous case. We do not adopt the rotating wave approximation in this section. We show how coherent electric fields E 1 are affected by Z = 2|ψ 1 | 2 − |ψ 0 | 2 .
In Equation (101), the second derivatives of coherent fields on the right-hand side are written by, where we have used Equation (100). As a result, we arrive at, Re∆ g,11,R (X, p)F 00 (X, p) + ∆ g,11,F (X, p)Re∆ 00,R (X, p) with the x 1 direction of the dipole moment (density) given by µ 1 = 2ed e ψ * 1ψ 0 +ψ * 0ψ 1 , F 11 (X, p) = ∆ 21 11 (X,p)+∆ 12 11 (X,p) 2 , F 00 (X, p) = ∆ 21 00 (X,p)+∆ 12 00 (X,p) 2 and ∆ g,11,F (X, p) = ∆ 21 g,11 (X,p)+∆ 12 g,11 (X,p) 2 . In the Appendix A we have shown the detailed derivation for the second, third, fourth and fifth terms in the above equations. We have assumed the self-energy Σ 00 = Σ 11 = 0 in deriving the time derivatives of ∆ 10 and ∆ 01 in Equation (101). Even if we include contributions of self-energy in Equation (121), they are higher order O (ed e ) 4 in the coupling expansion. We have neglected higher order terms in the gradient expansion for quantum fluctuations. In Equation (121), we leave the −(∂ 2 ) 2 E 1 term on the left-hand side in the above equation to compare with the sign of − 2(ed e ) 2 Z I E 1 term. We find the Higgs mechanism with the mass squared − 2(ed e ) 2 Z I in the case of the normal population Z = 2|ψ 1 | 2 − |ψ 0 | 2 < 0. On the other hand, the tachyonic instability appears in the inverted population Z > 0 in the above equation. Then the electric field E 1 will increase exponentially until Z becomes negative. In Equation (121), the second term on the right-hand side is proportional to 2F 11 − F 00 − ∆ g,11,F . Near equilibrium states, we might find F 00 > 2F 11 − ∆ g,11,F , where statistical functions F 11 , F 00 and ∆ g,11,F are proportional to the Bose-Einstein distribution 1 e p 0 /T −1 plus 1 2 (with the Kadanoff-Baym ansatz) with different dispersion relations p 0 ∼ p 2 2m for F 00 and p 0 ∼ p 2 2m + 1 2I for F 11 and ∆ g,11,F , due to the energy difference 1 2I − 0 2I between the ground state and first excited states. So the 2F 11 − F 00 − ∆ g,11,F in the second term is negative near the equilibrium states, which might mean no tachyonic unstable terms appear from quantum fluctuations near equilibrium states. The contributions of quantum fluctuations on the right-hand side written by statistical functions (second, third, fourth and fifth terms) vanish at zero temperature T = 0. Quantum fluctuations represent finite temperature effects at equilibrium states, although we need not restrict ourselves to only the equilibrium case. We have shown general contributions of quantum fluctuations in both equilibrium and non-equilibrium case in this paper.
We can solve Equations (121)-(124) with real functions in this section and the Kadanoff-Baym equations with real statistical functions and pure imaginary spectral functions in Section 4, simultaneously.

Discussion
In this paper, we have derived time evolution equations, namely the Klein-Gordon equations for coherent photon fields, the Schrödinger-like equations for coherent electric dipole fields and the Kadanoff-Baym equations for quantum fluctuations, starting with the Lagrangian in quantum electrodynamics with electric dipoles in 2 + 1 dimensions. We have adopted the two-particle-irreducible effective action technique with the leading-order self-energy of the coupling expansion. We find that electric dipoles change their angular momenta due to coherent electric fields E 1 ± iαE 2 with α = ±1. They also change momenta and angular momenta by scattering with incoherent photons. The proof of H-theorem is possible for these processes as shown in Section 3. Our analysis provides the dynamics of both the order parameters with coherent fields and quantum fluctuations for incoherent particles.
In Section 2, we adopt two-energy level approximation for the angular momenta of dipoles. Then, we find that the i∆ −1 0 is written by 3 × 3 matrix with zero (−1, 1) and (1, −1) components. The form of the matrix is similar to 3 × 3 matrix in the analysis in open systems, the central region, left and right reservoirs as in [59,[61][62][63]. Hence we can simplify the Kadanoff-Baym equations for dipole fields in an isolated system with the same procedures as those in open systems. The difference between QED with dipoles and φ 4 theory in open systems is that the coherent electric field changes the momenta of dipoles when the phase αζ in E 1 ± iαE 2 with α = ±1 is dependent on space-time. The space dependence of coherent electric fields might disappear in the time evolution due to the instability by the lower entropy of the system, then electric fields will change angular momenta of dipoles but not change momenta p due to ∂ζ = 0. We can also trace the dynamics with ∂ζ = 0. By setting the initial conditions with the symmetry α → −α, namelyψ ( * ) α =ψ ( * ) −α , ∆ α0 = ∆ −α0 and ∆ 0α = ∆ 0−α , with initial conditions E 2 = 0 and ∂ 0 E 2 = 0 in spatially homogeneous systems in ∂ ν F ν2 = J 2 in Equation (20), we can show E 2 = 0 at any time. Then we can use ∂ζ = 0. This condition simplifies numerical simulations in the Kadanoff-Baym equations since we need not estimate the momentum shift p → p ± α∂ζ in the finite-size lattice for the momentum space. As a result, the simulations for Kadanoff-Baym equations for dipoles and photons will be similar to those in QED with charged bosons in [72].
In Section 3, we have introduced a kinetic entropy current and shown the H-theorem in the leading-order of the coupling expansion with ed e . This entropy approaches the Boltzmann entropy in the limit of zero spectral width as in [58]. The mode-coupling processes between dipoles and photons produce entropy. When there are deviations between (00) and (αα) components of distribution functions, entropy production occurs. Entropy production stops when the Bose-Einstein distribution is realized in the dynamics of Kadanoff-Baym equations.
We can also derive the energy shifts in dispersion relations due to nonzero electric fields by using the retarded Green's functions in Section 3. The 0th order equations for retarded Green's functions are given by, with ∆ g,11,R = −1 . Multiply p 0 − p 2 2m − 1 2I , take the imaginary parts in the above equations and remember the imaginary parts of retarded Green's functions are the spectral functions, then we find, W ρ 00 ρ 11 = 0, with, By setting determinant |W| to be zero, we find the following solutions for dispersion relations, Here we assumed the symmetry for α = ±1 for Green's functions and zero self-energy Σ 00 = Σ 11 = 0. We find how electric fields shift two energy levels 0 and 1 2I . The above energy shift is similar to the energy shift given in [27] in 3 + 1 dimensions due to nonzero electric fields.
In Section 5.1, we have derived the super-radiance from time evolution equations for coherent fields. We find that it is possible to derive the Maxwell-Bloch equations from our Lagrangian with the probability conservation law and the energy conservation law. Super-radiant decay with intensity of the order ∝ N 2 (N: The number of dipoles) appears in a similar way to [70,71]. It is possible to derive the maximum energy of electric fields by use of Equation (113). We know that the moment of inertia of water molecule is I = 2m H R 2 with m H = 940 MeV with R = 0.96 × 10 −10 m. Hence the k 0 = 1 2I = 1.1 × 10 −3 eV. Since B = N V = 3.3 × 10 28 /m 3 for liquid water, we find When we multiply the volume of all microtubules (MTs) in a brain, If we maintain our brain 100 s without energy supply, we need at least 0.1 × 10 −2 cal/s or 86 cal/day to maintain the ordered states of memory.
We can compare 86 cal/day with 4000 cal/day = 2000 kcal/day × 0.2 (energy consumption rate of brain) × 0.01 (energy rate to maintain the ordered system). The 86 cal/day is within the 4000 cal/day, which is consistent with our experiences. In this derivation, we have used coefficients in 2 + 1 dimensions and the number density of water molecules in 3 + 1 dimensions.
In Section 5.2, we have derived time evolution equations for electric field E 1 . The Higgs mechanism appears in this equation in normal population Z < 0. As a result, the dynamical mass generation occurs with the maximum mass Ω Higgs = 2ed e √ k 0 B = 30k 0 where the number density of dipoles is B = 2|ψ 1 | 2 + |ψ 0 | 2 = N V . The period is 2π/Ω Higgs = 1.3 × 10 −13 s. In normal population Z < 0, the Meissner effect appears with the penetrating length 1/Ω Higgs = 6.3 µm. On the other hand, the tachyonic instability occurs in inverted population Z > 0. The electric field E 1 increases exponentially with exp(ΩX 0 ) (with Ω ≤ Ω max ) where the time scale is 1/Ω max = 2.1 × 10 −14 s with Ω max = Ω Higgs . Due to energy conservation, since Z decreases as the absolute value of the electric field increases, tachyonic instability stops in Z < 0.
We have prepared for numerical simulations with time evolution equations, namely the Schödinger-like equations for coherent electric dipole fields, the Klein-Gordon equations for coherent electric fields and the Kadanoff-Baym equations for quantum fluctuations. Our simulations might describe the dynamics towards equilibrium states for quantum fluctuations and the dynamics of super-radiant states for coherent fields. Our analysis is also extended to simulations in open systems by preparing the left and the right reservoirs like those in [59] or networks [73].

Conclusions
We have derived the Schrödinger equations for coherent electric dipole fields, the Klein-Gordon equations for coherent electric fields and the Kadanoff-Baym equations for quantum fluctuations in QED with electric dipoles in 2 + 1 dimensions. It is possible to derive equilibration for quantum fluctuations and super-radiance for background coherent fields simultaneously. Total energy consumption to maintain super-radiance in microtubules is consistent with energy consumption in our experiences. We can describe dynamical information transfer with super-radiance via microtubules without violation of the second law in thermodynamics. We have also derived the Higgs mechanism in normal population and the tachyonic instability in inverted population. These dynamical properties might be significant to form and maintain coherent domains composed of dipoles and photons. We are ready to describe memory formation processes towards equilibrium states in 2 + 1 dimensions with equations in this paper. Furthermore, our approach might pave the way to understand the dynamical thinking processes with memory recalling in QBD by investigating the case in open systems with the Kadanoff-Baym equations. This work will be extended to the 3 + 1 dimensional analysis to describe memory formation processes in numerical simulations. We should derive the Schödinger-like equations, the Klein-Gordon equations and the Kadanoff-Baym equations by starting with the single Lagrangian in QED with electric dipoles in 3 + 1 dimensions in the future study. These equations in 3 + 1 dimensions will describe more realistic and practical dynamics in QBD. in Equation (101) with the symmetry ∆ 10 = ∆ −10 and ∆ 01 = ∆ 0−1 . It appears in taking the time derivative of J 1 (given by Equation (21) ) in Equation (20). Here ∆ 10 (x, x) and ∆ 01 (x, x) can be rewritten by, ∆ 01 (x, x) = − ed e i w ∆ 00 (x, w)E 1 (w)∆ g,11 (w, x), where we have used Equations (31) and (34) by setting E 2 = 0. We rewrite second time derivatives of ∆ 10 (x, x) and ∆ 01 (x, x). We shall rewrite Equation (30) without self-energy Σ αα as, where we have multiplied ∆ g,11 from the right and left of Equation (30). By using the above equations and Equations (32) and (33) with Equations (A1) and (A2) and ∆ −1 0,00 (x, y) = i ∂ ∂x 0 + ∆ 00 + 2∆ g,11 ed e E 1 ∆ 01 E 1 − ∆ g,11 E 1 iδ C = ed e 1 2I ∆ g,11 + iδ C E 1 ∆ 00 + 2∆ g,11 ed e E 1 ∆ 01 E 1 − ∆ g,11 E 1 iδ C , ∂ ∂x 0 ∆ 01 (x, x) = ed e (−2ed e E 1 ∆ 10 + iδ C ) E 1 ∆ g,11 + ∆ 00 E 1 − where δ C represents the delta function in the closed-time path. Here the terms proportional to ∇ 2 i are cancelled in spatially homogeneous systems. By use of the above two equations, we can show and, ∂ 2 ∂(x 0 ) 2 (∆ 10 + ∆ 01 ) = ed e 2iI ∆ g,11