Extended Harmonic Map Equations and the Chaotic Soliton Solutions Gang Ren,* and Yi-Shi Duan,* The Molecular Foundry, Lawrence erkeley National Laboratory, erkeley C 9470, US The Department of Physics, Lanzhou University, Lanzhou, Gansu 730000, China * To whom correspondence should be addressed: G. Ren, Lawrence erkeley National Laboratory, the Molecular Foundry, Cyclotron Road, erkeley, C 9470, US, E-mail: gren@lbl.gov or Y.S. Duan, The Department of Physics, Lanzhou University, Lanzhou, Gansu 730000, China, E-mail: ysduan@lzu.edu.cn Key Words: Extended harmonic maps, Chaotic soliton solutions bstract In this paper, the theory of harmonic maps is extended. The soliton or travelling wave solutions of Euler's equations of the extended harmonic maps are studied. In certain cases, the chaotic behaviours of these partial equations can be found for the particular case of the metrics and the potential functions of the extended harmonic equations. Introduction ll physical processes are described by differential equations. Over the last several decades, many nonlinear ordinary differential equations have been found to exhibit chaotic behaviours, such as the nonlinear equation for a parametrically excited pendulum given by d x dx k [ Cos( t] Sin( x 0 ( dt dt where k, α and β are control parameters, and the equation for a harmonic system in periodic fields given by d x dx 3 k x x bcos( t ( dt dt where k, α and β are control parameters. There is also the Lorentz model 3, dx dy dz ( x y; xy x y; xy kx (3 dt dt dt where k, α and β are control parameters. Over the last several decades, many nonlinear ordinary differential equations possessing chaotic behaviours have been discovered 4. To date, the study of chaotic behaviours using the space distribution of nonlinear partial differential equations has been rare and is currently only in its initial stage. In this paper, we use the extended theory of harmonic maps to show that chaotic behaviours exist for a certain class of nonlinear partial differential equations. Three decades ago, Duan proposed a method to solve the Euler's equations of harmonic maps (HM 5, in which the solution of the Euler's equations of HM was found by solving the
Laplace-eltrami s equation and the geodesic equations. Using this method, the general solutions of Einstein equations in general relativity have been discussed. From the viewpoint of physics, the traditional theory of harmonic maps contains only the kinetic energy term. For a wider application of this method, the Lagrangian of the traditional HM theory must be supplemented with a potential energy term 6. This extended HM theory is helpful for studying the travelling wave or soliton solutions of some types of nonlinear partial differential equations. Here, we first establish an extended theory of harmonic maps and discuss the travelling wave or soliton solution of the Euler's equations of the extended HM. We then use the extended theory of HM and study two special cases for which the chaotic behaviours of some types of partial differential equations are found. The theory and method could be helpful for studying the existence of soliton solutions and the chaotic behaviours of more types of partial equations, possibly providing a further link between the extended HM theory and quantum entanglement via common classical chaos 7. The Extended Theory of Harmonic Maps The theory of harmonic maps 8-0 became an important branch of mathematical physics many years ago and has been applied to a wide variety of problems in mathematics and theoretical physics. In this section, we first describe the formulation of the extended theory of harmonic maps. Let M and N be two Riemannian manifolds with local coordinates x µ (µ =,,, m on M and local coordinates Φ (=l,,, n on N. The metrics on M and N are denoted by respectively. mapping dl g ( x dx dx dl G ( d d dim(m=m : M N x dim(n=n is called an extended harmonic map if it satisfies the Euler's equation resulting from the variational principle I 0, using the action I n d x g [ g G ( V ( ] (4 where x, V( Φ is the potential function of Φ, and n V ( V (,,...,. It is obvious that the traditional HM corresponds to the case of V ( 0. The conditions for a map to be harmonic are given by the Euler's equations L L 0 = l,,, n (5
where L g G ( g V ( g. (6 y substituting (6 into (5, we can obtain the Euler's equations of extended HM given by C V ( ( gg g C G =0 (7 g where C are the Christoffels symbols on manifold N given by D GD GCD GC C G [ ]. (8 C D To study a special type of solution of partial differential equation (7, following reference 5, we investigate the case where (=l,,, n are functions solely of the argument σ, and σ is a function of x on the manifold M: Φ = Φ (σ (9 σ = σ (x; 0 y substituting (9 into (7, equation (7 can be written as g ( g g d G d V ( d ( d C d d C d g d 0 ( It is easy to see that if the function ( x satisfies the Laplace-eltrami equations g ( g g 0 ( equation ( will take the form d d C d d V ( C G d d f (3 where f g (4 is a scalar function on the manifold M. ecause f is not the function of V ( fu (, equation (3 can also be written as
d d C d d U( C G d d (5 Equation (5 can viewed as the geodesic equation of a particle on the Riemannian manifold N subjected to an external force F G U ( (6 This will become equivalent to the usual geodesic equations on the manifold N if the external force F 0 (=,,..., n. more general form of equation (5 can be written because the existence and uniqueness of the solution of the above equation has been discussed. Therefore, generally, a solution of (5 exists when the initial conditions are given by 0 and d d. 0 From (, ( and (4, we see that if Φ satisfies the geodesic equations with an external force on the manifold N, i.e., (3, and σ = σ (x satisfies the Laplace-eltrami equations on the manifold M (, the functions ( ( x (=l,,, n will be the special solutions of the Euler's equations of the extended HM (7. In the following, we use instead of g if the manifold M is the pseudo-riemannian space-time. Now, we show that there exist travelling wave or soliton solutions for partial differential equations (7. It can be proven that a solution of the Laplace-eltrami equation ( will be kx (7 0 where k (µ =0,,, m- are components of a constant vector and x t is time. Substituting (7 into (, we find g g ( gg k 0 (8 and choosing Fock's harmonic coordinates that always exist in pseudo-riemannian M, i.e., g ( g g 0, (9 equation (8 is evidently satisfied. This means that ( and that the scalar function f defined by (3 is of the form k x is indeed a solution of equation f g kk. (0 Therefore, the solution of the equations of motion (5 on the manifold N with given initial conditions can be formally expressed as
( ( k x ( k x k0t (=,,,n ( The above expression shows that the partial differential equations of the extended HM (7 in pseudo-riemannian space-time possess solutions. where In the case that the manifold M is the pseudo-euclidean space-time g is a Lorentz metric, equation ( is simplified to 0 ( which has an obvious solution k x, and the function f defined by (0 is a constant f kk (3 In this case, the extended HM equation (7 becomes C V ( G 0 (4 C Chaotic soliton solutions of the extended HM equations In the previous section, it has been noted that there exist travelling wave or soliton solutions of the partial differential equations of the extended HM (7 in the pseudo-riemannian space-time. The form (shape of travelling wave or soliton ( is determined by the solution of nonlinear ordinary differential equations (4. In the following, we study the simple case that the M is the pseudo-euclidean space-time and N is a -dimensional manifold with coordinates {, }. s in the above section, we suppose that Φ and Φ are functions of an argument, i.e., x k x ( and further assume that ( ; ; k x (5 Then, equation (4 can be expressed as d d ( d d d ( d ( G f V (, G V (, 0 (6
d d V (, V (, ( ( ( G G 0 d d f. (7 Here, we have denoted the potential function as V (, V(, d Eliminating the term ( from the differential equations (6 and (7, we obtain d d d ( ( d d [ G f V (, G V (, ( G V (, G V (, ] 0 (8 To make the equation have the same form as equation (, we assume that the metrics on manifold N are diagonal and take the following form: G e ; G ( k e ; G G 0 (9 where k is a constant, and the potential function V(Φ, σ takes the form V (, fe u(, (30 3 3 3 3 u(, [ ] ( 4 4 b b Cos( Sin( 4 4 (3 Substituting (8, (9, (30 and (3 into equation (8, we obtain d k d d d 3 bcos( 0; k x (3 Comparing the above equation for with equation (, we find that they are of the same form. It has been shown that equation (, i.e., the equation of a harmonic system in a periodic field, possesses chaotic states that are characterized by the existence of a strange attractor in the phase space. This implies that equation (3 for ( should possess the same chaotic states as (. ecause equation (3 is deduced from the special case of equation (4, it determines the form (shape of soliton ( ( kx, as discussed in the previous section. Therefore, the chaotic states of equation (3 mean that the form (shape of solution determined by this equation exhibits chaotic behaviours, and we call this type of solution of (3 the chaotic soliton.
Substituting (8, (9, (30, (3 and ( ; k x two partial differential equation for. Eliminating the quadric term two equations, we have the partial differential equation into (4, we find from these 3 k k f [ bcos( k x ] 0 (33 If we let ( and k x and substitute these equations into (33, it is easy to verify that equation (3 can be directly derived from (33. This means that equation (33 has the soliton solution ( k x and that equation (3 determines the form of the soliton (. From the above discussion, it is obvious that the partial differential equation (33 possesses a special type of a chaotic soliton solution. In the following, we study another example, for which the function u (, in the potential function (9 is chosen as u(, Sin( Cos( [ 5 5 4 ( [ 4 ( [ 4 ( [ 4 ( ] Sin( Cos( 4 ( ] Cos( Sin( 4 ( ] Cos( Cos( 4 ( ] Sin( Sin( ] 4 ( (34 Substituting (34 into (8, we obtain the following equation corresponding to equation (3: d d k [ Cos( ]sin( 0 d d (35 This is merely the nonlinear equation of a parametrically excited damped pendulum (. Reference 4 has explored in detail the chaotic behaviours of this equation. The partial differential equation corresponding to (35 can be obtained using (4, (9, (30 and (34. That is, k k f [ Cos( k x ] Sin( 0 (36 s discussed in the first example, we know that this partial differential equation also has chaotic soliton solutions. The chaotic behaviour of such solitons of equation (36 is determined by the time-evolution curve of the amplitude of the parametric pendulum model.
The two examples discussed above are both nonlinear damped oscillator systems that are the special cases of the Duffing equation 4. d d k [ p( g(, ] 0 d d. (37 Here, g (, is a function of and, which is a periodic function of. p ( is a nonlinear function of. In this general case, the function u (, in (3 is the solution of the following equation u(, u(, u (, p( g(, (38 Using (38, we can derive the following partial differential equation with chaotic soliton solutions k k [ p( g(, ] 0 (39 with the form (shape of the solitons determined by the time-evolution curve of the amplitude described by the diffusion equation (37 Conclusion The theory presented in this paper is helpful for studying the existence of the soliton and chaotic soliton solutions for certain types of nonlinear partial differential equations. Using the Lorentz model 3, which is the standard model for the study of chaos theory, we have deduced a system of nonlinear partial differential equations. These equations possess chaotic travelling wave or soliton solutions with the space distribution. Reference Leven, R. W. & Koch,. P. Chaotic ehavior of a Parametrically Excited Damped Pendulum. Phys Lett 86, 7 74, doi:doi 0.06/0375 960(89067 5 (98. Huberman,.. & Crutchfield, J. P. Chaotic States of nharmonic Systems in Periodic Fields. Phys Rev Lett 43, 743 747, doi:doi 0.03/PhysRevLett.43.743 (979. 3 Lorentz, E. N. Deterministic nonperiodic flow. Journal of the. tmospheric Sciences 0 (963. 4 Hao,. L. ifurcation, chaos, strange attractor, turbulence and all that. Progress in Physics 3 (983. 5 Duan, Y. S. Solution of the Euler's equation of Harmonic Maps and its pplication to General Relativity. cta Phys. Sinica. 33 (984. 6 Duan, Y. S. Extended Harmonic Maps Equations and the Soliton Solutions. Journal of Lanzhou University (Natural Sciences S (99. 7 Neill, C. et al. Ergodic dynamics and thermalization in an isolated quantum system. Nature Physics, doi:0.038/nphys3830 (06. 8 Misner, C. W. Harmonic Maps as Models for Physical Theories. Physical Review D 8 (978.
9 Nutku, Y. & Halil, M. Colliding Impulsive Gravitational Waves. Phys Rev Lett 39, 379 38, doi:doi 0.03/PhysRevLett.39.379 (977. 0 Eris,. & Nutku, Y. Harmonic Mappings of Riemannian Manifolds and Stationary Vacuum Space Times with Whole Cylinder Symmetry. J Math Phys 6, 43 434 (975. Milnor, J. Morse Theory. (Princeton University, 963. Weinberg, S. Gravitation and Cosmology: Principles and pplications of the General Theory of Relativity. (John. Wiley, 97. cknowledgments: We thank Mr. Wencheng Ren to preserve the original results and manuscript. Work at the Molecular Foundry was supported by the Office of Science, the Office of asic Energy Sciences, and the U.S. Department of Energy, under Contract No. DE-C0-05CH3. uthor Contributions: This project was initiated by GR and YD., GR calculated the solution, and YS validated the solution. GR drafted the initial manuscript, which was revised by YD. dditional Information Competing financial interests The author(s declare no competing financial interests.