Renormalization in Quantum Brain Dynamics

: We show renormalization in Quantum Brain Dynamics (QBD) in 3 + 1 dimensions, namely Quantum Electrodynamics with water rotational dipole ﬁelds. First, we introduce the Lagrangian density for QBD involving terms of water rotational dipole ﬁelds, photon ﬁelds and their interactions. Next, we show Feynman diagrams with 1-loop self-energy and vertex function in dipole coupling expansion in QBD. The counter-terms are derived from the coupling expansion of the water dipole moment. Our approach will be applied to numerical simulations of Kadanoff–Baym equations for water dipoles and photons to describe the breakdown of the rotational symmetry of dipoles, namely memory formation processes. It will also be extended to the renormalization group method for QBD with running parameters in multi-scales.


Introduction
The human brain is one of the most complex and poorly understood biological structures studied by science. This is especially true regarding the molecular mechanisms of how the human brain works, especially its information processing and memory encoding functions. While much is known about the architecture of brain cells and their individual activities in terms of the action potential, ion channels and ionic currents [1], much less is known about such issues as where our memories are stored or which molecular mechanisms are involved in information processing and higher cognitive functions. Much has been speculated recently about the possibility of at least some cognitive functions requiring the operation at a level of quantum physics or even quantum field theory [2,3]. Experimental determination or even corroboration of such ideas is currently beginning to be achievable in view of the modern experimental methods and instrumentation available to researchers [4].
Electromagnetic (EM) field theories of consciousness propose that consciousness results when a brain produces an electromagnetic field with specific characteristics. Pockett [5] and McFadden [6] have proposed EM field theories; Some electromagnetic theories are also quantum mind theories of consciousness; examples include the Quantum Brain Dynamics (QBD) approaches of Jibu and Yasue [7] and Vitiello [8], who proposed a quantum field theory explanation of how memory is created. In general, quantum mind theories other than these QBD approaches do not treat consciousness as an electromagnetic phenomenon. The basis for McFadden and Pockett's CEMI theory is the neuron's firing, which generates an action potential that is then the trigger for a postsynaptic potential in the neighboring neuron. Moreover, this ionic current also perturbs the surrounding electromagnetic field. McFadden has proposed that the brain's electromagnetic field represents an informationrich signal in the neurons. This hypothesis is supported by demonstrated effects of weak EM such as ecosystems, neuron net and cellular metabolism. In a cortical system, neuronal avalanches indicate power-law behaviors and emerge as a scale-free phenomenon [41]. Moreover, magnetoencephalographic data suggest that fractal-like brain functional networks emerge [42]. Thus the analysis of fractal-like phenomena in a brain has been a prospective approach.
Fractals are isomorphic to squeezed coherent states, as suggested by Vitiello [43,44]. Self-similarity patterns in fractals represent scale-free phenomena in physical systems. Scale-free properties emerge as fixed points in the Renormalization Group (RG). The RG method was argued by Stüeckelberg and Petermann [45] and by Gell-Mann and Low [46] and applied to statistical physics by Wilson and Kogut [47]. In calculating quantum corrections in QFT, we encounter ultra-violet (UV) divergences. The divergences are renormalized by a finite number of counter-terms in each QFT model if the model is a renormalizable theory. In renormalization, we then encounter the arbitrariness of finite parts in counter-terms in renormalization prescription. To achieve the invariance of theory under renormalization prescription, physical parameters such as mass and coupling can evolve in changing physical scales as coarse graining procedures, which are described by RG equations [47][48][49]. In RG, the contributions in smaller scales are renormalized into parameters such as mass and coupling as coarse graining procedures, and we then encounter an effective theory in macroscopic scales (from microscopic scales) adopted to describe multiscale brain functions in the beginning with the QBD Lagrangian and provide diversity in the effective theory for multi-scales in QBD. To derive RG equations for parameters, we need to analyze ultra-violet (UV) divergences appearing in self-energy and vertex corrections for cancellation of divergences by counter-terms in each QFT model. We might encounter fixed points where physical parameters become independent of renormalization scales for further coarse graining. The fixed points in RG represent scale-free fractal-like phenomena.
The aim of this paper is to cancel ultra-violet divergences in self-energy and vertex functions by counter-terms in the renormalization in QBD in 3 + 1 dimensions. We begin with the Lagrangian density in QBD and expand terms involving counter-terms for renormalization. We next calculate self-energy and vertex functions and cancel ultra-violet divergences by counter-terms in renormalizations. We extend our analysis to the 2-Particle-Irreducible (2PI) effective action technique [50,51] in the derivation of Kadanoff-Baym (KB) equations [52][53][54]. Counter-terms will be adopted to achieve the cancellation of UV divergences in lattice simulations for the time-evolution of KB equations in QBD. Properties of divergences will be used in deriving RG equations in future studies, where fixed point solutions might describe the fractality of water-photon systems. Recently, we have shown a super-radiant solution in water-photon systems in 3 + 1 dimensions, cooperative spontaneous photon emission [13]. Next, we have proposed the integration of QBD and holography using super-radiance and investigated the time-evolution of holograms [14]. Furthermore, we have provided control theory of coherent fields in the φ 4 model toward manipulating holograms or our subjective experiences [15]. The present work provides renormalization of UV divergences in quantum corrections in an interacting water-photon system, which will open up a new way to describe time-evolving holograms involving incoherent dipoles and photons, provide a control theory of holograms and introduce a RG method for QBD.
Our paper is organized as follows. In Section 2, we introduce the Lagrangian density for QBD with counter-terms. In Section 3, we calculate ultra-violet divergences in selfenergy and vertex functions. In Section 4, we extend our approach to 2PI effective action technique. In Section 5, we discuss our results. In Section 6, we provide concluding remarks and perspectives. We adopt the natural unit with the light speed and the Planck constant divided by 2π set to unity. The metric tensor is set to be g µν = diag(1, −1, −1, −1) with space-time subscript µ, ν = 0, 1, 2, 3 and spatial subscript i, j, n = 1, 2, 3.

