
Conservation laws are formulated for systems of differential equations by
using symmetries and adjoint symmetries, and an application to systems of
evolution equations is made, together with illustrative examples. The
formulation does not require the existence of a Lagrangian for a given system,
and the presented examples include computations of conserved densities for the
heat equation, Burgers' equation and the Kortewegde Vries equation.

In this paper, a new generalized $5\times5$ matrix spectral problem of
AblowitzKaupNewellSegur(AKNS) type associated with the enlarged matrix Lie
super algebra is proposed and its corresponding super soliton hierarchy is
established. The super variational identities is used to furnish
superHamiltonian structures for the resulting super soliton hierarchy.

Integrable couplings are associated with nonsemisimple Lie algebras. In this
paper, we propose a new method to generate new integrable systems through
making perturbation in matrix spectral problems for integrable couplings, which
is called the `completion process of integrable couplings'. As an example, the
idea of construction is applied to the AblowitzKaupNewellSegur integrable
coupling. Each equation in the resulting hierarchy has a biHamiltonian
structure furnished by the componenttrace identity.

Using linear combinations of Lax matrices of soliton hierarchies, we
introduce trigonal curves by their characteristic equations, and determine
Dubrovin type equations for zeros and poles of meromorphic functions defined as
ratios of the BakerAkhiezer functions. We straighten out all flows in soliton
hierarchies under the AbelJacobi coordinates associated with Lax pairs, and
generate algebrogeometric solutions to soliton hierarchies in terms of the
Riemann theta functions, through observing asymptotic behaviors of the
BakerAkhiezer functions. We analyze the fourcomponent AKNS soliton hierarchy
in such a way that it leads to a general theory of trigonal curves applicable
to construction of algebrogeometric solutions of an arbitrary soliton
hierarchy.

Lump solutions are analytical rational function solutions localized in all
directions in space. We analyze a class of lump solutions, generated from
quadratic functions, to nonlinear partial differential equations. The basis of
success is the Hirota bilinear formulation and the primary object is the class
of positive multivariate quadratic functions. A complete determination of
quadratic functions positive in space and time is given, and positive quadratic
functions are characterized as sums of squares of linear functions. Necessary
and sufficient conditions for positive quadratic functions to solve Hirota
bilinear equations are presented, and such polynomial solutions yield lump
solutions to nonlinear partial differential equations under the dependent
variable transformations u=2(ln f)_x and u=2(ln f)_{xx}, where x is one spatial
variable. Applications are made for a few generalized KP and BKP equations.

Two multicomponent generalizations of the AKNStype spectral problems
associated with $sl(2,\mathbb{R})$ and $so(3,\mathbb{R})$ are introduced and
the corresponding two hierarchies of generalized multicomponent AKNStype
soliton equations are presented by the standard procedure, respectively. By
virtue of the trace identity, biHamiltonian structures which lead to a common
recursion operator are established for each of the two resulting soliton
hierarchies. And thus the Liouville integrability is shown for all systems in
each of the two new generalized soliton hierarchies, seperately.

The two matrix spectral problems of AblowitzKaupNewellSegur (AKNS) and
KaupNewell (KN) types associated with so(3,R) are generalized. The
corresponding hierarchies of generalized soliton equations are derived by the
standard procedure using the zero curvature formulation. Recursion operators
and biHamiltonian structures are explicitly constructed for the resulting two
generalized soliton hierarchies of AKNS and KN types, which shows their
Liouville integrability.

A counterpart of the WadatiKonnoIchikawa (WKI) soliton hierarchy,
associated with so(3,R), is presented through the zero curvature formulation.
Its spectral matrix is defined by the same linear combination of basis vectors
as the WKI one, and its Hamiltonian structures yielding Liouville integrability
are furnished by the trace identity.

The invariant subspace method is refined to present more unity and more
diversity of exact solutions to evolution equations. The key idea is to take
subspaces of solutions to linear ordinary differential equations as invariant
subspaces that evolution equations admit. A twocomponent nonlinear system of
dissipative equations was analyzed to shed light on the resulting theory, and
two concrete examples are given to find invariant subspaces associated with
2ndorder and 3rdorder linear ordinary differential equations and their
corresponding exact solutions with generalized separated variables.

A multiple expfunction method to exact multiple wave solutions of nonlinear
partial differential equations is proposed. The method is oriented towards ease
of use and capability of computer algebra systems, and provides a direct and
systematical solution procedure which generalizes Hirota's perturbation scheme.
With help of Maple, an application of the approach to the $3+1$ dimensional
potentialYuTodaSasaFukuyama equation yields exact explicit 1wave and
2wave and 3wave solutions, which include 1soliton, 2soliton and 3soliton
type solutions. Two cases with specific values of the involved parameters are
plotted for each of 2wave and 3wave solutions.

An explicit Bargmann symmetry constraint is computed and its associated
binary nonlinearization of Lax pairs is carried out for the super Dirac
systems. Under the obtained symmetry constraint, the nth flow of the super
Dirac hierarchy is decomposed into two super finitedimensional integrable
Hamiltonian systems, defined over the supersymmetry manifold $R^{4N2N}$ with
the corresponding dynamical variables $x$ and $t_n$. The integrals of motion
required for Liouville integrability are explicitly given.

A fivedimensional symmetry algebra consisting of Lie point symmetries is
firstly computed for the nonlinear Schroedinger equation, which, together with
a reflection invariance, generates two fiveparameter solution groups. Three
ansaetze of transformations are secondly analyzed and used to construct exact
solutions to the nonlinear Schroedinger equation. Various examples of exact
solutions with constant, trigonometric function type, exponential function type
and rational function amplitude are given upon careful analysis. A bifurcation
phenomenon in the nonlinear Schroedinger equation is clearly exhibited during
the solution process.

A direct approach to exact solutions of nonlinear partial differential
equations is proposed, by using rational function transformations. The new
method provides a more systematical and convenient handling of the solution
process of nonlinear equations, unifying the tanhfunction type methods, the
homogeneous balance method, the expfunction method, the mapping method, and
the Fexpansion type methods. Its key point is to search for rational solutions
to variablecoefficient ordinary differential equations transformed from given
partial differential equations. As an application, the construction problem of
exact solutions to the 3+1 dimensional JimboMiwa equation is treated, together
with a B\"acklund transformation.

Riemann theta functions are used to construct oneperiodic and twoperiodic
wave solutions to a class of (2+1)dimensional Hirota bilinear equations. The
basis for the involved solution analysis is the Hirota bilinear formulation,
and the particular dependence of the equations on independent variables
guarantees the existence of oneperiodic and twoperiodic wave solutions
involving an arbitrary purely imaginary Riemann matrix. The resulting theory is
applied to two nonlinear equations possessing Hirota bilinear forms:
$u_t+u_{xxy}3uu_y3u_xv=0$ and $u_t+u_{xxxxy}(5u_{xx}v+10u_{xy}u15u^2v)_x=0$
where $v_x=u_y$, thereby yielding their oneperiodic and twoperiodic wave
solutions describing one dimensional propagation of waves.

A general Casoratian formulation is proposed for the 2D Toda lattice
equation, which involves coupled eigenfunction systems. Various Casoratian type
solutions are generated, through solving the resulting linear conditions and
using a Baecklund transformation.

The discrete variational identity under general bilinear forms on semidirect
sums of Lie algebras is established. The constant $\gamma$ involved in the
variational identity is determined through the corresponding solution to the
stationary discrete zero curvature equation. An application of the resulting
variational identity to a class of semidirect sums of Lie algebras in the
Volterra lattice case furnishes Hamiltonian structures for the associated
integrable couplings of the Volterra lattice hierarchy.

Based on the Orlov and Shulman's M operator, the additional symmetries and
the string equation of the CKP hierarchy are established, and then the higher
order constraints on $L^l$ are obtained. In addition, the generating function
and some properties are also given. In particular, the additional symmetry
flows form a new infinite dimensional algebra $W^C_{1+\infty}$, which is a
subalgebra of $W_{1+\infty}$.

A relation between semidirect sums of Lie algebras and integrable couplings
of lattice equations is established, and a practicable way to construct
integrable couplings is further proposed. An application of the resulting
general theory to the generalized Toda spectral problem yields two classes of
integrable couplings for the generalized Toda hierarchy of lattice equations.
The construction of integrable couplings using semidirect sums of Lie algebras
provides a good source of information on complete classification of integrable
lattice equations.

A broad set of sufficient conditions consisting of systems of linear partial
differential equations is presented which guarantees that the Wronskian
determinant solves the Kortewegde Vries equation in the bilinear form. A
systematical analysis is made for solving the resultant linear systems of
secondorder and thirdorder partial differential equations, along with
solution formulas for their representative systems. The key technique is to
apply variation of parameters in solving the involved nonhomogeneous partial
differential equations. The obtained solution formulas provide us with a
comprehensive approach to construct the existing solutions and many new
solutions including rational solutions, solitons, positons, negatons,
breathers, complexitons and interaction solutions of the Kortewegde Vries
equation.

Complexiton solutions (or complexitons for short) are exact solutions newly
introduced to integrable equations. Starting with the solution classification
for a linear differential equation, the Kortewegde Vries equation and the Toda
lattice equation are considered as examples to exhibit complexiton structures
of nonlinear integrable equations. The crucial step in the solution process is
to apply the Wronskian and Casoratian techniques for Hirota's bilinear
equations. Correspondence between complexitons of the Kortewegde Vries
equation and complexitons of the Toda lattice equation is provided.

A set of coupled conditions consisting of differentialdifference equations
is presented for Casorati determinants to solve the Toda lattice equation. One
class of the resulting conditions leads to an approach for constructing
complexiton solutions to the Toda lattice equation through the Casoratian
formulation. An analysis is made for solving the resulting system of
differentialdifference equations, thereby providing the general solution
yielding eigenfunctions required for forming complexitons. Moreover, a feasible
way is presented to compute the required eigenfunctions, along with examples of
real complexitons of lower order.

A bridge going from Wronskian solutions to generalized Wronskian solutions of
the Kortewegde Vries equation is built. It is then shown that generalized
Wronskian solutions can be viewed as Wronskian solutions. The idea is used to
generate positons, negatons and their interaction solutions to the Kortewegde
Vries equation. Moreover, general positons and negatons are constructed through
the Wronskian formulation. A few new exact solutions to the KdV equation are
explicitly presented as examples of Wronskian solutions.

A systematic construction of St\"{a}ckel systems in separated coordinates and
its relation to biHamiltonian formalism are considered. A general form of
related hydrodynamic systems, integrable by the HamiltonJacobi method, is
derived. One Casimir biHamiltonian case is studed in details and in this case,
a systematic construction of related hydrodynamic systems in arbitrary
coordinates is presented, using a cofactor method and soliton symmetry
constraints.

Binary symmetry constraints are applied to the nonlinearization of spectral
problems and adjoint spectral problems into socalled binary constrained flows,
which provide candidates for finitedimensional Liouville integrable
Hamiltonian systems. The resulting constraints on the potentials of spectral
problems give rise to a kind of involutive solutions to zero curvature
equations, and thus the integrability by quadratures can be shown for zero
curvature equations once the corresponding binary constrained flows are proved
to be integrable. The whole process to carry out binary symmetry constraints is
called binary nonlinearization. The principal task of binary nonlinearization
is to expose the Liouville integrability for the resulting binary constrained
flows, which can usually be achieved as a consequence of the existence of
hereditary recursion operators. The theory of binary nonlinearization is
applied to the multiwave interaction equations associated with a 4x4 matrix
spectral problem as an illustrative example. The Backlund transformations
resulted from symmetry constraints are given for the multiwave interaction
equations, and thus a kind of involutive solutions is presented and the
integrability by quadratures is shown for the multiwave interaction equations.

A biHamiltonian formulation is proposed for triangular systems resulted by
perturbations around solutions, from which infinitely many symmetries and
conserved functionals of triangular systems can be explicitly constructed,
provided that one operator of the Hamiltonian pair is invertible. Through our
formulation, four examples of triangular systems are exhibited, which also show
that biHamiltonian systems in both lower dimensions and higher dimensions are
many and varied. Two of four examples give local 2+1 dimensional biHamiltonian
systems and illustrate that multiscale perturbations can lead to
higherdimensional biHamiltonian systems.