Next Article in Journal
Symmetries of the Large Scale Structures of the Universe as a Phenomenology of a Fractal Turbulence: The Role of the Plasma Component
Previous Article in Journal
Computational Design of an In-Line Coaxial-to-Circular Waveguide Adapter with More Than an Octave Bandwidth
Previous Article in Special Issue
Highly Dispersive Optical Solitons in the Absence of Self-Phase Modulation by Lie Symmetry
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Canonical Equations of Hamilton with Symmetry and Their Applications

1
Guangdong Provincial Key Laboratory of Nanophotonic Functional Materials and Devices, South China Normal University, Guangzhou 510631, China
2
School of Electrical & Electronic Engineering, Shangqiu Normal University, Shangqiu 476000, China
*
Author to whom correspondence should be addressed.
Symmetry 2024, 16(3), 305; https://doi.org/10.3390/sym16030305
Submission received: 3 January 2024 / Revised: 22 January 2024 / Accepted: 8 February 2024 / Published: 5 March 2024
(This article belongs to the Special Issue Nonlinear Schrödinger Equations and Symmetry)

Abstract

:
Two systems of mathematical physics are defined by us, which are the first-order differential system (FODS) and the second-order differential system (SODS). Basing on the conventional Legendre transformation, we obtain a new kind of canonical equations of Hamilton (CEH) with some kind of symmetry. We show that the FODS can only be expressed by the new CEH, but do not by the conventional CEH, while the SODS can be done by both the new and the conventional CEHs, on basis of the same conventional Legendre transformation. As an example, we prove that the nonlinear Schrödinger equation can be expressed with the new CEH in a consistent way. Based on the new CEH, the approximate soliton solution of the nonlocal nonlinear Schrödinger equation is obtained, and the soliton stability is analysed analytically as well. Furthermore, because the symmetry of a system is closely connected with certain conservation theorem of the system, the new CEH may be useful in some complicated systems when the symmetry considerations are used.

1. Introduction

The Hamiltonian viewpoint offers a theoretical framework in lots of physical fields. In classical mechanics, it constitutes the basis for further developments, including the Hamilton-Jacobi theory, the perturbation approaches and the chaos [1]. The canonical equations of Hamilton (CEH) in the field of classical mechanics are expressed as [1]
q i ˙ = H p i , p i ˙ = H q i ( i = 1 , , n )
where q i and p i are respectively the generalized coordinate and momentum, while q ˙ i = d q i / d t is the generalized velocity. The generalized momentum p i is defined as p i = L q ˙ i with L being the Lagrangian. By the Legendre transformation H = i = 1 n q ˙ i p i L , the Hamiltonian H is then obtained.
The CEH (1) can be extended to the continuous system [1]
q ˙ s = δ h δ p s , p ˙ s = δ h δ q s , ( s = 1 , , N )
where the subscript s indicates the components of the quantity of the continuous system [1], δ h δ q s = h q s x h q s , x and δ h δ p s = h p s x h p s , x denote the functional derivatives of h with respect to q s and p s with q s , x = q s x and p s , x = p s x , q s and p s are respectively the generalized coordinate and momentum, and h is the Hamiltonian density. The generalized momentum p s is defined as
p s = l q ˙ s .
The Hamiltonian density h is acquired by the following Legendre transformation
h = s = 1 N q ˙ s p s l ,
where l is the Lagrangian density. But it is significantly different for the continuous system that q s and p s are now the functions of both the time coordinate t and the spatial coordinate x. It should be noted that the spatial coordinate x is not the generalized coordinate, but is only the continuous index replacing the discrete i in Equation (1). To avoid confusion, we refer to time t as the evolution coordinate. h is a function of q s , p s and q s , x but not p s , x [1], so δ h δ p s = h p s , then the first equation of Equation (2) can be also expressed by
q ˙ s = h p s .
To our knowledge, the CEH in all current literatures are of the form (2), which are constructed on the basis of the second-order differential system (SODS). Besides the SODS, there are a lot of the first-order differential systems (FODS) to describe physical phenomena. Among them, the nonlinear Schrödinger equation (NLSE) is just the universal FODS. A question is raised in the nature of things: what is the form of the CEH applicable for the FODS? Are the conventional CEH, Equation (2), still valid for the FODS? In this paper, we gain a new CEH of the formal symmetry valid for the FODS, from which the NLSE can be expressed in a consistent manner. We also prove that the symmetric CEH is equivalent to the conventional CEH for the SODS. However, the conventional CEH can not model the FODS.
It is well known that the symmetry plays an important role in theoretical physics [2]. The search for and the discovery of new symmetries promote the exploration of fundamental laws of physics. Based on the idea, the CEH with the symmetry found by us, although this symmetry is only formal, might find their appropriate position in modern theoretical physics.