Lagrangian Density
In this section, we introduce the Lagrangian density for Quantum Brain Dynamics in 3 + 1 dimensions and show Feynman diagrams for propagators of renormalized fields and counter terms for mass and vertices.
First, we begin with the Lagrangian density for QBD in 3 + 1 dimensions based on [13,30]. It is written by, with F µν = ∂ µ A ν − ∂ ν A µ with electric scalar potential A 0 and magnetic vector potential , the gauge fixing parameter ξ = 1, mass of dipoles m, moment of inertia for dipoles I with 1 I = 4 meV, chemical potential µ c and dipole moment for a water molecule 2ed e with elementary charge e and d e = 0.2 Å. We adopt the temporal axial gauge A 0 = 0. We assume the isotropic moment of inertia I for water dipoles. We adopt the two-energy-level approximation for angular momentum squared l 2 = − 1 for rotational degrees of freedom of water dipoles. We expand water electric dipole fields by the dipole field for the ground state ψ s (x) and the dipole fields for the 1st excited states ψ α (x) with α = 0, ±1 as, using the spherical harmonics for the ground state Y 00 (θ, ϕ) corresponding to the eigenvalue l 2 = 0 and the 1st excited states Y 1α (θ, ϕ) with α = 0, ±1 corresponding to eigenvalues l 2 = 2 shown by, The 1st and 2nd terms in Equation (1) are rewritten as, The 3rd term in Equation (1) is then written as, The 4th term in Equation (1) is written by, The 5th term in Equation (1) with the electric fields with the dipole moment density M i (i = 1, 2, 3) defined as, We can then rewrite as, Next, we rewrite the Lagrangian density with several counter terms. Using the bare (α = 0, ±1) involving the renormalization coefficient Z ph with respect to wave functions for photons, coefficient Z s for dipoles in the ground state and coefficient Z 1 for dipoles in the 1st excited state, we write the Lagrangian density in QBD by the bare fields as, where we adopt the bare mass m s and m 1 , the bare moment of inertia I 0 , the bare chemical potential µ 0 and the bare dipole moment 2(ed e ) 0 . The above Lagrangian density is rewritten as, Next, we introduce the parameters as follows, Then the Lagrangian density is rewritten as, The 5th, 6th, 7th and 8th terms in the above Lagrangian represent counter-terms. Finally, we show Feynman diagrams for propagators, vertices and counter-terms in Fourier transformation. Using the above Lagrangian density, we show the diagrams in Figure 1. The propagator for photons D ij (k), the propagator for dipoles in the ground state G ss (p) and the propagators for dipoles in the 1st excited states G αα (p) with α = 0, ±1 are written by, respectively. The vertices represent, with matrix element α| · |k, s of z (iM i (z)E i (z)) (shown in Equation (13)) describing processes of incoming photons with factor e −ik·u 3 , incoming dipoles in the ground state s with factor e −ip·u 1 and outgoing dipoles in 1st excited states α with factor e il·u 2 involving integration u 1 ,u 2 ,u 3 , and with k, s| · |α representing processes of outgoing photons with factor e ik·w 3 , outgoing dipoles in the ground state with factor e ip·w 1 and incoming dipoles in the 1st excited states with factor e −il·w 2 involving integration w 1 ,w 2 ,w 3 . We show vertices in Figure 1. Here we find that the vertex with incoming dipoles in 1st excited states involves the factor with (ik 0 ) × · · · and that the vertex with outgoing dipoles in 1st excited states involves the factor with (−ik 0 ) × · · ·. Using Figure 1, we perform renormalization of self-energy and vertices in the next section. Figure 1. Diagrams for propagators, vertices and counter-terms. The wavy line represents a photon propagator with subscripts i, j = 1, 2, 3 for the polarization of photons, the solid line represents a dipole propagator for the ground state, and the dotted line represents a dipole propagator for the 1st excited states α = 0, ±1. The gray-circle represents the counter-terms.

