Discretized Classical Model of the New Equation
- SourceCodeOf_HumanGenome > Theory of Quantum History Entangled in Time-like Direction @ 2008/12/6 11:52
- SourceCodeOf_HumanGenome > Probability Formula in New Grammar @ 2013/5/18 13:32
- SourceCodeOf_HumanGenome > Examine Probability Formula @ 2013/5/28 16:29
- SourceCodeOf_HumanGenome > Twist Remove Normalization @ 2013/5/27 11:23
- SourceCodeOf_HumanGenome > Line Integral in Function Space @ 2013/5/26 16:32
- SourceCodeOf_HumanGenome > New Grammar Version of Schrödinger Equation @ 2011/8/14 14:51
- SourceCodeOf_HumanGenome > Circular Time Regularization @ 2011/8/24 10:54
- SourceCodeOf_HumanGenome > Solution by using circular time in the case V=0 @ 2011/10/12 10:38
- SourceCodeOf_HumanGenome > Similar problem in the old theories @ 2011/10/26 15:44
- SourceCodeOf_HumanGenome > Unentangled Fundamental Solution @ 2011/8/22 15:32
- SourceCodeOf_HumanGenome > Reduced Form of the New Equation @ 2011/8/16 14:58
- » SourceCodeOf_HumanGenome > Discretized Classical Model of the New Equation @ 2014/5/13 17:31
- SourceCodeOf_HumanGenome > Cartoon of the Equation @ 2013/4/3 14:39
- SourceCodeOf_HumanGenome > Meta Gauge Field @ 2011/8/21 10:32
SourceCodeOf_HumanGenome > Discretized Classical Model of the New Equation @ 2014/5/13 17:31 |
---|
In this post, I try to make n=∞
straight time version of Cartoon of the Equation. Below I solve classically the equation of motion of the mechanical system defined by the following Hamiltonian. H(p; x) = Qφ(x) + [1／(2M)] Σ_{k=－∞}^{∞} [p^{k} － QA^{k}(x)]^{2}. M = mεα^{2}, Q = m. A^{k}(x) = α(x^{k+1} － x^{k}), φ(x) = －[1／(2ε)] Σ_{k=－∞}^{∞} (x^{k+1} － x^{k})^{2}. Derivation of this Hamiltonian is the same as one of Cartoon of the Equation. Hamilton's canonical equations of motion of this system is as follows. dx^{j}／dt = ∂H／∂p^{j} = (p^{j} － QA^{j})／M, dp^{j}／dt = －∂H／∂x^{j} = (Q／M)Σ_{k} (p^{k} － QA^{k})∂A^{k}／∂x^{j} － Q∂φ／∂x^{j}. From the first of these equations, it is derived that p^{j} = Mdx^{j}／dt + QA^{j}. Substituting this into the second of the equations, it is derived that M(d／dt)^{2} x^{j} + QdA^{j}／dt = (Q／M)Σ_{k} (p^{k} － QA^{k})∂A^{k}／∂x^{j} － Q∂φ／∂x^{j} = (Q／M)Σ_{k} (Mdx^{k}／dt)∂A^{k}／∂x^{j} － Q∂φ／∂x^{j}. For the term dA^{j}／dt in the left hand side of this equation, substituting that dA^{j}／dt = Σ_{k} (∂A^{j}／∂x^{k})dx^{k}／dt, it is derived that M(d／dt)^{2} x^{j} = (Q／M)Σ_{k} (Mdx^{k}／dt)∂A^{k}／∂x^{j} － QΣ_{k} (∂A^{j}／∂x^{k})dx^{k}／dt － Q∂φ／∂x^{j}. ∴ M(d／dt)^{2} x^{j} = QΣ_{k} (dx^{k}／dt)(∂A^{k}／∂x^{j} －∂A^{j}／∂x^{k}) － Q∂φ／∂x^{j}. Here I use the concrete expressions of A and φ. ∂A^{k}／∂x^{j} = α(δ_{k+1,j} － δ_{k,j}), ∂A^{j}／∂x^{k} = α(δ_{j+1,k} － δ_{j,k}) = α(δ_{k-1,j} － δ_{k,j}) ∴ ∂A^{k}／∂x^{j} － ∂A^{j}／∂x^{k} = α(δ_{k+1,j} － δ_{k-1,j}) = α(δ_{j-1,k} － δ_{j+1,k}). ∂φ／∂x^{j} = －(1／ε)[(x^{j} － x^{j-1}) － (x^{j+1} － x^{j})] = －(2／ε)[x^{j} － (1／2)(x^{j-1} + x^{j+1})]. Substituting these expressions into the above result, we get the concrete form of the equation of motion of the main problem of this post as follows. M(d／dt)^{2} x^{j} = Qα(dx^{j-1}／dt － dx^{j+1}／dt) + Q(2／ε)[x^{j} － (1/2)(x^{j-1} + x^{j+1})]. Substituting mεα^{2} and m for M and Q in this equation, we get the following expression as the final form of the equation of motion. (d／dt)^{2} x^{j} = [1／(εα)](dx^{j-1}／dt － dx^{j+1}／dt) + [2／(ε^{2}α^{2})][x^{j} － (1/2)(x^{j-1} + x^{j+1})]. This equation is the opposite of the following equation representing, as it were, an automatic collision avoidance system. M(d／dt)^{2} x^{j} = Qα[(dx^{j-1}／dt － dx^{j}／dt) + (dx^{j+1}／dt － dx^{j}／dt)] － Q(2／ε)[x^{j} － (1/2)(x^{j-1} + x^{j+1})]. --- --- Multiplying both sides of the equation by dx^{j}／dt and summing over j, we get the following equation representing the energy conservation law. (d／dt)Σ_{a≦j≦b}(M／2)(dx^{j}／dt)^{2} = Qα[(dx^{a-1}／dt)(dx^{a}／dt) － (dx^{b}／dt)(dx^{b+1}／dt)] + (Q／ε)(d／dt)Σ_{a≦j≦b}[(x^{j})^{2} － x^{j}x^{j+1}] － (Q／ε)(x^{a-1} dx^{a}／dt － x^{b} dx^{b+1}／dt). ∴ (d／dt)Σ_{a≦j≦b}[(M／2)(dx^{j}／dt)^{2} － (Q／ε)x^{j}(x^{j} － x^{j+1})] = Qα[(dx^{a-1}／dt)(dx^{a}／dt) － (dx^{b}／dt)(dx^{b+1}／dt)] － (Q／ε)(x^{a-1} dx^{a}／dt － x^{b} dx^{b+1}／dt). ∴ (d／dt)Σ_{a≦j≦b}[(M／2)(dx^{j}／dt)^{2} － (Q／ε)x^{j}(x^{j} － x^{j+1})] = Q[α(dx^{a-1}／dt)(dx^{a}／dt) － (1／ε)x^{a-1}(dx^{a}／dt)] － Q[α(dx^{b}／dt)(dx^{b+1}／dt) － (1／ε)x^{b}(dx^{b+1}／dt)]. The right hand side of this equation represents the time derivative of energy possessed by the part such that a≦j≦b, and the left hand side of this equation represents (energy flow at j=a) - (energy flow at j=b). Therefore this equation represents the energy conservation law. --- --- To investigate the character of this problem briefly, now I omit the scalar potential term and try solving the following equation. (d／dt)^{2} x^{j} = [1／(εα)](dx^{j-1}／dt － dx^{j+1}／dt). Introducing a Matrix A defined as A_{jk} = δ_{j-1,k} － δ_{j+1,k}, the equation is rewritten as (d／dt)^{2} x^{j} = [1／(εα)]Σ_{k} A_{jk} dx^{k}／dt. Then the solution is writen as follows. dx／dt = a exp[tA/(εα)] v, x = aεα exp[tA／(εα)] u + w, where a is a constant real number, and u, v, w are infinite dimensional constant real column vectors such that v = Au. --- --- From the point of view of energy conservation, velocity does not increase or decrease as a whole. This situation is expected to be represented well by dx／dt = a exp(it × somematrix) v, but the above solution does not include the imaginary unit at all. Perhaps this superficial contradiction is cancelled by the fact that the matrix A has a property analogous to the property of a pure imaginary number. As for the case of a 2×2 matrix, the matrix A defined by the equation: A_{jk} = δ_{j-1,k} - δ_{j+1,k} is actually a representation of the imaginary unit. We can fully express this understanding by redefining the matrix A' as A'_{jk} = －iA_{jk} and rewriting as dx／dt = a exp[itA'／(εα)] v. Here A' is Hermitian. We know well from quantum mechanics that this represents a rotation in an infinite dimensional space. --- --- If A^{-1} exists, u = A^{-1} v. However, perhaps A^{-1} does not exist, because there is an invariant vector under the transformation exp[tA／(εα)] so zero is an eigenvalue of A so the determinant of A is zero. This is understood by considering about A' = -iA. A' is Hermitian so it is diagonalized to the form whose all diagonal elements are eigenvalues of A' so the determinant of A' is the product of all its eigenvalues. Therefore, if zero is an eigenvalue of A' then the determinant of A' is zero. The eigenvalues of A are completely same as the eigenvalues of A'. The determinant of A is different from one of A' only by a factor, so if zero is an eigenvalue of A then the determinant of A is zero. It is shown below that there is an invariant vector under the transformation exp[tA／(εα)]. The space of all classical solutions of the main equation of this post which is an equation of motion is { x | ∃ξ,η; ∀j; x^{j} = jξ + η }. Every vector in this space is written in the form: a ∂x/∂ξ + b ∂x/∂η. The j-th component of this vector is (a ∂x/∂ξ + b ∂x/∂η)^{j} = aj + b. This individual vector is not necessarily invariant under exp[tA/(εα)], but the space which consists of all such vectors is invariant under exp[tA/(εα)]. I prove it below. [A(a ∂x/∂ξ + b ∂x/∂η)]_{j} = Σ_{k} A_{jk} (a ∂x/∂ξ + b ∂x/∂η)^{k} = Σ_{k} (δ_{j-1,k} - δ_{j+1,k})(ak + b) = [(j-1)a + b] - [(j+1)a + b] = -2a. (So, the vector is invariant if a = 0.) [A^{2} (a ∂x/∂ξ + b ∂x/∂η)]_{j} = Σ_{k} A_{jk} (-2a) = Σ_{k} (δ_{j-1,k} - δ_{j+1,k})(-2a) = 0. A^{n} (a ∂x/∂ξ + b ∂x/∂η) = A^{n-2} A^{2} (a ∂x/∂ξ + b ∂x/∂η) = A^{n-2} 0 = 0 (n≧3). {exp[tA/(εα)](a ∂x/∂ξ + b ∂x/∂η)}_{j} = {[1 + tA/(εα)](a ∂x/∂ξ + b ∂x/∂η)}_{j} = (aj + b) + [t/(εα)](-2a) = aj + [b - 2at/(εα)]. ∴ exp[tA/(εα)](a ∂x/∂ξ + b ∂x/∂η) = a ∂x/∂ξ + [b - 2at/(εα)]∂x/∂η ∈ the original space. This justifies our thinking of exp[tA/(εα)] as a rotation around the space of all classical solutions as an invariant axis. --- --- To write down exp[tA/(εα)] explicitly, below I try to get explicit expressions of all powers of A. (A^{2})_{jk} = Σ_{s}A_{j,s}A_{s,k} = Σ_{s}(δ_{j-1,s} - δ_{j+1,s})(δ_{s-1,k} - δ_{s+1,k}) = δ_{j-2,k} - 2δ_{j,k} + δ_{j+2,k} (A^{3})_{jk} = Σ_{s}(δ_{j-2,s} - 2δ_{j,s} + δ_{j+2,s})(δ_{s-1,k} - δ_{s+1,k}) = δ_{j-3,k} - 3δ_{j-1,k} + 3δ_{j+1,k} - δ_{j+3,k} (A^{4})_{jk} = Σ_{s}(δ_{j-3,s} - 3δ_{j-1,s} + 3δ_{j+1,s} - δ_{j+3,s})(δ_{s-1,k} - δ_{s+1,k}) = δ_{j-4,k} - 4δ_{j-2,k} + 6δ_{j,k} - 4δ_{j+2,k} + δ_{j+4,k} (A^{5})_{jk} = Σ_{s}(δ_{j-4,s} - 4δ_{j-2,s} + 6δ_{j,s} - 4δ_{j+2,s} + δ_{j+4,s})(δ_{s-1,k} - δ_{s+1,k}) = δ_{j-5,k} - 5δ_{j-3,k} + 10δ_{j-1,k} - 10δ_{j+1,k} + 5δ_{j+3,k} - δ_{j+5,k} (A^{n})_{jk} = Σ_{0≦s≦n/2}(-1)^{s}c_{1}(s,n)δ_{j-n+2s,k} + Σ_{0≦s＜n/2}(-1)^{n+s}c_{2}(s,n)δ_{j+n-2s,k}, c_{r}(s,n+1) = c_{r}(s,n) + c_{r}(s-1,n), c_{r}(0,n) = 1. Using these results in the expression: dx／dt = a exp[tA/(εα)] v, it is derived that dx^{j}／dt = Σ_{n=0}^{∞} (1/n!)[t／(εα)]^{n}[Σ_{0≦s≦n/2} (-1)^{s} c_{1}(s,n) δ_{j-n+2s,0} - Σ_{0≦s＜n/2} (-1)^{s} c_{2}(s,n) δ_{j+n-2s,0}] in the case that a = 1 and v_{j} = δ_{j0}. I can not go further because the recurrence relation for c can not be solved fully. It seems to be impossible to write down fully c_{1}(s,n) and c_{2}(s,n) as functions of s and n. --- --- To take the scalar potential term into account, I rewrite the main equation of this post to the following form. (d／dt + F)(d／dt + G)x = 0. By comparing this expression with the original equation, the condition to determine F and G is found to be as follows. (F + G)_{jk} = [1／(εα)](δ_{j+1,k} － δ_{j－1,k}), (FG)_{jk} = [1／(ε^{2}α^{2})](δ_{j－1,k} + δ_{j+1,k} － 2δ_{jk}). Then the solution is written as follows. (d／dt + G)x = b exp(－tF) v_{1}, x = exp(－tG)[v_{1} + b ∫_{c}^{t}dt' exp(t'G)exp(－t'F)v_{2}]. Below I try to know F and G concretely. F_{jk} = Σ_{s=－∞}^{∞} f_{s} δ_{j+s,k}, G_{jk} = Σ_{s=－∞}^{∞} g_{s} δ_{j+s,k}. (FG)_{jk} = Σ_{r} F_{jr} G_{rk} = Σ_{r,s1,s2} f_{s1} δ_{j+s1,r} g_{s2} δ_{r+s2,k} = Σ_{s1,s2} f_{s1} g_{s2} δ_{j+s1+s2,k} The condition for FG means that (FG)_{jk} = 0 if |j－k| ≧ 2. So f and g must satisfy the following condition. |r| ≧ 2 ⇒ Σ_{s1+s2=r} f_{s1} g_{s2} = 0, ∴ |r| ≧ 2 ⇒ Σ_{s} f_{s} g_{r－s} = 0. The condition for FG means that (FG)_{jk} = －2／(ε^{2}α^{2}) if j = k. So f and g must satisfy the condition: Σ_{s} f_{s} g_{－s} = －2／(ε^{2}α^{2}) derived as the case that r = 0. The condition for FG means that (FG)_{jk} = 1／(ε^{2}α^{2}) if j+1 = k. So f and g must satisfy the condition: Σ_{s} f_{s} g_{1－s} = 1／(ε^{2}α^{2}) derived as the case that r = 1. The condition for FG means that (FG)_{jk} = 1／(ε^{2}α^{2}) if j－1 = k. So f and g must satisfy the condition: Σ_{s} f_{s} g_{－1－s} = 1／(ε^{2}α^{2}) derived as the case that r = －1. From the condition for F + G, the following conditions for f and g are derived. f_{1} + g_{1} = 1／(εα), f_{－1} + g_{－1} = －1／(εα), not(s = 1 or s = －1) ⇒ f_{s} + g_{s} = 0. Using the condition derived from the condition for F + G, now I eliminate f in the condition derived from the condition for FG. |r| ≧ 2 ⇒ f_{1} g_{r－1} + f_{－1} g_{r+1} + Σ_{1≠s≠－1} f_{s} g_{r－s} = 0. ∴ |r| ≧ 2 ⇒ [1／(εα) － g_{1}]g_{r－1} + [－1／(εα) － g_{－1}]g_{r+1} + Σ_{1≠s≠－1} (－g_{s})g_{r－s} = 0. ∴ |r| ≧ 2 ⇒ [1／(εα)]g_{r－1} － [1／(εα)]g_{r+1} = Σ_{s} g_{s} g_{r－s}. As for the part as r = 0, f_{1} g_{－1} + f_{－1} g_{1} + Σ_{1≠s≠－1} f_{s} g_{－s} = －2／(ε^{2}α^{2}). ∴ [1／(εα) － g_{1}]g_{－1} + [－1／(εα) － g_{－1}]g_{1} + Σ_{1≠s≠－1} (－g_{s})g_{－s} = －2／(ε^{2}α^{2}). ∴ [1／(εα)]g_{－1} － [1／(εα)]g_{1} + 2／(ε^{2}α^{2}) = Σ_{s} g_{s} g_{－s}. As for the part as r = 1, f_{1} g_{1－1} + f_{－1} g_{1+1} + Σ_{1≠s≠－1} f_{s} g_{1－s} = 1／(ε^{2}α^{2}). ∴ [1／(εα) － g_{1}]g_{0} + [－1／(εα) － g_{－1}]g_{2} + Σ_{1≠s≠－1} (－g_{s})g_{1－s} = 1／(ε^{2}α^{2}). ∴ [1／(εα)]g_{0} － [1／(εα)]g_{2} － 1／(ε^{2}α^{2}) = Σ_{s} g_{s} g_{1-s}. As for the part as r = －1, f_{1} g_{－1－1} + f_{－1} g_{－1+1} + Σ_{1≠s≠－1} f_{s} g_{－1－s} = 1／(ε^{2}α^{2}). [1／(εα) － g_{1}]g_{－2} + [－1／(εα) － g_{－1}]g_{0} + Σ_{1≠s≠－1} (－g_{s})g_{－1－s} = 1／(ε^{2}α^{2}). [1／(εα)]g_{－2} － [1／(εα)]g_{0} － 1／(ε^{2}α^{2}) = Σ_{s} g_{s} g_{－1－s}. Therefore, the complete set of conditions for determine g are as follows. |r| ≧ 2 ⇒ [1／(εα)]g_{r－1} － [1／(εα)]g_{r+1} = Σ_{s} g_{s} g_{r－s}, [1／(εα)]g_{－1} － [1／(εα)]g_{1} + 2／(ε^{2}α^{2}) = Σ_{s} g_{s} g_{－s}, [1／(εα)]g_{0} － [1／(εα)]g_{2} － 1／(ε^{2}α^{2}) = Σ_{s} g_{s} g_{1－s}, [1／(εα)]g_{－2} － [1／(εα)]g_{0} － 1／(ε^{2}α^{2}) = Σ_{s} g_{s} g_{－1－s}. After determining g from these conditions, f is determined from the condition derived from the condition for F + G. The number of equations for determining g is equal to the number of elements of g, and so perhaps the solution of the set of equations exists, but it is too difficult to be solved. --- --- Though I can not write the exact solution because of such a difficulty, perturbation method saves us from giving up knowing further. In convenience, here I paste again the main equation to be solved. (d／dt)^{2} x^{j} = [1／(εα)](dx^{j-1}／dt － dx^{j+1}／dt) + [1／(ε^{2}α^{2})](2x^{j} － x^{j-1} － x^{j+1}). To apply perturbation method to this equation, I replace the coefficient of the scalar potential part by a new constant k. (d／dt)^{2} x^{j} = [1／(εα)](dx^{j-1}／dt － dx^{j+1}／dt) + k(2x^{j} － x^{j-1} － x^{j+1}). Then I substitute a power expansion of x in terms of k: x^{j} = Σ_{n=0}^{∞} k^{n} x_{n}^{j} into the equation, and I equate both sides order by order as usual. Below I try calculating the functions corresponding to the first two orders. (d／dt)^{2} x_{0}^{j} = [1／(εα)](dx_{0}^{j-1}／dt － dx_{0}^{j+1}／dt). (d／dt)^{2} x_{1}^{j} = [1／(εα)](dx_{1}^{j-1}／dt － dx_{1}^{j+1}／dt) + (2x_{0}^{j} － x_{0}^{j-1} － x_{0}^{j+1}). x_{0} has already been derived above when I omitted the scalar potential term. From the second equation, it is derived that (d／dt)^{2} x_{1} = [1／(εα)] A dx_{1}／dt + B x_{0}, where A_{jk} = δ_{j-1,k} － δ_{j+1,k}, B_{jk} = 2δ_{jk} － δ_{j-1,k} － δ_{j+1,k}. Solving this equation, we get the expression for dx_{1}／dt: dx_{1}／dt = exp[tA／(εα)]∫_{c}^{t}dt' exp[－t'A／(εα)]Bx_{0}(t'). Substituting x_{0} into this expression, we get the following expression. dx_{1}／dt = exp[tA／(εα)]∫_{c}^{t}dt' exp[－t'A／(εα)]B{aεα exp[t'A／(εα)] u + w}. To reduce this expression, here I calculate the commutation relation between A and B. [A, B]_{jk} = Σ_{s} A_{js}B_{sk} － Σ_{s} B_{js}A_{sk} = Σ_{s}(δ_{j-1,s} － δ_{j+1,s})(2δ_{sk} － δ_{s-1,k} － δ_{s+1,k}) － Σ_{s}(2δ_{js} － δ_{j-1,s} － δ_{j+1,s})(δ_{s-1,k} － δ_{s+1,k}) = (2δ_{j-1,k} － δ_{j-2,k} － δ_{j,k} － 2δ_{j+1,k} + δ_{j,k} + δ_{j+2,k}) － (2δ_{j-1,k} － 2δ_{j+1,k} － δ_{j-2,k} + δ_{j,k} － δ_{j,k} + δ_{j+2,k}) = 0. Because of this vanishment, exp[－t'A／(εα)] B exp[t'A／(εα)]= B, so the expression for dx_{1}／dt is reduced to the following form. dx_{1}／dt = exp[tA／(εα)]∫_{c}^{t}dt' {aεαBu + exp[－t'A／(εα)]Bw} = exp[tA／(εα)]{aεα(t － c)Bu － εα exp[－tA／(εα)]w' + εα exp[－cA／(εα)]w'} = εα{a(t － c)exp[tA／(εα)]Bu － w' + exp[(t － c)A／(εα)]w'}. where w' and w are chosen such that Aw' = Bw. Let's consider the special case such that w = w' = 0 to sketch roughly the effect of scalar potential term. (d／dt)(x_{0} + kx_{1}) = a exp[tA／(εα)]v + kεαa(t － c)exp[tA／(εα)]Bu = a exp[tA／(εα)][v + kεα(t － c)Bu]. This means that a gradually varying vector, which is a linear function of t, is transformed by rotation with some constant angular velocity. --- Last edited at 2014/06/09/17:20JST |