2. The CEH for the FODS

The Newton’s second law, the base of the Hamiltonian formulation, is modeled by a second-order differential equation of the evolution coordinate (the time coordinate). Here, we define the system governed by second-order partial differential equations of evolution coordinates as the SODS. Similarly, the FODS is the system governed by the first-order partial differential equations of evolution coordinates. The Lagrangian density of the SODS of the continuous systems is expressed in general as [1]
l = s = 1 N k = 1 N A s k q ˙ s q ˙ k + s = 1 N B s q ˙ s + C ,
where A s k , B s , C depend on not only q s but also q s , x in general. The generalized momentum can be obtained by the definition (3) as
p s = k = 1 N ( A s k + A k s ) q ˙ k + B s ,
which is the function of q s , q ˙ s and q s , x . Equation (7) in fact have N equations and contain 4 N variables, which are q s , q ˙ s , p s and q s , x . So the degree of freedom of Equation (7) is 3 N . Then we take q s , p s and q s , x as independent variables, and express the generalized velocities q ˙ s by these independent variables.
Besides the SODSs, there are a number of the FODSs. The nonlinear Schrödinger equation (NLSE)
i φ t + 1 2 2 φ x 2 + | φ | 2 φ = 0
is just a universal model that can be applied to hydrodynamics [3], nonlinear optics [4,5,6], nonlinear acoustics [7], Bose-Einstein condensates [8]. In nonlinear optics [4,5,6], the NLSE (8) governs the propagation of the slowly-varying light-envelope with the evolution coordinate t being the propagation direction coordinate. The light-envelope φ is a cw paraxial beam in a planar waveguide [5] or a narrow spectral-width pulse in optical fibers [4,6]. x is a transverse space coordinate for the beam and a frame moving at the group velocity (the so-called retarded frame) for the pulse, respectively.
For the FODS, the Lagrangian density should be a linear function of the generalized velocities q ˙ s . If the Lagrangian density is a quadratic function of the generalized velocities like Equation (6), the equation of motion, i.e., the Euler-Lagrange equation
t l q ˙ s δ l δ q s = 0
will be the second-order partial differential equation of the evolution coordinate t, which is in contradiction with the definition of the FODS. Therefore, the Lagrangian density of the FODS can only be expressed as
l = s = 1 N F s ( q s ) q ˙ s + Q ( q s , q s , x ) .
Besides, F s in Equation (10) is not the function of q s , x . If F s is also the function of q s , x , there will be such terms as q s , x q ˙ s appearing in Equation (10). Substitution of Equation (10) into Equation (9) leads to the appearance of the mixed partial derivative terms 2 q s x t . Via the coordinate rotation transform, terms 2 q s t 2 and 2 q s x 2 will appear instead. Therefore, the Euler-Lagrange equation (9) expressed with the canonical form of the second-order partial differential equation [9] is the second-order partial differential equation about the evolution coordinate t. According to our definition, the system is the SODS. Consequently, the generalized momentum p s , obtained by definition (3)
p s = F s ( q s ) ,
is only a function of q s . This is of significant difference from the case of the SODS, where the generalized momentum p s is the function of not only q s , but also q ˙ s and q s , x , as shown in Equation (7). Equation (11) have N equations, but contain 2 N variables, q s and p s . Therefore, the degree of freedom of the system described by Equation (11) is N. We can take q 1 , , q ν and p 1 , , p μ as the independent variables, where ν + μ = N . The rest of generalized coordinates and momenta can be expressed by the independent variables q α = f α ( q 1 , , q ν , p 1 , , p μ ) ( α = ν + 1 , , N ) , and p β = g β ( q 1 , , q ν , p 1 , , p μ ) ( β = μ + 1 , , N ) .
We now derive the CEH for the FODS. The total differential of the Hamiltonian density h can be obtained by using Equation (4)
d h = s = 1 N p s d q ˙ s + η = 1 μ q ˙ η d p η + β = μ + 1 N q ˙ β λ = 1 ν g β q λ d q λ + η = 1 μ g β p η d p η d l ,
while the total differential of the Lagrangian density l ( q s , q ˙ s , q s , x ) with respect to its arguments is
d l = λ = 1 ν l q λ d q λ + α = ν + 1 N l q α λ = 1 ν f α q λ d q λ + η = 1 μ f α p η d p η + s = 1 N l q ˙ s d q ˙ s + s = 1 N l q s , x d q s , x .
Substitution of Equation (13) into Equation (12) yields
d h = λ = 1 ν β = μ + 1 N q ˙ β g β q λ α = ν + 1 N l q α f α q λ l q λ d q λ + η = 1 μ q ˙ η + β = μ + 1 N q ˙ β g β p η α = ν + 1 N l q α f α p η d p η s = 1 N l q s , x d q s , x .
Since the total differential of h ( q 1 , , q ν , p 1 , , p μ , q s , x ) with respect to its arguments can be written as d h = λ = 1 ν h q λ d q λ + η = 1 μ h p η d p η + s = 1 N h q s , x d q s , x , by comparing this equation with Equation (14), we obtain 2 N equations
h q λ = β = μ + 1 N q ˙ β g β q λ α = ν + 1 N l q α f α q λ l q λ ,
h p η = q ˙ η + β = μ + 1 N q ˙ β g β p η α = ν + 1 N l q α f α p η ,
h q s , x = l q s , x ,
where λ = 1 , , ν , η = 1 , , μ , s = 1 , , N . From the Euler-Lagrange Equation (9), we obtain l q s = t l q ˙ s + x l q s , x , substitution of which into Equations (15) and (16) yields
h q λ = p ˙ λ + β = μ + 1 N q ˙ β g β q λ α = ν + 1 N p ˙ α f α q λ α = ν + 1 N x l q α , x f α q λ x l q λ , x ,
h p η = q ˙ η + β = μ + 1 N q ˙ β g β p η α = ν + 1 N p ˙ α f α p η α = ν + 1 N x l q α , x f α p η .
Then substituting Equation (17) into Equations (18) and (19), we obtain N CEHs for the FODS
δ h δ q λ = p ˙ λ + β = μ + 1 N q ˙ β g β q λ α = ν + 1 N p ˙ α f α q λ + α = ν + 1 N x h q α , x f α q λ ,
δ h δ p η = q ˙ η + β = μ + 1 N q ˙ β g β p η α = ν + 1 N p ˙ α f α p η + α = ν + 1 N x h q α , x f α p η .
To obtain Equation (21), we have used δ h δ p η = h p η , because h is not a function of p η , x . The CEH, Equations (20) and (21), can be expressed in a symmetric form as
δ h δ q λ = s = 1 N q ˙ s p s q λ p ˙ s q s q λ + α = ν + 1 N x h q α , x f α q λ ,
δ h δ p η = s = 1 N q ˙ s p s p η p ˙ s q s p η + α = ν + 1 N x h q α , x f α p η
( λ = 1 , , ν , η = 1 , , μ , and ν + μ = N ), because p ˙ λ = λ = 1 ν p ˙ λ q λ q λ , q ˙ η = η = 1 μ q ˙ η p η p η , η = 1 μ q ˙ η p η q λ = 0 , and λ = 1 ν p ˙ λ q λ p η = 0 . The CEH above can be extended to the discrete system
H q λ = s = 1 N q ˙ s p s q λ p ˙ s q s q λ , H p η = s = 1 N q ˙ s p s p η p ˙ s q s p η ,
where λ = 1 , , ν , η = 1 , , μ , and ν + μ = N .
We will prove that the symmetric CEH (22) and (23) can also describe the SODS. In Equations (22) and (23), if all the generalized coordinates q s and momenta p s are independent, we can obtain that g β q λ = f α q λ = g β p η = f α p η = 0 ( α = ν + 1 , , N , β = μ + 1 , , N ), then the CEH will be reduced to Equation (2). It does be the case of the SODS, where all the generalized coordinates and momenta are independent. So, the symmetric CEH obtained by us can express both the FODS and the SODS. In other words, the new CEH and the conventional CEH are equivalent when describing the SODS, but the former are of some formally symmetry. The conventional CEH, Equation (2), can only be used to expresses the SODS.