Renormalization
In this section, we write Feynman diagrams for self-energy and vertices and cancel ultra-violet (UV) divergences by counter-terms. We have referred to [49] for calculations.

1-Loop Self-Energy
In this section, we show a 1-loop self-energy for photons and dipoles, and we find that UV divergences are canceled by counter-terms.
First, we shall see the self-energy for photons Π ij (k) in 1-loop order shown in Figure 2. We find that 1-loop self-energy for photons contains the factor, 4 , since the 1st derivative of residue −1 by p 0 is zero. Here we have used the Feynman's formula, Similarly, we also find the relation, As a result, the self-energy for photons in 1-loop order is zero. We then find δ Z ph = 0 or Z ph = 1 in 1-loop order. Next, we shall see self-energy for dipoles in 1st excited states, as shown in Figure 3. Using the Feynman rules in Figure 1, the self-energy Σ αα with α = 0, ±1 can be written by, The detailed calculations are shown in Appendix A. The result is given by, where the subscript E represents the Euclidean 4-dimensional momenta, and [numerator 1] represents, We find Σ αα (p) involves cubic, quadratic, linear and logarithmic divergences in Equation (27); for example, the term (l 0 (28) induces cubic divergence. They are canceled by counter-terms in Figures 1 and 3. We then find the relation, We find that δ Z 1 and δ m −1 1 are independent of the chemical potential µ c . Finally, we show counter-terms for Σ ss . Since we can write, the counter-terms can be derived by changing coefficient 4 3 to 4 and factor µ c to µ c − 1 I for Σ αα (p) in Equation (26). We can then derive, with,

1-Loop 3-Point Vertex
In this section, we show a 1-loop 3-point vertex and find that UV divergences are canceled by a counter-term. We shall see a 1-loop 3-point vertex shown in Figure 4. The UV divergences in this figure are canceled by counter-terms involving δ ed e for vertices in Figure 1. The diagram in Figure 4a represents the term given by, The detailed calculations are given in Appendix B. We encounter the counter-term δ ed e as, The vertex correction involves quadratic, linear and logarithmic divergences canceled by counter-terms with δ ed e . We have derived the counter-term for 1-loop vertex correction. We can derive the same result from Figure 4b.