3. Application of the Symmetric CEH for Continuous Systems to the NLSE

In this section, we will apply the new CEH with symmetry, Equations (22) and (23), to the NLSE. The Lagrangian density for the NLSE is stated as [10] l = i 2 ( φ * φ t φ φ * t ) + 1 2 | φ x | 2 1 2 | φ | 4 . The NLSE is complex, so it is in fact an equation of two real functions, one is of its real part and the other is of its imaginary part. Alternatively, we can take the fields φ and φ * as two independent functions. In this sense, the components of the quantity N for the NLSE is equal to two. In other words, for the NLSE there are two generalized coordinates, q 1 = φ * and q 2 = φ , and two generalized momenta
p 1 = i 2 φ , p 2 = i 2 φ * .
By using Equation (4) [4], we can obtain the Hamiltonian density
h = 1 2 φ x 2 + 1 2 | φ | 4 .
If we take q 1 and p 1 as the independent variables, q 2 and p 2 can be expressed by the relations (25) as q 2 = 2 i p 1 and p 2 = i 2 q 1 , respectively. It also should be noted that h is also the function of q s , x . Then, we can express the Hamiltonian density (26) with independent variables q 1 , p 1 , q 1 , x and q 2 , x as
h = 1 2 q 1 , x q 2 , x 2 q 1 2 p 1 2 .
For the NLSE, ν = μ = 1 and N = 2 , therefore the Equation (22) in fact have only one equation. It is the same case for Equation (23). Consequently, for the NLSE, the CEH (22) and (23) can produce two equations. From the left side of Equation (22), we can obtain δ h δ q 1 = 1 2 2 φ x 2 + | φ | 2 φ . From its right side, we can have p ˙ 1 + q ˙ 2 p 2 q 1 = i φ ˙ . Then the NLSE (8) is obtained. While, for the other CEH (23), the left side is
δ h δ p 2 = 4 q 1 2 p 1 = 2 i | φ | 2 φ * ,
and the right side is q ˙ 1 p ˙ 2 q 2 p 1 + x h q 2 , x q 2 p 1 = 2 φ ˙ * + i 2 φ * x 2 , which results in the generation of the complex conjugate of the NLSE. As a result, the CEHs (22) and (23) are consistent. From one of the two CEHs, the NLSE is expressed; from the other, the complex conjugate of the NLSE is expressed.
We now demonstrate that the conventional CEH (2) can not be used to express the NLSE. By Equation (3), we have p φ = l / φ ˙ = i / 2 φ * . Substituting the Hamiltonian density (26) into the second equation of Equation (2) only yields i 2 φ t + 1 2 2 φ + | φ | 2 φ = 0 , which does not be the NLSE (8). By substituting the Hamiltonian density (26) into Equation (5), we can obtain the left side φ * t , and the right side h p φ * = h φ φ p φ * = 2 i | φ | 2 φ * , where p φ * = l φ ˙ * . Then the equation i 2 φ * t + φ 2 φ * = 0 can be obtained. It is surely not the complex conjugate of the NLSE.

4. Application of the Symmetric CEH for Discrete Systems to Light-Envelope Propagations

In this section, we use the symmetric CEH for discrete systems, i.e., Equation (24), to discuss light-envelope propagations in nonlocally nonlinear media, which is modeled by the following (1+D)-dimensional nonlocal nonlinear Schrödinger equation (NNLSE) [11,12,13,14]
i φ z + 2 φ + Δ n φ = 0 .
Δ n ( r , z ) is the nonlinear refractive index, which can phenomenologically be expressed as a convolution between the response function R ( r ) and the light intensity
Δ n ( r , z ) = R ( r r ) | φ ( r , z ) | 2 d D r .
When R is the Dirac delta function, the NNLSE (29) is reduced to the NLSE (8).
We assume the trial solution of Equation (29)
φ ( r , z ) = q A ( z ) exp r 2 q w 2 ( z ) exp i q c ( z ) r 2 + i q θ ( z ) ,
where q A , q θ are the amplitude and phase of the complex amplitude φ , respectively, q w is the width, q c is the phase-front curvature, and they all vary with z. We consider the response function
R ( r ) = 1 ( π w m ) D exp r 2 w m 2 .
Substituting the trial solution (31) into the Lagrangian density l = i 2 φ * φ z φ φ * z φ 2 + 1 2 | φ ( r , z ) | 2 Δ n ( r , z ) , and performing the integration L = l d D r , we have
L = 2 2 D π D / 2 q A 2 q w 2 + D ( w m 2 + q w 2 ) D / 2 [ 2 q A 2 q w 2 + D + 2 D / 2 ( w m 2 + q w 2 ) D / 2 ( 4 D + 4 D q c 2 q w 4 + D q w 4 q ˙ c + 4 q w 2 q ˙ θ ) ] ,
which is the function of generalized coordinates q A , q w , q c and velocities q ˙ c , q ˙ θ .
Then the generalized momenta can be obtained
p A = p w = 0 , p c = 2 2 D 2 D π D / 2 q A 2 q w 2 + D , p θ = π 2 D / 2 q A 2 q w D .
By Legendre transformation, the Hamiltonian can be determined
H = 2 1 D π D / 2 q A 2 q w 2 + D ( w m 2 + q w 2 ) D / 2 [ q A 2 q w 2 + D + 2 1 + D 2 D ( w m 2 + q w 2 ) D / 2 ( 1 + q c 2 q w 4 ) ] ,
which will be proved to be a constant, that is H ˙ = 0 .
Four generalized coordinates and four generalized momenta are contained in the four equations (34). It indicates that Equation (34) have four degrees of freedom. Here, we can take q c , q θ , p c and p θ as independent variables. By solving Equation (34), we can express the generalized coordinates q A and q w by generalized momenta p c and p θ as q A = ( p θ ) 1 / 2 [ D p θ / ( 2 π p c ) ] D / 4 and q w = [ 4 p c / ( D p θ ) ] 1 / 2 . Then, the Hamiltonian (35) is rewritten as
H = D 2 p θ 2   +   16 p c 2 q c 2 4 p c 1 2 π D / 2 ( 4 p c D p θ + w m 2 ) D / 2 .
By using the CEH (24) with μ = ν = 2 and n = 4 , we can obtain four equations as
q ˙ c = D 2 p θ 2 4 p c 2 4 q c 2   +   D π D / 2 p θ 2 ( 4 p c D p θ   +   w m 2 ) D / 2 4 p c   +   D p θ w m 2 ,
q ˙ θ = ( 4   +   D ) π D / 2 p c p θ ( 4 p c D p θ   +   w m 2 ) D / 2 4 p c   +   D p θ w m 2 D 2 p θ 2 p c D π D / 2 p θ 2 w m 2 ( 4 p c D p θ   +   w m 2 ) D / 2 4 p c   +   D p θ w m 2 ,
p ˙ c = 8 p c q c ,
p ˙ θ = 0 .
One can find that the generalized coordinate q θ does not appear in the Hamiltonian (36), therefore q θ is a cyclic coordinate. Because the generalized momentum conjugate to a cyclic coordinate is conserved [1], the generalized momentum p θ conjugate to the generalized coordinate q θ is a constant, which can be confirmed by Equation (40). In fact, this represents that the power P 0 = φ 2 d D r = q A 2 ( π / 2 q w ) D is conservative. Then we can obtain
q A 2 = P 0 ( π / 2 q w ) D .
Taking derivative with respect to z on two sides of the third equation of (34), then comparing it with Equation (39) we obtain q c = q ˙ w 4 q w , the substitution of which into the Hamiltonian (35) yields H = T + V , where
T = 1 16 D P 0 q ˙ w 2 , V = D P 0 q w 2 1 2 π D / 2 P 0 2 w m 2 + q w 2 D / 2
are the generalized kinetic energy and potential, respectively.
From the Hamiltonian point of view, the dynamics of light-envelopes in nonlinear media can be regarded as a problem of small oscillations of a Hamiltonian system about its equilibrium position. The equilibrium state of the system described by the Hamiltonian H corresponds to the soliton solutions of the NNLSE, which can be gained as the extremum points of generalized potential V. The equilibrium position is stable if a small disturbance from equilibrium leads to small bounded motion about the rest position. While, if an infinitesimal disturbance produces unbounded motion, the equilibrium is unstable [1]. In other words, the equilibrium must be stable when the extremum of the generalized potential is a minimum, and unstable otherwise. In this sense, therefore, the viewpoint in a few literatures [15,16,17,18], where solitons were taken for the extremum of the Hamiltonian rather than the generalized potential, might be somewhat ambiguous. In these literatures [15,16,17,18] the trial solution has an invariable profile (soliton profile), and the soliton state is the static state in fact. In this case, the kinetic energy is zero, and the Hamiltonian is just equal to the potential. In this connection, for the static system the extrema of the Hamiltonian and the generalized potential are equal only in value. Although in such literatures [15,16,17,18] the obtained soliton solutions are correct, it is more reasonable to regard the soliton solutions as the extremum points of the generalized potentials rather than the Hamiltonian.
To find the equilibrium state (soliton solution), we let V / q w = 0 , then we obtain 32 q w 3 + 8 π D / 2 P 0 q w w m 2 + q w 2 1 D 2 = 0 . The critical power is then obtained as
P c = 4 π D / 2 w m 2 + q w 2 1 + D 2 q w 4 ,
with which the light can propagate with a changeless profile. Besides, P 0 = P c also leads to q ˙ c = q c = 0 , which indicates that the soliton wavefront is a plane.
Then we address the soliton stability by discussing the generalized potential V. Performing the second-order derivative of V with respect to q w , and substituting P c into it, we have
Υ 2 V q w 2 P 0 = P c = 64 q w 4 2 2 + D 2 1 + σ 2 ,
where σ = w m / q w represents the degree of nonlocality. When Υ > 0 , generalized potential V has a minimum, then the soliton is stable. From Equation (44) we can get the stability criterion of solitons
σ 2 > 1 4 ( D 2 ) ,
which also agrees with the Vakhitov-Kolokolov criterion [19].