2-Loop Self-Energy
In this section, we estimate divergences in self-energy in 2-loop order.
First, we shall investigate the self-energy 1 for photons in Figure 5. This diagram contains the following factor, (35) where we have used the Feynman's formula. The function f is independent of p 0 . Since the 3rd derivative of residue one by p 0 is zero, the contribution of this diagram is zero. Next, we shall estimate the self-energy 2 on the left of Figure 6. This diagram contains a divergent part of the propagator of dipoles in the ground state, which can be canceled by the right diagram with counter-term. Even if the finite part, such as Cp 0 with finite constant C remains after the cancellation, by using Feynman's formula, we encounter the form of where g is a function independent of p 0 . Since the 2nd derivative of the residue Cp 0 by p 0 is zero, the remaining part in the above equation is zero. The same estimation is possible for the case when the propagator of dipoles in the 1st excited states contains divergent self-energy in the inner loop. We find δ Z ph is zero or Z ph = 1 in 2-loop order. Next, we shall estimate the self-energy for dipoles in the ground state in 2-loop order in the left diagram in Figure 7a. This diagram contains quadratic divergences from vertex correction times cubic divergences from the remaining part. The quadratic divergences from vertex corrections are canceled by the two middle diagrams in Figure 7a with counterterms. The sum of the left diagram and these middle diagrams contains cubic divergence in their loops. The cubic divergence can be canceled by the right counter-term in a similar way to 1-loop self-energy. We estimate the self-energy in the left diagram in Figure 7b. This diagram involves cubic divergence in the inner loop. The cubic divergence is canceled by the middle diagram with a counter-term. The sum of the left and middle diagrams contains cubic divergence at most. The cubic divergence can be canceled by the right counter-term in a similar way to 1-loop self-energy. We can also estimate the self-energy in Figure 7c. This diagram contains zero contributions in inner loops so the contribution is zero. Finally, we shall show the self-energy for dipoles in the 1st excited states in Figure 8. The estimation of divergences is similar to the cases of self-energy for dipoles in the ground state. We just need to show the labels of external lines for dipoles in the 1st excited states α and β will be the same one. The labels α, β = 0, ±1 represent the three 1st excited states in our model originating from spherical harmonics Y 1α for rotational degrees of freedom of dipoles. We can consider the three cases shown in Figure 8d, where there are two incoming dipoles in 1st excited states to vertices connected by photon propagator (left), two outgoing dipoles in 1st excited states to vertices connected by propagator (middle), and one incoming and one outgoing dipole in 1st excited states connected by propagator (right). In the right diagram, we encounter the factor, where we have used the Feynman rules in Figure 1. In α = 1, this factor is nonzero only when β = 1. (In α = 0, it is nonzero only when β = 0). We find that this factor is nonzero only in the cases α = β in the right diagram in Figure 8d. We shall consider the middle diagram. It contains the factor This factor is nonzero only in α = −β. Similarly, we find that the left diagram is nonzero only in α = −β. Using these rules, we find that nonzero contributions are only in α = −γ and γ = −β in the left diagram in Figure 8a, so we find α = β. In the left diagram in Figure 8b, we also find α = β. In self-energy for dipoles in 1st excited states, the label of incoming dipoles α is equal to that of outgoing dipoles β at least in 1-loop and 2-loop order.

Renormalization in 2-Particle-Irreducible Effective Action Technique
In this section, we write 2-Particle-Irreducible (2PI) effective action [50,51,55] with counter-terms in Quantum Brain Dynamics and describe how to cancel ultra-violet diagrams in self-energy by counter-terms in the Kadanoff-Baym equations.
Beginning with the Lagrangian density in QBD in Equation (17), we can derive 2PI effective action. The action Γ 2PI is given by, where we have used the full propagator of photons D ij (x, y) = TA i (x)A j (y) (the brackets representing expectation values), that of dipoles in the ground state G ss (x, y) = Tψ s (x)ψ * s (y) and that of dipoles in 1st excited states G αα (x, y) = Tψ α (x)ψ * α (y) and shown the case of vanishing background fields A i = 0, ψ ( * ) s = 0 and ψ ( * ) α = 0. Here, iD 0,ij , iG −1 0,ss and iG −1 0,αα are written by, We find that counter-terms are in iD 0,ij , iG −1 0,ss and iG −1 0,αα . In 2PI effective action, the term Γ 2 2 represents all the 2-Particle-Irreducible loop diagrams in loop-expansion technique [53]. This term is expressed in Figure 9 for up to O (ed e ) 4 . We find diagrams involving counter-terms. Differentiating 2PI action by propagators, we can derive the Kadanoff-Baym (KB) equations [52][53][54]. By the relation δΓ 2PI δD = 0, we can derive the KB equations for photons as, with Π ij (x, y) ≡ i δΓ 2 δD ji (y,x) . ( Π is not renormalized due to δ Z ph = 0.) By the relation δΓ 2PI δG ss = 0, we derive the KB equations for dipoles in the ground state as, with Σ ss ≡ i 2 δΓ 2 δG ss . In Fourier transformation in this equation, we can use the renormalized self-energy Σ ss,ren as, where we adopt counter-terms derived in perturbation theory in the previous section. In the relation δΓ 2PI δG αα = 0, we derive the KB equations for dipoles in the 1st excited states as, with Σ αα ≡ i 2 δΓ 2 δG αα . Here we can use the renormalized self-energy Σ αα,ren in Fourier transformation as, The right-hand side in this equation is depicted in Feynman diagrams in Figure 10. We adopt the same counter-terms as those in perturbation theory in the previous section.

Discussion
In this paper, we have introduced the Lagrangian density for Quantum Brain Dynamics (QBD) in 3 + 1 dimensions and shown the Feynman rules for propagators, vertices and their counter-terms. We have investigated how to cancel ultra-violet (UV) divergences for 1-loop and 2-loop self-energy and 1-loop 3-point vertex corrections in the perturbation theory in the coupling expansion of ed e . The photon self-energy has no UV divergent terms, at least in 2-loop order. The dipole self-energy involves, at most, cubic-divergent terms in 1-loop order canceled by counter-terms. The vertex corrections for 1-loop 3-point vertices involve, at most, quadratic divergent terms canceled by counter-terms. We can adopt these counter-terms derived in the perturbation theory to cancel UV divergences in the selfenergy emerging in the Kadanoff-Baym (KB) equations derived in the 2-Particle-Irreducible effective action technique.
Using the counter-terms derived in this paper, we can cancel UV divergences in selfenergy for dipoles in lattice simulations of the KB equations in QBD. We set the momenta l µ → 2πn µ 2a s N s with n µ = −N s , −N s + 1, · · ·N s − 1, N s on the lattice with lattice spacing a s . The finite UV cutoff Λ in the integral dl 0 is then ± π a s . We can then calculate finite counterterms in self-energy with a finite cutoff on the lattice. Using them in numerical simulations of time-evolution of KB equations, we might find convergent behaviors in increasing the cutoff on the lattice π a s (decreasing lattice spacing a s ). Juchem et al. [56] have reported that the convergent behaviors appear in increasing cutoff by renormalizing the non-local sunset self-energy in the φ 4 model in 2 + 1 dimensions. Without renormalization, the numerical results of KB equations are dependent on the cutoff in φ 4 theory. The convergent behaviors are also expected to appear in numerical simulations in KB equations in QBD using the counter-terms derived in this paper, although we require numerical costs involving parallel computations to check the convergence since we need to trace the time-evolution of KB equations in 3 + 1 dimensions during equilibration on the lattice with several cutoffs and discuss numerical results in cases with and without renormalization.
We can consider the Ward-Takahashi identity in QBD. We shall consider the case when we do not adopt axial gauge fixing conditions, such as A 0 = 0, but adopt Lorentz invariant gauge fixing conditions. Using electric fields written by E i = −∂ i A 0 + ∂ 0 A 1 and the 4th term in the Lagrangian density in Equation (17), the Feynman rules for vertices in Figure 1 are changed by, and, with subscripts of spatial coordinate i = 1, 2, 3 and subscripts of space-time coordinate µ = 0, 1, 2, 3. We find that inner products with external momenta k µ with the above vertices are zero. Using this property, Feynman diagrams involving inner products of vertices and momenta k µ are zero even if the condition k 2 = 0 is not satisfied so we find the Ward-Takahashi identity in QBD. Using this identity, we find that the longitudinal component with k µ k ν k 2 in propagators for photons does not contribute to physical quantities, such as the S-matrix.
We shall discuss the renormalizability in QBD in 3 + 1 dimensions. We have investigated UV divergences in self-energy in the expansion of dipole moment 2ed e . Up to 2-loop order self-energy, we need not introduce new counter-terms in our Lagrangian density. However, we might need to introduce terms such as ψ * s ψ s A i A i , ψ * s ψ s ψ * α ψ α (α = 0, ±1), and so on, in higher loop order self-energy. Similar problems will occur in Quantum Electrodynamics for non-relativistic charged bosons labeled by field ϕ ( * ) with interaction term ∂Aϕ * ϕ, where we need to introduce additional ϕ 4 interaction terms to renormalize 4-point vertices of ϕ ( * ) appearing in (∂Aϕ * ϕ) 4 . In QBD, the terms such as ψ * s ψ s A i A i and ψ 4 might be introduced. If further additional terms need not emerge, we will arrive at renormalizable theory. Otherwise, the QBD Lagrangian in this paper is non-renormalizable and will represent a low-energy effective theory. We have introduced the hierarchy k ∼ 1 I m Λ with 1 I = 4 meV, mass m = 18 × 940 MeV for water molecules and cutoff Λ. The dipole moment 2ed e times 1 I = 4 meV (or k) is the order 10 −7 due to elementary charge e = 0.3 and d e = 0.2 × 10 −10 m for water molecules. Then the perturbation is achieved for expansion with sufficiently small (2ed e /I) n (n = 2, 3, 4, · · ·). However, if we expand by (2ed e Λ) n and our QBD Lagrangian is non-renormalizable, we find that the order of Λ is 10 4 eV in 2ed e Λ ∼ 1 and that QBD is effective in typical frequency p 0 Λ ∼ 10 4 eV. Then we need to adopt the hierarchy 1 Λ m with expansions by Λ m .) Our results will be applied to trace our decision-making or thinking processes in brain dynamics. Human decision-making is mathematically described in quantum-like theory in quantum cognition adopting quantum superposition and entanglement, which cannot be explained classically [57][58][59][60][61][62], although no physical degrees of freedom are provided in quantum cognition. QBD might provide concrete quantum degrees of freedom in quantum cognition. We adopt the theory of QBD and holography [14]. A single photon has a superposition state of two or more paths, as given by, The photon superposition state propagating through two or more paths will provide a quantum state involving entanglement with holograms. Optical parallel information processing through holograms will provide quantum decision-making involving entanglement, as shown in the above equation. Increasing the number of photons, statistical properties of photons propagating through holograms in water-photon system will emerge. We then encounter quantum many-body problems in a water-photon system to describe our decision-making or thinking processes. Quantum many-body properties of water dipoles and photons can be described by the KB equation in quantum field theory. We can then adopt renormalization in an interacting water-photon system to trace our thinking processes in brain dynamics using KB equations.
In canceling UV divergences emerging in self-energy and vertex corrections by renormalization, we encounter arbitrariness in finite parts in counter-terms in renormalization prescription [48]. In achieving the invariance of the theory in renormalization prescriptions, parameters (m, 1/I, µ c , 2ed e ) will become dependent on the renormalization scale a s (lattice spacing). The contributions from smaller scales than a s are renormalized into those parameters, so that the coarse graining procedures cause group, namely Renormalization Group (RG). The changes in parameters for scales a s are described by the RG equations in QBD. The RG equations describe running parameters of (m, 1/I, µ c , 2ed e ) dependent on the renormalization scale a s . In an increasing scale a s with renormalizing contributions from smaller scales, we might encounter the fixed point in RG equations independent of further coarse graining procedures. The fixed point represents scale-free physical phenomena in QBD. We then find the fractality in water-photon systems representing scale-invariant physics by solving fixed-point solutions in RG equations. Furthermore, when properties larger than the momentum cutoff Λ or those smaller scales than a s (with the momentum cutoff Λ ∼ 1 a s ) are renormalized in QBD theory, we might encounter a new effective theory in macroscopic scales from microscopic water-photon dynamics. The effective theory might provide macroscopic conscious phenomena in brain functions. Renormalization in this work will be extended to describe the diversity in multi-scale brain dynamics involving running parameters, fractality in water-photon systems and so on.
Our paper, following on from a series of previous studies [13][14][15], provides a basis for the proper construction of the Hamiltonian structure of the Quantum Brain Dynamics with respect to a renormalization group application that shows the cancellation of divergent terms. Importantly, the resultant Hamiltonian exhibits fractal self-similar properties that could lead to a better understanding of the quantum degrees of freedom in the human brain in a hierarchical manner. We will explore these repercussions in the context of dipole-dipole interactions within neurons in future work.

Concluding Remarks and Perspectives
We have introduced the Lagrangian density in QBD in 3 + 1 dimensions and shown the counter-terms in the Lagrangian. We have derived counter-terms for the cancellation of ultra-violet (UV) divergences in self-energy and vertex corrections in perturbation theory with the coupling expansion of the dipole moment. Counter-terms will be adopted to cancel UV divergences in numerical simulations of Kadanoff-Baym equations during equilibration. Our analysis will be extended to the renormalization group (RG) method in QBD to investigate running parameters of mass, dipole moment, moment of inertia and chemical potentials and show fractal-like properties for water-photon systems in fixed points in RG equations.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A. Calculations of 1-Loop Self-Energy of Dipoles in 1st Excited States
We shall calculate self-energy in Equation (26). For calculations of the coefficient for self-energy, we have used the relation, When we use p 2 2m p 2 + m 2 − m and the relation, in the denominator in Equation (26), we can rewrite Σ αα by, where we have used Equation (24). In the above equation, the denominator is where we have set, The next step is to expand the numerator in Equation (A3). Using the hierarchy given by, we find, where we have used the relation, We can use, where we have omitted odd power in l 0 . We also use, where we have omitted the odd power of l µ and used the relation l i l j → l 2 3 δ ij satisfied in the integration. Using Equations (A11) and (A12), we find the numerator in Equation (A3) is written by, We use the Wick rotation shown in Figure A1 where the integrand in the complex l 0 curve falls off sufficiently rapidly. Using the rotation with l 0 = il 0 E and l = l E for the Euclidean variable l ν E , we can derive, with, where we neglected terms without UV divergences. Adopting the hierarchy in nonrelativistic cases given by, we can expand ∆ α (p) in Equation (A7) as, and use the relation, Using the above relation and the hierarchy in Equation (A16) for Equation (A15), we arrive at, with, where we set, and, and we have used Feynman's formula given by, The inside of the bracket of denominator 3 is expanded as, where we set q 0 and q as, and define ∆ as, (A32) Next, we calculate numerator 3. Using Equation (A29), we can write, Using Equation (A30), we find the relations, where we have used relation |p|, |k| m |q| and Equation (A10). Similarly, we find the relation, (p + l + k) 2 + m 2 = (q − xp − y(p + k) + p + k) 2 We shall set the renormalization point as, We then expand Equation (A32) in denominator 3 at the point as, ∆ = (x + y)m 2 + y(1 − y) 1 We also expand numerator 3 in Equation (A38) at the renormalization point as, (A42) Due to k 2 = 1 I 2 , we find, (A43) Neglecting the odd power of q 0 and q and removing finite contributions in the integral q , we find, Neglecting the anti-symmetric term for interchange of x and y and using q i q j → |q| 2 δ ij 3 , we find, [numerator 3] → (q 0 ) 2 (q 0 ) 2 + (x + y) 2 m 2 + (y − x) 2 1 I 2 −4(q 0 ) 2 (x + y)m + (y − x) 1 I m − (x + y)m + (x − y) 1 I + |q| +(q 0 ) 2 |q| 2 + 2|q|(1 − x − y)m + y 2 + (1 − y) 2 3I 2 + m 2 +|q| 2 (x + y) 2 m 2 + (y − x) 2 1 I 2 As a result, we arrive at, [numerator 3] → (q 0 ) 2 (q 0 ) 2 + |q| 2 + 6m 2 (x + y)(x + y − 1) + 2m|q|(1 − 3(x + y)) + 2m 2 Finally, we perform the Wick rotation q 0 = iq E , q = q E and adopt the following procedure for Equation (A41), (A47) Using the above relation and Equation (A46), Equation (A24) can be rewritten as,