4.1. The Local Case

When w m 0 , R ( r ) δ ( r ) , and the NNLSE is reduced to the NLSE (8). Then, Equations (43) and (44) will be reduced to P c = 4 π D / 2 q w D 2 , Υ = 32 q w 4 ( 2 D ) . In the case of D = 1 , P c = 4 π / q w , it is the same as Equation (42) of Ref. [10]. In the case of D = 2 , P c = 4 π , it is the same as Equation (16a) of Ref. [20]. We can obtain Υ > 0 if D < 2 , Υ < 0 if D > 2 , and Υ = 0 if D = 2 . Therefore, in the local case, the soliton is stable if D = 1 , but unstable if D > 2 . In the case of D = 2 , it needs further analysis because Υ = 0 . In the case of D = 2 , the generalized potential V = 4 π P 0 P 0 / 2 π q w 2 , which does not have extreme when P 0 4 π . When P 0 = P c = 4 π , we obtain that V = 0 . It is the extreme rather than the minimum. Hence, (1+2)-dimensional local solitons are always unstable. When P 0 = P c , the potential V = 0 is constant, the light-envelope without the external disturbance will stay in its initial state. But the ideal condition can not occur in experiment. If the external disturbance makes P 0 > P c , the beam will become narrower and narrower, and the optical beam will eventually collapse. If the external disturbance makes P 0 < P c , the optical beam will diffract at last. These conclusions all agree with those of Refs. [21,22,23].

4.2. The Nonlocal Case

When w m 0 and D 2 , the condition (45) is satisfied automatically. It means that the (1+1)-dimensional and the (1+2)-dimensional solitons are always stable in nonlocally nonlinear media of Gaussian response. It is consistent with the conclusion of Ref. [24]. When D > 2 the stabe solitons can exist when the degree of nonlocality should be strong enough that satisfies the criterion (45), which is also the same as Ref. [24].

5. Conclusions

We obtain a new CEH with symmetry in form. It can express both the FODS and the SODS, while the conventional CEH can only express the SODS but impossibly express the FODS on basis of the same conventional Legendre transformation. By using the new CEH we analytically obtain the approximate solution of the NNLSE and analytically discuss the stability of the solitons. We will use the new CEH to analytically deal with two important problems related to solitons in diverse nonlinear systems, that is, the existence and the stability, which will save amounts of time for the numerical explorations.

Author Contributions

Writing—original draft preparation, writing—review and editing, G.L.; conceptualization, methodology, X.C. and Z.R.; formal analysis, G.L.; conceptualization, methodology, supervision, project administration, Q.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by “Natural Science Foundation of Guangdong Province of China, grant number 2021A1515012214” and “The APC was funded by Natural Science Foundation of Guangdong Province of China, grant number 2021A1515012214”.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Goldstein, H.; Poole, C.; Safko, J. Classical Mechanics; Addison-Wesley: Boston, MA, USA, 2001. [Google Scholar]
  2. Gross, D.J. The role of symmetry in fundamental physics. Proc. Natl. Acad. Sci. USA 1996, 93, 14256. [Google Scholar] [CrossRef]
  3. Nore, C.; Brachet, M.E.; Fauve, S. Numerical study of hydrodynamics using the nonlinear Schrödinger equation. Phys. D 1993, 65, 154. [Google Scholar] [CrossRef]
  4. Hasegawa, A.; Kodama, Y. Solitons in Optical Communications; Clarenoon Press: Oxford, UK, 1995. [Google Scholar]
  5. Haus, H.A. Waves and Fields in Optoelectronics; Prentice-Hall: Englewood Cliffs, NJ, USA, 1984; pp. 277–283. [Google Scholar]
  6. Agrawal, G. Nonlinear Fiber Optics; Academic Press: San Diego, CA, USA, 2001. [Google Scholar]
  7. Bisyarin, M.A. Weak-nonlinear acoustic pulse dynamics in a waveguide channel with longitudinal inhomogeneity. AIP Conf. Proc. 2008, 1022, 38. [Google Scholar]
  8. Seaman, B.T.; Carr, L.D.; Holland, M.J. Nonlinear band structure in Bose-Einstein condensates: Nonlinear Schrödinger equation with a Kronig-Penney potential. Phys. Rev. A 2005, 71, 033622. [Google Scholar] [CrossRef]
  9. Arfken, G.B.; Weber, H.J. Mathematical Methods for Physicists; Academic Press: Cambridge, MA, USA, 2005; p. 538. [Google Scholar]
  10. Anderson, D. Variational approach to nonlinear pulse propagation in optical fibers. Phys. Rev. A 1983, 27, 3135. [Google Scholar] [CrossRef]
  11. Guo, Q.; Lu, D.; Deng, D. Nonlocal spatial optical solitons. In Advances in Nonlinear Optics; Chen, X., Zhang, G., Zeng, H., Guo, Q., She, W., Eds.; De Gruyter: Berlin, Germany, 2015; Chapter 4; pp. 227–305. [Google Scholar]
  12. Snyder, A.W.; Mitchell, D.J. Accessible Solitons. Science 1997, 276, 1538. [Google Scholar] [CrossRef]
  13. Krolikowski, W.; Bang, O.; Rasmussen, J.J.; Wyller, J. Modulational instability in nonlocal nonlinear Kerr media. Phys. Rev. E 2001, 64, 016612. [Google Scholar] [CrossRef] [PubMed]
  14. Guo, Q.; Hu, W.; Deng, D.; Lu, D.; Ouyang, S. Features of strongly nonlocal spatial solitons. In Nematicons: Spatial Optical Solitons in Nematic Liquid Crystals; Assanto, G., Ed.; John Wiley & Sons: New York, NY, USA, 2012; Chapter 2. [Google Scholar]
  15. Seghete, V.; Menyuk, C.R.; Marks, B.S. Solitons in the midst of chaos. Phys. Rev. A 2007, 76, 043803. [Google Scholar] [CrossRef]
  16. Picozzi, A.; Garnier, J. Incoherent soliton turbulence in nonlocal nonlinear media. Phys. Rev. Lett. 2011, 107, 233901. [Google Scholar] [CrossRef] [PubMed]
  17. Lashkin, V.M.; Yakimenkoa, A.I.; Prikhodko, O.O. Two-dimensional nonlocal multisolitons. Phys. Lett. A 2007, 366, 422. [Google Scholar] [CrossRef]
  18. Petroski, M.M.; Petrović, M.S.; Belić, M.R. Quasi-stable propagation of vortices and soliton clusters in saturable Kerr media with square-root nonlinearity. Opt. Commun. 2007, 279, 196. [Google Scholar] [CrossRef]
  19. Vakhitov, N.G.; Kolokolov, A.A. Stationary solutions of the wave equation in a medium with nonlinearity saturation. Radiophys. Quantum Electron. 1973, 16, 783. [Google Scholar] [CrossRef]
  20. Desaix, M.; Anderson, D.; Lisak, M. Variational approach to collapse of optical pulses. J. Opt. Soc. Am. B 1991, 8, 2082. [Google Scholar] [CrossRef]
  21. Berge, L. Wave collapse in physics: Principles and applications to light and plasma waves. Phys. Rep. 1998, 303, 259. [Google Scholar] [CrossRef]
  22. Moll, K.D.; Gaeta, A.L.; Fibich, G. Self-similar optical wave collapse: Observation of the Townes profile. Phys. Rev. Lett. 2003, 90, 203902. [Google Scholar] [CrossRef] [PubMed]
  23. Sun, C.; Barsi, C.; Fleischer, J.W. Peakon profiles and collapse-bounce cycles in self-focusing spatial beams. Opt. Express 2008, 16, 20676. [Google Scholar] [CrossRef] [PubMed]
  24. Bang, O.; Krolikowski, W.; Wyller, J.; Rasmussen, J.J. Collapse arrest and soliton stabilization in nonlocal nonlinear media. Phys. Rev. E 2002, 66, 046619. [Google Scholar] [CrossRef] [PubMed]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Liang, G.; Chen, X.; Ren, Z.; Guo, Q. Canonical Equations of Hamilton with Symmetry and Their Applications. Symmetry 2024, 16, 305. https://doi.org/10.3390/sym16030305

AMA Style

Liang G, Chen X, Ren Z, Guo Q. Canonical Equations of Hamilton with Symmetry and Their Applications. Symmetry. 2024; 16(3):305. https://doi.org/10.3390/sym16030305

Chicago/Turabian Style

Liang, Guo, Xiangwei Chen, Zhanmei Ren, and Qi Guo. 2024. "Canonical Equations of Hamilton with Symmetry and Their Applications" Symmetry 16, no. 3: 305. https://doi.org/10.3390/sym16030305

APA Style

Liang, G., Chen, X., Ren, Z., & Guo, Q. (2024). Canonical Equations of Hamilton with Symmetry and Their Applications. Symmetry, 16(3), 305. https://doi.org/10.3390/sym16030305

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop